Notes de Cours-Gravimetrie PDF

GLQ2200 - Géophysique Appliquée I
Bernard Giroux, Michel Chouteau
Notes de cours - Gravimétrie
Laboratoire de géophysique appliquée

Automne 2008
Table des matières
1 Théorie 1
1.1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1
1.2 Principes de base . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1
1.2.1 Lois de l’attraction universelle . . . . . . . . . . . . . . . . . . . . . . . . 1
1.2.2 Champ gravitationnel . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2
1.2.3 Potentiel gravitationnel . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2
1.2.4 Équations du potentiel . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5
1.3 Une référence pour la terre . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5
1.3.1 Un ellipsoı̈de de révolution : le sphéroı̈de . . . . . . . . . . . . . . . . . . 5
1.3.2 Le géoı̈de . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7
1.4 Densité des roches . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7
2 Les données gravimétriques 11

2.1 Corrections et références . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11
2.1.1 Correction de dérive . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11
2.1.2 Correction de latitude . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 14
2.1.3 Correction d’altitude . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15
2.1.4 Correction de plateau . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 16
2.1.5 Correction de relief . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 17
2.1.6 Méthode de Nettleton . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 23
2.1.7 Anomalie Bouguer . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 25
2.2 Levé gravimétrique . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 26
2.2.1 Numérotation des stations . . . . . . . . . . . . . . . . . . . . . . . . . . . 26
2.2.2 Nivellement . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 27
i
ii Table des matières
2.2.3 Résumé pour faire un levé gravimétrique . . . . . . . . . . . . . . . . . . 29

2.3 Instrumentation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 30
2.3.1 Mesures absolues . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 30
2.3.2 Mesures relatives . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 32
2.4 Traitements des données . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 36
2.4.1 Séparation régionale - résiduelle . . . . . . . . . . . . . . . . . . . . . . . 36
2.4.2 Prolongement vers le haut . . . . . . . . . . . . . . . . . . . . . . . . . . . 38
2.4.3 Exemples de superposition d’une anomalie avec une régional . . . . . . 39
2.4.4 Cône de sources . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 43
3 Interprétation 45
3.1 Modèles simples . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 45
3.1.1 La sphère . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 45
3.1.2 Le cylindre horizontal . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 47
3.1.3 Le cylindre vertical . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 48
3.1.4 Le feuillet vertical . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 49
3.1.5 La plaque mince horizontale infinie . . . . . . . . . . . . . . . . . . . . . 49
3.1.6 Le prisme rectangulaire . . . . . . . . . . . . . . . . . . . . . . . . . . . . 51
3.2 Modèles complexes . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 51
3.2.1 Les méthodes graphiques . . . . . . . . . . . . . . . . . . . . . . . . . . . 51
3.2.2 Méthode analytique . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 53
3.2.3 Gravité 3-D . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 53
3.3 Excès de masse . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 56
3.3.1 Calcul de l’excès de masse . . . . . . . . . . . . . . . . . . . . . . . . . . . 56
3.3.2 Unicité de la solution . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 57
3.3.3 Considérations pratique . . . . . . . . . . . . . . . . . . . . . . . . . . . . 58
3.3.4 Exemple de calcul de tonnage . . . . . . . . . . . . . . . . . . . . . . . . . 59
3.3.5 Exemple de calcul avec G-2 Marmora . . . . . . . . . . . . . . . . . . . . 59
4 Signature des structures géologiques en gravimétrie 67

4.1 Levés régionaux et tectoniques . . . . . . . . . . . . . . . . . . . . . . . . . . . . 67
4.1.1 Étude structurale à grande échelle . . . . . . . . . . . . . . . . . . . . . . 67
4.1.2 Études régionales . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 70
Table des matières iii
4.2 Pétrole . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 71
4.2.1 Formations récifales . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 72
4.2.2 Dômes de sel . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 74
4.2.3 Anticlinaux . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 78
4.3 Vallée alluvionnaire . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 83
4.4 Batholite gravimétrique . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 85
4.5 Gisements métallifères . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 86
4.6 Archéologie, travaux publics . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 88
4.7 Autres Exemples . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 90
4.7.1 Prolongement et filtrage ; modélisation . . . . . . . . . . . . . . . . . . . 90
4.7.2 Modélisation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 93
4.7.3 Applications : Étude de cas . . . . . . . . . . . . . . . . . . . . . . . . . . 98
5 Références 99
A Correction de latitude 101
B Obtenir la latitude géocentrique par rapport à la latitude géographique 103
C Système de coordonnées et Théorèmes fondamentaux 105

C.1 Coordonnées cartésiennes . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 105
C.2 Coordonnées cylindriques . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 105
C.3 Coordonnées sphériques . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 106
C.4 Théorèmes fondamentaux . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 107
D Exercices 109
D.1 Correction de données gravimétriques . . . . . . . . . . . . . . . . . . . . . . . . 109
D.2 Corrections gravimétriques . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 111
D.3 Calcul modélisation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 111
D.3.1 Collecteur d’égouts . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 111
D.3.2 Produits toxiques . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 111
iv Table des matières
1 Théorie
1.1 Introduction
En géophysique appliquée, la gravimétrie sert à détecter les contrastes de densité du sous-

sol. Pour ce faire, on mesure en plusieurs points de l’espace les variations de l’accélération
gravitationnelle. Examinons d’abord la nature du champ gravitationnel terrestre.
1.2 Principes de base
1.2.1 Lois de l’attraction universelle
1ère loi de Newton
Deux particules de masse m1 et m2 séparées par une distance r sont attirées l’une vers
l’autre par une force F telle que
Gm1 m2
F=− r (1.1)
r12
où F est la force appliquée sur la masse m2 , r le vecteur unitaire (voir figure 1.1), r1 la distance
entre les masses m1 et m2 , et G, la constante universelle de la gravité. Les termes r1 et G sont
données par
q
r1 = ( x2 − x1 )2 + ( y2 − y1 )2 + ( z2 − z1 )2
G = 6.6742 × 10−8 dyne cm2 /g2 [CGS]
−11 2 2
= 6.6742 × 10 N m /kg [SI] (1.2)
2ème Loi de Newton
Il faut appliquer une force F à une masse m pour lui faire subir une accélération a. Ceci se
traduit par la relation
F = ma. (1.3)
1
2 1. Théorie
z
m2
r1
x
r
m1
y
F IG . 1.1: Deux masses dans un référentiel cartésien.
En utilisant les équations (1.1) et (1.3), on trouve que l’accélération d’une masse m à la surface
du sol s’exprime par
GMT
a=− 2 r=g (1.4)
RT
où MT est la masse de la terre (5.9736 × 1024 kg) et R T le rayon moyen de la terre (6370 km). g
est dite accélération de la gravité, ou simplement gravité, et vaut en moyenne 9.797 m/s2 .
En l’honneur de Galilée, on a nommé l’unité d’accélération gravitationnelle le gal avec :
1 gal = 1 cm/s2 = 10−2 m/s2

1 mgal = 10−3 gal = 10−5 m/s2
La précision d’un gravimètre d’exploration est de l’ordre de 0.01 mgal (10−7 m/s2 ). Les
gravimètres pour les études géodynamiques ou géotechniques sont sensibles au µgal, soit
10−8 m/s2 , environ le milliardième de g.
1.2.2 Champ gravitationnel
Soit une particule immobile en un point A de l’espace. Toutes les particules se trouvant
autour de la masse m du point A subissent une accélération (voir figure 1.2). Chaque point
de l’espace est alors caractérisé par un vecteur accélération qui point vers A et qui est pro-
portionnel à l’inverse de la distance au carré. L’ensemble de ces vecteurs constitue le Champ
Gravitationnel de la masse m.
1.2.3 Potentiel gravitationnel
Le champ gravitationnel est un champ conservatif, c’est-à-dire que le travail fournit pour
déplacer une masse dans ce champ est indépendant du chemin parcouru. Il n’est fonction que
des points de départ et d’arrivée. Donc, si on revient au point de départ, le bilan énergétique
est nul.
1.2 Principes de base 3
a
a
A a
m
F IG . 1.2:
La force qui engendre un champ conservatif peut être obtenue par le gradient de la fonction
du potentiel scalaire U par
F
∇U = g = , (1.5)
m2
où l’opérateur ∇ est donné par :
∂U ∂U ∂U
∇U = i+ j+ k. (1.6)
∂x ∂y ∂z
L’équation du potentiel nous donne donc (m2 : masse unité) :

" #
1 R
Z R Z R
dr Gm1
U= g · r dr = − Gm1 2
= Gm 1 = . (1.7)
∞ ∞ r r ∞ R

Ainsi, on peut définir le Potentiel Gravitationnel en un point du champ gravitationnel comme

le travail requis pour déplacer une masse unitaire de l’infini jusqu’à ce point (figure 1.3).
S’il y a plusieurs masses, on somme la contribution de chacune des masses
N N
∑ Ui = G ∑
mi
U= . (1.8)
i =1 i =1
Ri
Si l’on a une distribution continue de masse dans un volume V extérieur au point (voir fi-
gure 1.4), le potentiel U au point P est
Z
ρ
U=G dv, (1.9)
V r
où ρ est la densité (g/cm3 ) et dv l’élément de volume (cm3 ).

Lors de l’interprétation des données, il est parfois plus simple de travailler avec le potentiel
U, pour ensuite obtenir la gravité g par la relation (1.5).
4 1. Théorie
ds dr
O r
m2
R P
m1
F IG . 1.3: Potentiel gravitationnel.
z P(x,y,z)
V
dv
S
x
y
F IG . 1.4:
1.3 Une référence pour la terre 5
1.2.4 Équations du potentiel
Partant du théorème de Gauss, nous avons

Z Z Z
g · ds = gn ds = ∇ · g dv, (1.10)
S S V
où S est la surface qui enveloppe le volume V, et gn est l’accélération normale à S, positive
vers l’extérieur. Si le volume V est vide de masse, ∇ · g = 0 et, en vertu de l’équation (1.5),
nous avons
∇2 U = 0, (1.11)
qui est l’équation de Laplace. Si le volume contient une masse m, et que la surface S est une
sphère centrée sur m, alors

Gm
Z
4πr2

gn ds = − 2
S r
= −4πGm.
Le signe négatif signifie que l’attraction est dirigée vers le centre de la sphère. Si la surface S
contient plusieurs masses et que leur total vaut M, alors on peut montrer que
Z Z
∇ · g dv = gn ds = −4πGM.
V S
Pour un élément de volume infinitésimal, l’intégrale tombe et nous avons
∇ · g = −4πGρ,
où ρ est la densité de cet élément de volume. Toujours en vertu de l’équation (1.5), nous trou-
vons
∇2 U = 4πGρ, (1.12)
qui est l’équation de Poisson.
Des équations (1.11) et (1.12), il découle que le potentiel sur une surface donnée peut être
produit par des distributions de masses différentes. Ce fait peut être la source d’ambiguı̈tés
lors de l’interprétation des données.
1.3 Une référence pour la terre
1.3.1 Un ellipsoı̈de de révolution : le sphéroı̈de
Pour prédire le champ gravitationnel de la terre en tout point, sa forme et ses variations
de densité doivent être connus. À cause de sa rotation, la terre n’est pas sphérique. Sa forme
peut être approximée par une ellipsoı̈de de révolution quelques fois appelée sphéroı̈de et ca-
ractérisé par son coefficient d’aplatissement
Req − R po 1
f = = , (1.13)
Req 298.257
6 1. Théorie
Rpo
R
Req
Sphère Ellipsoïde
F IG . 1.5:
os
φ
2r c
ω
r ω2r
φ
Req
F IG . 1.6: Paramètres liés au calcul de l’accélération centrifuge : ω est la vitesse angulaire ; ω 2 r

l’accélération centrifuge ; ω 2 r cos ϕ la composante dans la direction de g
où Req est le rayon de la terre à l’équateur (6378.137 km) et R po le rayon de la terre au pôle
(voir figure 1.5).
Sur l’ellipsoı̈de, la gravité de référence go pour un point de latitude géodésique φ est (for-
mule WGS-84) :
1 + 0.00193185265241 sin2 φ
gth (φ) = 9.7803253359 q (1.14)
1 − 0.00669437999014 sin2 φ
La valeur de la gravité ainsi obtenue est celle qui serait observée au niveau de la mer sur une
terre de forme sphéroı̈dale (approximant de près sa forme réelle) et dont la densité ne varie
qu’en profondeur et non pas latéralement.
La différence de 5186 mgals entre la valeur aux pôles et à l’équateur est causée par :
1. L’effet de la rotation de la terre : plus on approche du pôle, plus l’accélération centrifuge

est faible, donc g est maximum (voir figure 1.6).
2. La différence entre le rayon équatorial et le rayon polaire, i.e. le rayon étant plus petit
au pôle, la gravité y est plus forte. Cet effet est cependant atténué par la distribution de
masse plus importante à l’équateur.
La différence de 5186 mgals se répartie environ 2/3 pour l’accélération centrifuge et 1/3 pour
l’aplatissement.
1.4 Densité des roches 7
Sphéroïde Sphéroïde Continent

Géoïde de référence
de référence
Géoïde
Excès Océan
de masse
F IG . 1.7: Sphéroı̈de et géoı̈de
1.3.2 Le géoı̈de
La formule de gth donnée précédemment suppose (1) que le niveau des océan est lisse et
(2) que la densité ne varie qu’en profondeur. Or, il n’en est rien dans la nature. On sait que
cette surface présente des rehaussements et des dépressions de plusieurs dizaines de mètres
en certains endroits, et que la densité peut varier suivant toutes les directions. Ceci nous amène
alors à définir le concept de géoı̈de que l’on définit par la surface équipotentielle correspondant
à la surface des océans aux repos. Sur les continents, le géoı̈de correspond à la surface définie
par l’eau contenue dans un canal étroit reliant les océans de part et d’autre du continent. Par
définition, le géoı̈de est partout perpendiculaire à la verticale telle qu’indiquée par le fil à plomb. Le
géoı̈de et le sphéroı̈de ne coı̈ncident pas en tout point. Il existe des cartes de la hauteur de
géoı̈de par rapport au sphéroı̈de. Les deux plus grandes variations sont au sud de l’Inde (-
105m) et en Nouvelle-Guinée (+73m).
Jusqu’ici, aucune interprétation reliant les lignes de contour du géoı̈de à la surface du globe
(frontières océan-continents, rides mid-océaniques, etc) ne s’est avérée possible. On a émis
l’hypothèse qu’elle pourrait être expliquée par des hétérogénéités du manteau inférieur.
1.4 Densité des roches
La densité des roches est fonction de la nature des minéraux les composants, et de la po-
rosité. Les tableaux suivants donnent des valeurs de densité pour un grand nombre d’entre
elles.
8 1. Théorie
TAB . 1.1: Densités des roches ignées (g/cm3 )
Type de roche Intervalle Moyenne Type de roche Intervalle Moyenne

Rhyolite vitreuse 2.20-2.28 2.24 Diorite quartzeuse 2.62-2.96 2.79
Obsidienne 2.20-2.40 2.30 Diorite 2.72-2.99 2.85
Vitrophyre 2.36-2.53 2.44 Laves 2.80-3.00 2.90
Rhyolite 2.35-2.70 2.52 Diabase 2.50-3.20 2.91
Dacite 2.35-2.80 2.58 Essexite 2.69-3.14 2.91
Phonolite 2.45-2.71 2.59 Norite 2.70-3.24 2.92
Trachyte 2.42-2.80 2.60 Basalte 2.70-3.30 2.99
Andésite 2.40-2.80 2.61 Gabbro 2.70-3.50 3.03
Néphéline- Syénite 2.53-2.70 2.61 Hornblende- Gabbro 2.98-3.18 3.08
Granite 2.50-2.81 2.64 Péridotite 2.78-3.37 3.15
Granodiorite 2.67-2.79 2.73 Pyroxénite 2.93-3.34 3.17
Porphyre 2.60-2.89 2.74 Ignées acides 2.30-3.11 2.61
Syénite 2.60-2.95 2.77 Ignées basique 2.09-3.17 2.79
Anorthosite 2.64-2.94 2.78
TAB . 1.2: Densités des roches métamorphiques (g/cm3 )
Type de roche Intervalle Moyenne Type de roche Intervalle Moyenne

