See discussions, stats, and author profiles for this publication at: [Link]

net/publication/347948671

LE GEOIDE ET LES SYSTEMES D'ALTITUDES

Preprint · December 2020

DOI: 10.13140/RG.2.2.22282.90560

CITATIONS READS

0 179

1 author:

Abdelmajid Ben Hadj Salem

ST2I, Tunisia
144 PUBLICATIONS 15 CITATIONS

SEE PROFILE

Some of the authors of this publication are also working on these related projects:

Modernization of the Cadaster of Cameroon View project

Resolution of Number's Conjectures View project

All content following this page was uploaded by Abdelmajid Ben Hadj Salem on 28 December 2020.

The user has requested enhancement of the downloaded file.

 LE GEOIDE ET LES SYSTEMES
D’ALTITUDES

Abdelmajid BEN HADJ SALEM

Ancien Ingénieur Général Géographe

à l’Office de la Topographie et du Cadastre

Abstract : In this paper, we give the different definitions of heights and

the corrected terms to take in consideration. An example is presented on
how to correct a line of the spirit leveling by introducing gravity observa-
tions.

Résumé : Dans cet article, nous donnons les différentes définitions des
altitudes et les termes corrigés à prendre en considération. Un exemple est
présenté sur la manière de corriger une ligne de nivellement de précision
en introduisant des observations de gravité.

Décembre 2020

Version 2.

abenhadjsalem@[Link]
Table des matières 1

Table des matières

1 Le Champ de Pesanteur . . . . . . . . . . . . . . . . . . . . . 2
1.1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2
1.2 Le géoïde . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3
2 Systèmes d’Altitudes et Déviation de la Verticale . . . 4
2.1 Les Altitudes Dynamiques . . . . . . . . . . . . . . . . . . . . . . 5
2.2 Les Altitudes Orthométriques . . . . . . . . . . . . . . . . . . . . 5
2.3 Les Altitudes Normales ou de Molodensky . . . . . . . . . . . . . 6
3 Références . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6

A Annexes . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7
A.1 Le Champ du Potentiel . . . . . . . . . . . . . . . . . . . . . . . 7
A.2 Gradient . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7
A.3 Le Champ Réel ou Champ du Potentiel de la Pesanteur . . . . . 8
A.4 Calculs de corrections des altitudes . . . . . . . . . . . . . . . . . 9
Liste des Figures12Liste des Tables13
2 1 Le Champ de Pesanteur

LE GEOIDE ET LES SYSTEMES D’ALTITUDES

Abdelmajid BEN HADJ SALEM, Ancien Ing.

Général
A l’Office de la Topographie et de la Cartographie (OTC)

Résumé

Ce papier représente les notes rédigées à partir de l’action de formation donnée à

l’OTC, du 11 au 20 juin 1984, sur le thème " Les Satellites Doppler ".

1 Le Champ de Pesanteur

1.1 Introduction
Le phénomène fondamental qui régit la forme de la Terre est la pesanteur.
Elle est le résultat de l’attraction newtonienne du corps terrestre et de la force
centrifuge due à la rotation de la Terre autour de son axe. Voyons ce-ci en détail.

La force newtonienne, pour la masse ponctuelle est écrite sous la forme :

Mm
F =G (1)
r2
où :
- G la constante de gravitation universelle,
- M, m les masses des corps,
- r la distance des masses ponctuelles.

La force centrifuge est donnée par :

f = mω 2 rcosϕ (2)

où ω est la vitesse angulaire de rotation de la Terre. On peut calculer la com-

posante verticale de la force centrifuge avec la formule :

fv = f cosϕ = mω 2 rcos2 ϕ (3)

Sous l’action de la pesanteur, un fil à plomb au repos se met en équilibre suivant

la verticale physique : c’est la définition même de la verticale en un point.

Si l’on considérait un ensemble de fils à plomb disposés à la suite des uns des
autres, ils dessinaient dans l’espace une courbe gauche que l’on appelle ligne
1.2 Le géoïde 3

de force du champs de pesanteur.

La courbure et la torsion de cette courbe sont très petites de l’ordre d’une

fraction de seconde par kilomètre, mais elles ne sont pas tout à fait négligeables
dans certains cas.

Nous pouvons exprimer la pesanteur par les expressions suivantes :

∂W
g=− (4)
∂r
g = −gradW (5)

La différence entre la force newtonienne et la force centrifuge c’est la force de

la pesanteur :
P = F − fv (6)
On appelle l’accélération de la pesanteur ou la pesanteur :
dm
Z
g=k − ω 2 rcosϕ (7)
T r2
Le déplacement dr de la masse homogène sphérique m dans le champ de la
pesanteur d’une masse M produit le travail élémentaire :