Quartzite 2.50-2.70 2.90 Serpentine 2.40-3.10 2.78
Schiste 2.39-2.90 2.64 Ardoise 2.70-2.90 2.79
Grauwacke 2.60-2.70 2.65 Gneiss 2.59-3.10 2.80
Granulite 2.52-2.73 2.65 Schiste à chlorite 2.75-2.98 2.87
Phyllite 2.68-2.80 2.74 Amphibolite 2.90-3.04 2.96
Marbre 2.60-2.90 2.75 Éclogite 3.20-3.54 3.37
Ardoise quartzique 2.63-2.91 2.77 Métamorphique 2.40-3.10 2.74
F IG . 1.8: Densités moyennes d’échantillons de surface et de carottes

1.4 Densité des roches 9
TAB . 1.3: Densités des minéraux (g/cm3 )
Minéral Intervalle Moyenne Minéral Intervalle Moyenne

Cuivre - 8.7 Sulphures,Arséniures
Argent - 10.5 -Sphalérite 3.5-4.0 3.75
Or 15.6-16.4 - -Covellite - 3.8
-Malachite 3.9-4.03 4.0
-Charcopyrite 4.1-4.3 4.2
Oxydes,carbonates -Stannite 4.3-4.52 4.4
-Limonite 3.5-4.0 3.78 -Stibnite 4.5-4.6 4.6
-Sidérite 3.7-3.9 3.83 -Pyrrhotine 4.5-4.8 4.65
-Rutile 4.18-4.3 4.25 -Molybdénite 4.4-4.8 4.7
-Manganite 4.2-4.4 4.32 -Marcassite 4.7-4.9 4.85
-Chromite 4.3-4.6 4.36 -Pyrite 4.9-5.2 5.0
-Illmenite 4.3-5.0 4.67 -Bornite 4.9-5.4 5.1
-Pyrolusite 4.7-5.0 4.82 -Millérite 5.3-5.65 5.4
-Magnétite 4.9-5.2 5.12 -Charcocite 5.5-5.8 5.65
-Franklinite 5.0-5.22 5.12 -Cobaltite 5.8-6.3 6.1
-Hématite 4.9-5.3 5.18 -Arsénopyrite 5.9-6.2 6.1
-Cuprite 5.7-6.15 5.92 -Smaltite 6.4-6.6 6.5
-Cassitérite 5.8-7.1 6.92 -Bismuthinite 6.5-6.7 6.57
-Wolframite 7.1-7.5 7.32 -Argentite 7.2-7.36 7.25
-Uraninite 8.0-9.97 9.17 -Niccolite 7.3-7.67 7.5
-Galène 7.4-7.6 7.5
-Cinabre 8.0-8.2 8.1
TAB . 1.4: Densités des minéraux non-métalliques et des minéraux divers (g/cm3 )
Type Intervalle Moyenne Type Intervalle Moyenne

Neige - 0.125 Gypse 2.20-2.60 2.35
Pétrole 0.60-0.90 - Bauxite 2.30-2.55 2.45
Glace 2.88-0.92 - Kaolinite 2.20-2.63 2.53
Eau de mer 1.01-1.05 - Orthoclase 2.50-2.60 -
Tourbe - 1.05 Quartz 2.50-2.70 2.65
Asphalte 1.10-1.20 - Calcite 2.60-2.70 -
Lignite 1.10-1.25 1.19 talc 2.70-2.80 2.71
Houille grasse 1.20-1.50 1.32 Anhydrite 2.90-3.00 2.93
Anthracite 1.34-1.80 1.50 Biotite 2.70-3.20 2.92
Brique - 1.50 Magnésite 2.90-3.12 3.03
Carnallite 1.60-1.70 - Fluorine 3.01-3.25 3.14
Soufre 1.90-2.10 - Épidote 3.25-3.50 -
Craie 1.53-2.60 2.01 Diamant - 3.52
Graphite 1.90-2.30 2.15 Corindon 3.90-4.10 4.0
Sel gemme 2.10-2.60 2.22 Barite 4.30-4.70 4.47
Zircon 4.00-4.90 4.57
10 1. Théorie
TAB . 1.5: Densités des matériaux rocheux typiques (g/cm3 )
Roches Ignées
Roche Nb échan. Intervalle
Granite 155 2.516-2.809
Granodiorite 11 2.668-2.785
Syénite 24 2.630-2.899
Diorite 13 2.721-2.960
Norite 11 2.720-3.020
Gabbro 27 2.850-3.120
Diabase 40 2.804-3.110
Péridotite 3 3.152-3.276
Dunite 1 3.289
Pyroxénite 8 3.10-3.318
Anorthosite 12 2.640-2.920
Valeurs tirées de :
“Handbook of Pysical Constants”,
édité par Francis Birch,
Geological Society of America, Special Paper 36, 1942.
2 Les données gravimétriques
De façon générale en géophysique appliquée, on mesure la composante verticale de l’accélération

gravitationnelle, que l’on note gz . Étant donnée que l’on cherche à mesurer une variation très
faible du champ, la mesure est délicate à réaliser et requiert également une grande précision
lors du positionnement des stations de mesure.
2.1 Corrections et références
Afin d’obtenir les variations du champ gravitationnel dues à des causes géologiques, il est
nécessaire de corriger nos lectures de toutes les autres causes extérieures pouvant les influen-
cer (dérive de l’appareil, marée, ellipticité de la terre, . . .).
2.1.1 Correction de dérive
Par cette correction, on tente d’éliminer l’influence apportée sur les mesures par les marées
(figure 2.1) et la fatigue de l’instrument.
Dans ce but, il est nécessaire de suivre un certain cheminement entre les stations de lec-
tures. Dans la pratique, on fait une série de mesures en suivant un cheminement en boucle : la
série débute habituellement en un point donnée et se termine à ce même point (figure 2.2). Le
point de départ de la boucle est normalement relié à une station de base.
En général, les mesures du début et de la fin à la station de base ne sont pas semblables.
Cette différence, appelée dérive, est due en partie au gravimètre, en partie au marée lunaire.
Les valeurs mesurées sont donc entravées d’erreurs puisqu’une de leurs composantes provient
de la dérive et ne reflète pas un changement dans la valeurs de g dû à des hétérogénéités du
sous-sol.
La correction est faite en supposant que la dérive est linéaire dans le temps. Donc, si on est
passé à la station de base, aux temps t1 et t2 et que les valeurs mesurées étaient respectivement
v1 et v2 , le taux de dérive ∆0d est défini par
v2 − v1
∆0d = . (2.1)
t2 − t1
Lorsque la dérive est positive, c’est que les mesures ont été surestimées, il faut donc les dimi-
nuer. La correction de dérive sera négative. Inversement, dans le cas où la dérive est négative,
11
12 2. Les données gravimétriques
F IG . 2.1: Exemple de l’influence des marées sur la gravité.
F IG . 2.2: Illustration d’une boucle de mesure simple.
les mesures sont sous-estimées et la correction devra être positive. Ainsi, toute valeur v prise
au temps t (où t1 ≤ t ≤ t2 ) est corrigée par la formule suivante :
vcor (t) = vlu (t) − ∆0d (t − t1 ). (2.2)
Exemple
Station Lecture Temps

1 1032.1 12h15
2 12h20
3 12h25
4 12h31
5 12h35
6 12h39
1 1031.0 13h05
Le taux de dérive est
1031.0 − 1032.1 1.1

∆0d = =− = −0.022 div./minutes.
13h05 − 12h15 50
2.1 Corrections et références 13
Lecture
5 6 v2
7
Différence
4
observée
3 8 m2
Dérive
réelle
2
Diff.
9 v1
1 10
m1
Stn1 Stn2
Base Base
Temps
t1 t2
F IG . 2.3: Exemple d’une boucle avec plusieurs références.
Donc, pour la lecture de la station 4, prise 16 minutes après la 1re lecture de la station 1, la
correction est de :
16 × (−0.022) = −0.352 div. ' −0.4 div.
Le principe demeure le même si au lieu de boucler sur la station de départ, la dernière

mesure se fait sur une autre station de base (figure 2.3). Évidemment, les stations initiales et
finales doivent auparavant avoir été reliées entre elles.
Si lors de l’établissement des deux stations de bases on a trouvé des valeurs égales à ce
que la différence entre les nouvelles valeurs observées v1 et v2 soit semblable à celle qui existe
entre m1 et m1 , obtenus antérieurement à v1 et v2 . La dérive est égale à
dérive = différence réelle − différence observée

= ( m2 − m1 ) − ( v2 − v1 ). (2.3)
De la même manière qu’auparavant, la formule de correction est
( m2 − m1 ) − ( v2 − v1 )
vcor (t) = vlu (t) + × ( t − t1 ) + ( m1 − v1 ). (2.4)
t2 − t1
Notons le terme supplémentaire à la fin de l’équation, m1 − v1 . Il a pour but de ramener les

valeurs à un niveau de référence semblable pour chaque partie du levé.
Exemple
Soit la ligne de base suivante :
Station lecture
BL1 1030.1
BL2 1032.0
BL3 1031.7
BL4 1032.4
et les mesures suivantes prises le lendemain :
Station Temps Lecture

BL1 8h50 1027.9
BL2 8h53
BL3 8h56
BL4 9h00
... ...
BL10 10h00
BL2 10h30 1028.7
La dérive est de 1.1 division (1032.0-1030.1-1028.7+1027.9), et le taux de correction est de

0.011 div./minute (1.1/1h40). La correction de niveau vaut 2.2 div. (1030.1-1027.9). Ainsi, pour
la station BL1 (2ème journée) :
vcor = 1027.9 + 0.011 × 0 + 2.2 = 1030.1 ≡ m1
et pour la station BL2 (2ème journée) :
vcor = 1028.7 + 0.011 × 100 + 2.2 = 1032.0 ≡ m2
2.1.2 Correction de latitude
Cette correction tient compte des variations de g avec la latitude dues à la rotation de la
terre et à son aplatissement. À partir des mesures géodésiques mondiales, on sait que la terre
est un ellipsoı̈de de révolution presque parfait. Sur cette surface, le champ gravitationnel peut
être décrit par l’équation suivante (WGS-84) :
1 + 0.00193185265241 sin2 φ
gth (φ) = 9.7803253359 q (2.5)
1 − 0.00669437999014 sin2 φ
où gth (φ) est la valeur du champ au point de latitude géodésique φ. La correction ∆ L pour un
déplacement l suivant un méridien est donc
∂gth
∆L = · l, (2.6)
∂l
avec
dl = R(φ)dφ ≈ Re dφ, (2.7)
où Re est le rayon équatorial de la terre (6378.136 km). Finalement, si l est en mètres, on peut
simplifier l’expression à
∆ L ≈ 8.1669 × 10−4 · l · sin 2φ [mgal] : (N → S) (2.8)
Notez que la latitude géodésique se distingue de la latitude géocentrique. La seconde est

l’angle entre la direction du centre de la terre vers le point considéré en surface et le plan
N
-0.030 mgal
25 m
-0.015 mgal
0 (48°44'N)
0.015 mgal
0.030 mgal
F IG . 2.4:
équatorial. La première, qui est la plus courantes sur les cartes, est l’angle entre la normale à
l’ellipsoı̈de de référence et le plan équatorial. On montre en annexe B la relation entre ces deux
latitudes.
L’équation est linéaire (i.e. φ = cte) sur une distance d’environ 1.5 km. Puisque gth est
plus élevée aux pôles qu’à l’équateur, il faut additionner ∆ L (correction positive) pour un
déplacement vers l’équateur. Notez que pour obtenir une précision acceptable, on doit cher-
cher à positionner les différentes stations avec une précision d’une dizaine de mètres (par
exemple à partir de photos-aériennes). Pour une précision recherchée de 0.01 mgal, il faut
connaı̂tre à environ ± 10 m la distance entre deux stations séparées de 100 m si φ =45˚. Il est
à noter que les corrections sont positives lorsque les stations se localisent au sud de la ligne
de référence et négative pour celles se situant au nord. Aucune correction de latitude n’est
apportée pour un cheminement est-ouest.
Dans un levé local, les corrections ne sont pas calculées pour chacune des stations à par-
tir de la formule générale, mais sont plutôt déterminées à partir d’une grille proprement
graduée. Par exemple, supposons qu’un levé de gravimétrie a été effectué autour de la latitude
géodésique 48˚44’N. L’échelle des cartes de travail est de 1 : 2000 (20m/cm) et nos stations de
mesure sont espacées de 25 m. On a φ =48˚44’, ce qui donne en décimales 48.7333˚. On trouve
alors la correction de latitude correspondante, soit :
∆ L = 8.1669 × 10−4 sin(2φ)l = 6.14 × 10−4 l mgal (N → S)
Ainsi, chaque déplacement de 1.25 cm du nord vers le sud (N→S) sur la carte au 1 : 2000
entraı̂nera une correction de 0.015 mgal (0.000614 × 25). La grille peut donc être graduée en
multiples de 0.015, la correction zéro étant affectée aux stations se trouvant à la latitude 48˚44’N
(voir figure 2.4).
2.1.3 Correction d’altitude
Afin de pouvoir interpréter les données en fonction du sous-sol, celles-çi doivent être rat-
tachées à une référence équipotentielle unique pour le levé. Or, les lectures d’un levé gra-
vimétrique ne sont pas forcément prises au-dessus d’un terrain plat. Ainsi, plus on se rap-
proche du niveau de référence, plus g augmente car on se rapproche du centre de la terre. Les
mesures obtenues présentent donc des variations qui ne sont dues qu’à la position de la station
de mesure et non pas à des hétérogénéités du sous-sol. Il faut donc corriger les mesures.
On sait que
Gm
gr = , (2.9)
r2
où r est le rayon de la terre au niveau de référence. Si on se déplace d’une hauteur h par rapport
à ce niveau de référence, alors
Gm Gm
gh = 2
= 2 . (2.10)
(r + h ) r (1 + 2(h/r ) + (h/r )2 )
Puisque l’on a r h, alors
Gm(1 − 2h/r )
gh = = gr − 2hgr /r (2.11)
r2
et donc
gh − gr = −2hgr /r. (2.12)
En prenant r comme rayon moyen de la terre, la correction d’altitude (aussi nommée cor-
rection d’air libre) est donnée par (h en mètres, positif vers le haut)
∆h = 0.3086 h mgal ; (h > 0) (2.13)
Donc ∆h est positif si on est au-dessus du référentiel et négatif si on est en-dessous. Pour une
précision d’environ 0.01 mgal, il faut connaı̂tre à ± 3 cm la hauteur de la station par rapport
au référentiel.
2.1.4 Correction de plateau
La correction de plateau vise à tenir compte de la masse comprise entre le référentiel et

la station de mesure. On considère que cette masse peut être approximée par une tranche
horizontale homogène et d’extension infinie. Dans bien des cas, cette approximation n’est pas
juste, et il faudra apporter une correction supplémentaire, dite correction de terrain.
Pour une tranche de hauteur h, l’attraction est donnée par
∆ p = 2πGρ B h (2.14)
où G = constante universelle de la gravitation et ρ B est la densité présumée de la masse de

la tranche (ρ B = 2.67 g/cm3 en moyenne). Comme ∆ p augmente lorsque h augmente, il faut
soustraire ∆ p lorsque h > 0 et donc :
∆ p = −0.04191 ρ B h mgal ; (h > 0), (2.15)
où ρ B est en g/cm3 et h en mètres. Comme pour la correction d’altitude, il faut connaı̂tre
précisément l’élévation du gravimètre à chaque station (h = ±10 cm si la précision recherchée
est de ± 0.01 mgal).
z station P en trop
manquant h tranche de Bouguer
référence
F IG . 2.5: Tranche de Bouguer et correction de relief.
Si la mesure gravimétrique est faite sous terre (dans une mine ou un tunnel), la tranche
fictive ajoutée entre la station de référence (à z1 ) et la station souterraine (à z2 , plus bas que z1 )
exerce une attraction vers le bas à z1 et vers le haut à z2 . La correction est doublée, et devient
4πGρ(z2 − z1 ).
De façon courante, on combine la correction d’altitude et la correction de plateau pour
obtenir ce que l’on appelle alors la correction de Bouguer (attention, ceci n’est pas l’anomalie de
Bouguer) :
∆hB = (0.3086 − 0.04191 ρ B )h mgal ; (h > 0). (2.16)
Si l’on choisit ρ B = 2.67 g/cm3 , on obtient
∆hB = 0.1966h mgal ; (h > 0). (2.17)
2.1.5 Correction de relief
Lors de la correction de Bouguer, on enlève l’attraction d’une tranche de terrain d’épaisseur

h. Si les variations topographiques sont telles qu’on ne peut approximer le terrain par une
tranche uniforme, il faut intégrer numériquement d’une part les parties qui dépassent et d’autre
part les parties qui manquent à la tranche de Bouguer.
Au point P (voir figure 2.5), la gravité verticale causée par un morceau de volume V, en
trop ou manquant, est donnée par
zo
Z
∆gt = G∆ρ dv. (2.18)
V ( xo2 + y2o + z2o )3/2
Pour les morceaux en trop, zo et ∆ρ sont positifs et pour les morceaux manquant, zo et ∆ρ
sont négatifs. Ainsi, la correction de relief est toujours positive. Par exemple, considérons un
volume de roche au dessus de la station. Ce volume génère une attraction vers le haut, de sens
contraire à l’attraction terrestre. Il a donc pour effet de réduire g et il faut additionner le terme
de correction pour éliminer son effet.
En général, l’intégration se fait numériquement au moyen d’un ordinateur, en utilisant des
cartes topographiques numérisées. Pour simplifier les calculs, on peut discrétiser le terrain en
cylindres creux concentriques, ces mêmes cylindres eux-mêmes scindés en morceaux dont les
sommets sont ajustés à la topographie moyenne (voir figure 2.6).
Topo
Station
hi
h
Référence
r1
r2
F IG . 2.6: Définition du cylindre creux utilisé pour accomplir la correction de relief.
L’expression donnant l’attraction gravitationnelle g, sur l’axe d’un morceau de cylindre

creux et d’épaisseur r2 − r1 est la suivante
q q
2πGρ 2 2
∆ti = 2 2
r2 − r1 + r1 − h i − r2 − h i (2.19)
N
où ∆ti est l’attraction du ie secteurs du cylindre, hi la hauteur du morceau de cylindre par
rapport à la station, ρ la densité du cylindre et N le nombre de secteurs par lequel le cylindre
a été divisé. La correction totale pour le cylindre entier est :
N
∆T = ∑ ∆ti . (2.20)
Par ailleurs, au lieu d’utiliser la formule donnée précédemment, il est possible d’utiliser un
réticule (voir figure 2.7) que l’on supperpose aux cartes topographiques et des tables préparées
par Hammer (et complétées par Bible). Ces tableaux (2.1 et 2.2) nous donnent, pour différentes
valeurs de h, les corrections en mgal qu’il nous faut apporter pour chacun des secteurs du
réticule.
Les différentes étapes à exécuter pour effectuer la correction sont :
1. Localiser la station à corriger sur une carte topographique assez précise.
2. Centrer le réticule sur ce point.
3. Déterminer l’élévation moyenne de chacun des secteurs de l’abaque.
4. Prendre la différence positive ou négative entre l’élévation de la station considérée et
chacun des secteurs.
5. Pour chacun des secteurs, à l’aide de la table de Hammer, calculer la correction à apporter
(la correction est toujours positive).
TAB . 2.1: Tables de Hammer pour la correction de relief ; densité : ρ = 2 g/cm3 ; Zones B à H ;
Correction : ∑ ti ∗ 0.01 mgal
Zone B C D E F G H
Secteurs 4 6 6 8 8 12 12
Rayon 2 à 16.6 à 53.3 à 170.1 à 390 à 895 à 1529 à
16.5 m 53.3 m 170.1 m 390 m 895 m 1529 m 2615 m
ti ±h en m ±h en m ±h en m ±h en m ±h en m ±h en m ±h en m
0 0.0 à 0.3 0.0 à 1.3 0.0 à 2.4 0.0 à 5.5 0.0 à 8.2 0.0 à 17.6 0.0 à 22.9
0.1 0.3 à 0.6 1.3 à 2.3 2.4 à 4.1 5.5 à 9.1 8.2 à 14.0 17.6 à 30.5 22.9 à 40.0
0.2 0.6 à 0.8 2.3 à 3.0 4.1 à 5.3 9.1 à 12.0 14.0 à 18.3 30.5 à 39.3 40.0 à 51.5
0.3 0.8 à 0.9 3.0 à 3.5 5.3 à 6.3 12.0 à 14.3 18.3 à 21.6 39.3 à 46.6 51.5 à 61.0
0.4 0.9 à 1.0 3.5 à 4.0 6.3 à 7.1 14.3 à 16.2 21.6 à 24.4 46.6 à 52.7 61.0 à 68.9
0.5 1.0 à 1.1 4.0 à 4.4 7.1 à 7.8 16.2 à 17.7 24.4 à 27.0 52.7 à 58.0 68.9 à 76.0
1 1.1 à 2.0 4.4 à 7.3 7.8 à 13.1 17.7 à 29.6 27.0 à 45.0 58.0 à 97.0 76.0 à 126
2 2.0 à 2.7 7.3 à 9.8 13.1 à 17.0 29.6 à 38.3 45.0 à 58.0 97.0 à 125. 126. à 163
3 2.7 à 3.5 9.8 à 11.9 17.0 à 20.2 38.3 à 45.4 58.0 à 68.0 125 à 148 163 à 193
4 3.5 à 4.2 11.9 à 13.8 20.2 à 23.1 45.4 à 51.8 68.0 à 77.0 148 à 158 193 à 219
5 4.2 à 5.0 13.8 à 15.6 23.1 à 25.7 51.8 à 57.6 77.0 à 86.0 158 à 186 219 à 242
6 5.0 à 5.7 15.6 à 17.4 25.7 à 28.1 57.6 à 62.9 86.0 à 94.0 186 à 202 242 à 263
7 5.7 à 6.5 17.4 à 19.1 28.1 à 30.4 62.9 à 67.8 94.0 à 101 202 à 213 263 à 283
8 6.5 à 7.3 19.1 à 20.8 30.4 à 32.6 67.8 à 72.4 101 à 108 213 à 233 283 à 302
9 7.3 à 8.2 20.8 à 22.6 32.6 à 34.7 72.4 à 76.8 108 à 114 233 à 247 302 à 320
10 8.2 à 9.1 22.6 à 24.4 34.7 à 36.7 76.8 à 81.1 114 à 120 247 à 260 320 à 337
11 ... 24.4 à 26.1 36.7 à 38.7 81.1 à 85.3 120 à 126 260 à 272 337 à 353
12 ... 26.1 à 27.9 38.7 à 40.6 85.3 à 89.3 126 à 131 272 à 284 353 à 368
13 ... 27.9 à 29.7 40.6 à 42.6 89.3 à 93.2 131 à 137 284 à 296 368 à 383
14 ... 29.7 à 31.6 42.6 à 44.5 93.2 à 97.0 137 à 142 296 à 308 383 à 397
15 ... 31.6 à 33.5 44.5 à 46.4 97.0 à 100.8 142 à 147 308 à 319 397 à 411
F IG . 2.7: Exemple de réticule de Hammer.
TAB . 2.2: Tables de Hammer (suite) ; Zones I à M.

Zones I J K L M
Secteur 12 16 16 16 16
Rayon 2615 à 4470 à 6650 à 9900 à 14750 à
4470 m 6650 m 9900 m 14750 m 21950 m
ti ±h en m ±h en m ±h en m ±h en m ±h en m
0 0 à 30.2 0 à 51 0 à 62 0 à 76 0 à 93
0.1 30.2 à 52.1 51 à 88 62 à 108 76 à 131 93 à 160

0.2 52.1 à 67.1 88 à 114 108 à 139 131 à 170 160 à 207
0.3 67.1 à 79.6 114 à 135 139 à 165 170 à 201 207 à 245
0.4 79.6 à 90.2 135 à 153 165 à 187 201 à 228 245 à 278
0.5 90.2 à 100 153 à 169 187 à 206 228 à 252 278 à 307
1 100 à 164 169 à 280 206 à 341 252 à 416 307 à 507
2 164 à 213 280 à 361 341 à 441 416 à 537 507 à 655
3 213 à 253 361 à 427 441 à 521 537 à 636 655 à 776
4 253 à 287 427 à 485 521 à 591 636 à 721 776 à 880
5 287 à 317 485 à 537 591 à 654 721 à 797 880 à 978
6 317 à 344 537 à 584 654 à 711 797 à 867 973 à 1058
7 344 à 369 584 à 628 711 à 764 867 à 932 1058 à 1136
8 369 à 393 628 à 669 764 à 814 932 à 993 1136 à 1210
9 393 à 416 669 à 708 814 à 861 993 à 1050 1210 à 1280
10 416 à 438 708 à 745 861 à 906 1050 à 1104 1280 à 1346
11 438 à 459 745 à 780 ... ... ...

12 459 à 479 780 à 813 ... ... ...
13 479 à 498 813 à 845 ... ... ...
14 498 à 516 845 à 877 ... ... ...
15 516 à 534 877 à 908 ... ... ...
6. Additionner la contribution de chacun des secteurs et convertir à la densité moyenne

utilisée pour le levé (car les tables ont été calculées en supposant une densité de 2 g/cm3 ).
Si l’on doit adopter une autre valeur ρ’, la correction évaluée d’après les tables doit être
multipliée par ρ’/2.
Exemple de correction de relief
Un réticule de Hammer a été superposé à la carte topographique (figure 2.8) et l’altitude

moyenne de chacun des secteurs a été notée dans le réticule. La table de Hammer a alors
été utilisée afin de calculer les corrections de terrain en centième de milligal. En additionnant
les corrections de chacun des secteurs (voir tableau 2.8), on obtient une correction totale de
0.032 mgal pour les zones B et C.
Il est évident que l’on ne peut songer à utiliser la même carte topographique pour évaluer
les altitudes moyennes des zones B et M ; à l’échelle 1 : 1000, une distance de 2 mètres est
représentée par 2 mm alors qu’une distance de 22 km conduirait à un abaque de 22 m de
rayon ; au 1 : 100, 000, ces mêmes longueurs deviennent respectivement 0.02 mm et 22 cm. La
première échelle convient pour évaluer l’influence des zones centrales, alors que la seconde
n’est utilisable que pour les zones éloignées. Il sera donc nécessaire, dans tous les cas où l’on
calculera complètement la correction de relief, de la scinder en plusieurs parties.
Choix du niveau de référence
Jusqu’à présent, nous avons admis que la surface de référence est le géoı̈de (cote 0). En
réalité, pour réduire le caractère imparfait de la correction de Bouguer, on sera souvent amené
à donner à la surface de référence la cote zo de la station la plus basse de l’étude. Dès lors,
l’épaisseur des tranches de terrain intéressées par la correction de Bouguer sera réduite de la
quantité zo , et le rayon à donner à la zone extérieure pour effectuer les corrections de terrain
se trouvera également diminué. Ce procédé est d’ailleurs plus justifié, pour une prospection
d’étendue limitée, que celui consistant à adopter systématiquement le géoı̈de comme surface
de référence.
Considérations pratiques
Nous avons vu que l’influence des défauts et des excès de masse était de même sens ; il
convient, en conséquence, de traiter de manière particulière les compartiments du réticule
présentant à la fois des altitudes plus élevées et plus faibles que celle de la station. Le plus
souvent, il sera possible d’éviter une telle disposition par une rotation appropriée du réticule
(bien entendu, ce dernier doit conserver le même azimut pour tout le calcul relatif à une même
zone), conduisant par exemple à placer la courbe de niveau de cote Z le long d’un rayon vec-
teur. Si la chose est impossible, on devra (pour chaque compartiment présentant cette par-
ticularité) déterminer séparément l’altitude moyenne des terrains de cote supérieure à z et
celle des terrains de cote inférieure, puis calculer la moyenne des deux valeurs absolues des
différences avec z, pondérée en fonction des surfaces occupées dans chaque cas.
Remarques : Le calcul de la correction d’une station est assez long et peut prendre entre
Zone Secteurs Total

B 0.3 / 1.0 / 0.2 / 0.3 1.8
C 0.1 / 0.5 / 0.3 / 0.1 / 0.3 / 0.1 1.4
F IG . 2.8: Exemple d’utilisation du réticule de Hammer

Zones D E F G H I J K
h en mètres ± 13 ± 30 ± 45 ± 97 ± 126 ± 164 ± 280 ± 341
TAB . 2.3: Élévations d’influence négligeable, par zone du réticule de Hammer
30 minutes et une heure dépendant de la topographie. Ce temps peut être réduit à 15 minutes
grâce aux remarques qui suivent :
1. De faibles irrégularités dans la topographie n’ont que peu d’effet sur la gravité et peuvent
être négligées. (Moins de 1/20 de la distance à la station, par exemple moins de 0.1 m
d’élévation à l’intérieur d’un rayon de 2 m autour de la station).
2. Afin de simplifier les corrections de terrain, l’emplacement des différentes stations doit
être judicieux : éviter les ravins, les pieds de montagne, etc, car ce sont les dénivellations
les plus près de la station qui ont le plus d’influence.
Le tableau 2.3 résume, pour les différentes zones, les élévations de la topographie qui en-
traı̂nent une correction inférieure à 0.01 mgal lorsque la densité du terrain vaut 2 g/cm3 :
2.1.6 Méthode de Nettleton
Pour la correction de Bouguer, il est important d’essayer de définir la densité ρ B avec une
précision de 0.1 g/cm3 . En ce qui concerne la correction de relief, le rôle joué par ρ B est moins
considérable. En effet, on observe rarement une différence d’un milligal entre stations voi-
sines ; une variation de 0.2 sur la densité n’introduira dans ces conditions qu’un écart de l’ordre
de 0.10 mgal sur une valeur qui n’est souvent définie qu’avec une moins bonne approximation.
La densité ρ B peut être déterminée à partir d’échantillons prélevés sur le terrain et repré-
sentatifs de la géologie, pour lesquels une mesure est faite en laboratoire. Elle peut également
être déterminée à partir des données gravimétriques même. Une fois la correction d’altitude
appliquée, l’anomalie offre une corrélation très forte avec la topographie du terrain. La méthode
de Nettleton consiste donc à représenter, sur une même figure, un profil topographique et les
profils de l’anomalie de Bouguer qui lui correspondent, calculés pour plusieurs densités (voir
figure 2.9 pour un exemple). On choisit une région au relief assez accidenté pour que le rôle
de la correction d’altitude soit déterminant vis-à-vis de la forme des profils de gravité. Parmi
ceux-ci, une partie reflétera assez fidèlement les irrégularités topographiques (groupe des den-
sités de 2.4 à 2.8) et une autre donnera une image inversée du relief (groupe des densités de
1.8 à 2.2).
Dans les deux cas, il existera une certaine corrélation entre les formes topographiques et
gravimétriques ; le profil le plus satisfaisant sera évidemment celui qui se situera entre les
deux groupes précédemment définis, et pour lequel, visiblement, la corrélation entre le relief
et l’anomalie de Bouguer sera la moins nette. Un certain nombre de profils de ce genre, ju-
dicieusement distribués sur toute l’étendue prospectée, conduiront à différentes valeurs de
densité. Si ces valeurs ne sont pas trop dispersées, leur moyenne pourra être retenue comme
caractérisant au mieux la région, et les corrections d’altitude et de relief seront calculées en
utilisant cette valeur.
L’emploi de cette méthode repose sur une hypothèse rarement rencontrée dans la nature,
en tout cas invérifiable : définir comme le plus probable le profil gravimétrique reflétant le
1.8
1.9
2.0
2.1
2.8 2.2
2.3
2.4
1.8 2.5
2.6
2.7
2.8
F IG . 2.9: Exemple d’application de la méthode de Nettleton

moins l’allure de la topographie, c’est admettre implicitement qu’il n’y a pas à l’aplomb de ce
profil d’anomalies d’origine géologique dont la forme soit en relation avec le relief.
Il y a donc lieu de préciser dans quel cas cette méthode offre le plus de chances de conduire
à un résultat convenable. Ce sera évidemment lorsque les conditions géologiques permettent
d’affirmer qu’il ne peut exister aucune relation entre la forme du relief et celle des couches pro-
fondes. Un très bon exemple de réalisation d’une telle condition est fourni par les dunes de
sable dont le relief, d’ailleurs variable, est déterminé avant tout par l’orientation et le régime
des vents. Ou encore, si le relief est assez accidenté pour que l’on observe un creux relatif au
voisinage du sommet d’une colline ou une ondulation au centre d’une colline ou une ondula-
tion au centre d’une vallée, il est assez probable que de tels mouvements secondaires, résultant
en général de l’action de l’érosion, n’ont pas une origine géologique profonde.
Le profil gravimétrique le plus probable sera donc celui qui “effacera” le mieux de tels
mouvements, même si son dessin doit offrir une certaine parenté avec la forme générale du
relief.
2.1.7 Anomalie Bouguer
L’anomalie de Bouguer constitue la réponse gravimétrique causée par les hétérogénéités

de densité du sous-sol. L’interprétation se fait donc à partir de celle-ci. Elle est donnée par
∆gB = ∆gobservée ± les 5 corrections
où les corrections sont :

1. Correction de dérive de l’appareil et des marées ;
2. Correction de latitude ∆ L = 8.1669 × 10−4 sin 2φ mgal ;
3. Correction d’altitude ∆h = 0.3086 h mgal
4. Correction de plateau ∆ B = −0.04191 ρ B h mgal
5. Correction de terrain ∆ T
et où h est en mètres et positif si la station est au-dessus du référentiel et négatif dessous, et
∆gobservée = gobservée − gre f .
Remarque :
Le gravimètre ne donne pas une valeur absolue de g, mais bien une valeur relative.
∆gB = gobservée ± ∆g − gre f (2.21)
Si pour un niveau z de référence, gobservée en x, y et z a été rendue comparable à gre f calculée

sur l’ellipsoı̈de de référence en ( x, y, z = 0) grâce aux corrections. On en conclut donc que
gobservée a été réduite au niveau de l’ellipsoı̈de.
C’est inexact.
ligne 15N
Stn : 10N-05W
ligne 10N
ligne 05N
Stn : 00-00
ligne 00
ligne 05S
ligne 10S
Stn : 15S-10E
ligne 15S
F IG . 2.10: Exemple de numérotation des stations
En fait, il faudrait plutôt écrire

∆gB = gobservée − ( gre f ± ∆g) (2.22)
Au point ( x, y, z), on peut faire correspondre un point ( x, y, 0) sur l’ellipsoı̈de où la gravité
normale vaut gre f . En effet, gobservée et gre f ne sont pas comparables car ( x, y, z) et ( x, y, 0) ne
sont pas à la même altitude ni affectés du même relief. On fait donc subir à gre f les corrections
nécessaires pour l’emmener dans la position désirée et nous permettre de disposer d’une va-
leur théorique convenable de g en ( x, y, z). L’anomalie de Bouguer est donc attachée au point
( x, y, z) et non pas au point ( x, y, 0) comme on tend à le croire.
2.2 Levé gravimétrique
2.2.1 Numérotation des stations
Suivant le genre de levé, le numérotation des stations pose plus ou moins de difficultés.
Dans un levé local, le long de traverses bien définies, les points de mesure portent habituel-
lement le numéro de la traverse et un numéro indiquant son ordre sur celle-ci. On aurait par
exemple la station 02 + 09 qui serait la neuvième station sur la traverse deux ou, comme in-
diqué sur la figure 2.10, la station ”10N-05W” qui serait la station située à 5 unité vers l’ouest
par rapport à la station de référence (station 00-00) sur la ligne à 10 unité par rapport à la ligne
de référence (ligne 00).
Dans un levé régional (par exemple pour toute la Gaspésie) les points de mesure seraient
principalement déterminés par leur longitude et leur latitude. Habituellement, on utilise de
plus une numérotation arbitraire pour les stations.
Recommandations :
– La représentativité de la station : il faut à tout prix éviter l’usage d’un même nom pour
deux stations différentes.
2.2 Levé gravimétrique 27
– Le nom de la station doit fournir une idée approximative de la position géographique

du point considéré. Pour les principaux réseaux, il convient de choisir un seul nom ou
numéro pour identifier un réseau en particulier. On aurait, par exemple, 4 pour chemin
Lévesque, etc. Un nom (ou l’initial) d’un lac, d’une rivière, d’une route ou d’une munici-
palité convient parfaitement. Ainsi par la connaissance du début du nom d’une station,
on peut déterminer la zone géographique où il est le plus probable de la trouver.
– La concision : les noms trop longs doivent être évités parce qu’ils sont souvent écrits sous
forme abrégée dans les carnets de terrain, ce qui constitue une source d’erreur poten-
tielle. On devrait au maximum accepter trois groupes différents de lettres ou de chiffres
comportant au maximum deux caractères. Seraient acceptables des noms tels que : S-2,
TH-33-A, TL-6-10, BK-7-3, etc. Les noms : TL-6-24-2, TL-14-13-2, 323-28-2, et autres sem-
blables devraient être rejetés. Le nom des stations pourrait aussi compter un indice de la
qualité du levé au point considéré. Ainsi, les points dont l’élévation est mesurée à l’aide
d’un baromètre pourraient voir leur nom précédé d’un B, on aurait : BK-5-8, BTH-8-2,
etc. Lorsqu’un réseau secondaire se déploie à partir d’un point donné, ce réseau peut
prendre le nom du point de départ. Ainsi, le point BK-5-8 qui a son début au point K-5,
etc.
2.2.2 Nivellement
Le nivellement est extrêmement important en gravimétrie. Comme nous l’avons vu dans la

correction d’altitude, une erreur de quelques cm entraı̂ne une erreur importante dans la valeur
corrigée de la gravité. Si on considère, la correction d’altitude et la correction de plateau avec
une densité de 2.76 g/cm3 , nous avons une correction de 0.1966h. Ainsi, une erreur d’estima-
tion de h de 10 cm correspond à une erreur 0.02 mgal sur g, 50 cm correspond à 0.1 mgal et 1 m
à 0.2 mgal. Donc, une erreur de 50 cm dans la mesure des élévations introduit une erreur dans
la valeur corrigée de la gravité qui est de l’ordre des anomalies recherchées (géotechnique,
archéologie, minier). Heureusement, le nivellement à l’aide d’un niveau ou d’une station to-
tale nous donne une précision meilleure que 5 cm, soit 0.01 mgal. Cette méthode de mesure
requiert toutefois beaucoup de personnel (une personne pour la mire et une pour l’instrument)
et du temps.
Dans la pratique, avant chaque mesure gravimétrique, il faut déterminer la hauteur relative
de la station de mesure par rapport à la station de base. À partir des deux visées obtenues avec
le théodolite, la différence de hauteur entre les deux stations est donnée par la différence entre
les deux lectures (voir figure 2.11). Attention cependant, la lecture sur la mire nous donne la
hauteur relative entre le niveau de la station et la ligne de visée. Ainsi, une ambiguı̈té peut
provenir du fait que plus la valeur lue sur la mire est importante, plus l’altitude absolue du
point de mesure est faible. Par exemple, si on obtient des lectures de 1.53 m et 1.75 m pour les
stations A et B respectivement, alors le niveau de la station B est dessous celui de la station A.
Sur le terrain, il est très rare de pouvoir faire un levé entier avec une seule position pour
la station totale ou le théodolite (terrain accidenté, visée trop lointaine). Il est alors nécessaire
de changer la position de l’instrument pour poursuivre les mesures de nivellement. Dans ce
cas, il ne faut surtout pas oublier de relier entre eux les différents segments. Cette opération
se réalise en prenant pour une même station deux mesures de niveau : une avant et l’autre
après avoir déplacé le théodolite (voir figure 2.11). Il sera alors possible de relier les mesures
théodolite ligne de visée
h h’
S2 surface
S1 h - h'
T2
T1 h[S3/T2]
S3
S1 h[S2/T1] h[S2/T2]
S2
F IG . 2.11:
faites suivant le second segment avec celles du premier. Un exemple est présenté sur la figure
2.11. En considérant la station S1 comme référence pour notre levé, le niveau de la station S2
est tout simplement donné par ∆hS2/S1 = h[S2/T1] − h[S1/T1]. Pour la station S3, il faut faire
l’opération en deux étapes. D’abord, calculer le niveau de la station S3 par rapport à la station
S2 : ∆hS3/S2 = h[S3/T2] − h[S2/T2], puis ramener cette mesure par rapport à la station S1 :
∆hS3/S1 = ∆hS3/S2 + ∆hS2/S1 .
Après avoir complété le levé et avant de ranger les instruments, il faut s’assurer d’avoir fait
un bon levé du niveau relatif des stations. Ceci signifie qu’il faut fermer la boucle du nivellement.
Pour ce faire, on réalisera un second levé topographique depuis la dernière station du levé
jusqu’à la première et ceci en utilisant un nombre reduit de position pour le théodolite (voir
figure 2.12). On obtiendra ainsi deux mesures distinctes de la différence de niveau entre la
première et la dernière station : l’une à l’aller et la seconde au retour. Un bon levé est bien
entendu un levé où les deux mesures sont identiques (∆e = 0). Néanmoins, il est fréquent
d’avoir une différence notable entre ces deux mesures. Il faut alors redistribuer cette erreur
sur les N stations de mesures, en supposant que le chemin retour est juste. Pour ce faire, on
corrigera chaque mesure de niveau de la quantité ∆e /N.
Pour les levés régionaux, il est possible d’utiliser des baromètres pour mesurer l’élévation
des stations gravimétriques. Le nivellement au baromètre peut avoir une précision meilleure
que 1.5 m aux conditions suivantes :
1. Utiliser des instruments de précision.
2. Faire usage d’un baromètre témoin (enregistreur).
3. Opérer les instruments seulement dans le cas de température favorable.
4. Revenir le plus souvent possible au point de base ou à un point dont l’élévation est
connue. (Boucle d’une heure au maximum).
5. Procéder au levé en évitant les longues boucles (grande distance entre les points).
2.2 Levé gravimétrique 29
aller
Stations
1 2 3 4 5 6 7 8 9
retour
F IG . 2.12:
6. Prendre un levé précis de la température de l’air en même temps que celui de la pression.
Les causes d’erreur les plus communes dans ce genre de levé sont dues principalement à :
1. Les variations périodiques et non périodiques de la pression atmosphérique.
2. La turbulence de l’air.
3. L’effet du vent.
4. Les variations dans la température de l’air.
Suggestions :
1. Dans tout levé gravimétrique, il faudrait minimiser au maximum l’usage du baromètre.

2. Si le levé couvre un grand territoire, on peut utiliser un GPS de haute résolution ; sinon,
faire un levé topo au niveau optique.
3. Dans le cas de levé au niveau, il faudrait faire usage d’un niveau automatique afin de
diminuer les temps d’opération.
4. Si le territoire est vaste il serait souhaitable d’avoir deux équipes de niveau sur le terrain.
2.2.3 Résumé pour faire un levé gravimétrique

– Nivellement : ne pas oublier que pour une précision de 0.01 mgal, il faut connaı̂tre
l’élévation à ± 3 cm pour la correction d’air libre et à ± 9 cm pour la correction de
plateau.
– Manipulation du gravimètre : le gravimètre est une appareil extrêmement sensible et ne
doit pas subir de chocs. Un soin particulier doit être apporté lors de son transport et de
son utilisation.
F IG . 2.13:
F IG . 2.14:
– Dérive : il faut repasser à un point de contrôle à toutes les 3-4 heures. Le cheminement
employé dépend du terrain sur lequel les mesures sont prises et le temps requis pour
faire ces mesures. Le plus important est d’établir un bon réseau de stations de base. Sur
une grille traditionnelle d’exploration, on établiera les stations de base sur la ligne de
base ou sur une des lignes de rattachement (voir figures 2.13 et 2.14).
– Positionnement des stations : il est souhaitable d’éviter le faire des mesures à proximité
des accidents topographiques, de façon à faciliter les corrections de relief.
– Personnel : une équipe de gravimétrie devrait être composée d’au moins deux per-
sonnes : l’opérateur du gravimètre et son aide. Ce dernier pourrait, lorsque les condi-
tions de terrain ne sont pas trop difficiles, faire les calculs nécessaires à l’obtention de la
carte de Bouguer.
– Véhicule : selon les terrains étudiés, il devrait avoir les qualités suivantes :
1. être tout terrain et avec treuil ;
2. être fiable ;
3. être ni trop gros, ni trop pesant ;
4. posséder une bonne manoeuvrabilité.
2.3 Instrumentation
Il y a deux types de mesures : absolues et relatives.
2.3.1 Mesures absolues
Pendule
La mesure est obtenue à partir de la relation :
4π 2 I
g= (2.23)
T 2 mh
2.3 Instrumentation 31
F IG . 2.15:
F IG . 2.16: Principe des appareils à corde vibrante.
où I est le moment d’inertie, T la période d’oscillation, m la masse et h la distance du pivot

au centre de masse du pendule (voir figure 2.15). Cette équation est remplacée par la formule
pour le pendule idéal (la connexion entre le pivot et la masse est parfaitement rigide et sans
poids) pour devenir
4π 2 I
g= . (2.24)
T2
On peut mesurer la valeur absolue de g avec le pendule, mais la précision n’est pas très
élevée : on n’a pas encore réussi à en construire qui soit précis à 0.1 mgal.
Cordre vibrante
Le principe de la corde vibrante (voir figure 2.16) est de déterminer la fréquence de résonnance
entre d’une part la corde soutenant la masse m et le circuit (solénoı̈des) électronique, cette
fréquence étant proportionelle à g.
Les appareils reposant sur ce principe sont encore au stade expérimental.
F IG . 2.17:
Chute libre
Suivant le principe de la chute libre, on a :

1 2
z= gt (2.25)
2
et donc, si un corps de vitesse initiale inconnue, tombe de distance z1 et z2 dans des temps
respectifs de t1 et t2 , alors
( z2 − z1 )
g=2 . (2.26)
( t2 − t1 )2
Pour une précision de 1 mgal sur une chute de 1 à 2 m, le temps doit être connu à 10−8 s
et la distance à 0.5 µm. On a réussi à construire des appareils reposant sur ce principe grâce à
l’emploi de lasers (voir figure 2.17).
2.3.2 Mesures relatives
Trois instruments : pendule, balance à torsion et gravimètre.
Pendule
Suivant le principe exposé précedement, on a gT 2 =constante, et donc que

d( gT 2 ) = 2gT∆T + T 2 ∆g = 0. (2.27)
lecture
lampe
miroir
plaques
reflechissantes
m m
h 2l
F IG . 2.18: Principe de la balance à torsion.
Alors,
∆T t2 − t1
∆g = −2g = −2g . (2.28)
T t1
Si la période T peut être mesurée à 1 µs, la précision sur ∆g sera de l’ordre de 1 µgal. Dans les
années 30, les appareils offraient une sensibilité de ∼ 0.25 µgal. Avec ces appareils, le temps
requis pour effectuer une mesure est d’environ 30 minutes.
Balance à torsion
La balance de torsion est l’ancêtre du gravimètre. L’apareil est formée par deux masses
égales séparées par une barre rigide de longueur 2l (horizontale) et une d’une hauteur h (ver-
ticale). Le système est suspendue en son centre par une fibre de torsion à laquelle est attachée
un petit miroir afin de mesurer la rotation d’un rayon lumineux fourni par une lampe (voir
figure 2.18).
A partir des variations du faisceau lumineux lues sur l’écran, mesurera le gradient (varia-
tion) horizontal de la gravité. On ne mesure pas gz car le mouvement n’est que rotationnel et
causé par de petites différences dans la composante horizontale de g agissant sur deux masses.
Les mesures sont données en Eötvös, égales à 10−6 mgal/cm.
Les gravimètres
Développés pour mesurer ∆gz sur le terrain (∼ 1930), le système correspond essentielle-
ment à une balance extrêmement sensible dans laquelle une masse est reliée à un ressort. Les
variations de la gravité se traduisent alors par une élongation du ressort amplifiée mécanique-
ment ou électriquement. Les gravimètres modernes utilisent deux ressorts : un dont la tension
correspond à une valeur moyenne pour la région et un autre plus sensible relié à une vis mi-
crométrique qui sert à faire la lecture.
On retrouve deux types de gravimètres : les gravimètres stables et les instables.
F IG . 2.19:
Les gravimètre stables
Ces gravimètres sont établis suivant le principe de la loi de Hooke. La variation de la gra-
vité est égale à la variation de la force exercée par le ressort (voir figure 2.19) :
k
∆g = ∆x (2.29)
m
où k est la constante du ressort, m la masse et ∆x l’élongation du ressort. Pour un ressort ayant
une période d’oscillation T, nous obtenons
4π 2
g= ∆x (2.30)
T2
où la période d’oscillation est donnée par
r
m
T = 2π . (2.31)
k
Pour détecter ∆g de 0.1 µgal, ∆x doit être de l’ordre de 10−7 cm. Il faut dont des ressort
extrêmement sensibles.
Un exemple, le Gravimètre Gulf :
– Il mesure la rotation d’un ressort construit à partir d’un ruban métallique plat en s’étirant,
il engendre un mouvement relatif transmis à un miroir qui lui est fixé, la déflexion du
rayon lumineux est amplifiée puis lue ;
– sensibilité > 0.1µgal ;
– initialement 100 lbs, réduit à 25.
F IG . 2.20:
F IG . 2.21:
Les gravimètre astables
Les Gravimètre astables sont plus précis que les gravimètres stables.
– La Coste - Romberg : développé en 1934 par J.B. LaCoste, il est basé sur le principe du
ressort de longueur zéro : la tension ∝ longueur du ressort (voir figures 2.20 et 2.21).
Son degré de précision est de ±0.01 mgal, voire mieux. Ils sont fabriqués en métal avec
faible extension thermique et sont isolés et thermostatés ±0.002˚C. Les premiers pesaient
environ 80 lbs (1940), maintenant, seulement 6 lbs.
– Worden (Sodin) : développé en 1948, il possède un mécanisme en quartz (très léger, le
mécanisme est gros comme un poing). Sa sensibilité aux ∆T et ∆P est réduite parce que le
mécanisme est sous vide. Il possède également un système de compensation thermique.
Le mouvement est similaire au LaCoste - Romberg (voir figures 2.22). Ses caractéristiques
techniques sont : ∼ 1000 de haut, 500 de diamètre, et environ 6 lbs. Le prospector de Sodin
F IG . 2.22:
a une précision de ∼ 10 µgals ; Coût ∼ 29,000 Can $ (1998).

– Scintrex CG-3 : gravimètre “électronique” : les mesures sont automatiques et multiples
(nivellement automatique, correction de marées, interface avec ordinateur) Sensibilité
∼ 1µgal ; Précision ∼ 3µgals ; Coût ∼ 45,000 Can $ (1998). Version autonivellante
héliportée pour les régions d’accès difficiles : l’Héligrav de Scintrex.
2.4 Traitements des données
2.4.1 Séparation régionale - résiduelle
L’anomalie de Bouguer peut provenir de plusieurs niveaux :

1. grande profondeur : ex. : variations du socle métamorphique ;
2. profondeur moyenne : ex. : lentille de sel à l’intérieur d’une colonne sédimentaire ;
3. faible profondeur : variations de l’épaisseur du mort-terrain.
Plus la source est profonde, plus l’anomalie est évasée (voir figure 2.24).
Une fois toutes les corrections appliquées, on obtient une carte de l’anomalie de Bouguer
qui démontre en général deux caractéristiques (l’anomalie de Bouguer représente la somme
de tous les corps sous la surface) :
2.4 Traitements des données 37
F IG . 2.23:
F IG . 2.24:
1. Des variations du champ gravitationnel régulières et continues sur de grande distance

appelées variations régionales. Elles sont produites par les hétérogénéités à grandes pro-
fondeurs. Composante qui varie lentement en ( x, y).
2. Superposées à ces variations régionales, et souvent masquées par celles-ci, on observe
de petites perturbations locales du champ gravitationnel qui sont secondaires en dimen-
sions, mais d’intérêts potentiellement primordiales selon l’application.
Selon le but du levé, il faut :
1. Lisser et enlever les effets de surface pour ne retenir que les effets de profondeur (régionale).
2. Lisser les effets de sources profondes et les soustraire pour obtenir les anomalies de sur-
face (résiduelle).
Les anomalies dites résiduelles sont surtout produites par des hétérogénéités situées dans
la partie supérieure de l’écorce terrestre. Ce sont souvent le résultat de minéralisation ou de
réservoirs. Afin de pouvoir interpréter ces anomalies de façon quantitative, il est nécessaire
de soustraire l’anomalie régionale de nos données. Pour séparer la régionale et la résiduelle,
plusieurs possibilités s’offrent au géophysicien :
– faire une lissage graphique sur le profil ;
– faire un lissage graphique sur les cartes de contours ;
– calculer la régionale analytiquement ou appliquer un filtre (effectué sur ordinateur) ;
– calculer l’effet de la source à éliminer si sa géométrie et sa densité sont connues afin de
le soustraire à l’anomalie Bouguer (modélisation).
2.4.2 Prolongement vers le haut
Le fait que le potentiel gravimétrique obéisse à l’équation de Laplace nous permet de

déterminer sa valeur sur une surface arbitraire si, d’une part, on le connaı̂t complètement
sur une autre surface, et si, d’autre part, aucune masse ne se situe entre les deux surfaces.
Considérons les deux demi-espaces de la figure 2.25. En se référant aux équations de La-
place et de Poisson (équations (1.11) et (1.12)), on a
∇2 UQ = −4πGρ (2.32)
Z
ρ
UP = G dv. (2.33)
V R
En isolant ρ, on obtient
Z
1 1
UP = − ∇2 UQ dv (2.34)
4π V R
qui devient, en utilisant le théorème de Green
Z Z
1 g
UP = dxdy (2.35)
2π x y Rs
où Rs est donné par

R2s = ( x − xO )2 + (y − yO )2 + h2 (2.36)
F IG . 2.25: Prolongement vers le haut.
Finalement, la gravité au point P vaut

Z Z
∂UP h g
gP = = dxdy (2.37)
∂z 2π x y R3s
2.4.3 Exemples de superposition d’une anomalie avec une régional
La sphère
Sur la figure 2.26, on a superposé sur la moitié inférieure l’effet d’une régionale avec gra-
dient posistif vers l’ouest de 0.3 µgal/100m. La régionale a une valeur de 0.2µgal à 475 m au
nord du centre de la sphère.
Cylindre horizontal
Sur la figure 2.27, on a superposé sur la moitié inférieure l’effet d’une régionale avec gra-
dient positif vers l’ouest de 0.3 µgal/100m. La régionale a une valeur de 0.2 µgal à 1000 m à
l’ouest du centre du cylindre.
Faille verticale
Sur la figure 2.28, on a superposé sur la partie supérieure une régionale de gradient 0.3 µ
gal/100m vers l’est et dans la partie inférieure, une régionale de gradient 0.3 µgal/100m vers
l’ouest.
Autres exemples
Figure 2.29
La régional est estimée par une droite de pente négative. La résiduelle est calculée par
simple différence entre les valeurs mesurée (corrigées) et cette droite.
F IG . 2.26: Contours aux 0.2 µgal Rayon de la sphère 100 m ; Profondeur au centre de la sphère
85 m
F IG . 2.27: Contours aux 0.2 µgal Rayon de cylindre 100 m ; Profondeur au centre du cylindre
105 m
F IG . 2.28: Contours aux 0.2 µgal Déplacement vertical de la faille 90 m
F IG . 2.29:
F IG . 2.30:
F IG . 2.31:
Figure 2.30
La régionale est estimée par un plan (cas 2D) déterminer par le prologement de chacunes
des lignes de contours depuis les bordures vers le centre de la carte.
Figure 2.31
Différentes courbes sont calculées en calculant des moyennes sur 5 ou 11 points. Noter ici
que plus on prend un nombre important de points, plus on a tendance à éliminer les petites
variations pour obtenir presque une droite (cas 11 pts).
F IG . 2.32:
2.4.4 Cône de sources
Sur la figure 2.32, la sphère (1) est le corps le plus profond qui peut produire approxima-
tivement l’anomalie graviétrique présentée. Des corps plus superficiels et plus larges, tels que
(2) et (3), pourraient aussi produire des anomalies semblables. Tous auraient la même anoma-
lie de masse totale.
3 Interprétation
En géophysique appliquée, les applications classiques de la gravimétrie sont la délinéation

des corps sous-terrains (la calcul de la position des limites des interfaces entre les corps), et
le calcul des tonnages. Souvent, des modèles géométriques simples ayant une réponse analy-
tique connue donnent un résultat satisfaisant. Si la géométrie est trop complexe, la modélisation
numérique vient à la rescousse.
3.1 Modèles simples
3.1.1 La sphère
L’anomalie d’une sphère peut s’écrire sous la forme (voir figure 3.1) :
Gm 4 3
gr = où M= πa ρ. (3.1)
r2 3
La composante verticale s’exprime donc par
z 4 z
gz = gr cos θ = gr = πa3 Gρ 2 . (3.2)
r 3 ( x + z2 )3/2
Si il existe un contraste de densité entre la sphère et le matériel encaissant, ∆ρ = ρ1 − ρ2 ,

alors l’anomalie gravimétrique de la sphère est
4 z
∆gz = ( πa3 ) G∆ρ 2 . (3.3)
3 ( x + z2 )3/2
L’anomalie maximale se trouve à x = 0 et est égale à
4 3 ∆ρ
∆gmax = πa G 2 . (3.4)
3 z
Examinons maintenant l’endroit particulier de la courbe où ∆g = ∆gmax /2 (voir figure 3.1)
définissant le point x = x1/2 (x1/2 est appelé la «demi-largeur» à la «mi-hauteur»). On a, en
45
46 3. Interprétation
F IG . 3.1: Anomalie gravimétrique et paramètres géométriques de la sphère.
posant V = 43 πa3 ,
∆g = ∆gmax /2
z 1
VG∆ρ 2 2 3/2
= VG∆ρ 2
( x1/2 + z ) 2z
z3 1
2
=
( x1/2 + z2 )3/2 2
" #3
z 1
2
=
( x1/2 + z2 )1/2 2
 3
 1  1
=
 
1/2 
2

 2
x1/2 
1+ z2
x 2 3/2
1/2
1+ = 2
z
x 2 3
1/2
1+ = 4
z
d’où on trouve finalement

z = 1.306 x1/2 (3.5)
Il est donc possible de connaı̂tre la profondeur z du centre de la sphère à partir de l’anomalie

gravimétrique qu’elle produit. Losque z est connu, on peut alors calculer l’excès de masse de
3.1 Modèles simples 47
F IG . 3.2: Anomalie gravimétrique et paramètres géométriques d’un cylindre horizontal.
la spère (∆M)
∆gmax · z2
∆M = (3.6)
G
et, si les densités de milieu encaissant et de la sphère sont connues, la mase réelle de la sphère
(M) est alors obtenue par règle de trois
ρsphère
M = ∆M (3.7)
∆ρ
où ∆ρ = ρsphère − ρencaissant
3.1.2 Le cylindre horizontal
L’anomalie gravimétrique d’un cylindre de longueur L, dont la profondeur du centre est z

et dont le contraste de densité est ∆ρ = ρ1 − ρ2 (voir figure 3.2), est donné par
 
πGR2 ∆ρ  1 1
∆gz = 1/2 − (3.8)

z(1 + ( x/z)2 )
 1/2 
( x 2 + z2 ) ( x 2 + z2 )
1+ ( y + L )2
1+ ( y − L )2
Si le cylindre est infiniment long (L > 10 z), alors l’équation se simplifie pour donner
2πGR2 ∆ρ
∆gz = (3.9)
z(1 + ( x/z)2 )
avec un maximum en x = 0 donné par

2πGR2 ∆ρ
∆gmax = . (3.10)
z
F IG . 3.3: Différence entre l’anomalie gravimétrique causée par une sphère et un cylindre hori-
zontal.
En utilisant ces derniers résultats et, de la même manière que pour la sphère, on a, au point
x = x1/2 , en posant C = 2πGR2
∆gmax
∆g =
2
C∆ρ 1 C∆ρ
=
z(1 + ( x1/2 /z)2 ) 2 z
1 1
x1/2 2 =
(1 + z ) 2
h x i2
(1 + 1/2 ) = 2
z
d’oz on tire
z = x1/2 (3.11)
La profondeur du cylindre est trouvée directement par la valeur de x1/2 . De plus, le cylindre
donne une amomalie plus large que celle d’une sphère (voir figure 3.3).
3.1.3 Le cylindre vertical
Dans le cas du cylindre vertical, on doit intégrer un petit élément dgz donné par
dgz = 2πGρdl sin φdφ (3.12)
pour φ allant de 0 à arctan( R/l ) et puis pour dl variant de z à z + l. Finalement, après quelques
calculs, la valeur maximum de l’anomalie ∆gmax est donnée par
h 1/2 1/2 i
∆gmax = 2πG∆ρ L + z2 + R2 − ( z + L )2 + R2 . (3.13)
Il est à noter ici que la correction de terrain fait par le réticule est donnée par intégration
sur une partie du cylindre
h i
∆gt = Gρθ (r2 − r1 ) + (r12 + L2 )1/2 − (r22 + L2 )1/2 . (3.14)
3.1 Modèles simples 49
F IG . 3.4: Paramètres géométriques du feuillet vertical
3.1.4 Le feuillet vertical
Pour un feuillet vertical (voir figure 3.4), l’anomalie est donnée par
( h + l )2 + x 2

∆gz = 2Gt∆ρ ln . (3.15)
x 2 + h2
Si h ≈ l, le maximum est donné par
∆gmax = 2Gt∆ρ ln(4) (3.16)
et en x = x1/2 (∆gz = ∆gmax /2)

" #
(2h)2 + x1/2
2
1
2Gt∆ρ ln 2 + h2
= Gt∆ρ ln(4)
x1/2 2
4h2 + x1/2
2
2 + h2
= 2
x1/2
2
x1/2
h2 =
z
ce qui donne
h = 0.7x1/2 (3.17)
3.1.5 La plaque mince horizontale infinie
L’anomalie pour la plaque mince horizontale infinie (voir figure 3.5) est
hπ x i
∆gz = 2Gt∆ρ + arctan (3.18)
2 h
F IG . 3.5: Anomalie gravimétrique et paramètres géométriques de la plaque mince horizontale

infinie
où h est la profondeur au plan médian de la faille et non pas la profondeur au toit.
Certains paramètres de la plaque mince peuvent être trouvés à partir de l’anomalie qu’elle
produit. Ainsi, lorsque :
x → −∞ ∆g → 0
x → +∞ ∆g → 2πGt∆ρ
et donc, l’anomalie maximum ∆gmax est égale à
∆gmax = |∆g+∞ − ∆g+∞ | = 2πGt∆ρ. (3.19)
Ceci permet de déterminer un premier paramètre caratérisant la plaque, c’est à dire t∆ρ
∆gmax
t∆ρ = . (3.20)
2πG
De plus, la dérivée (i.e. la pente au point d’inflexion de la courbe) en x = 0 nous donne
∆gmax

∂g h
= 2πGt∆ρ 2 2
= . (3.21)
∂x x=0
h + x x =0 πh
qui nous donne donc le second paramètre pour la plaque, h
∆gmax
h= (3.22)
∂g
π ∂x
0
Si la plaque a une longueur finie L, alors
L−x
x
∆gz = 2Gt∆ρ arctan − arctan . (3.23)
h h
3.2 Modèles complexes 51
F IG . 3.6: Paramètres géométriques du prisme rectangulaire.
3.1.6 Le prisme rectangulaire
L’anomalie pour le prisme rectangulaire (voir la figure 3.6 pour les différents paramètres)
est donné par
r1 r4 r4
∆gz = 2G∆ρ z2 θ2 − z1 θ1 + x ln + b ln (3.24)
r2 r3 r3
Si b → ∞, alors θ1 → ϕ1 , θ2 → ϕ2 et r3 → r4 . L’équation devient alors

r1
∆gz = 2G∆ρ z2 ϕ2 − z1 ϕ1 + x ln (3.25)
r2
Pour un contact (z1 = 0, b → ∞), alors :
(1 + x2 /t2 )1/2
x
x π
∆gz = 2G∆ρ ln + + arctan (3.26)
t x/t 2 t
3.2 Modèles complexes
Lorsque les corps étudiés ne peuvent raisonnablement être approximés par les formes
simples dont on connaı̂t analytiquement les réponses, il est nécessaire de recourir à d’autres
outils. Pour calculer l’anomalie, il existe les méthodes graphiques et les méthodes analytiques.
3.2.1 Les méthodes graphiques
On peut utiliser des graticules pour calculer l’effet des corps. Un graticule est une série
de cellules de différentes formes et grosseur, chacune couvrant une surface correspondant à
une contribution connue (et généralement uniforme) à la valeur de gz au point de mesure en
surface. Il y a deux types.
F IG . 3.7: Graticule à compartiments trapézoı̈daux.
F IG . 3.8: Graticule “Dot chart”
Premier type
Les compartiments trapézoı̈daux (voir figure 3.7) sont formés par des lignées horizontales
équidistantes coupant une série de lignes radiales séparées d’un même angle. L’effet des cel-
lules est
gz = 2Gρ(θ2 − θ1 )(z2 − z1 ). (3.27)
Puisqu’en général θ2 − θ1 = cst et z2 − z1 = cst, l’effet de chacune des cellules est le même. Il
s’agit alors de somme la contribution de chacun des compartiments pour obtenir l’anomalie.
Deuxième type : “Dots charts”
Les compartiments (voir figure 3.8) sont formés par la rencontre de lignes radiales et d’arcs
de cercle. Leur contribution n’est pas uniforme mais est dépendante du nombre de points
qu’ils contiennent. Chacun des points représente une contribution constante à gz au point de
mesure.
Pour calculer l’effet gravimétrique d’un corps enfouis, le “Dot chart”, imprimé sur un
transparent, est superposé sur une section transversale, le vertex est placé sur la position de la
surface où l’effet gravimétrique est désiré.
Remarques :
Si les structures ne sont pas réellement 2-D, on peut appliquer des corrections pour les
effets de bordure. Dans le cas des structures 3-D, la méthode graphique s’applique aussi. Les
corps sont divisés en une série de plaques horizontales dont les effets sont calculés à l’aide de
graticules. Cette approche est compliquée puisqu’on doit introduire un facteur d’échelle pour
tenir compte de la profondeur de chacune des plaques. On voit donc que cette méthode s’ap-
plique bien mieux à l’aide d’un ordinateur (proposé par Hubbert en 1948, repris par Talwani
sur ordinateur en 1959 (2D) et en 1960 (3D)).
3.2.2 Méthode analytique
On peut montrer que g produit par un corps 2D de section quelconque est égale à une
intégrale de ligne autour de cette section
I
g = 2G∆ρ zdθ. (3.28)
Si la section est approximée par un polygone de n côtés (voir figure 3.9), l’intégrale devient
une sommation I n
zdθ = ∑ Zi . (3.29)
i =1
où Zi est l’intégrale de ligne pour le ie côté, qui est égale à
cos θi (tan θi − tan φi )

Zi = ai sin θi cos θi (θi − θi+1 ) + tan θi × (3.30)
cos θi+1 (tan θi+1 − tan φi+1 )
et où les différents paramètres sont définis par

zi
θi = arctan ;
xi
zi+1 − Zi

φi = arctan ;
xi+1 − xi
ai = xi+1 − zi+1 cot φi
xi+1 − xi

= x i +1 − z i +1 .
z i +1 − z i
3.2.3 Gravité 3-D
En trois dimensions, on peut procéder de façon similaire, mais cette fois on coupe le volume
en tranches horizontales (figures 3.10 et 3.11). L’expression pour anomalie revêt la forme
rdrdθ
Z ZZ
∆g(O) = G∆ρ zdz
S (r 2
+ z2 )3/2
m
rdrdθ
Z Z
= G∆ρ ∑ Wj z j 2
. (3.31)
j =1 S j (r + z j )3/2
2
F IG . 3.9: Paramètres géométriques d’un polygône.
F IG . 3.10: Paramètres géométriques du volume découpé en tranches.

F IG . 3.11: Paramètres géométriques du volume découpé en tranches (suite).

z
S1
R
y
x
S2 M
F IG . 3.12: Définition de la demi-sphère utilisé pour le calcul de l’excès de masse.
3.3 Excès de masse
3.3.1 Calcul de l’excès de masse
Supposons une masse M causant un champ gravitationnel g. Entourons cette masse d’une
demi-sphère présentant sa surface plane à z=0 (voir figure 3.12), et appliquons le théorème de
Gauss : Z Z
(∇ · g) dv = g · n ds. (3.32)
V S
où n est normal à la surface de la sphère. Or, g = −∇U. D’où :
Z Z
− ∇2 U dv = g · n ds (3.33)
V S
ou encore, en utlisant la formule de Poisson : ∇2 U = 4πGρ, où ρ est la densité de la masse (ρ

étant égal à zéro partout ailleurs) :
Z Z
4πG ρ dv = 4πGM = − g · n ds (3.34)
V S
où M est la masse du corps causant l’anomalie. Le flux de la force d’attraction à travers toute
surface fermée située dans un champ gravitationnel est égal à 4πG fois la masse située à
l’intérieur de la surface. Notons que ceci est valide pour toute forme de surface et toute dis-
tribution de masse (volume, position, forme). On peut donc calculer la masse produisant une
anomalie gravimétrique à partir du levé sans connaı̂tre la forme ou la profondeur de la masse.
Dans le cas d’un levé gravimétrique, on connaı̂t ∆gz et on peut supposer que la surface de
mesure était plane, ce qui correspond à la surface z = 0 de la demi-sphère.
Sur le plan z = 0, on a (le signe (-) intervient parce que n pointe vers l’extérieur de la
surface) :
∆g · n = −∆gz (3.35)
Pour la surface r = R de la demi-sphère, on aura, en coordonnées sphériques :
∂U
∆g · n = −∇U · r = − (3.36)
∂r
3.3 Excès de masse 57
et donc Z Z Z
∂U
4πG∆M = − ∆g · n ds = ∆gz ds + ds. (3.37)
S S1 S2 ∂r
Si R → ∞, on aura
Z +∞ Z +∞ Z π Z 2π
∂U 2
4πG∆M = ∆gz dx dy + lim R sin θ dθ dϕ (3.38)
−∞ −∞ R→∞ θ =π/2 ϕ =0 ∂r
Pour R très grand, le potentiel n’est plus vu que comme celui d’une masse ponctuelle M à
l’origine, ce qui conduit à
G∆M
UR → ∞ = . (3.39)
R
Alors
Z +∞ Z +∞
G∆M
Z π
∆gz dx dy = 4πG∆M − 2π lim sin θ dθ
−∞ −∞ R→∞ π/2 R2
= 4πG∆M − 2πG∆M
= 2πG∆M (3.40)
L’excès de masse est ainsi

Z +∞ Z +∞
1
∆M = ∆gz dx dy (3.41)
2πG −∞ −∞
Puisque nous travaillons avec ∆gz , ∆M est l’excès de masse du corps causant l’anomalie
et non pas sa masse réelle. Pour connaı̂tre la masse réelle du corps, il faut utiliser un facteur
correctif
ρM
M = ∆M (3.42)
ρ M − ρE
où ρ M est la densité de la masse et ρ E est la densité du milieu encaissant.
3.3.2 Unicité de la solution
Nous savons que des corps de forme, volume et positions spatiales différentes peuvent
donner des anomalies gravimétriques identiques. On peut donc se demander si la détermina-
tion de la masse causant l’anomalie est unique.
Supposons deux masses M1 et M2 , produisant des anomalies ∆gz1 et ∆gz2 partout égales
sur le plan de référence. On aura donc :
Z Z
∆gz1 ds = 2πG∆M1
Z Z
∆gz2 ds = 2πG∆M2 (3.43)
qui est indépendant de la géométrie (forme, volume, position spatiale) des masses puisque
par hypothèse ∆gz1 = ∆gz2 . Les deux intégrales sont donc identiques et M1 doit égaler M2 . La
solution au calcul de l’excès de masse est donc unique.
3.3.3 Considérations pratique
Dans le cas d’intégration numérique, l’équation(3.42) devient

S
ρM
Tmin = 23.9 ×
ρ M − ρE ∑ ∆g∆s (3.44)
où Tmin est le tonnage minimum en tonnes métriques, ρ M la densité de la masse (g/cm3 ), ρ E la
densité de la roche encaissante (g/cm3 ), ∆g l’intensité de l’anomalie (mgal), ∆S la surface (m2 )
où l’anomalie a la valeur ∆g et ∑ la sommation sur la superficie S du produit ∆g∆s. Puisqu’en
général la surface d’intégration S est finie, on doit appliquer un facteur correctif dépendant de
la géométrie de l’anomalie.
1. Si l’anomalie est ± circulaire

1
T = Tmin × (3.45)
1 − h/R
où T est le tonnage total, R le rayon de la superficie circulaire sur laquelle la sommation
est faite et h la profondeur du centre de gravité
√ de la masse. Si S n’est pas parfaitement
circulaire, on prendra le rayon moyen R = S/π.
2. Si l’anomalie est rectangulaire
 
 π/2 
T = Tmin ×   (3.46)
xy
 
arctan 1h √ 2
x + y2
où 2x2y est la superficie rectangulaire de sommation (x > y) et h la profondeur du centre

de gravité de la masse.
3. Si S n’est pas parfaitement rectangulaire et qu’il est difficile d’obtenir des valeurs précises
de x et y, il est préférable de trouver le rapport approximatif x/y et d’utiliser
 
π/2
T = Tmin ×  q  . (3.47)
1 S x/y
arctan 2h 1+( x/y)2
Détermination approximative de h
1. Si l’anomalie est circulaire en plan et relativement aiguë, h = 1.3x1/2 , où x1/2 est la demi-
largeur à la demi-hauteur.
2. Si l’anomalie est au moins cinq fois plus longue que large en plan et relativement aigüe,
h = x1/2 .
3. Pour les longueurs intermédiaires, choisir un coefficient entre 1 et 1.3.
4. Si l’anomalie présente un maximum très plat, la valeur de h peut être beaucoup trop
grande et T sera une valeur maximum.
3.3.4 Exemple de calcul de tonnage
• G1
Paramètres pour G1 (figure 3.13)

– Echelles : 1 pouce carré = 14 900 m2
– Intensités = milligal
– Minerai : densité ρ1 = 2.92 g/cm3
– Roche encaissante : densité ρ2 = 2.8 g/cm3
• G2

– Echelles : 1 pouce carré = 14900 m2
• G3

– Echelles : 1 pouce carré = 21850 m2
3.3.5 Exemple de calcul avec G-2 Marmora
Données :
– 1 po = 1 000 pi ; 1 po2 = 92903 m2
– 1 careau = 1/400 po2
– ρ1 = 4.00 g/cm3 ; ρ2 = 2.80 g/cm3
Mesures des surfaces : on compte les carreaux !

– Contours 8.0-10.0 milligals : 28 carreaux
ligne Nb de carr.
1 8
2 13
3 7
F IG . 3.13:
F IG . 3.14:
ligne Nb de carr. ligne Nb de carr.

1 2 7 7
2 19 8 1
3 11 9 7
4 10 10 10
5 6 11 13
6 12
ligne Nb de carr. ligne Nb de carr.

1 10 10 3
2 3 11 6
3 10 12 complète
4 4 13 1
5 14 14 7
6 7 15 6
7 6 16 5
8 6 17 10
9 14 18 6
F IG . 3.15:
ligne Nb de carr. ligne Nb de carr. ligne Nb de carr.

1 8 11 8 21 12
2 17 12 2 22 6
3 7 13 5 23 1
4 8 14 24 24 10
5 6 15 6 25 6
6 5 16 16 26 5
7 5 17 3 27 1
8 10 18 4 28 5
9 19 19 15 29 13
10 2 20 17
ligne Nb de carr. ligne Nb de carr. ligne Nb de carr.

01 07 13 12 24 13
02 05 14 20 25 11
03 06 15 1 26 5
04 06 16 13 27 4
05 06 17 19 28 8
06 05 18 3 29 5
07 06 19 13 30 2
08 07 20 21 31 9
09 06 21 2 32 5
10 02 22 9 33 3
11 09 23 8 34 2
12 08
Calcul de T min
Contours ∆g Nb ∆S ∆S ∆g × ∆S
(mgal) (mgal) div. po2 m2 (mgal m2 )
8.0-10. 9.0 28 0.0700 6 503 58 527
6.0-8.0 7.0 88 0.2000 20 439 143 073
4.0-6.0 5.0 118 0.2950 27 406 137 030
2.0-4.0 3.0 146 0.6150 59 135 171 405
0.0-2.0 1.0 262 0.6550 60 852 60 852
TOTAL : 742 1.8550 172 335 570 887
F IG . 3.16:
S
ρ1
ρ1 − ρ2 ∑
Tmin = 23.9 ∆g × ∆s

4.00
= 23.9 570 880 ≈ 45.5 × 106 (3.48)
4.00 − 2.80
Calcul de T
Par la meilleure formule pour G-2
 
π/2
T = Tmin ×  h q i (3.49)
1 S x/y
arctan 2h 1+( x/y)2
1. Calcul de x, y (voir figure 3.16). On obtient :
2x ≈ 2.40 po 2y ≈ 0.75 po
( x/y) = 3.20 2x × 2y = 1.80 po2 ≈ S = 1.855 po2

(3.50)
2. Calcul de h (voir figure 3.17) : On a :
2x1/2 = 3.40 po x1/2 = h = 1.7 po
Nous avons alors

p √
1 S x/y 1 1.855 × 3.20
p = p = 2.137 (3.51)
2h 1 + ( x/y) 2 0.34 1 + (3.20)2
F IG . 3.17:
Ainsi
π/2
T = Tmin ×
arctan(2.137)
= 63.1 × 106 tonnes métriques (3.52)
Par l’autre formule pour S rectangulaire
 
 π/2 
T = Tmin ×   (3.53)
xy
 
arctan 1h √ 2
x + y2
avec
1 xy 1 1.2 × 0.375
p = √ = 2.105 (3.54)
h x +y
2 2 0.17 1.22 + 0.3752
Ainsi
π/2
T = Tmin ×
arctan(2.105)
Ce résultat n’est pas tellement différent du précédent. Ceci provient du fait que 2x × 2y ≈
S.
Par la formule pour S circulaire

1
T = Tmin × (3.56)
1 − h/R
Comme l’anomalie n’est pas vraiment circulaire, il faut prendre, pour calculer T, la valeur
de R donné par : r r
S 1.855
R= = = 0.7684 (3.57)
π π
Ainsi
" #
1
T = Tmin × 0.17
1− 0.7684
Ce dernier résultat est assez diférrent des deux précédent, mais beaucoup mieux que Tmin .
Dans ce cas, il faut accepter TM = 63.1 × 106 tonnes métrique.
4 Signature des structures
géologiques en gravimétrie
4.1 Levés régionaux et tectoniques
4.1.1 Étude structurale à grande échelle
Canada : Fig.4.1
F IG . 4.1: Carte d’anomalie gravimétrique du Canada. L’intervalle des lignes de contour est de
20 ; l’intervalle des lignes de contour de couleurs est de 40 mgal.
67
68 4. Signature des structures géologiques en gravimétrie
États-Unis d’Amérique : Fig.4.2 et 4.3
F IG . 4.2: Carte de la gravimétrie régionale de Bouguer des États-Unis d’Amérique composée

d’ondes de longueur plus grande que 250 km. L’intervalle des lignes de contour est de 20 mgal.
F IG . 4.3: Carte montrant les éléments géologiques majeurs observables sur la carte d’anomalie
gravimétrique.
4.1 Levés régionaux et tectoniques 69
Amérique du sud : Fig.4.4 et 4.5
F IG . 4.4: Amérique du sud.
F IG . 4.5: Carte géologique de l’Amérique du sud.

4.1.2 Études régionales
L’information gravimétrique permet de caractériser les provinces géologiques et de suivre

leurs tracés si l’opposition est faible.
État de Washington : Fig.4.6
F IG . 4.6: Gravimétrie régionale de Bouguer de l’ouest de l’État de Washington (partie ouest).

L’intervalle de lignes de contour est de 5 mgal. Elle indique un minimum coı̈ncidant avecla
chaı̂ne des Cascades, un maximum avec la Coast Range et un minimum avec Puget Sound et
Olympic Mountains.
4.2 Pétrole 71
Vallée Owens en Californie : Fig.4.7
F IG . 4.7: Gravimétrie de Bouguer du sud de la Vallée Owens en Californie ainsi que la relation
à la géologie de cette région. Les contours montrent un minimum fermé asymétrique avec son
axe au centre de la vallée. La gravité est faible dans les régions riches en aluminium, dont l’ori-
gine a été interprétée comme provenant d’une épaissse dépression sédimentaire (sédiments
clastiques). L’interprétation indique que le socle est faillé et est couvert au centre par presque
10,000 de sédiments.
4.2 Pétrole
Dans ce cas, on cherche des structures pièges : dômes de sel, anticlinaux, formations récifales
(“reef”) ou à connaı̂tre l’épaisseur des sédiments ou d’un bassin. Le pétrole a été découvert
dans des anciennes formations récifales de calcaire et on s’est intéressé à connaı̂tre leur signa-
ture en gravité. Ceci n’est possible que s’il y a contraste de densité ρ entre la formation récifale
et les roches sédimentaires qui l’entourent. Malheureusement, la densité et porosité du cal-
caire et du matériel de remplacement sont si variables qu’il n’y a pas de règles générales ;
seulement, dépendant du contexte et de la situation géographique, on peut développer des
signatures particulières.
4.2.1 Formations récifales
Dawn no 156, Ontario Fig.4.8 à 4.10
Ce champ est situé dans le bassin de l’Illinois, recouvert de dépôts glaciaires lourds, d’al-
lure irréguliers, pouvant atteindre 30m d’épaisseur. Il n’y a pas de correction de relief, le pays
est plat. C’est un bassin primaire avec des couches de sel d’épaisseur variable dans le silu-
rien. Les formations recherchées récifales sont d’âge silurien. Le Précambrien donne de larges
anomalies de types type régional, dues à des variations de densité du socle.
Stations aux 300m ; Erreurs de nivellement < 15 cm ; Erreurs de planimétrie < 8 m
F IG . 4.8: Carte géologique du sud-ouest de l’Ontario.

4.2 Pétrole 73
F IG . 4.9: Gravimétrie de Bouguer (intervalle - 0.1 mgal).
F IG . 4.10: Carte résiduelle (intervalle - 0.1 mgal).
Anomalie positive due au calcaire récifal entouré de sel. Le recif silurien, noyé dans le
sel (Dawn,156) apparaı̂t décalé largement vers l’ouest. La résiduelle (fig 4.10) permet de le
localiser correctement. On savait aussi à l’avance que les récifs siluriens sont souvent alignés
et donnent des anomalies circulaires en Ontario. L’anomalie de gravité dans le cas précis peut
résulter :
1. Soit d’une accumulation de galets glaciaires.
2. Soit d’une variation de l’épaisseur de sel.
3. Soit de la présence d’un relief. La porosité d’une formation récifale pouvant beaucoup
varier, il ne faut pas s’attendre à une localisation précise de l’anomalie recherchée.

La formation récifale recherchée a été découverte à une profondeur de 300m.
Champ de Cement (Oklahoma) Fig.4.11 et 4.12
Difficile à localiser sur l’anomalie de Bougeur à cause d’une forte régionale (fig 4.11), la
carte de la dérivée seconde (fig 4.12) donne une bonne image de l’anticlinal.
F IG . 4.11: Champ de Cement (Oklaoma). Bouguer avec un intervalle de 0.5 mgal.
F IG . 4.12: Champ de Cement (Oklaoma). Dérivée seconde obtenue avec un pas de 1 km.
4.2.2 Dômes de sel
Grand Saline à l’Est du Texas Fig.4.13
Le dôme de sel du Grand Saline à l’Est du Texas à son sommet à 250’ de la surface. Les
corrections d’air libre, Bouguer et de latitude ont été effectuées, mais il n’y a aucun ajustement
pour la direction de la régionale. Chaque unité gravimétrique est de 0.1 mgal, alors l’intervalle
entre le lignes de contour est de 1 mgal. L’anomalie est négative (fig 4.13). La zone pointillée
représente la position de la masse de sel.
4.2 Pétrole 75
F IG . 4.13: Grand Saline
Dômes (4) de sel sur la côte du Golfe du Mexique au Texas Fig.4.14 à 4.17
L’anomalie de Bouguer (fig 4.14) donne une image floue et la résiduelle (fig 4.15) qui
montre une anomalie négative ne permet pas de séparer les effets des diverses sources. Sur la
carte de la dérivée seconde (fig 4.16), calculée avec un pas de 1.8km, les 4 sources apparaissent
clairement. Ce traitement est à comparer avec l’exemple théorique des 3 sphères (figure 4.17).
F IG . 4.14: Bouguer. Intervalle des courbes 0.5 mgal.
F IG . 4.15: Résiduelle.
4.2 Pétrole 77
F IG . 4.16: Dérivé seconde.
F IG . 4.17: Exemple téorique de la représentation par le calcul de la dérivé Baranov d’une ano-
malie provoquée par trois sphères pesantes.
4.2.3 Anticlinaux
Anticlinal de Cole Creek (Wyonning 1968) Fig.4.18 et 4.19
La régionale est régulièrement pentée, positive d’est en ouest, brusquement accidentée

d’un replat important (fig 4.18). L’anticlinal assez régulier présente une image bien reflétée
par la résiduelle gravimétrique (fig 4.19).
F IG . 4.18: Bouguer de la région de cole Creek. Graduation tous les 0.5 mgal.
F IG . 4.19: Carte résiduelle de Cole Creek.

4.2 Pétrole 79
Anticlinal de Wellington (Colorado) Fig.4.20 à 4.22
L’anticlinal à des flancs assez raides. Son expression est peu visible sur la carte de Bouguer
(fig 4.20) à cause de la régionale plane (gradient de 2.2 mgal/mile). La résidulle obtenue par
soustraction de la régionale N-S de l’anomalie de Bouguer marque très nettement l’anticlinal
(fig 4.21). Il s’agit d’un cas idéal : fort gradient des isogals, rentrant des courbes très accusé et
régulier. Les mesures ont été prises avec un gravimètre sensible à 0.03 mgal.
F IG . 4.20: Bouguer de Wellington avec isogrammes tous les 0.2 mgal.

F IG . 4.21: Résiduelle graphique de Wellinton par la méthode des profils.
F IG . 4.22: Coupe selon A-B
Chapelet d’anticlinaux de Kettleman Hills, Californie (Boyd, 1946) : Fig.4.23 à 4.25
Le chapelet d’anticlinaux est bien dessinés par la résiduelle analytique (fig 4.24) sur une
Bouguer à fort gradient (fig 4.23). Il s’agit d’un chapelet d’anticlinaux, disposés NO-SE de-
4.2 Pétrole 81
puis Kettleman Nord, au nord jusqu’à Hills, au sud. Deux structures sont à l’aplomb d’un
maximum gravimétrique (+), une autre structure, qui appartient au même alignement cor-
respond au contraire à un minimum (-). Cette inversion est due à un changement de faciès
stratigraphique du nord au sud. à l’aplomb d’un minimum gravimétrique, on a des argiles à
diatomées, très légères, et vers le nord, des marnes plus denses. L’anticlinal occupe le volume
situé à l’aplomb d’un maximum de g. Des isobathes (obtenus par forage) (fig 4.25) tracés sur
la formation Temblor tous les 500’ (180m) montrent la très bonne corrélation entre la forme de
la structure et la carte de Bouguer.
F IG . 4.23: Bouguer du Chapelet d’anticlinaux de Kettleman Hills (Californie), contour tous les
0.5 mgal.
F IG . 4.24: Résiduelle du Chapelet d’anticlinaux de Kettleman Hills (Californie), contour tous

les 0.5 mgal.
4.3 Vallée alluvionnaire 83
F IG . 4.25: Isobathes obtenus par forage, sur des horizons présentant des variations brusque de
densité. Chapelet d’anticlinaux de Kettleman Hills (Californie).
4.3 Vallée alluvionnaire
L’objectif ici est de déterminer la géométrie et l’épaisseur du remplissage d’une vallée par
les alluvions.
Vallée de Portmadoc en Galles du Nord : Fig.4.26 et 4.27
L’épaisseur d’alluvion est dépendante du contraste de densité considéré. On trouve 240 m

avec ∆ρ = 0.7g/m3 et 360 m avec ∆ρ = 0.5g/m3 .
F IG . 4.26: Anomalie de Bouguer, entre Portmadoc et Harlech.

4.4 Batholite gravimétrique 85
F IG . 4.27: Anomalie résiduelle, entre Portmadoc et Harlech.
4.4 Batholite gravimétrique
Batholite granitique : Fig 4.28
Au milieu de volcaniques ou métamorphiques, l’anomalie circulaire est négative alors

qu’au milieu de sédimentaires, elle est positive.
F IG . 4.28: Effet gravimétrique d’un batholite granitique.
4.5 Gisements métallifères
Il faut s’attendre à des anomalies aussi petites que 0.05 mgal économiquement intéressantes.
La précision dans les relevés de terrain doit être très élevée : c’est-à-dire des lectures précises
au gravimètre et au nivellement. Une erreur de 0.016 mgal est probable sur une simple obser-
vation.
Cuba : Fig 4.29
Recherche de minerais de chromite (1945) de densité de ρ ∼ 4.0g/cm3 particulièrement

favorable à l’exploration par gravimétrie.
4.5 Gisements métallifères 87
F IG . 4.29: Cuba : Anomalie gravimétrique au dessus d’un dépot connu de chromite. L’in-
tervalle entre les lignes de contour est de 0.05 mgal. Les cercles blanc sont des stations gra-
vimétriques et les noirs sont des trous de forages.
New Hosco, Abitibi (1950) (ouest de Matagami) Fig 4.30
Détermination de l’anomalie en gravité pour mieux évaluer le corps (cuivre, zinc). Il s’agit
d’un gisement de 2400’ de long, 20 à 150’ de large et de pendage de 65˚ à 70˚, avec 1, 800, 000 tonnes
à 2.5% de Cu et 780, 000 t à 8.15% de Zn.
F IG . 4.30: New Hosco (Abitibi)
4.6 Archéologie, travaux publics
Il s’agit bien souvent de caractériser d’anciennes galeries ou cavités (ex : port de Montréal,
Louis-Riel). Les mesures sont prises avec une très grandes précision. Il faut une précision
au cm sur l’altitude (2 à 3 µgal) et d’environ 10 cm sur la position. Le gravimètre doit être
également extrèment précis (5 µgals = 0.005 mgal) et sensible (1 µgal). Les contrastes de den-
sités sont forts (roche/air) et sont de l’ordre de 1.5 à 2.5 pour des cavités en travaux publics.
4.6 Archéologie, travaux publics 89
Carriere : Fig.4.31
F IG . 4.31: Anomalie gravimétique au-desus d’une carrière. Les courbes sont exprimées en 0.01
mgal.
4.7 Autres Exemples
Consulter : Hinze, W.J. (1985). The utility of regional Gravity and Magnetic anomaly maps.
SEG
4.7.1 Prolongement et filtrage ; modélisation
Texas (Ouest) : Fig.4.32 à 4.38
F IG . 4.32: Localisation : Partie ouest du Texas

4.7 Autres Exemples 91
F IG . 4.33: Anomalie de Bouguer de la partie ouest du Texas. Contours aux 5 mgals
F IG . 4.34: Anomalie de Bouguer ; Traitement : Prolongement vers le haut 5 km ; Contours aux

5 mgals
F IG . 4.35: Anomalie de Bouguer ; Traitement : Filtre Passe-Haut et prolongement vers le haut

2 km (longueurs d’ondes de 80 à 4 km) ; Contours aux 2 mgals
F IG . 4.36: Anomalie de Bouguer ; Traitement : Filtre Passe-Bas et prolongement vers le haut

2 km (longueurs d’ondes de ∞ à 80 km) ; Contours aux 5 mgals
F IG . 4.37: Interprétation d’une coupe EW passant dans le “Salt Basin Grabben” (voir fig 4.36)
F IG . 4.38: Interprétation d’une coupe SW-NE passant dans le “Delaware basin” et le “Central
basin platform” (voir fig 4.36)
4.7.2 Modélisation
Grennville Fig.4.39 et 4.40
Une des anomalie négative majeure du Canada poursuit le front du Grennville sur environ
1200 km (Voir article de Thomas).
F IG . 4.39:
F IG . 4.40:
Domaine volcanique Fig.4.41 à 4.48
F IG . 4.41:
F IG . 4.42: F IG . 4.43:
F IG . 4.44:
F IG . 4.45: F IG . 4.46:
F IG . 4.47: F IG . 4.48:
F IG . 4.49:
4.7.3 Applications : Étude de cas
Articles à consulter
– V.K. Gupta and F.S. Grant, Mineral-exploration aspects of gravity and aeromagnetic sur-
vey in the Sudbury-Cobalt area, Ontario
– R.J. Whitelay. Geophysical case study of the Woodlawn Orebody New South Wales, Aus-
tralia.
5 Références
• Adller, J.L. [ 1941 ]. Simplification of Tidal Corrections for gravity meter-surveys. An-
nual Meeting, Houston, Texas, Avril 1941.
• Blaricom, R.V. , [ 1980 ]. Practical Geophysics for the Exploration Geologist. Northwest-
Mining Association, Spokane.
• Bible, G.L. [ 1962 ]. Terrain Correction Tables for gravity. Geophysics, vol. 27, p. 716-718.
• Compagnie Générale de Géophysique , Prospection gravimétrique. Manuel. 2e édition,
tome II, Paris.
• Dobrin M.B. [ 1988 ]. Intoduction to geophysical propecting. McGraw-Hill.
• Grant, F.S. et West, G.F. [ 1965 ]. Interpretation theory in Applied Geophysics. McGraw-
Hill.
• Griffin, W.R. [ 1949 ]. Residuel Gravity in Theory and Practice. Geophysics, vol. 14, pp.
39-56.
• Hammer, S. [ 1945 ]. Estimating ore masses in gravity prospecting, Geophysics, Vol.
10(1), p. 50-62
• Hinze, W.J. [ 1985 ]. The utility of regional Gravity and Magnetic anomaly maps. Society
of Exploration Geophysicists
• Ivanhoe, L.F. [ 1957 ]. Chart to check elevation factor effects on gravity anomalies. Geo-
physics, vol. 22, no. 3, pp. 643-645.
• Keary, P. et Brooks, M. , [ 1991 ]. An introduction to Geophysical Exploration. Blackwell
Scientific.
• Nettleton, L.L. [ 1939 ]. Determination of density for reduction of gravimeter observa-
tions. Geophysics, vol. 4(3), p. 176-183.
• Nettleton, L.L. [ 1971 ] Elementary Gravity and Magnetics for Geologists and Seismolo-
gists. Society of Exploration Geophysicists, Monograph Series 1.
• Parasnis, D.S. [ 1962 ]. Principles of Applied Geophysics. Chapman & Hall.
• Schoeffler J. [ 1975 ]. Gravimétrie appliquée aux recherches structurales et à la prospec-
tion pétrolière et minière. Editions Technip.
• Reynolds, J. M. [ 1997 ]. An Introduction to Applied and Environmental Geophysics.
John Wiley & Sons.
• Telford, W.M., Geldart, L.P. et Sherif, R.E. [ 1990 ]. Applied Geophysics. Cambridge
University Press.
99
100 5. Références
A Correction de latitude
L’équation pour g( ϕ) s’écrit sous la forme
g( ϕ) = c1 1 + c2 sin2 ϕ + c3 sin2 2ϕ

(A.1)
Donc
dg
= c1 [2c2 sin ϕ cos ϕ + 4c3 sin 2ϕ cos 2ϕ]
dϕ
= c1 [c2 sin 2ϕ + 2c3 sin 4ϕ]] (A.2)
En considérant que l = Rϕ, alors