dL = [Link] (8)

Le travail total d’une force de masse m de la distance r à l’infini est donné par
l’intégrale : Z +∞
L= P dR (9)
r
Divisée par la masse m, nous obtenons la valeur du potentiel de la pesanteur :
Z +∞
1
W = P dr (10)
m r

Dans le cas d’une masse homogène sphérique m, nous obtenons :

m ω2 r2
W =k + (11)
r 2
L’intégration de cette formule n’est pas possible car nous ne connaissons pas la
distribution de la masse dans la Terre. La distance entre deux surfaces équipo-
tentielle est donnée par :
dW
dR = − (12)
g

1.2 Le géoïde
La surface des liquides et des fluides se met en équilibre perpendiculairement
à la verticale. Si on considérait un ensemble fluide recouvrant toute la Terre, il
définissait donc une surface de niveau de la pesanteur.
4 2 Systèmes d’Altitudes et Déviation de la Verticale

La surface moyenne du niveau des mers, abstraction faite des marées et

corrigée des variations, définissait une surface de niveau unique pour le monde
entier. On peut définir d’autres surfaces de niveau de proche en proche à par-
tir d’un point quelconque pris comme origine sur une verticale donnée, par la
condition que cette surface soit en tout point perpendiculaire à toutes les autres
verticales.

Ces surfaces de niveau successives que l’on peut numéroter sont des surfaces
fermées qui s’enveloppent les unes les autres. Par exemple, les surfaces de même
pression atmosphérique (isobares) sont théoriquement des surfaces de niveau.

Comme par définition, une surface de niveau est normale aux lignes de force
et que ces lignes sont des courbes gauches, les surfaces de niveau ne sont pas
parallèles entre elles, c’est-à-dire que la distance de deux surfaces de niveau
n’est pas constante, ce qui reste constant c’est le travail qu’il faut accomplir
contre la pesanteur pour déplacer une masse d’un point donné de ces surfaces
à un point quelconque d’une autre de celles-ci.

Définition 1.1. On appelle géoïde la surface de niveau qui coincide avec la

surface moyenne des mers et qui se prolonge sous les continents par la condition
d’y rester normale à toutes les lignes de force.

On peut dans ces conditions considérer que la Terre est consituée par le
géoïde, surmontée du relief dont l’altitude au dessus du niveau moyen de la
mer sera par définition égale à la distance qui le sépare du géoïde. L’expérience
prouve que le géoïde s’écarte très peu d’un ellipsoïde de révolution parceque
le géoïde a une forme mathématique très compliquée, alors nous utilisons en
géodésie la surface mathématique à savoir l’ellipsoïde de révolution.

2 Systèmes d’Altitudes et Déviation de la Verticale

Nous avons vu précédemment que la distance de deux surfaces équipoten-

tielles est donnée : par
dW
dR = −
g
c’est-à-dire que la différence de potentiel de deux points A et B s’exprime par
la formule suivante :
∆WAB = −[Link] (13)
Cette différence peut être évaluée pour une masse unité par :
B
X
∆WAB = [Link] (14)
A

où g est l’accélération de la pesanteur mesurée en chaque point du trajet et dh

l’élément vertical qui précisemment se mesure à l’aide du niveau.
2.1 Les Altitudes Dynamiques 5

2.1 Les Altitudes Dynamiques

Définition 2.1. L’altitude dynamique dérive des côtes géopotentielles par la

formule : Z B
d W0 − WB 1
HB = = [Link] (15)
γ0 (45◦ ) γ0 (45◦ ) 0
où γ0 (45°) est la pesanteur normale à l’altitude zéro sous la latitude géodésique
ϕ = 45°.

Pour la différence des altitudes des points A et B, nous obtenons :

g − γ0 (45◦ )
Z B Z B
g

d d
HB − HA = dh = 1+ dh
A γ0 (45◦ ) A γ0 (45◦ )
g − γ0 (45◦ )
Z B Z B
= dh + dh (16)
A A γ0 (45◦ )

En pratique, on utilise la formule suivante :

B B
d d
X X gi − γ0 (45◦ )
HB − HA = ∆hi + ∆hi (17)
A A
γ0 (45◦ )

où le deuxième terme du membre à droite est la correction dynamique.

2.2 Les Altitudes Orthométriques

Le concept de l’altitude orthométrique cherche à revenir à la définition géo-

métrique litérale de l’altitude en tant que distance point-surface du géoïde comp-
tée selon la normale. Comme la différence des potentiels d’un point B0 sur le
géoïde et d’un point B sur le terrain vaut :
Z B
m
W0 − WB = [Link] = gB HB (18)
B0

m la pesanteur moyenne entre les points B et B .

avec gB 0

Définition 2.2. On écrit l’équation qui définit l’altitude orthométrique par :

Z B
o 1
HB = m [Link] (19)
gB 0

Pour la différence des altitudes des points A et B, nous obtenons :

Z B m
gB − gB m
gA − gA 1
Z B
HB − HA = dh + HB − HA − dg (20)
A g g g A

avec g la valeur moyenne de la pesanteur dans la région. En pratique, nous

6 3 Références

utilisons la formule suivante :

B B
X χ 2 2 1X
HB − HA = ∆i + (HA − HB )− Hi ∆gi (21)
A
R g A
σ
avec χ = (22)
σm
- σ : densité de la croûte terrestre,
- σm : densité moyenne de la Terre,
- R : rayon moyen de la Terre.

Quand on ne dispose pas des valeurs gravimétriques, on remplace g et ∆g

par les valeurs γ et ∆γ de la pesanteur normale et la différence des altitudes
sera :
B
X
HB − HA = ∆hi − 0.0053sin2ϕm .Hm ∆ϕ (23)
A

Avec ∆ϕ l’amplitude en latitude géodésique.

2.3 Les Altitudes Normales ou de Molodensky

Elles sont définies de la façon suivante :

Définition 2.3. En un point B connu par sa position géodésique dont on

connait par conséquent sa latitude géodésique, on appellera altitude normale
la quantité :
B
n 1 X
HB = m [Link] (24)
γB B
0

La différence des altitudes entre deux points A et B est :

B B
n n
X 1 1 X
HB − HA = ∆hi + (γ0B − γ0A )HAB + (g0 − γ0 )i ∆hi (25)
A
γ̄AB γ̄AB A

avec :
- γ̄AB = γ0ϕm − [Link] ,
- HAB l’altitude moyenne des points A et B,
- γ0B , γ0A la pesanteur normale respectivement de B et de A à l’altitude
zéro,
- (g0 − γ0 )0 l’anomalie à l’air libre,
- g0 = g + 0.3086H

3 Références

1. Helmut Moritz & and Ivan I. Mueller. Earth Rotation :Theory and Ob-
servation. Ungar Publishing Compagny. New York. 1988.
 7

2. Bernhard Hofmann-Wellenhof & Helmut Moritz. Physical Geodesy, ISBN-

10 3-211-23584-1, Springer Wien New-York, 422 pages, 157–172 (2005).

A Annexes

Soit le repère OXY Z tel que O soit le centre de gravité de la Terre et OZ

son axe de rotation. Le plan OXY contient le méridien de Greenwich.

Fig. 1: Le Repère 3D

A.1 Le Champ du Potentiel

On appelle champ du potentiel la fonction scalaire V définie par :

Gmm0
V = = V (X, Y, Z) (26)
r

A.2 Gradient
On appelle gradient d’une fonction scalaire U (X, Y, Z) le vecteur noté gradU
et de composantes :  ∂U


 ∂x



∂U
gradU = ∂y (27)





 ∂U
∂z
Si on pose :
r = OM = Xi + Y j + Zk (28)
Alors :
1 -r
 
grad = (29)
r r3
8 A Annexes

Calculons le gradient de la fonction scalaire donnée par l’équation (26) c’est-à-

dire le champ du potentiel, on a :

Gmm0 1 r
   
gradV = grad = Gmm0 grad = −Gmm0 (30)
r r r3
Remarquons si on pose :
r
n= (31)
r
On a n est un vecteur unitaire porté par OM et dans la direction OM . L’ex-
pression de la force F s’écrit :

Gmm0
F = −F n = − n (32)
r2
Gmm0
gradV = − n (33)
r2
D’où :
F = gradV (34)
On dit que la force F dérive du champ de potentiel V .

A.3 Le Champ Réel ou Champ du Potentiel de la Pesanteur

Un point M (X, Y, Z) de masse unité est soumis au potentiel V de gravita-
tion et au potentiel Φ de la force centrifuge due à la rotation de la terre.

L’expression de V est :

dm0
ZZZ
V =G (35)
T erre r
Malheureusement, cette expression n’est pas calculable car nous ignorons la
distribution des masses à l’intérieur de la Terre.

L’expression du potentiel Φ de la force centrifuge est donnée par :

1 2 
Φ= x + y2 ω2 (36)
2
où ω est la vitesse de la rotation de la Terre.