dl
=R (A.3)
dϕ
et donc
dg dg dϕ c
= = 1 [c2 sin 2ϕ + 2c3 sin 4ϕ] (A.4)
dl dϕ dl R
En prenant c1 =978.0524 gals, c2 = 5.297 × 10−3 , c3 = −5.9 × 10−6 et R = 6367 km (I.U.G.G.

1967), on obtient :
dg h i
= 8.137 × 10−4 sin 2ϕ − 1.813 × 10−6 sin 4ϕ
dl
= 0.081 sin(2ϕ) × dl mgal/100m (A.5)
101
102 A. Correction de latitude
B Obtenir la latitude géocentrique
par rapport à la latitude géographique
Considérons l’équation de l’ellipse d’axes a et b :
x2 y2
+ = 1. (B.1)
a2 b2
En choisissant les changements de variable
x = r cos ϕ et y = r sin ϕ (B.2)
où ϕ est la latitude géocentrique, l’équation de l’ellipse peut s’écrire
r2 cos2 ϕ r2 sin2 ϕ
+ = 1. (B.3)
a2 b2
On peut en déduire que

a2 b2
r2 = (B.4)
b2 cos2 ϕ + a2 sin2 ϕ
et donc,
ab cos ϕ
x= q . (B.5)
b2 cos2 ϕ + a2 sin2 ϕ
De même, on a
bp 2
y= a − x2 (B.6)
a
Soit

dx b 1 −1/2 b x
= (−2x ) a2 − x2 =− √ (B.7)
dy a 2 a a2 − x 2
103
104 B. Obtenir la latitude géocentrique par rapport à la latitude géographique
Par conséquent
  vu   2  −1

dx b ab cos ϕ ab cos ϕ
u
= −  t a2 −  q
q  u  
dy a

b2 cos2 ϕ + a2 sin2 ϕ b2 cos2 ϕ + a2 sin2 ϕ
  v  −1
a4 sin2 ϕ

b  ab cos ϕ
u
= −  u 
q t q
a

2 2 2 2 2 2 2 2 2
b cos ϕ + a sin ϕ ( b cos ϕ + a sin ϕ)
  q 

b  ab cos ϕ b2 cos2 ϕ + a2 sin2 ϕ
= −  
a2 sin ϕ
q
a b2 cos2 ϕ + a2 sin2 ϕ
b2 cos ϕ
= − (B.8)
a2 sin ϕ
Il est donc maintenant possible de calculer ψ, qui est l’angle entre le plan équatorial et la
parallèle au plan de l’ellipse :
2
dx b cos ϕ
ψ = − arctan − = − arctan − 2 (B.9)
dy a sin ϕ
Puisque la latitude géographique θ vaut par définition π/2 − ψ, alors
b2 cos ϕ