Définition A.1. On appelle Potentiel du champ réel W ou potentiel de la pe-

santeur la somme du potentiel V et Φ :

W =V +Φ (37)

Définition A.2. On appelle vecteur de gravité le vecteur g tel que :

g = gradW (38)

g mesure la gravité ou la pesanteur, a la dimension d’une accélération et

A.4 Calculs de corrections des altitudes 9

Fig. 2: Le Potentiel

exprimée en m/s2 (Unité Système International) ou en cm/s2 (1cm/s2 = 1gal

en hommage à Galilée 1 ). g mesure 978 gals à l’équateur et 983 gals aux pôles.

A.4 Calculs de corrections des altitudes

On a observé une ligne de nivellement de précision où on a mesuré les valeurs

de l’accélération terrestre dans les points caractéristiques du terrain :

Numéro du point Latitude 49◦ Altitude approximative en m g en mgal

A 20.30 739.49 980 816.2
1 19.30 764.92 980 809.9
2 18.8 0 823.10 980 798.8
3 18.60 798.52 980 802.9
4 17.9 0 827.58 980 797.2
5 17.70 914.51 980 782.6
6 17.4 0 912.91 980 762.0
B 17.20 915.01 980 778.5

Tab. 1: Tableau des Calculs 1

En utilisant les formules des corrections (17-23-25), on a le tableau suivant :

1. Galileo Galilei, dit Galilée : Astronome et physicien italien (Pise 1564-Arcetri 1642).
10 A Annexes

◦
Numéro du point ∆hi gi gi − γ045 Hi ∆gi
A
+25.43 980 813.0 197.7 752.20 −6.3
1
+58.18 980 804.4 188.5 794.01 −11.1
2
−24.58 980 800.8 184.9 810.80 4.1
3
29.06 980 800.0 184.1 813.05 −5.7
4
86.93 980 789.9 174.0 871.05 −14.6
5
−1.6 980 772.3 156.4 913.71 −20.6
6
+2.1 980 770.2 154.3 913.96 +16.5
B

Tab. 2: Tableau des Calculs 2

On a les données et les formules :

σ = 2.7 ; σm = 5.52 ; R = 6.37 × 106 m

B
X gi − γ0 (45◦ )
Correction Dynamique = C.D = ∆hi (39)
A
γ0 (45◦ )

B
χ 2 2 1X
Correction Orthométrique = C.O = (HA − HB )− Hi ∆gi (40)
R g A
Pas de données gravimétriques, Correction Orthométriques (sans g) =
[Link] = −0.0053sin2ϕm Hm ∆ϕ (41)
B
γ0B − γ0A X (g0 − γ0 )i
Correction Normale ou de Molodensky = HAB + ∆hi
γ̄AB A
γ̄AB
(42)
Avec :
γ̄AB = γ0 (ϕm ) − 0.15mgal HAB
g0 = g + 0.3086H
γ0 (ϕ) = 978030(1 + 0.005302sin2 ϕ − 0.000007sin2 2ϕ)
(HA + HB )
HAB =
2
A.4 Calculs de corrections des altitudes 11

Numéro du point g0 γ0 (ϕ) g0 − γ0 (g0 − γ0 )i ∆hi

A 981 044.4 981 007.2 +37.2
+38.9 +25.43
1 981 046.4 981 005.8 +40.6
+44.2 +58.18
2 981 052.8 981 005.0 +47.8
+46.2 −24.58
3 981 049.3 981 004.7 +44.6
+46.8 +29.06
4 981 052.6 981 003.6 +49.0
+55.2 +86.93
5 981 064.8 981 003.4 +61.4
+51.1 −1.06
6 981 043.7 981 002.9 +40.8
+49.6 +2.10
B 981 060.9 981 002.6 +58.3

Tab. 3: Tableau des Calculs 3

D’où les valeurs numériques :

C.D = +32.6 mm C.O = +19.2 mm [Link] = 3.9 mm

γAB = γ0 (ϕm ) − 0.15HAB = 981005 − 124 = 980881 mgals

C.N = +12.7 mm
12 Table des figures

Table des figures

1 Le Repère 3D . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7
2 Le Potentiel . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 9
 Liste des tableaux 13

Liste des tableaux

1 Tableau des Calculs 1 . . . . . . . . . . . . . . . . . . . . . . . . 9

2 Tableau des Calculs 2 . . . . . . . . . . . . . . . . . . . . . . . . 10
3 Tableau des Calculs 3 . . . . . . . . . . . . . . . . . . . . . . . . 11

View publication stats