θ = π/2 + arctan − 2 (B.10)
a sin ϕ
Ainsi
b2 cos(θ − π/2)

ϕ = arctan − 2 (B.11)
a sin(θ − π/2)
Comme cos(θ − π/2) = sin θ et sin(θ − π/2) = − cos θ, alors

2
b
ϕ = arctan 2 tan θ (B.12)
a
avec a = re = 6378.139 km et b = r p = 6356.754 km.

C Système de coordonnées et
Théorèmes fondamentaux
C.1 Coordonnées cartésiennes

– Élément de déplacement :
dl = dxi + dyj + dzk (C.1)
– Élément de surface (plan parallèle à x0y seulement) :
ds = dxdy (C.2)
– Élément de volume :
dv = dxdydz (C.3)
– Gradient de φ :
∂φ ∂φ ∂φ
∇φ = i+ j+ k (C.4)
∂x ∂y ∂z
– Divergence de V :
∂Vx ∂Vy ∂Vz
∇V = + + (C.5)
∂x ∂y ∂z
– Rotationel de V :

∂Vz ∂Vy ∂Vx ∂Vz ∂Vy ∂Vx
∇×V = − i+ − j+ − k (C.6)
∂y ∂z ∂z ∂x ∂x ∂y
– Laplacien de φ :
∂2 φ ∂2 φ ∂2 φ
∇2 φ = + 2+ 2 (C.7)
∂x2 ∂y ∂z
C.2 Coordonnées cylindriques

– Équations de transformation : 
 x = r cos θ
y = r sin θ (C.8)
z = z

105
106 C. Système de coordonnées et Théorèmes fondamentaux
dl = drar + rdθa` + dzk (C.9)
– Élément de surface (cylindre d’axe Oz seulement) :
ds = rdθdz (C.10)
dv = rdrdθdz (C.11)
∂φ 1 ∂φ ∂φ
∇φ = ar + a` + k (C.12)
∂r r ∂θ ∂z
1 ∂ 1 ∂Vθ ∂Vz
∇V = (rVr ) + + (C.13)
r ∂r r ∂θ ∂z

1 ∂Vz ∂V ∂Vr ∂Vz 1 ∂ ∂Vr
∇×V = − θ ar + − a` + (rVθ ) − k (C.14)
r ∂θ ∂z ∂z ∂r r ∂r ∂θ
1 ∂2 φ ∂2 φ

2 1 ∂ ∂φ
∇ φ= r + 2 2 + 2 (C.15)
r ∂r ∂r r ∂θ ∂z
C.3 Coordonnées sphériques

– Équations de transformation :

 x = r sin θ cos ϕ
y = r sin θ sin ϕ (C.16)
z = r cos θ

dl = drar + rdθa` + r sin θdϕa’ (C.17)
– Élément de surface (sphère de centre O seulement) :
ds = r2 sin θdθdϕ (C.18)
dv = r2 sin θdrdθdϕ (C.19)
∂φ 1 ∂φ 1 ∂φ
∇φ = ar + a` + a’ (C.20)
∂r r ∂θ r sin θ ∂ϕ
1 ∂ 2 1 ∂ 1 ∂Vϕ
∇V = 2
(r Vr ) + (sin θVθ ) + (C.21)
r ∂r r sin θ ∂θ r sin θ ∂ϕ
C.4 Théorèmes fondamentaux 107

1 ∂ ∂Vθ
∇×V = (sin θVϕ ) − ar +
r sin θ ∂θ ∂ϕ

1 1 ∂Vr ∂
− (rVϕ ) a` +
r sin θ ∂ϕ ∂r

1 ∂ ∂Vr
(rVθ ) − aϕ (C.22)
r ∂r ∂θ
∂2 φ

1 ∂ 2 ∂φ 1 ∂ ∂φ 1
∇2 φ = 2
r + 2
sin θ + 2 2 (C.23)
r ∂r ∂r r sin θ ∂θ ∂θ r sin θ ∂ϕ2
C.4 Théorèmes fondamentaux

– Théorème de Gauss : I ZZZ
U · ds = (∇ · U)dv (C.24)
S V
– Théorème de Stokes : I ZZ
U · ds = (∇ × U) · ds (C.25)
C S
– Première identité de Green :
I ZZZ
(∇ ϕ · ∇ψ) + ( ϕ∇2 ψ) dv

( ϕ∇ψ) · ds = (C.26)
S V
– Deuxième identité de Green :

I ZZZ
ϕ∇2 ψ − ψ∇2 ϕ dv

[ ϕ∇ψ − ψ∇ ϕ] = (C.27)
S V
108 C. Système de coordonnées et Théorèmes fondamentaux
D Exercices
D.1 Correction de données gravimétriques
Les données ci-jointes furent obtenues au-dessus d’un gisement de sulfures massifs près
de Noranda. A partir des données du tableau A on obtient la dérive de l’appareil à chaque
station de base en comparant la lecture corrigée établie dans la partie A et la lecture obtenue
lors de mesures subséquentes. La dérive aux autres stations est obtenue en interpolant entre
les valeurs de la dérive aux stations de base.
[A]
Corr. Lecture
Station Heure Lecture diurne corrigée
(div.) (div.) (div.)
1 +00 E 13h41 34.8 0.0 34.8
6 +00 E 13h53 39.3 -0.2 39.1
1 +00 E 14h02 35.1 -0.3 34.8
6 +00 E 14h16 39.8 -0.5 39.3
39.8 -0.6 39.2
12+00 E 13h31 22.6 -0. 8 21.8
6 +00 E 14h42 40.1 -0.9 39.2
12+00 E 14h50 23.2 -1.0 22.2
Complétez le tableau B des données en corrigeant les lectures de l’erreur de dérive. Cal-
culez les valeurs relatives de la gravité observée en multipliant les lectures corrigées par la
constante de l’appareil (K = 0.1074 mgals/div.).
109
110 D. Exercices
[B]
Corr. Lecture
Station Heure Lecture diurne corrigée ∆gobs
(div.) (div.) (div.) (mgal)
12+00 E 15h50 23.7 -1.7 22.0 2.3
11+00 E 15h53 24.3
10+00 E 15h57 27.2
9 +00 E 15h59 36.5
8 +00 E 16h03 49.4
7 +00 E 16h05 45.9
6 +00 E 16h10 41.2 -2.0 39.2
5 +00 E 16h13 40.0
4 +00 E 16h15 39.7
3 +00 E 16h17 38.5
2 +00 E 16h20 38.6 -2.2 34.8
1 +00 E 16h22 37.0
On peut maintenant calculer l’anomalie Bouguer en appliquant la correction appropriée. Si

on suppose que le terrain possède une densité moyenne de 2.70 g/cm3 . La correction combinée
d’altitude et de Bouguer est égale à 0.0594 mgal/pi. Complétez le tableau C en appliquant cette
correction. Comme le profil est orienté E-O, la correction de latitude n’est pas nécessaire. Étant
donné que le terrain de prospection est peu accidenté, on peut aussi négliger la correction
topographique.
[C]
Corr. Anom. Anom.
Station ∆g Elev. Bouger Bouger résiduelle
(mgal) (pi) (mgal) (mgal) (mgal)
1 +00 59.3
2 +00 60.7
3 +00 60.2
4 +00 59.3
5 +00 56.5
6 +00 52.8
7 +00 50.5
8 +00 46.4
9 +00 34.0
10+00 21.5
11+00 17.3
12+00 15.0
Dessinez le profil de l’anomalie Bouguer. Séparez graphiquement l’anomalie résiduelle de

l’anomalie régionale. Dessinez le profil de l’anomalie résiduelle et estimez la profondeur du
gisement.
D.2 Corrections gravimétriques 111
D.2 Corrections gravimétriques
Faire les corrections de dérive, latitude, altitude et plateau des données présentées dans le
tableau. Reportez dans le tableau la valeur de la correction et non pas la valeur de l’anomalie
corrigée. Calculez ensuite l’anomalie de Bouguer pour chacune des stations.
La latitude de la station 0+00 est 48˚45’10” et le profil est orienté suivant la direction N30E.
On utilisera une densité de ρ B = 2.67 g/cm3 .
Station alt. ∆gmes heure ∆t ∆h ∆gder ∆glat ∆galt ∆g pla ∆gB

(m) (mgal) (hhmm) (mn) (m) (mgal) (mgal) (mgal) (mgal) (mgal)
2+50 N 100.64 106.89 10h30
2+25 N 100.46 106.99 10h35
2+00 N 102.01 106.80 10h39
1+75 N 100.37 107.23 10h46
1+50 N 99.57 107.47 10h54
1+25 N 101.52 107.18 10h58
1+00 N 104.76 106.53 11h03
0+75 N 103.52 106.90 11h14
0+50 N 104.48 106.78 11h31
0+25 N 105.06 106.60 11h35
0+00 104.79 106.65 11h40
0+25 S 104.01 106.82 11h47
0+50 S 103.77 106.95 11h50
0+75 S 103.79 107.05 11h59
1+00 S 104.38 107.10 12h07
1+25 S 104.73 107.19 12h10
1+50 S 105.35 107.16 12h15
1+75 S 106.48 107.11 12h18
2+00 S 107.71 107.05 12h28
2+25 S 108.46 107.19 12h39
2+50 N 100.64 107.65 12h59
D.3 Calcul modélisation
D.3.1 Collecteur d’égouts
A quelle profondeur maximale doit se trouver le toit un collecteur d’égouts de 4 m de rayon

si l’on veut le détecter à l’aide d’un gravimètre de sensibilité de 0.05 mgal. On supposera le
tuyau horizontal. La roche encaissante a une densité de ρ = 2.2 g/cm3 .
D.3.2 Produits toxiques
On suppose qu’un baril contenant des produits toxiques est enfouis verticalement dans le
sol à une profondeur d’environ 4 m (profondeur du toit). Celui-ci aurait un rayon de 0.2 m
112 D. Exercices
pour une hauteur de 1 m. La densité de l’encaissant et du contenant du baril sont respective-

ment ρ=2.2 g/cm3 et ρ = 0.8 g/cm3 .
Pouriez vous le détecter à l’aide d’un gravimètre de sensibilité de 0.01 mgal ?
On rappel que l’anomalie maximale d’un cylindre vertical est
p q
∆gmax = 2πG∆ρ L + z + R − (z + L) + R
2 2 2 2
où z est profondeur du toit ; R le rayon du cylindre et L la hauteur du cylindre.

Notes de Cours-Gravimetrie PDF

Transféré par

Droits d'auteur :

Formats disponibles

Notes de Cours-Gravimetrie PDF

Transféré par

Informations du document

Titre original

Copyright

Formats disponibles

Partager ce document

Partager ou intégrer le document

Options de partage

Avez-vous trouvé ce document utile ?

Ce contenu est-il inapproprié ?

Droits d'auteur :

Formats disponibles

Notes de Cours-Gravimetrie PDF

Transféré par

Droits d'auteur :

Formats disponibles

GLQ2200 - Géophysique Appliquée I

Bernard Giroux, Michel Chouteau

Notes de cours - Gravimétrie

Laboratoire de géophysique appliquée

2 Les données gravimétriques 11

2.2.3 Résumé pour faire un levé gravimétrique . . . . . . . . . . . . . . . . . . 29

4 Signature des structures géologiques en gravimétrie 67

A Correction de latitude 101

B Obtenir la latitude géocentrique par rapport à la latitude géographique 103

C Système de coordonnées et Théorèmes fondamentaux 105

En géophysique appliquée, la gravimétrie sert à détecter les contrastes de densité du sous-

1.2 Principes de base

1.2.1 Lois de l’attraction universelle

1ère loi de Newton

2ème Loi de Newton

1 gal = 1 cm/s2 = 10−2 m/s2

1.2.2 Champ gravitationnel

1.2.3 Potentiel gravitationnel

L’équation du potentiel nous donne donc (m2 : masse unité) :

Ainsi, on peut définir le Potentiel Gravitationnel en un point du champ gravitationnel comme

où ρ est la densité (g/cm3 ) et dv l’élément de volume (cm3 ).

F IG . 1.3: Potentiel gravitationnel.

1.2.4 Équations du potentiel

Partant du théorème de Gauss, nous avons

Pour un élément de volume infinitésimal, l’intégrale tombe et nous avons

1.3 Une référence pour la terre

1.3.1 Un ellipsoı̈de de révolution : le sphéroı̈de

F IG . 1.6: Paramètres liés au calcul de l’accélération centrifuge : ω est la vitesse angulaire ; ω 2 r

1. L’effet de la rotation de la terre : plus on approche du pôle, plus l’accélération centrifuge

Sphéroïde Sphéroïde Continent

F IG . 1.7: Sphéroı̈de et géoı̈de

1.4 Densité des roches

TAB . 1.1: Densités des roches ignées (g/cm3 )

Type de roche Intervalle Moyenne Type de roche Intervalle Moyenne

TAB . 1.2: Densités des roches métamorphiques (g/cm3 )

Type de roche Intervalle Moyenne Type de roche Intervalle Moyenne

F IG . 1.8: Densités moyennes d’échantillons de surface et de carottes

TAB . 1.3: Densités des minéraux (g/cm3 )

Minéral Intervalle Moyenne Minéral Intervalle Moyenne

Type Intervalle Moyenne Type Intervalle Moyenne

TAB . 1.5: Densités des matériaux rocheux typiques (g/cm3 )

De façon générale en géophysique appliquée, on mesure la composante verticale de l’accélération

2.1 Corrections et références

2.1.1 Correction de dérive

F IG . 2.1: Exemple de l’influence des marées sur la gravité.

F IG . 2.2: Illustration d’une boucle de mesure simple.

vcor (t) = vlu (t) − ∆0d (t − t1 ). (2.2)

Station Lecture Temps

Le taux de dérive est

1031.0 − 1032.1 1.1

F IG . 2.3: Exemple d’une boucle avec plusieurs références.

Le principe demeure le même si au lieu de boucler sur la station de départ, la dernière

dérive = différence réelle − différence observée

De la même manière qu’auparavant, la formule de correction est

Notons le terme supplémentaire à la fin de l’équation, m1 − v1 . Il a pour but de ramener les

Soit la ligne de base suivante :

et les mesures suivantes prises le lendemain :

Station Temps Lecture

La dérive est de 1.1 division (1032.0-1030.1-1028.7+1027.9), et le taux de correction est de

Puisque l’on a r h, alors