ddd

Hydrologie et Hydrogologie

Responsable: V. MERRIEN-SOUKATCHOFF

Anne Universitaire 2008/2009

AVERTISSEMENT

Ce polycopi est un document de travail des cours de lEcole Nationale Suprieure des Mines
de Nancy.
Il souffre d'imperfections et toutes remarques constructives concernant ce texte seront les
bienvenues. Vous pouvez les envoyer [email protected].
Il est parfois incomplet il existe la bibliothque de l'ENSMN bon nombre d'ouvrages qui
complteront ce polycopi.
Certains dessins et tableaux ont t repris d'ouvrages vous trouverez les rfrences en fin de
polycopi.

Hydrologie - Hydrogologie Septembre 2008 Vronique MERRIEN-SOUKATCHOFF

Hydrologie et Hydrogologie
A

HYDROLOGIE DE SURFACE........................................................................................ 7
A.1

Le cycle de l'eau et bilans .................................................................................................... 7

A.2

Le Bassin versant ............................................................................................................... 10

A.2.1 Notion de bassin versant ............................................................................................................... 10

A.2.1.2
Le bassin versant hydrogologique ..................................................................................... 11
A.2.2 Caractristiques morphomtriques ................................................................................................ 12
A.2.2.1
En plan................................................................................................................................. 12
A.2.2.1.1 Surface ........................................................................................................................... 12
A.2.2.1.2 Caractristiques de longueur .......................................................................................... 13
A.2.2.1.3 Caractristiques de forme............................................................................................... 14
A.2.2.2
Caractristiques des altitudes (hypsomtrie) ....................................................................... 15
A.2.2.3
Indices de pente ................................................................................................................... 16
A.2.2.3.1 Pente moyenne ............................................................................................................... 16
A.2.2.3.2 Indice de pente de Roche ............................................................................................... 17
A.2.2.3.3 Indice de pente globale .................................................................................................. 17
A.2.2.3.4 Dnivele spcifique ...................................................................................................... 18
A.2.3 Caractristiques du rseau hydrographique................................................................................... 18
A.2.3.1
Hirarchisation du rseau .................................................................................................... 19
A.2.3.2
Lois de Horton..................................................................................................................... 20
A.2.3.3
Autres caractristiques......................................................................................................... 21
A.2.3.4
Profil en long ....................................................................................................................... 22
A.2.4 Caractristiques gologiques......................................................................................................... 22
A.2.5 Le couvert vgtal ......................................................................................................................... 23
A.2.6 Caractristiques glaciologiques..................................................................................................... 23

A.3
A.3.1
A.3.2
A.3.3
A.3.4

A.4

La Pluie ............................................................................................................................... 24
Mesure des prcipitations.............................................................................................................. 24
Synthse des mesures de prcipitations ........................................................................................ 24
Calcul de la pluie moyenne sur un bassin versant......................................................................... 24
Statistiques sur les donnes pluviomtriques ................................................................................ 25

Le retour de l'eau l'atmosphre : vaporation et vapotranspiration........................ 26

A.4.1 Gnralits. Dfinition du pouvoir vaporant. .............................................................................. 26

A.4.1.1
Les formes du retour de l'eau l'atmosphre....................................................................... 26
A.4.1.2
Pouvoir vaporant de l'atmosphre...................................................................................... 26
A.4.1.3
Pouvoir vaporant des surfaces d'eau libre.......................................................................... 27
A.4.2 Mesures de paramtres physiques conditionnant l'vaporation. ................................................... 28
A.4.2.1
Tempratures ....................................................................................................................... 28
A.4.2.2
Humidit de l'air .................................................................................................................. 28
A.4.2.2.1 Psychromtre.................................................................................................................. 28
A.4.2.2.2 Hygromtres................................................................................................................... 28
A.4.2.3
Pression................................................................................................................................ 28
A.4.2.4
Rayonnement solaire ........................................................................................................... 28
A.4.2.5
Vent ..................................................................................................................................... 29
A.4.3 Mesure de l'vaporation ................................................................................................................ 29
A.4.3.1
Mesures de l'vaporation partir d'une surface libre : les bacs ........................................... 29
A.4.3.1.1 Les bacs placs au-dessus du niveau du sol ................................................................... 29
A.4.3.1.2 Les bacs enterrs ............................................................................................................ 29
A.4.3.1.3 Les bacs flottants :.......................................................................................................... 30
A.4.3.2
A partir des surfaces poreuses (les atmomtres).................................................................. 30
A.4.3.3
Formules empiriques d'estimation du pouvoir vaporant .................................................... 30
A.4.3.3.1 Formule de Lugeon ........................................................................................................ 31
A.4.3.3.2 Formule de Meyer .......................................................................................................... 31
A.4.3.4
Ordre de grandeur du pouvoir vaporant............................................................................. 31

Hydrologie - Hydrogologie Septembre 2008 Vronique MERRIEN-SOUKATCHOFF

A.4.4 Mesures et estimation de l'vapotranspiration relle et potentielle ............................................... 32

A.4.4.1
Notion d'vapotranspiration relle et potentielle ................................................................. 32
A.4.4.2
Mesures directes .................................................................................................................. 32
A.4.4.2.1 Les cases lysimtriques (mesure de ETR)...................................................................... 32
A.4.4.2.2 Les parcelles d'essai ....................................................................................................... 32
A.4.4.3
Estimation de l'vapotranspiration....................................................................................... 33
A.4.4.3.1 Evapotranspiration potentielle ....................................................................................... 33
A.4.4.3.2 Evapotranspiration relle ............................................................................................... 33
A.4.4.4
Bilan au niveau d'un bassin versant..................................................................................... 33

A.5

Les coulements de surface ............................................................................................... 35

A.5.1 mesure des dbits .......................................................................................................................... 35

A.5.2 Estimation des dbits de frquence rare ........................................................................................ 35
A.5.2.1
A partir des statistiques sur les dbits.................................................................................. 35
A.5.2.2
A partir de modles pluie-dbit ........................................................................................... 35

Annexe : Quelques rappels de statistiques.................................................................................... 36

a - Rappels sur la notion de population ....................................................................................................... 36
b - Rappels sur la notion d'chantillon. ....................................................................................................... 36
c - Ajustement ............................................................................................................................................. 37
d - Danger d'apparition ................................................................................................................................ 37

Bibliographie Hydrosystmes ........................................................................................................ 39

Quelques revues............................................................................................................................... 39
Quelques sites web HYDRO & Environnement........................................................................... 40

B HYDROGEOLOGIE ET HYDRODYNAMIQUE : Le milieu naturel, description et

fonctionnement. Nappe, types et circulations ........................................................................ 41
B.1

Introduction - gnralits .................................................................................................. 41

B.2

L'eau dans le rservoir : la porosit ................................................................................. 41

B.2.1 Caractrisation : la porosit........................................................................................................... 41

B.2.2 Diffrents types de classification de la porosit............................................................................ 42
B.2.3 Classification de la porosit par rapport la taille des pores ........................................................ 42
B.2.4 Classification par rapport l'origine ............................................................................................. 43
B.2.5 Classification morphologique ....................................................................................................... 43
B.2.5.1
La porosit d'interstices (intergranulaire) : les pores ........................................................... 43
B.2.5.1.1 La porosit d'interstices simple ou nette......................................................................... 43
B.2.5.1.2 La porosit d'interstices restreinte .................................................................................. 43
B.2.5.1.3 La porosit d'interstices rduite ...................................................................................... 43
B.2.5.2
La porosit de fissure........................................................................................................... 44
B.2.5.2.1 Porosit de joints ............................................................................................................ 44
B.2.5.2.2 Porosit de diaclases....................................................................................................... 44
B.2.5.2.3 Porosit de failles ........................................................................................................... 44
B.2.5.2.4 Porosit de schistosit .................................................................................................... 44
B.2.5.2.5 Porosit de retrait............................................................................................................ 44
B.2.6 Ordre de grandeur de la porosit................................................................................................... 44
B.2.7 Mobilit de l'eau dans le sol : eau lie, eau libre, l'gouttage des roches...................................... 45
B.2.7.1
Eau libre, eau lie ................................................................................................................ 45
B.2.7.2
L'gouttage des roches......................................................................................................... 46
B.2.7.3
Relations air-eau pour diffrents degrs de saturation......................................................... 47
B.2.7.4
Pression capillaire................................................................................................................ 47
B.2.8 Profil hydrique des sols................................................................................................................. 48

B.3

Les mouvements de l'eau dans les roches : permabilit................................................ 48

B.3.1 Rappel sur la charge hydraulique .................................................................................................. 48

B.3.1.1
Charge hydraulique.............................................................................................................. 48
B.3.1.2
Cas des sols.......................................................................................................................... 49
B.3.1.2.1 Charge Hydraulique ....................................................................................................... 49
B.3.1.2.2 Notion de hauteur pizomtrique ................................................................................... 49

Hydrologie - Hydrogologie Septembre 2008 Vronique MERRIEN-SOUKATCHOFF

B.3.2 Exprience de Darcy ..................................................................................................................... 50

B.3.3 Exprience de Reynolds (pour mmoire)...................................................................................... 52
B.3.4 Ecoulement dans les roches stratifies .......................................................................................... 53
B.3.4.1
Permabilit horizontale : .................................................................................................... 53
B.3.4.2
Permabilit verticale. ......................................................................................................... 53
B.3.5 Mesures et estimation de la permabilit au laboratoire ............................................................... 54
B.3.5.1
Problmes poss par l'chantillonnage ................................................................................ 54
B.3.5.2
Estimation de la permabilit .............................................................................................. 54
B.3.5.2.1 Relation de Hazen : ........................................................................................................ 55
B.3.5.2.2 Relation de Casagrande :................................................................................................ 55
B.3.5.3
Permamtres....................................................................................................................... 55
B.3.5.3.1 Permamtre charge constante..................................................................................... 55
B.3.5.3.2 Permamtre charge variable....................................................................................... 56
B.3.6 Mesures in situ .............................................................................................................................. 56
B.3.7 Ordre de grandeur de la permabilit............................................................................................ 56
B.3.8 Gnralisation en 3 dimensions..................................................................................................... 57

B.4

Les nappes........................................................................................................................... 59

B.4.1 Conditions d'existence des nappes ................................................................................................ 59

B.4.1.1
Processus gnral................................................................................................................. 59
B.4.1.2
Facteurs d'existence d'une nappe ......................................................................................... 59
B.4.1.2.1 Cas d'une alimentation et d'une lithologie favorables : .................................................. 59
B.4.1.2.2 Cas d'une lithologie et d'une structure favorables : ........................................................ 59
B.4.1.2.3 Cas d'une alimentation et d'une structure favorables :.................................................... 60
B.4.2 Alimentation.................................................................................................................................. 60
B.4.3 Exutoires ....................................................................................................................................... 60
B.4.3.1
Les sources .......................................................................................................................... 60
B.4.3.1.1 Classification des sources............................................................................................... 60
B.4.3.1.2 Fonctionnement des sources........................................................................................... 64
B.4.3.2
Les exutoires cachs ............................................................................................................ 64
B.4.4 Classification des nappes .............................................................................................................. 64
B.4.4.1
Critres gologiques ............................................................................................................ 64
B.4.4.1.1 Nappes de terrains sdimentaires stratifis..................................................................... 64
B.4.4.1.2 Nappes de terrains sdimentaires mal (ou non) stratifis ............................................... 64
B.4.4.1.3 Nappes de terrains cristallins ou ruptifs........................................................................ 64
B.4.4.2
Critres hydrodynamiques ................................................................................................... 65
B.4.4.2.1 Nappe libre ..................................................................................................................... 65
B.4.4.2.2 Nappe captive ................................................................................................................. 65

B.5

Etude de quelques coulement dans les nappes............................................................... 67

B.5.1 Paramtres dont dpend l'coulement. Rseau d'coulement........................................................ 67

B.5.1.1
Transmissivit...................................................................................................................... 67
B.5.1.2
Coefficient d'emmagasinement............................................................................................ 67
B.5.1.3
Surface pizomtrique et isopizes...................................................................................... 68
B.5.1.4
Lignes de courants ............................................................................................................... 68
B.5.2 Effet des variations de la transmissivit ........................................................................................ 68
B.5.3 Ecoulement entre deux tranches.................................................................................................. 69
B.5.3.1
Cas d'une nappe libre........................................................................................................... 69
B.5.3.2
Cas d'une nappe captive....................................................................................................... 70
B.5.4 Ecoulement radial circulaire en rgime permanent ....................................................................... 71
B.5.4.1
Nappe captive ...................................................................................................................... 71
B.5.4.2
Nappe libre .......................................................................................................................... 72

B.6

Cas Gnral : quation fondamentale de l'hydrodynamique ........................................ 74

B.6.1 Etablissement de l'quation de diffusivit ..................................................................................... 74

B.6.2 Proprits de l'quation de diffusivit ........................................................................................... 75
B.6.2.1
Unicit de la solution (Cf. Cours TCS11, Maths 1A et SG232 EDP))................................ 75
B.6.2.2
Principe de superposition..................................................................................................... 75

B.7
B.7.1

Solutions de l'quation de diffusivit lors de pompage dans un puits........................... 75

Equation de diffusivit en coordonnes radiales........................................................................... 76

Hydrologie - Hydrogologie Septembre 2008 Vronique MERRIEN-SOUKATCHOFF

B.7.2 Solution de Theis .......................................................................................................................... 76

B.7.3 Solution de Jacob .......................................................................................................................... 76
B.7.4 Pompage d'essai ............................................................................................................................ 77
B.7.4.1
Interprtation graphique de la mthode de Theis................................................................. 77
B.7.4.2
Interprtation graphique de la mthode de Jacob ................................................................ 78
B.7.5 Remonte de la nappe ................................................................................................................... 79
B.7.6 Mthode des paliers enchans...................................................................................................... 79
B.7.7 Quelques cas particuliers de nappe ............................................................................................... 79
B.7.7.1
Drainance............................................................................................................................. 80
B.7.7.1.1 Cas des aquifres superposs semi-permable : Schma de Hantush ............................ 80
B.7.7.1.2 Solution de Boulton........................................................................................................ 80
B.7.7.2
Ralimentation latrale ........................................................................................................ 80
B.7.7.2.1 Mthode des puits images .............................................................................................. 80
B.7.7.2.2 Application l'interprtation d'un pompage avec ralimentation latrale du puits ........ 82
B.7.8 Essais Lugeon et Lefranc .............................................................................................................. 82

B.8
B.8.1
B.8.2

Les essais spcifiques : ....................................................................................................... 82

Essais par chocs hydrauliques : Pulse test..................................................................................... 82
Slug test (dbit variable partir d'une injection unique)............................................................... 83

B.9
Mthodes analytiques de rsolution de l'quation de diffusivit : fonctions
harmoniques .................................................................................................................................... 83
B.10

Mthodes numriques de rsolution de l'quation de diffusivit................................... 83

B.10.1
Equations rsoudre ................................................................................................................ 84
B.10.2
Conditions aux limites :............................................................................................................ 85
B.10.3
Principe des diffrences finis : ................................................................................................. 85
B.10.3.1 En rgime permanent :......................................................................................................... 86
B.10.3.2 En rgime transitoire : ......................................................................................................... 86
B.10.4
Principe des lments finis : ..................................................................................................... 86

B.11

Ecoulement dans les milieux fracturs............................................................................. 88

B.11.1
Ecoulement dans une fracture .................................................................................................. 88
B.11.1.1 Diffrents rgimes d'coulement ......................................................................................... 88
B.11.1.2 Conductivit hydraulique, permabilit, transmissivit d'une fracture................................ 93
B.11.1.3 Coefficient d'emmagasinement............................................................................................ 94
B.11.2
Ecoulement dans un rseau de fractures................................................................................... 94
B.11.2.1 Analyse de la dimension de l'coulement ............................................................................ 94
B.11.2.2 Modle d'coulement radial gnralis de Barker ............................................................... 95
B.11.3
Milieu continu quivalent......................................................................................................... 95

B.12

Transport de soluts en milieu poreux ............................................................................. 95

B.12.1
Position du problme et mcanisme de transport ..................................................................... 95
B.12.2
La convection (ou advection) ................................................................................................... 97
B.12.3
La diffusion .............................................................................................................................. 98
B.12.3.1 La diffusion molculaire...................................................................................................... 98
B.12.3.2 La diffusion en milieu poreux ............................................................................................. 99
B.12.4
La dispersion cinmatique........................................................................................................ 99
B.12.5
Les phnomnes d'change avec le milieu solide................................................................... 102
B.12.6
Les processus de dgradation ................................................................................................. 102
B.12.7
Expression de l'ensemble des termes du transport.................................................................. 102

Proprits physiques de l'eau : symboles, valeurs et units ................................................ 103

Bibliographie Eaux Souterraines................................................................................................. 105
Quelques revues............................................................................................................................. 106
Quelques sites web HYDRO & Environnement......................................................................... 107

Hydrologie - Hydrogologie Septembre 2008 Vronique MERRIEN-SOUKATCHOFF

A HYDROLOGIE DE SURFACE
A.1 Le cycle de l'eau et bilans
L'alimentation de l'eau souterraine provient de l'infiltration la surface du sol d'une fraction
de l'eau provenant des prcipitations, mais seule une partie des prcipitations s'infiltrera.

LE CYCLE DE L'EAU
Transpiration

Prcipitation

Evaporation
zone non-sature

ruissellement
Transpiration

recharge de la nappe

Evaporation

nappe

Figure 1 : Le cycle de l'eau

Examinons ce qu'il va advenir de l'eau issue des prcipitations :
- une partie de cette eau va ruisseler et alimenter les cours d'eau superficiels (problmes
d'hydrologie de surface) ;
- une partie de l'eau va retourner l'atmosphre, sous forme vapeur : c'est l'vapotranspiration
qui est la somme de deux phnomnes :
* l'vaporation (phnomne physique) qui intervient la surface des lacs, des cours
d'eau, mais aussi sur le sol.
* la transpiration (phnomne biologique) qui est le fait de la couverture vgtale.
- enfin une partie de l'eau issue des prcipitations va s'infiltrer.
Pour pouvoir effectuer un bilan sur le cycle de l'eau, il va falloir se dfinir une surface de
bilan et une unit de temps de bilan.
- espace : notion de bassin versant

Hydrologie - Hydrogologie Septembre 2008 Vronique MERRIEN-SOUKATCHOFF

- temps : anne hydrologique, mais on peut galement faire des bilans un pas de temps plus
faible.
Pourquoi faire un bilan ?
Il est intressant de connatre, les dbits qui transitent, le % d'eau qui va s'infiltrer pour tous
les projets qui vont utiliser l'eau et galement pour prvoir, les risques de pnurie,
d'inondations
Les eaux de surface peuvent tre rparties en 2 ensembles:

les eaux courantes se dplaant par coulement sur le sol

les eaux non courantes animes de mouvements (lacs, mers et ocans.)

Volumes (106 km3)

Part %

Ocans

1320

97,2

Neiges et glaces

30

2,15

Eaux souterraines (- de 800 m)

0,31

Eaux souterraines (+ de 800 m)

0,31

Zone non sature

0,07

5 10-3

Lacs en eau douce

0,12

9 10 -3

Lacs en eaux sale

0,1

8 10-3

Rivires

0,001

10-4

Atmosphre

0,0013

10-4

Tableau 1 : estimation des volumes d'eau disponibles dans le monde

Vous trouverez des lments sur le cycle de l'eau et l'origine de l'eau sur Terre l'adresse
http://www.u-picardie.fr/beauchamp/

Hydrologie - Hydrogologie Septembre 2008 Vronique MERRIEN-SOUKATCHOFF

Figure 2 : Rpartition ingale de l'eau dans le monde d'aprs

http://www.sagep.fr/html/alimenter/fr_potable.shtml

Pour l'ensemble de cette partie du cours vous trouverez des complments intressants
l'adresse http://hydram.epfl.ch/e-drologie/

10

Hydrologie - Hydrogologie Septembre 2008 Vronique MERRIEN-SOUKATCHOFF

A.2 Le Bassin versant

A.2.1 Notion de bassin versant
Le bassin versant, en une section d'un cours d'eau, est dfini comme la surface draine par ce
cours d'eau et ses affluents en amont de la section.
Tout coulement prenant naissance l'intrieur de cette surface doit traverser la section
considre, appele exutoire, pour poursuivre son trajet vers l'aval.
A.2.1.1.1.1 Le bassin versant topographique
Ligne de crte

Si le sous-sol est
impermable,
le
cheminement de l'eau
ne sera dtermin que
par la topographie ; le
bassin versant sera
limit par les lignes de
Ligne de crtes et les lignes de
crte
plus grande pente.

Sens de lcoulement

Le bassin versant est

l'unit spatiale qui va
tre
utilise
pour
effectuer
un
bilan
hydrologique.
Exutoire
ligne de plus
grande pente

Figure 3 : bassin versant topographique

Hydrologie - Hydrogologie Septembre 2008 Vronique MERRIEN-SOUKATCHOFF

11

A.2.1.2 Le bassin versant hydrogologique

Bassin versant topographique
Bassin versant hydrogologique

ruissellement

Station de jaugeage

Figure 4 : bassin versant hydrogologique

Si la rgion est permable, une partie des eaux tombes l'intrieur du bassin versant
topographique s'infiltrent et sortent souterrainement du bassin. Inversement, on peut avoir
l'entre d'eaux souterraines dans le bassin versant. Pour dterminer les limites du bassin
versant hydrogologique, il est donc ncessaire de prendre en compte les limites gologiques.
Une des limites du bassin versant hydrogologique est la ligne de partage des eaux. Cette
limite peut varier en fonction du niveau de l'eau dans la nappe

Ligne de partage des eaux

Ligne de partage des eaux

Zone synclinale

Zone synclinale

Si la pluie diminue, la nappe peut se vidanger

et dans ce cas la position de la ligne de
partage des eaux peut varier.
Figure 5 : variation des limites du bassin versant hydrogologique
Les limites d'un bassin versant sont donc variables en fonction du temps, ce qui entranera
quelques complications lorsque l'on voudra effectuer un bilan.
La diffrentiation entre bassin versant topographique et bassin versant hydrogologique est
valable pour des petits bassins versants ; quand la taille augmente :
- les apports et les pertes ont plus de chance de se compenser ;
- le dbit des cours d'eau augmente en fonction de la surface du bassin versant, par contre les

Hydrologie - Hydrogologie Septembre 2008 Vronique MERRIEN-SOUKATCHOFF

12

changes souterrains varient en fonction du primtre du bassin versant (changes aux

frontires). Comme la surface augmente plus rapidement que le primtre, les changes
souterrains diminuent en pourcentage par rapport aux dbits superficiels.
A.2.2 Caractristiques morphomtriques
Il est difficile de dcrire entirement un bassin versant. Pour diffrencier deux bassins
versants, il faudrait pouvoir comparer les cartes topographiques, gologiques et celles du
rseau hydrographique.
Pour faciliter la description, on dfinit un certain nombre de paramtres qui permettent de
caractriser la forme du bassin versant. Ces paramtres sont de 3 types qui caractrisent :
- la disposition en plan
- l'altitude
- l'indice de pente
A.2.2.1 En plan
A.2.2.1.1 Surface
Aire de la surface
Les limites topographiques et ventuellement hydrogologiques peuvent tre traces. En
planimtrant cette surface on obtient sa valeur qui est exprime en km2 .
Si le contour du bassin versant est schmatis
par un contour polygonal dfini par n points
de coordonnes xi et yi, l'aire de la surface A
peut tre calcule par la relation suivante :
A=

(x i x i 1 ) (y i + y i 1 )

i =1

n=0
n'
1

Figure 6
Remarque (pas trs importante, mais certains
aiment savoir d'o sortent les formules) : soit un
polygone form de n points numrots de 1 n
dans le sens anti-horaire (le point 0 = le point n).
si A1 est l'aire sous la courbe infrieure :
n'

A1 = ( xi xi 1 ) y i 1 + 1
i =1

n'

A1 =
i =1
n'

A1 =
i =1

( xi xi 1 )( yi yi 1 )

(x i x i 1 ) (2 y + y y )
i 1
i
i 1
2
(x i x i 1 ) ( y + y )
i 1
i
2

si A2 est l'aire sous la courbe suprieure :

n=0

n'
i -1

yi -1
yi

A1

x
xi -1

Figure 7

xi

Hydrologie - Hydrogologie Septembre 2008 Vronique MERRIEN-SOUKATCHOFF

A2 =
i =n'

(x i x i 1 ) ( y + y )
i 1
i
2

yi

Ici le signe est ngatif car les x sont dcroissants.

13

n=0

i -1

yi -1
A2

n'

x
xi

xi -1

Figure 8
n

A = A2 A1 =
i =n'

A=
i =1

(xi xi 1 ) ( y + y ) n' (xi xi 1 ) ( y + y ) donc

i 1
i
i 1
i
2
2
i =1

(x i 1 x i ) ( y + y )
i 1
i
2

Si les points sont numrots dans le sens horaire la formule devient :

n

A=
i =1

(x i xi 1 ) ( y + y )
i 1
i
2

On peut donc retenir que A =

( xi xi 1 ) ( yi + yi 1 )

i =1

Si on inverse les coordonnes

x et y le problme est
similaire.

A2
A1

Figure 9

Figure 10

A.2.2.1.2 Caractristiques de longueur

Il existe diffrentes caractristiques de longueur pouvant tre utilises.
- Primtre :
Pour le mesurer, on curvimtre le bassin versant. Souvent, les rsultats dpendent de l'chelle
de la carte (les dtails sont plus ou moins nombreux).
Si le contour du bassin versant est schmatis par un contour polygonal dfini par n points de
coordonnes xi et yi, le primtre P peut tre calcul par la relation suivante :
n

P=
i =1

( xi xi 1 )2 + ( yi yi 1 )2

14

Hydrologie - Hydrogologie Septembre 2008 Vronique MERRIEN-SOUKATCHOFF

- longueur du plus long talweg

Figure 11 : plus long Talweg

- Distance entre l'exutoire et le centre de

gravit

G
x

Figure 12 : distance exutoire centre de gravit

- L : plus grande longueur

l

- l : plus grande largeur (perpendiculaire

la plus grande longueur)

Figure 13 : plus grande longueur et largeur

A.2.2.1.3 Caractristiques de forme

A.2.2.1.3.1.1 coefficient K de Gravelius = coefficient de compacit

Kc =

P
2 A

0,28

P
A

Ce paramtre est suprieur ou gal 1 ; il est

proche de 1 si le bassin versant est proche
d'une forme circulaire.

Figure 14 : Kc proche de 1
A.2.2.1.3.1.2 Rapport de la plus grande longueur la plus grande largeur
A.2.2.1.3.1.3 Rectangle quivalent

Dimension du rectangle ayant mme primtre et mme surface

Hydrologie - Hydrogologie Septembre 2008 Vronique MERRIEN-SOUKATCHOFF

15

A.2.2.1.3.1.4 Courbe aire-distance

Cette courbe donne la surface S (en km2 ou en pourcentage du bassin versant) en fonction de
la distance hydraulique (point se trouvant une distance hydraulique suprieure d). La
distance hydraulique est la distance parcourue par une particule d'eau qui ruisselle d'un point
jusqu' l'exutoire.
courbe aire-distance
14 000

12 000

Distance

10 000
8 000

d50%

4 000
2 000
0%

20%

40% 50% 60%

80%

100%

Surface en pourcentage

Figure 15 : dtermination de la courbe aire-distance

Pour calculer la distance hydraulique, on recherche la ligne de plus grande pente par calculs
successifs (informatiss)
A.2.2.2 Caractristiques des altitudes (hypsomtrie)
700

h5%
600
500
Altitude

Nous nous intresserons la

dispersion des altitudes
plutt
qu'
l'altitude
moyenne. Pour cela on trace
la courbe hypsomtrique,
qui est la courbe de la
surface (en km2 ou en
pourcentage
du
bassin
versant) o les altitudes sont
suprieures une altitude h
donne.

400
300
200

h95%
0%

10%

20%

30%

40%

50%

60%

70%

80%

90% 100%

Surface en pourcentage

Figure 16 : courbe hypsomtrique

Deux mthodes permettent d'obtenir cette courbe :

-

on planimtre entre les courbes de niveau,

on chantillonne les altitudes selon un maillage et l'on considre que l'altitude au centre de
la maille correspond l'altitude de la maille.

A partir de cette courbe il est possible de dfinir la dnivele D par : D = h 5% h 95%

Hydrologie - Hydrogologie Septembre 2008 Vronique MERRIEN-SOUKATCHOFF

16

A.2.2.3 Indices de pente

A.2.2.3.1 Pente moyenne

j+1
La premire manire de
courbes de niveaux
caractriser les pentes est de
h
calculer la valeur moyenne I
j
pondre par les surfaces.
lj
hj+1
y
j-1
Soient
D
l'quidistance
D
(diffrence d'altitude) entre 2 h
j
courbes de niveau (sur les
limites du
cartes au 1/50 000, souvent D
dj
bassin versant
= 20 m), dj la largeur
moyenne de la bande j
comprise entre les lignes de
niveau j et j+1, nj la pente
x
moyenne de cette bande. La
pente moyenne sur cette
Figure 17 : dtermination de la pente moyenne
D
bande est donc : n j =
dj

L'aire de la surface comprise entre les courbes j et j+1 (aj ) est : a j = d j l j

La pente moyenne I est la moyenne pondre par les surfaces, des pentes entre les courbes de
niveau. Si A est la surface totale du bassin versant :

n a
I=
a
j

dj lj
=D

Si on note : L = l j la longueur de l'ensemble des courbes de niveau quidistantes de D :

I=

DL
A

Pour calculer la pente moyenne, il est donc ncessaire de connatre l'quidistance entre les
courbes de niveau, l'aire du bassin versant et la longueur des courbes de niveau. L'estimation
de cette expression simple est cependant laborieuse puisqu'il faut curvimtrer toutes les
courbes de niveau. Ceci explique que cet indice soit peu utilis dans la pratique.

Hydrologie - Hydrogologie Septembre 2008 Vronique MERRIEN-SOUKATCHOFF

17

A.2.2.3.2 Indice de pente de Roche

M. Roche a propos un indice de

pente plus facile calculer que
h
le prcdent partir de la courbe
hypsomtrique.
Ip
est
la
moyenne de la racine carre des
pentes, mesures sur le rectangle
quivalent, et pondre par les hj+1
hj
surfaces.

La pente moyenne nj sur la

D
bande j est : n j =
xj
xj

La surface de cette bande est :

aj = lxj

l
L

Figure 18 : Indice de pente de Roche

D'o l'expression de Ip : I p =

a n
a
j

x
j

l
A

D
xj

l D x j
Ll

D
xj
L

en posant j le pourcentage de la surface totale se trouvant entre h j+1 et h j :

j =
Ip =

xj l
Ll

xj
L

D
L j =
L

D
L

L'estimation de Ip est plus simple que celle de I puisque l'on travaille sur le rectangle
quivalent. Par ailleurs, la valeur de Ip est peu affecte par le choix de D (une dizaine de
classes suffit pour bien estimer Ip ).
A.2.2.3.3 Indice de pente globale
L'indice de Roche tant trop long valuer pour des tudes rapides, on a propos un indice
D
encore plus simple : la pente globale : I g =
L
D tant la dnivele : D = h 5% h 95% dfinie sur la courbe hypsomtrique ou mme
directement l'il sur la carte topographique.
L tant la longueur du rectangle quivalent.

18

Hydrologie - Hydrogologie Septembre 2008 Vronique MERRIEN-SOUKATCHOFF

R1

Relief trs faible

Ig < 0,002

R2

Relief faible

0,002 < Ig < 0,005

R3

Relief assez faible

0,005 < Ig < 0,01

R4

Relief modr

0,01 < Ig < 0,02

R5

Relief assez fort

0,02 < Ig < 0,05

R6

Relief fort

0,05 < Ig < 0,1

R7

Relief trs fort

0,1 < Ig

Tableau 2 : Classement ORSTOM du relief partir de l'indice de pente globale (pour

des bassins versants < 25 km2)

A.2.2.3.4 Dnivele spcifique

Ig diminue quand L augmente (donc en gnral quand la surface du bassin versant
augmente) ; il sera donc difficile de comparer deux bassins versants de taille diffrente. On
dfinit donc une grandeur qui ne prsente pas cet inconvnient : Ds : dnivele spcifique :
DS = Ig A
Ds est proportionnelle la longueur : D S =

D
l
Ll = D
L
L

Ds dpend de l'hypsomtrie D et de la forme du bassin.

R1

Relief trs faible

Ds < 10 m

R2

Relief faible

10 m < Ds < 25 m

R3

Relief assez faible

25 m < Ds < 50 m

R4

Relief modr

50 m < Ds < 100 m

R5

Relief assez fort

100 m < Ds < 250 m

R6

Relief fort

250 m < Ds < 500 m

R7

Relief trs fort

500 m < Ds

Tableau 3 : Classement ORSTOM du relief partir de la dnivele spcifique

(indpendant de la surface)
A.2.3 Caractristiques du rseau hydrographique

- Forme : par exemple, on veut sparer :

Hydrologie - Hydrogologie Septembre 2008 Vronique MERRIEN-SOUKATCHOFF

Figure 19 : Rseau en artes de poisson

19

Figure 20 : Rseau dendritique

- Densit du rseau

Figure 22 : mme rseau mais beaucoup

plus dense

Figure 21 : rseau faible densit

A.2.3.1 Hirarchisation du rseau

Il s'agit de numroter les cours d'eau en fonction de leur importance. Strahler propose de les
classer de la manire suivante :
-

les cours d'eau sans affluents sont d'ordre 1

au confluent de deux cours d'eau de mme ordre rsulte un cours d'eau d'ordre n+1

un cours d'eau recevant un affluent d'ordre infrieur garde son ordre ( par exemple si un
cours d'eau d'ordre 4 reoit un cours d'eau d'ordre 1 => toujours ordre 4)

Problme : dfinition de l'ordre 1 dpend de l'chelle laquelle on travaille. On peut

cependant tablir une correspondance entre l'chelle et l'ordre rel (rvl par photographie
arienne)
Ordre rel

Ordre lu sur la carte

Echelle de la carte

1/20 000

1/50 000

1/100 000

1/200 000

20

Hydrologie - Hydrogologie Septembre 2008 Vronique MERRIEN-SOUKATCHOFF

ordre 1
ordre 2
ordre 3
ordre 4

Figure 23 : hirarchisation du rseau hydrographique

Quand la hirarchisation du rseau hydrographique est tablie, il peut tre intressant de

caractriser son dveloppement, c'est--dire le nombre de cours d'eau et la longueur des cours
d'eau.
A.2.3.2 Lois de Horton

Les lois de Horton sont des "lois" empiriques qui relient le nombre, la longueur moyenne et
l'ordre des cours d'eau ;
Rapport de confluence : c'est le rapport du nombre de cours d'eau d'ordre n au rapport du
nombre de cours d'eau d'ordre n+1 :

R c(n) =

nombre de cours d' eau d' ordre n

nombre de cours d' eau d' ordre n + 1

On constate que, quel que soit n, ce rapport Rc est peu prs constant.
Rc va permettre de diffrencier des rseaux en arte de poisson pour lesquels Rc est important
et des rseaux dendritiques pour lesquels Rc est faible :

Hydrologie - Hydrogologie Septembre 2008 Vronique MERRIEN-SOUKATCHOFF

21

Figure 25 : Rc = 2

Figure 24 : Rc = 10

Si Rc est constant le nombre de cours d'eau d'ordre n suit une progression gomtrique.
Rapport de longueur : Rl (n) =

ln
l n 1

La dfinition est analogue la prcdente, mais avec la longueur moyenne des cours
d'eau ( ).
L encore on constate que les longueurs moyennes suivent aussi une progression gomtrique
et que Rl(n) est peu prs constant
Pour dterminer Rc et Rl on trace le nombre de cours d'eau d'ordre i (ou la longueur des cours
d'eau d'ordre i) en fonction de l'ordre sur du papier semi-log et on peut procder un
ajustement graphique.
A.2.3.3 Autres caractristiques

Il existe d'autres paramtres que Rl et Rc qui caractrisent le chevelu du rseau. On peut

citer :
- La densit de drainage : D d =

l
A

l : longueur totale des cours d'eau

A : surface du bassin versant
- Frquence des talweg d'ordre 1 F1 =

n1
A

Problme d'chelle encore

- Courbe aire-distance
Cette courbe a dj t cite prcdemment dans les caractristiques du bassin versant.
Pour pouvoir faire des comparaisons entre ces caractristiques il est important de toujours
prendre le mme support.

22

Hydrologie - Hydrogologie Septembre 2008 Vronique MERRIEN-SOUKATCHOFF

A.2.3.4 Profil en long

Ces profils sont tablis en portant sur un graphique, en abscisse les longueurs partir d'un
point de rfrence, en ordonne la cte (altitude) de la surface de l'eau ou parfois celle du
fond de la rivire. Le point de rfrence des longueurs est gnralement la source, mais le
profil peut tre aussi effectu sur un tronon de rivire.
Ces profils sont parfois tablis quand il existe
des problmes de navigation ou lors d'tudes
pour l'hydrolectricit.
tronon j

Les profils en long permettent d'estimer la

pente moyenne du cours d'eau ; cette dernire
est calcule en dcomposant le profil en long
en tronons de pente constante.

ij

ij : pente du tronon j
lj

Figure 26 : Profil en long

Cette estimation de pente moyenne est

gnralement utilise pour le calcul des temps
de concentration dans un bassin versant, or la
vitesse de transfert de l'eau varie comme
(I tant la pente).
En faisant une moyenne arithmtique ; les
deux profils ci-contre ont la mme pente
moyenne, alors qu'ils n'auront pas le mme
comportement.

1
2

Figure 27 : pente moyenne d'un profil en

long

On va donc rechercher la pente moyenne comme tant la pente qui provoquerait la mme
vitesse de propagation V. Comme V varie en fonction de la racine carre de la pente du bassin
versant I , le temps T varie comme 1 donc comme 1 . La pente moyenne est donc
V
I
n
lj
1
1
=
estime par :
I L j=1 i j
lj : longueur des tronons de pente constante
ij : pente du tronon j
A.2.4 Caractristiques gologiques
La gologie d'un bassin versant conditionne le rgime de ses cours d'eau :

en priode de crue, les dbits seront d'autant plus importants que le bassin sera
impermable,

par contre en priode de basses eaux, les dbits seront d'autant plus importants que les
nappes sont nombreuses (donc le sous-sol permable).

Hydrologie - Hydrogologie Septembre 2008 Vronique MERRIEN-SOUKATCHOFF

23

La gologie influe galement sur l'vapotranspiration :

-

par l'effet thermique d la couleur du sol,

par la vgtation qui varie en fonction de la nature du sol.

L'ORSTOM a propos une classification en 5 types, base sur la Permabilit :

Classe

Intitul

Exemple

P1

Zone
permable

aquifre Formations grseuses dont les exutoires

l'extrieur
du
drainant ou non drain (le terrain sont
qui constitue l'aquifre est important ; bassin versant.
l'aquifre n'alimente pas ou peu le
rseau hydrographique)

P2

Zone permable aquifre drain

Aquifre important alimentant une

partie du rseau hydrographique

Formations grseuses dont les sources

alimentent le rseau

P3

Permabilit moyenne ou faible.

Alternances marne-calcaire

P4

Zone karstique.

Formations calcaires avec permabilit

de fissures et dveloppement d'un rseau
souterrain.

P5

Impermable.

Terrains marneux, argileux, cristallins

(granite, schistes)

A.2.5 Le couvert vgtal

Le couvert vgtal :
- influe sur l'vapotranspiration,
- retarde le ruissellement,
Il est possible d'valuer le pourcentage de la surface du bassin versant occup par les
diffrents types de vgtation.
A.2.6 Caractristiques glaciologiques

C'est le pourcentage de la surface occupe par les glaciers.

Hydrologie - Hydrogologie Septembre 2008 Vronique MERRIEN-SOUKATCHOFF

24

A.3 La Pluie
A.3.1 Mesure des prcipitations

Le principe de la mesure de la prcipitation est de poser un rcipient, d'aire horizontale A au

sol et de mesurer le volume V de pluie tombe pendant un temps t. La hauteur de pluie Ht,
tombe pendant l'intervalle de temps t est :

H t =

V
A

Le rcipient de mesure est, en fait, normalis et s'appelle un pluviomtre s'il s'agit du rcipient
seul, un pluviographe, si l'appareil est reli un systme de mesure automatique. L'appareil
doit tre install de telle matire que la surface rceptrice soit bien horizontale et 1 m de la
surface du sol. Il existe des pluviographes ayant une surface de 400 cm2 et 1000 cm2.
Suivant le systme d'acquisition associ, le pas de temps minimum de cumul du volume de
pluie est plus ou moins fin.
A.3.2 Synthse des mesures de prcipitations

Les donnes sont enregistres avec un pas de temps minimum t en un point. On peut faire
des cumuls horaires, journaliers, mensuels (Nancy environ 700 mm/an).
Il faut les critiquer immdiatement :
-

Corrlation des postes 2 2 ;

Corrlation de la somme des prcipitations (mthode de la double masse) pour se rendre

compte des drives systmatiques ;

Frquence des chiffres utiliss ;

Probabilit d'une observation ;

A.3.3 Calcul de la pluie moyenne sur un bassin versant

Plusieurs mthodes peuvent tre utilises pour dterminer la pluie moyenne sur un bassin
versant :
-

moyenne des postes ;

mthode des polygones de Thiessen;

mthodes des isohytes ;

mthodes d'interpolation diverses, en particulier, utilisation de la gostatistique.

Notion de coefficient d'abattement

On constate que pour une mme frquence d'apparition, donc un mme temps de retour
l'intensit d'une pluie est d'autant plus forte que sa dure est courte. Ou, dure de pluie
gale, une prcipitation sera d'autant plus intense que sa frquence d'apparition sera petite
(donc que son temps de retour sera grand).
Les relations entre les intensits, la dure et la frquence d'apparition des pluies peuvent tre
reprsentes selon des courbes caractristiques : on parle gnralement de courbes IntensitDure-Frquence (IDF). La notion de frquence est en faite exprime par la notion de temps

Hydrologie - Hydrogologie Septembre 2008 Vronique MERRIEN-SOUKATCHOFF

25

de retour.

Figure 28 : Reprsentation schmatique des courbes IDF

(d'aprs Andr Musy http://hydram.epfl.ch/e-drologie/)
A.3.4 Statistiques sur les donnes pluviomtriques

On veut connatre :
- la probabilit pour qu'il pleuve x mm de pluie par an ;
- qu'il tombe plus de 200 mm en une journe ;
- moins de 2 mm pendant un mois
Notion de risques directs, indirects. --> cf. cours Risques

D'aprs
http://www.eatlantide.com/thematiques/shadok/index.htm#dessin
s

A.4 Le retour de l'eau l'atmosphre : vaporation et

vapotranspiration
Une partie de l'eau de pluie (mme pendant la pluie) est immdiatement r-vapore. En effet
l'humidit de l'atmosphre est rarement sature.
Beaucoup d'instruments cits dans ces pages sont trs bien dcrit sur le site
http://users.skynet.be/bk342309/meteo/index.html.
A.4.1 Gnralits. Dfinition du pouvoir vaporant.
A.4.1.1 Les formes du retour de l'eau l'atmosphre

Il y a deux aspects dans le retour de l'eau l'atmosphre :

- Evaporation-sublimation
C'est l'ensemble des processus physiques de transformation de l'eau liquide en vapeur.
C'est l'vaporation directe partir d'une surface d'eau libre (mer, lac).
Dans ce cas, la quantit d'eau qui repart l'atmosphre est lie des paramtres physiques
tels que temprature de l'air, de l'eau, vitesse du vent, degr de saturation (degr
hygromtrique), ensoleillement...
- Transpiration (phnomne biologique)
Elle dpend du couvert vgtal, de son stade de dveloppement...
Elle est trs difficile mesurer
On regroupe ces deux phnomnes, sans les distinguer sous le terme "vapotranspiration"
Le pouvoir vaporant d'un sol dpend de plusieurs paramtres :
- atmosphre (degr de saturation de l'atmosphre),
- terrain : aptitude du terrain accepter l'vapotranspiration,
- surface terrain/atmosphre. Cette surface va rgler la vitesse de l'vaporation et donc le taux
de l'vapotranspiration.
A.4.1.2 Pouvoir vaporant de l'atmosphre.

L'atmosphre contient une certaine quantit de vapeur d'eau. La pression partielle de vapeur
d'eau dans ce mlange est appele "tension de vapeur". Elle ne peut thoriquement pas
dpasser un certain seuil appel "tension de vapeur saturante". La tension de vapeur saturante
augmente avec la temprature.
Si l'on appelle :
est la tension de vapeur saturante
et la tension relle de vapeur

le pouvoir vaporant de l'atmosphre est fonction de est - et. Il est donc fonction des
paramtres suivant :
- Temprature ; quand la temprature augmente, la tension de vapeur saturante augmente,
donc le pouvoir vaporant augmente

Hydrologie - Hydrogologie Septembre 2008 Vronique MERRIEN-SOUKATCHOFF

27

- Pression totale : si la pression totale augmente, la pression partielle augmente c'est--dire

que es augmente et par consquent le pouvoir vaporant diminue. En effet est ne dpend pas
de la pression mais seulement de la temprature.
Remarque :
Pouvoir vaporant = K (est -e t ) = K' (1 - h r ) est
h r = 100

et

e st

: humidit relative

A.4.1.3 Pouvoir vaporant des surfaces d'eau libre.

Thoriquement le pouvoir vaporant d'une surface

d'eau libre est la quantit d'eau contenue dans les
plans d'eau. En pratique il faut se rfrer au pouvoir
vaporant de l'atmosphre.

t'

air

eau

Figure 29

pression de vapeur saturante (en mm Hg)

60

50

40

30

20

10

0
0

10

t'1

15

20

et'

condensation

et'

t'2

30

35

40

temprature

Figure 30 : pression de vapeur saturante en fonction de la temprature

t' < t l'atmosphre limite le pouvoir vaporant de l'eau est' (et' est forcment infrieur
est')

t' > t

et' < est : possibilit d'vaporation jusqu' est

et' > est : il va y avoir condensation de l'atmosphre vers le plan d'eau, au

voisinage de ce dernier

Quand la salinit de l'eau augmente cela limite le pouvoir vaporant de l'eau car est eau sale
infrieur est de l'eau douce

Hydrologie - Hydrogologie Septembre 2008 Vronique MERRIEN-SOUKATCHOFF

28

A.4.2 Mesures de paramtres physiques conditionnant l'vaporation.

L'vapotranspiration est fonction de diffrents paramtres mtorologiques, dont nous

indiquons brivement les moyens de mesure.
A.4.2.1 Tempratures

- thermomtre
- thermomtre minima et maxima
A.4.2.2 Humidit de l'air

A.4.2.2.1 Psychromtre

th

mousse
+ eau

Un psychromtre se compose de deux

thermomtres : un thermomtre sec (en fait
l'humidit de l'air ambiant) et un humide.
L'vaporation consommant de l'nergie, la
temprature du thermomtre humide est infrieure
celle du thermomtre en contact avec l'air
atmosphrique.

Figure 31 : Psychromtre

Il existe des tables psychromtriques qui permettent connaissant tsec, thumide et la pression
atmosphrique d'valuer le pouvoir vaporant de l'atmosphre.
L'cart psychromtrique ts - th permet donc d'obtenir le degr de saturation de l'air

A.4.2.2.2 Hygromtres
Le psychromtre ne permet pas de faire une mesure en continu. On utilise dans un autre type
d'appareil la proprit qu'ont certains corps de s'allonger quand l'humidit augmente. Souvent
une mche de cheveux va servir de capteur. La dilatation de cette mche est amplifie,
transmise un stylet et les variations de longueur des cheveux sont enregistres sur un
diagramme entran par un mouvement d'horloge.
Problmes : tarage, vieillissement du capteur.
A.4.2.3 Pression

Il existe plusieurs types de baromtres :

- mercure,
- ressort.
Ces baromtres sont ventuellement associs un enregistrement.
A.4.2.4 Rayonnement solaire

- pyrhliomtre (mesure du rayonnement solaire direct)

- pyranomtre (mesure du rayonnement solaire global ou du rayonnement solaire diffus)

Hydrologie - Hydrogologie Septembre 2008 Vronique MERRIEN-SOUKATCHOFF

29

A.4.2.5 Vent

- anmomtres : permettent de mesurer la vitesse du vent et la quantit de "vent pass" en

hectomtres.
- girouettes. Elles permettent de mesurer la direction du vent.
A.4.3 Mesure de l'vaporation

Les mesures de l'vaporation peuvent se faire de diffrentes manires selon les buts
poursuivis
A.4.3.1 Mesures de l'vaporation partir d'une surface libre : les bacs

Il parat vident qu'une des faons la plus simple de mesurer l'vaporation est de mettre une
cuvette d'eau et de mesurer au bout d'un certain temps la quantit d'eau vapore.
Le problme vient de la normalisation de cette mesure ; il faut tre capable de rpter la
mme mesure dans les mmes conditions, en diffrents sites.
Diffrents types d'appareils ont t conus pour ce genre de mesure. Ils sont installs dans la
nature et en gnral regroups avec d'autres instruments : pluviomtre, thermomtre flottant
et anmomtre.
L'vaporation dans ces diffrents appareils dpend des conditions d'installation.
Les appareils peuvent tre classs en 3 catgories:

A.4.3.1.1 Les bacs placs au-dessus du niveau du sol

Ils ont l'avantage d'avoir des conditions d'installation simple et de ne pas recevoir de gouttes
de pluie provenant du rejaillissement des gouttes de pluie du terrain avoisinant.
Ces appareils sont par contre sensibles aux variations de tempratures et l'insolation. Pour
remdier en partie au problme il est possible d'isoler les parois extrieures du bac
exemple de bac : le bac dit de classe A
Ce bac est un cylindre de 121,9 cm et de 25,4 cm de haut. On maintient dans le bac une
paisseur d'eau de 17,5 20 cm. L'appareil est pos sur un caillebotis 15 cm du sol. Le
caillebotis permet l'aration du bac.

A.4.3.1.2 Les bacs enterrs

L'eau y est maintenue au niveau du sol.

Les fluctuations thermiques de l'appareil correspondent plutt celles du sol et ces bacs sont
donc moins sensibles la temprature ambiante.
Ils ont l'inconvnient de recueillir les gouttes de pluies ayant rebondi sur le sol mais aussi les
dtritus.
Leur installation et leur entretien sont dlicats en particulier en cas de fuite cette dernire est
difficilement dtectable.
Exemple : le bac Colorado

30

Hydrologie - Hydrogologie Septembre 2008 Vronique MERRIEN-SOUKATCHOFF

A.4.3.1.3 Les bacs flottants :

Ils sont utiliss quand on veut tudier l'vaporation de grandes surfaces. Leur installation est
difficile de mme que les mesures (en particulier en cas de vent).
Les deux principaux problmes poss par ces bacs sont donc les conditions d'installation
(reprsentativit du milieu et protection du bac) et les dimensions du bac. Pour corriger
l'vaporation par rapport la pluie, un pluviomtre de mme diamtre est associ au bac. Ce
pluviomtre doit tre install dans les mmes conditions (au sol si le bac est au sol)
En gnral les mesures sont effectues deux fois par jour 6h et 18h .
Les bacs permettent une estimation de l'vaporation d'une surface d'eau libre, mais il est
ncessaire d'apporter une correction la mesure qui tient compte des dimensions du bac. On
multiplie donc l'vaporation mesure par un coefficient de bac qui varie entre 0,6 et 0,99 et
qui augmente avec les dimensions du bac (un grand bac est assimilable une surface d'eau
libre alors que dans un petit bac l'vaporation est plus importante).
A.4.3.2 A partir des surfaces poreuses (les atmomtres)

Ces appareils sont destins mesurer une grandeur caractristique du pouvoir vaporant de
l'air ambiant.
Pour mesurer le pouvoir vaporant de l'atmosphre les agronomes et les mtorologues
utilisent des surfaces poreuses : sphres poreuses ou plaques de porcelaine blanche.
Le taux d'vaporation d'un tel systme est considr comme tant voisin de celui des plantes
(stomates).
L'atmomtre type Livingston est une sphre poreuse de 5 cm de diamtre et de 1 cm
d'paisseur remplie d'eau distille provenant d'un rservoir gradu qui sert l'alimentation et
la mesure du volume vapor.
L'vaporomtre de Piche est un type d'atmomtre. C'est un appareil simple d'emploi et de
faible cot. Il est constitu d'un tube de verre en forme de U gradu. A l'extrmit de ce tube
est fixe une capsule de buvard blanc. Le tube assure la fois l'alimentation et la mesure ; il
est gradu et ferm la partie suprieure. L'ouverture infrieure est donc obture par une
feuille circulaire de papier buvard. L'appareil ayant t rempli d'eau distille, la diminution du
niveau de l'eau dans le tube permet de calculer le taux de l'vapotranspiration.
Ces deux types d'appareils permettent de mesurer l'vaporation travers une surface poreuse.
Ils doivent avoir :
- une faible inertie thermique,
- ne pas tre perturbs par le vent,
- ne pas modifier l'humidit de l'air ambiant ( donc tre de petite taille)
A.4.3.3 Formules empiriques d'estimation du pouvoir vaporant

En l'absence de mesures directes, il existe un certain nombre de formules empiriques, bases

sur des mesures climatologiques qui tentent d'estimer l'vapotranspiration. Elles sont bases
sur la loi de Dalton :
E (vaporation) = k (est-es) / P

Hydrologie - Hydrogologie Septembre 2008 Vronique MERRIEN-SOUKATCHOFF

31

A.4.3.3.1 Formule de Lugeon

Cette formule donne l'vapotranspiration potentielle au cours d'un mois de n jours
E = 0,398 n (e sm em )

(273 + t m ) 760
273 (P e sm )

esm : tension de vapeur saturante pour la temprature maximale mensuelle moyenne

em : tension relle moyenne mensuelle
P : pression en mm de mercure.

A.4.3.3.2 Formule de Meyer

V
E = C (esm em ) 1 +
10
V : vitesse moyenne mensuelle du vent en milles/heure
Em : vaporation moyenne mensuelle
C : coefficient empirique : C = 15 pour les bacs Colorado et C = 11 pour les rservoirs, les
lacs profonds.
- en URSS il existe une formule quivalente la formule de Meyer
- Formule de Coutagne
Cette formule est base sur l'cart psychromtrique
Remarque : ces formules correspondent des rservoirs de grandes dimensions et donc une
vaporation potentielle.
A.4.3.4 Ordre de grandeur du pouvoir vaporant

En rgion parisienne :

1mm/j en hiver

3 mm/j en t
Dans les Bouches du Rhne : 3 mm/j en hiver
9 mm/j en t
mer Morte 2400 mm/an
lac Tchad 2260 mm/an
Sud de la France 1000 1500 mm/an
l'vaporation moyenne en France est de 600 mm/an

32

Hydrologie - Hydrogologie Septembre 2008 Vronique MERRIEN-SOUKATCHOFF

A.4.4 Mesures et estimation de l'vapotranspiration relle et potentielle

A.4.4.1 Notion d'vapotranspiration relle et potentielle

On appelle vapotranspiration relle note ETR la quantit d'eau (en gnral exprime en
mm) vapore ou transpire par le sol, les vgtaux et les surfaces d'eau libre d'un bassin
versant.
On appelle vapotranspiration potentielle (ETP), la quantit d'eau susceptible d'tre vapore,
si la quantit d'eau disponible tait illimite.
A.4.4.2 Mesures directes

A.4.4.2.1 Les cases lysimtriques (mesure de ETR)

Cet appareil est une cuve tanche enterre et remplie d'un bloc de sol de quelques mtres
carrs et d'environ 2 mtres d'paisseur. Le sol est drain par un lit de cailloux la base de la
case. L'eau d'infiltration (I) est recueillie. Un collecteur qui fait le tour de la cuve permet
d'valuer le ruissellement.
Les informations suivantes peuvent tre mesure :
- P : pluie (grce un pluviomtre)
- Q : ruissellement
- I : infiltration (vers la nappe)
- R : variation des rserves
(ces variations sont mesurs soit par pese soit au moyen de sondes)
Il est donc possible d'crire le bilan suivant :
P = E + Q + I + R
dans lequel seule E est inconnue.
Problmes poss par les cases :
les cases perturbent les relations sol/atmosphre et sol/sous-sol.
La vgtation que l'on peut installer au niveau d'une case n'est pas reprsentative de celle d'un
bassin versant (il n'est pas possible de planter un arbre sur une case lysimtrique !)

A.4.4.2.2 Les parcelles d'essai

Pour essayer d'avoir un chantillon plus reprsentatif, il est possible de choisir une parcelle
sur laquelle on tentera d'effectuer un bilan.
Au niveau de cette parcelle les mesures suivantes seront effectues : la pluie, le degr de
saturation du sol. En effet grce des forages le profil hydrique du sol en diffrents points
pourra tre tabli grce une sonde neutrons.
Pour contrler les apports amont et les dparts aval, la parcelle est entoure de bton.

Hydrologie - Hydrogologie Septembre 2008 Vronique MERRIEN-SOUKATCHOFF

33

A.4.4.3 Estimation de l'vapotranspiration

Il existe diffrentes formules permettant d'estimer l'ETP et l'ETR (cf. TD)

A.4.4.3.1 Evapotranspiration potentielle

- Formule de Turc
- Formule de Thorwaithe

A.4.4.3.2 Evapotranspiration relle

- Formule de Turc
- Mthode du bilan de Thorwaithe
A.4.4.4 Bilan au niveau d'un bassin versant

P, Q et E sont estimes
P = Q + E +/- S +/- SS +/- Qs
S : variation des rserves du sol
SS : variation des rserves du sous-sol
Qs : changes souterrains cachs ou invisibles
Quand on fait le bilan au niveau de l'anne, en moyenne S et SS sont pratiquement nuls,
par contre Qs n'est pas forcment ngligeable; il l'est si l'on considre un bassin versant de
grande dimension. Nous avons donc :
P=Q+E
On appelle dficit d'coulement D : D = P - Q et l'on compare cette valeur l'vaporation E .
On a observ que ce dficit d'coulement varie peu sous nos climat. Il est de l'ordre de 400
600 mm ; il varie avec la temprature l'altitude et la pression.
Il existe des formules pour approcher la valeur de ce dficit d'coulement. Il est estim
partir des paramtres mtrologiques mesurs.
-

Variations D en fonction de la temprature

Cette formule est de la forme : D = a t + b

D : valeur moyenne annuelle
t : temprature moyenne annuelle
a et b dpendent de la gographie
Sous nos latitudes pour P = 800 mm +/- 20% Coutagne propose la relation :
D (en mm) = 210 + 10 t
-

variations de D en fonction de la pluie

On a toujours D < P ( vaporation < pluie)

Le dficit atteint un maximum et dcrot ou n'augmente plus ensuite en fonction de la
pluviomtrie.
Dm correspond la tension de vapeur saturante

34
-

Hydrologie - Hydrogologie Septembre 2008 Vronique MERRIEN-SOUKATCHOFF

variations de D en fonction de la pluviomtrie et de la temprature

Dm = P - P2
= 1 / ( 0.8 + 0.14 t ) t: temprature en d C
La formule est applicable si 1 / 8 < P < 2
-

si P < 1 / 8

P = D ( faibles prcipitations)

si P > 1 / 2

D est indpendant des prcipitations

Figure 32 : ETP en France, d'aprs J. BESSEMOULIN (dir.), 1980. Atlas climatique de

la France. Ed. rduite. Red., Ministre des transports, Paris, n.p

Hydrologie - Hydrogologie Septembre 2008 Vronique MERRIEN-SOUKATCHOFF

35

Figure 33 : Rapport pluviomtrie ETP, d'aprs J. BESSEMOULIN (dir.), 1980. Atlas

climatique de la France. Ed. rduite. Red., Ministre des transports, Paris, n.p

A.5 Les coulements de surface

A.5.1 mesure des dbits

Problme : Etablir les variations de dbits Q l'exutoire d'un bassin versant en fonction du
temps. Le problme rside en l'tablissement de la courbe Q(t). Pour connatre ces dbits, des
stations de jaugeage sont installes sur les rivires.

A.5.2 Estimation des dbits de frquence rare

A.5.2.1 A partir des statistiques sur les dbits
A.5.2.2 A partir de modles pluie-dbit

Analyse d'un hydrogramme de crue

36

Hydrologie - Hydrogologie Septembre 2008 Vronique MERRIEN-SOUKATCHOFF

Annexe : Quelques rappels de statistiques

a - Rappels sur la notion de population

Soit X une variable alatoire

La probabilit que X soit infrieure une valeur x s'crit : P( X x) = F ( x)
F(x) est la fonction de rpartition de la variable alatoire X. 0 F ( x) 1

Si F est continue, elle est drivable et on peut dfinir sa densit de probabilit f(x) par :
xi

F(x i ) =

f (x)

dx

Paramtres d'une loi de probabilit :

Paramtres de tendance centrale

moyenne
mdiane
mode

Les moments :
xb

k
le moment d'ordre k est : mk = x f ( x) dx
xa

le moment centr d'ordre k est : k =

(x x )

xb

f ( x) dx

xa

Le moment d'ordre 1 est la moyenne, le moment centr d'ordre 2 est la variance, sa racine
carre l'cart type.

b - Rappels sur la notion d'chantillon.

Soit X une variable alatoire (par exemple en hydrologie : la pluie qu'il tombe en 24 heures,

l'altitude sur un bassin versant), on ne pourra rarement connatre totalement cette variable
(jamais en ce qui concerne la pluie, par exemple), ce que l'on a, c'est un chantillon.
si est un paramtre de la variable alatoire que suit la population (ex : la moyenne ), on ne
peut, partir de l'chantillon, que calculer une estimation l du paramtre (ex : m la
moyenne de l'chantillon est l'estimateur de ).
Soit un chantillon (x1, x2, xi , xn)
L'estimation de la moyenne de la population mre est : x =

1 n
xi
n i =1

Hydrologie - Hydrogologie Septembre 2008 Vronique MERRIEN-SOUKATCHOFF

L'estimation de la variance de la population mre est 2 =

1
n 1

(x

37

x)

i=1

c - Ajustement

On connat un chantillon d'une population et on va chercher ajuster les paramtres de la

population mre par rapport l'chantillon que l'on a. Plusieurs types d'ajustement sont
possibles (mthode des moments, mthode graphique, mthode du maximum de
vraisemblance)
-

mthode des moments : on calcule partir de l'chantillon un estimateur des premiers

moments (moyenne, variance et ventuellement moment d'ordre 3). On estime autant de
moments que la loi comporte de paramtres. En assimilant les moments calculs partir
de l'chantillon ceux de la population mre, il est possible de trouver les paramtres de
la loi.

mthode graphique.

A chaque valeur mesure on va associer une frquence. Pour cela on classe les valeurs
mesures par ordre croissant (x1, x2xixn). i est donc le nombre de fois o a t observ
une valeur de x infrieure ou gale xi parmi n observations. On peut montrer que dans ce
cas la frquence F(xi) associ xi est une variable alatoire dont une valeur approche de la
mdiane est :
)
i 0,5
Fi = F ( xi ) = Prob( x xi ) =
n

On estime donc la frquence de la valeur de rang i par : Fi =

i 0,5
n

ou encore Fi est un estimateur de la fonction de rpartition ou de la frquence au non

dpassement F(xi).
En reportant sur un graphique les valeurs associes leur frquence sur du papier gradu
selon la loi employe, il est possible de vrifier l'ajustement graphique et d'en dduire les
paramtres de la loi.
d - Danger d'apparition

Priode de retour :
La priode de retour T est l'inverse de la frquence au non dpassement F, pour des
frquences infrieures 0,5 et l'inverse de la frquence au dpassement, pour des frquence
suprieures 0,5.
F ou F(x) : Frquence au non dpassement (probabilit que la pluie soit infrieure x)
(1-F) est donc la frquence au dpassement (probabilit que la pluie soit suprieure x)
si F < 0,5
si F < 0,5 : T =

1
F

si F > 0,5 : T =

1
1 F

38

Hydrologie - Hydrogologie Septembre 2008 Vronique MERRIEN-SOUKATCHOFF

Les Frquences ont pour unit T-1, puisqu'une dure est associe chaque donne (pluie
annuelle, pluie maximale mensuelle). La priode de retour a donc la dimension d'un temps.

Hydrologie - Hydrogologie Septembre 2008 Vronique MERRIEN-SOUKATCHOFF

39

Bibliographie Hydrosystmes
[1] ALERY R., GRISOLLET H.et GUICHMET B., Climatologie. Mthodes et pratiques.
[2] BEAR J., TSANG C.-F., DE MARSILY G., 1993, Flow and Contaminant transport in
fractured Rock, Academic Press, Inc.
[3] BRAS R. L., 1990, Hydrology, an introduction tohydrologic Sciences, Addison-Wesley
Publishing Compagny
[4] HILLEL Daniel, 1988, L'eau et le Sol. Principes et processus physiques, collection
Pedasup 5, 2me dition.
[5] REMENIERAS G., 1965, Hydrologie de l'ingnieur, 2e dition, Eyrolles, Paris., ROCHE
M., 1963, Hydrologie de surface., Gauthier-Villars, Paris
[6] TRIPLET J.P. et ROCHE G. Mtorologie gnrale. Edition de l'Ecole Nationale de la
Mtorologie.

Quelques revues

La Houille Blanche

Les cahiers de l'ORSTOM

Journal of Hydrology

Hydrology

40

Hydrologie - Hydrogologie Septembre 2008 Vronique MERRIEN-SOUKATCHOFF

Quelques sites web HYDRO & Environnement

http://waternet.com/

Welcome to WaterNet

http://water.usgs.gov/

USGS Water Resources Information

news:sci.geo.hydrology
http://www.usc.edu/dept/fccchr/ccvlib/
http://dir.yahoo.com/Science/earth_sciences/h
ydrology/
http://hydram.epfl.ch/edrologie/general/index.html

Newsgroup: sci.geo.hydrology
The World Wide Web Virtual Library: CrossConnection Control/Backflow Prevention
Yahoo! - Science:Earth Sciences:Hydrology
Cours Hydrologie Gnrale (Prof. Andr
Musy) Site fortement conseill

http://www.ecologie.gouv.fr/rubrique.php3?id Ministre de l'cologie et du dveloppement _rubrique=4

Domaine de leau
http://www.cig.ensmp.fr/~thiry/courshydro/hydro00.htm

Expliquer l'hydrogologie par Mdard Thiry

http://www.epa.gov/

US Environmental Protection Agency

http://www.cig.ensmp.fr/~hubert/glu/aglo.htm

Glossaire International d'Hydrologie

et
http://www.cig.ensmp.fr/~hubert/glu/HINDF
R.HTM
http://www.u-picardie.fr/beauchamp/

Plein de choses intressantes

Hydrologie - Hydrogologie Septembre 2008 Vronique MERRIEN-SOUKATCHOFF

41

B HYDROGEOLOGIE ET HYDRODYNAMIQUE : Le milieu

naturel, description et fonctionnement. Nappe, types et
circulations
B.1 Introduction - gnralits
Tous les terrains contiennent un pourcentage plus ou moins important de vides.
L'eau peut gnralement pntrer dans ces vides, y circuler et parfois s'y accumuler.
La prsence d'eau dans les terrains est importante, d'une part parce qu'elle constitue des
rserves en eau potable et industrielle (qu'il est parfois possible d'exploiter), d'autre part parce
qu'elle influe sur les proprits mcaniques des terrains.
Nous ne nous intresserons ici qu' l'aspect rserve, mais l'aspect mcanique n'est pas pour
autant ngligeable et bien des problmes gotechniques surviennent parce que l'eau n'a pas
t prise en compte ou a t sous-estime dans les calculs.
L'eau dans les terrains sera tudie dans ce qui suit deux chelles et sous deux aspects :
- Nous nous intresserons dans un premier temps ce qui se passe au niveau d'un
chantillon : nous dcrirons les rapports entre les terrains et l'eau d'un point de vue statique ;
ce sera l'tude de la porosit, puis nous dfinirons la permabilit qui caractrise la
circulation de l'eau dans le sol.
- Dans une deuxime partie, nous tudierons les nappes et les terrains aquifres : leur
description et les circulations.

B.2 L'eau dans le rservoir : la porosit

B.2.1 Caractrisation : la porosit

Considrons un volume lmentaire de terrain reprsentatif (VER) i.e. un volume

suffisamment important pour que les proprits du terrain ne varient pas quand on augmente
le volume de l'chantillon, mais suffisamment petit pour tre caractristique des proprits
locales.

42

Hydrologie - Hydrogologie Septembre 2008 Vronique MERRIEN-SOUKATCHOFF

La porosit caractrise l'aptitude d'un sol contenir un fluide.

Si un volume V de terrain l'chelle macroscopique contient un
volume Vv de vides et un volume Vs de solide (V = Vv + Vs) la
porosit est le rapport n :
VV
V

n=

Vv
eau

(ce rapport est souvent exprim en %)

VS

Ce rapport est aussi not par les hydrogologues et d'autres

notations sont employes par les mcaniciens des fluides.
Un autre paramtre est galement utilis, plutt par les
mcaniciens des sols ; c'est l'indice des vides e :
e=

VV
VS

avec la relation :

Figure 34

n=

e
1+e

ou

e=

n
1n

Les vides peuvent contenir plus ou moins d'eau et le degr de saturation (S) caractrise le
pourcentage d'eau contenu dans les vides :
Sr =

Volume d ' eau contenu dans les vides du matriau

Volume total des vides

Sr =

Vw
VV

On utilise parfois la teneur en eau volumique

=

Volume d ' eau contenu dans les vides du matriau

Volume total du matriau

B.2.2 Diffrents types de classification de la porosit

Plusieurs critres peuvent tre utiliss pour diffrencier la porosit :

- la taille des pores ;
- la description des pores (i.e. le type de porosit) ;
- l'origine de la porosit.
B.2.3 Classification de la porosit par rapport la taille des pores

La taille des pores est trs variable. On parlera de :

- porosit rticulaire quand la taille des vides est de l'ordre de l'angstrm (1 = 10 -10
m). L'eau contenue dans ces vides est mobilisable par vaporisation par suite des variations du
degr hygromtrique de l'air;
- porosit collodale pour des vides d'environ 100 . Elle correspond aux vides des
agrgats collodaux. C'est une porosit qui peut tre importante et que l'on rencontrera
principalement dans les argiles. L'eau contenue dans ces vides est mobilisable par compaction
naturelle ou provoque (centrifugation, filtration sous presse ou sous vide);
- microporosit jusqu' 2 10-7 m;

Hydrologie - Hydrogologie Septembre 2008 Vronique MERRIEN-SOUKATCHOFF

43

- porosit capillaire entre 2 10-7 m et 2 10-3 m ;

- macroporosit au-del de 2 mm.
Les vides dont on pourra extraire l'eau et qui vont donc intresser l'hydrogologue
correspondent la macroporosit, la porosit capillaire et dans une moindre mesure la microporosit.
B.2.4 Classification par rapport l'origine

L'origine de la porosit peut tre primaire ou secondaire :

- la porosit primaire est forme par les pores crs au cours de la gense de la roche :
lors de la sdimentation, au cours de la cristallisation ou du refroidissement ;
- la porosit secondaire est acquise aprs la gense soit par fracturation, soit par
dissolution (ex: grs ciment calcaire ; la dissolution du ciment calcaire va entraner
l'acquisition d'une porosit secondaire).
B.2.5 Classification morphologique

On distingue deux grands types morphologiques de vides : les pores et les fissures.
B.2.5.1 La porosit d'interstices (intergranulaire) : les pores

C'est l'ensemble des vides compris entre les diffrentes particules

d'un terrain ; elle spare les "grains".
La porosit peut tre ouverte ou ferme (cas de certaines laves
volcaniques) selon que les vides communiquent ou non les uns
avec les autres.
Suivant la taille des pores, il sera possible de distinguer une
porosit d'interstices rticulaire (entre les cristaux des roches
magmatiques et mtamorphiques), collodale (argiles), une
microporosit et une macroporosit.

Figure 35 : porosit
d'interstice

B.2.5.1.1 La porosit d'interstices simple ou nette

Quand les grains sont bien classs, c'est--dire sont de taille quivalente et que les vides qu'ils
laissent ne sont pas remplis par des grains de plus petite taille, la porosit sera qualifie de
nette.
La structure peut tre plus ou moins compacte selon le tassement. Si on suppose que l'on a
des particules sphriques de mme diamtre, la porosit dpendra de l'arrangement des
sphres : la disposition pourra varier d'une disposition en carr (arrangement le plus lche)
une disposition losangique (ou rhombodrique en 3 dimensions) qui donneront des porosits
de 45 % (au maximum) 25 %.
B.2.5.1.2 La porosit d'interstices restreinte
Ce type de porosit provient d'un mauvais tri des grains qui entrane un remplissage par des
particules fines des vides laisss entre les gros grains.
B.2.5.1.3 La porosit d'interstices rduite

44

Hydrologie - Hydrogologie Septembre 2008 Vronique MERRIEN-SOUKATCHOFF

Le volume des vides peut tre "rduit" par un dpt (carbonate de chaux, hydroxyde de fer,
silice ) qui se fait sur la surface des grains et diminue la taille des pores.
En plus de ces trois types (simple, restreint, rduite), il est possible de trouver une porosit
double, quand les "gros" lments sont eux-mmes composs de grains et de pores plus petits
que les vides laisss par les "gros" lments.
B.2.5.2 La porosit de fissure

Il existe plusieurs sortes de "fissures". Nous

pourrons distinguer plusieurs types de porosits
de fissures en fonction de la nature de ces
dernires.

Figure 36 : porosit de fissure

B.2.5.2.1 Porosit de joints

Elle est due aux joints stratigraphiques. Cette porosit est primaire.
B.2.5.2.2 Porosit de diaclases
C'est une porosit secondaire lie aux diaclases donc des fissures sans rejet orthogonales ou
obliques par rapport la stratification.
B.2.5.2.3 Porosit de failles
Nous parlerons de porosit de faille lorsqu'il existe un rseau de fractures bien dvelopp li
la prsence d'une faille proximit.
B.2.5.2.4 Porosit de schistosit
Le long des plans de schistosit, si ces plans se dcollent plus ou moins, il peut se former
quelques vides.
B.2.5.2.5 Porosit de retrait
Ce type de porosit, relativement restreint, est li au refroidissement des roches ruptives.
B.2.6 Ordre de grandeur de la porosit

Sables, grs 15 % 25 %.

Argiles 40 % 90 % (la porosit des argiles peut parfois tre suprieure 100 % car le
volume total augmente). Cette porosit ne correspond pas de l'eau mobilisable et
n'intresse donc gnralement pas l'hydrogologue.

Marnes : 30 % 50 %, mais une partie de cette porosit est collodale.

Calcaires : Quelques % 25 % (dans le cas d'un calcaire dtritique fissur).

Roches cristallines : quelques %. Cette porosit peut augmenter du fait de la fracturation

et de l'altration.

Hydrologie - Hydrogologie Septembre 2008 Vronique MERRIEN-SOUKATCHOFF

45

Remarques sur l'ordre de grandeur de la porosit

La porosit la plus importante correspond la porosit d'interstices : une disposition en

carr de sphre rgulire donnerait une porosit de 45 %, mais un grs ciment peut avoir
une porosit d'environ 5 % seulement. La porosit de fissure est moins importante ; pour
se donner une ide, il est possible d'imaginer un bloc de 20 cm 20 cm 25 cm sur le
bord duquel se trouve une fracture de 1 mm de large ; ceci correspond une porosit de
0,4 %. Gnralement, la porosit de fissure est infrieure 5 %, mais cette porosit est
trs importante du point de vue de la circulation des eaux (trajet prfrentiel).
Les diffrentes porosits peuvent s'ajouter, par exemple une porosit de fissures et
d'interstices dans un grs. La porosit double est trs intressante (elle peut correspondre
une porosit de fissures et d'interstices combines) car le dbit traversant une section est
proportionnel au carr du diamtre des vides.
Dans le cadre de ce cours nous supposerons en gnral que la porosit en un point n'est pas
volutive, alors que les variations de l'tat de contrainte et de la pression de l'eau conduisent
des phnomnes coupls, qui engendre, entre autre, l'volution de cette porosit. La porosit
dcroit quand le module des contraintes auquel est soumis le terrain augmente. Ce
phnomne et l'volution ptrographique des roches de la surface vers la profondeur, fait que
gnralement la porosit dcroit quand la profondeur augmente.
B.2.7 Mobilit de l'eau dans le sol : eau lie, eau libre, l'gouttage des roches

Seule une partie de l'eau contenue dans un terrain est mobilisable. Les auteurs qui
s'intressent ce problme l'abordent de diffrentes manires (cf. notamment [5] et [8]). Nous
allons essayer de retranscrire les diffrentes observations qui sont ralises, mais les
diffrentes perceptions ne sont pas toujours faciles relier entre elles.
B.2.7.1 Eau libre, eau lie

Si on ne s'intresse dans un premier temps qu'au milieu satur (cf. De Marsily [8]), on peut
distinguer l'eau libre et l'eau lie.
L'eau lie est rattache la surface des grains par le jeu de forces d'attraction molculaire.
Ces forces dcroissent avec la distance au grain :

l'eau adsorbe constitue un film continu, une pellicule de 1/10 de microns (quelques
dizaines de molcules). Son volume relatif augmente en proportion inverse de la taille des
particules : c'est 2 5% du volume d'eau dans les sables grossiers et jusqu' 50% dans les
argiles;

entre 0,1 et 0,5 les molcules d'eau supportent une attraction non ngligeable et sont
immobiles;

au del de 0,5 1 (donner une limite est un peu arbitraire), les forces d'attraction sont
ngligeables et l'eau est dite libre.

46

Hydrologie - Hydrogologie Septembre 2008 Vronique MERRIEN-SOUKATCHOFF

L'eau libre est donc celle qui

est en dehors du champ
d'attraction
des
particules
solides et qui est susceptible de
se dplacer sous l'effet de la
gravit ou des gradients de
pressions.
particule
solide :
"grain"

attraction molculaire
eau adsorbe

molcules
libres

0,5 1

distance
au grain
eau libre

eau lie

Figure 37 : Schma de la structure de l'eau souterraine

au voisinage d'un grain, d'aprs Polubrina-Kochina
(1962)

Le volume o l'eau peut

circuler est infrieur la
porosit totale puisque l'eau
lie peut tre considre du
point de vue du dplacement
des fluides comme faisant
partie du solide. La porosit
efficace (ou cinmatique) ne
(ou c) exprime le pourcentage
d'eau qui peut circuler par
rapport au volume total du
terrain; elle est donc infrieure
la porosit totale. Des
phnomnes
autres
que
l'adsorption peuvent limiter la
porosit cinmatique : les pores
non connects, les pores culs de
sac.

B.2.7.2 L'gouttage des roches

Si on extrait l'eau d'un chantillon de petite dimension par gouttage [5], puis centrifugation,
on va mobiliser des volumes croissants d'eau :
-

Un chantillon (de volume V) initialement satur, que l'on laisse goutter sur une grille,
libre un volume Ve d'eau. Ce volume libr par la gravit est appel eau gravitaire. Au
bout d'un certain temps (environ une journe), l'chantillon ne libre pratiquement plus
d'eau. On peut observer que le rapport Ve/V augmente avec la taille de l'chantillon et la
granulomtrie du terrain.

Si le mme chantillon goutt est plac dans une centrifugeuse, il libre encore un
certain volume d'eau dite eau pelliculaire. L'eau pelliculaire reprsente un film de l'ordre
du micron qui peut se dplacer la surface des grains sous l'action des molcules d'eau
voisines.

L'chantillon renferme encore un certain volume d'eau : l'eau adsorbe.

L'eau pelliculaire et adsorbe constituent ce que l'on appelle l'eau de rtention (rtention la
gravit)
On ne peut pas mettre exactement en parallle ce que nous avons qualifi d'eau libre et d'eau
lie avec l'eau gravitaire et l'eau de rtention car dans un milieu non satur, il faut donc
prendre en compte la prsence de l'air en plus de l'eau et du solide. Nous avons prsent les
diffrents stades de la libration de l'eau par gouttage, puis centrifugation. Il est galement
possible de dcrire ce qui est observ en fonction du degr de saturation.

Hydrologie - Hydrogologie Septembre 2008 Vronique MERRIEN-SOUKATCHOFF

47

"On appelle passoire tout instrument sur

lequel on peut dfinir trois sousensembles: l'intrieur, l'extrieur et les
trous"
L'intrieur est gnralement plac audessus de l'extrieur et se compose le
plus souvent de nouilles et d'eau.
Les trous ne sont pas importants. En
effet, une exprience simple permet de
se rendre compte que l'on ne change pas
notablement les qualits de l'instrument
en rduisant de moiti le nombre des
trous, puis en rduisant cette moiti de
moiti... etc... etc... et la limite jusqu'
ce qu'il n'y ait plus de trous du tout.
D'o THEOREME:

Remarque :

La notion de passoire
est indpendante de la
notion de trous et
rciproquement.

On appelle passoires du Premier Ordre

les passoires qui ne laissent passer NI
les nouilles NI l'eau.
Pour prolonger vos connaissances sur le
sujet cf.
http://www.lesshadoks.com/index2.php?
page=11

B.2.7.3 Relations air-eau pour diffrents degrs de saturation

Dans un sol contenant la fois de l'air et de l'eau, on remarque que l'eau libre entoure les
grains tandis que l'air a tendance tre dispos au milieu des vides. Suivant le degr de
saturation, diffrents stades peuvent tre observs :
-

Dans un sol proche de la saturation, l'eau peut circuler sous l'influence de la gravit, c'est
l'eau gravifique, ou encore qualifie de funiculaire, car elle forme des fils continus. La
phase air, qui peut occuper jusqu' 10 15% des vides, est discontinue et ne circule pas.

Lorsque la phase eau est continue, mais ne circule pas sous l'action de la gravit et que la
phase air est galement continue mais ne circule pas, le sol est la saturation d'quilibre
ou encore la capacit de rtention capillaire.

Pour de faibles teneurs en eau, l'eau entoure les grains du terrain et occupe des espaces
qui ne sont continus qu'aux points o les grains du sol sont proches. On parle d'eau
pendulaire ou d'anneaux pendulaires pour ces fils d'eau qui entourent les grains. L'air est
continu.

Si la teneur en eau dcrot encore par des phnomnes d'vaporation ou de transpiration

des vgtaux, la teneur en eau va pouvoir encore dcrotre et il ne restera que l'eau lie ou
hygroscopique.

B.2.7.4 Pression capillaire

L'quilibre entre l'air et l'eau dpend des forces de tension superficielles. La diffrence de
pression entre les deux fluides est appele pression capillaire. Cette pression capillaire

48

Hydrologie - Hydrogologie Septembre 2008 Vronique MERRIEN-SOUKATCHOFF

explique que des zones puissent tre pratiquement satures alors que la pression de l'eau est
infrieure la pression atmosphrique.
On distingue l'eau capillaire continue dont la prsence est due l'ascension capillaire et qui
est donc en continuit avec une zone sature et l'eau capillaire suspendue qui existe mme
dans les zones non satures.
B.2.8 Profil hydrique des sols

Nous avons vu que suivant le degr de saturation l'eau pouvait se trouver diffrents stades.
Ces diffrents types d'eau se retrouvent dans un sol. Nous pouvons sparer le sol en une zone
non sature et une zone sature :
-

la zone non sature ou d'aration peut-tre subdivise en fonction des teneurs en eau
qui croissent vers le bas en trois sous-zones :
la zone d'vapotranspiration la surface. Cette partie du sol est soumise des
variations importantes de teneur en eau provoques par l'infiltration et
l'vapotranspiration. Sa profondeur est variable suivant le type de sol et le climat;
La zone de transition o la teneur en eau est voisine de la capacit de rtention
(cf. 0);
Au-dessus du niveau de la nappe se trouve une zone pratiquement sature 100%
o la pression de l'eau est infrieure la pression atmosphrique. La saturation du terrain
se fait par ascension capillaire. Cette zone est appel frange capillaire. Il peut y avoir une
certaine quantit d'air pige dans cette zone; la saturation est alors lgrement infrieure
100% (85 90%).

La zone sature ou nappe, dont le niveau peut-tre mesur par un tube pizomtrique.

B.3 Les mouvements de l'eau dans les roches : permabilit

Nous nous contenterons, dans cette partie, d'une tude essentiellement descriptive. Une vision
plus mcaniste est dveloppe dans une autre partie.
Pour que l'eau circule dans un terrain il est ncessaire que les vides (pores et fissures) soient
interconnects. L'aptitude d'un terrain se laisser traverser par les fluides est caractrise par
la permabilit de ce terrain par rapport au fluide.
B.3.1 Rappel sur la charge hydraulique
B.3.1.1 Charge hydraulique

Nous rappelons ici quelques notions de mcanique des fluides.

Considrons un fluide parfait i.e. incompressible et non visqueux. Si ce fluide est en
mouvement et que sa vitesse ne varie pas dans le temps, c'est--dire que son mouvement est
permanent, les particules suivent des trajectoires invariables dans le temps. Dans ce cas la
trajectoire = filet liquide = ligne de courant (nous rappelons que la ligne de courant est la
ligne tangente au vecteur vitesse en chacun de ces points l'instant considr).
On appelle charge hydraulique la quantit H :
V2
P
H=
+
+z
2 g g

Hydrologie - Hydrogologie Septembre 2008 Vronique MERRIEN-SOUKATCHOFF

49

z tant l'altitude du point, P la pression, w la masse volumique du fluide (w est le poids

volumique)
Remarque : On note parfois en mcanique des sols la pression u, donc :
V2
u
H=
+
+z
2g w
Remarque : Thorme de Bernouilli : Si le liquide est parfait, la charge hydraulique reste
constante.
En fait gnralement un fluide n'est pas parfait et il existe des forces de viscosit ou de
frottement visqueux. C'est le cas pour l'eau s'coulant sous l'action de la pesanteur travers
les vides d'un terrain : il existe des forces de viscosit entre les molcules ; ces frottements
vont dissiper de l'nergie et il y aura perte de charge.
Souvent nous nous intresserons plus aux variations de charge dans l'espace qu' la charge
elle-mme. Lorsqu'une particule parcourt la distance L, le gradient hydraulique I est dfini
par :
I=

H H 2 H 1 dH
=
=
= grad (H )
L
L
dl

B.3.1.2 Cas des sols

B.3.1.2.1 Charge Hydraulique

Les vitesses d'coulement dans le sol sont toujours faibles (mme dans un sol trs permable
l'ordre de grandeur est 0,1 m/s). Par consquent dans l'expression de la charge hydraulique, le
terme V2 / (2 g) est ngligeable par rapport aux autres. On notera donc :
H=

+z

On exprime souvent les charges par rapport au nivellement gnral (NGF), comptes comme
des altitudes topographiques. La pression atmosphrique est souvent omise, dans cette
expression, dans la mesure o on s'intresse principalement aux niveaux saturs et donc aux
surpressions par rapport la pression atmosphrique.
B.3.1.2.2 Notion de hauteur pizomtrique
Considrons un coulement d'eau dans un terrain et un point M la cote z. Faisons descendre
un tube plein jusqu' ce point M. Nous observons une remonte de l'eau dans ce tube jusqu'
la cote z'. La charge hydraulique au point M peut s'crire :
H=

+z=

w ( z ' z )
+ z = z'
w

La charge hydraulique au point M est donc gale l'altitude de la surface d'eau dans un tube
plein que l'on qualifie de tube pizomtrique ou pizomtre. La hauteur z' sera appele
hauteur pizomtrique.
La hauteur pizomtrique en un point est donc la hauteur d'quilibre de l'eau avec la pression
atmosphrique.

50

Hydrologie - Hydrogologie Septembre 2008 Vronique MERRIEN-SOUKATCHOFF

ATTENTION cette hauteur est gnralement diffrente de la surface libre de la nappe. En

effet, si dans un milieu satur, la nappe s'coule horizontalement et que la charge reste la
mme sur une verticale, la cote de la surface libre reste toujours celle mesure par le
piezomtre quelle que soit sa profondeur. Si par contre l'coulement n'est pas horizontal, la
charge varie avec la profondeur du pizomtre (la pression ne correspond plus la hauteur
d'eau au dessus du point) et la surface libre est dfinie par la cote obtenue quand le
pizomtre commence pntrer dans le milieu satur.

tube pizomtrique
surface libre
niveau pizomtrique
P
: hauteur pizomtrique
w g
x

Substratum
Figure 38 : niveau et hauteur pizomtrique
B.3.2 Exprience de Darcy

Exprimentalement, Le Chevalier
Henry Darcy (vers 1856) trouve la
relation suivante :

Q = K S

H
L

S : section du massif sableux

K est un coefficient qui dpend du
fluide et du terrain. Il a la dimension
d'une vitesse (L T-1). Ce coefficient
est le coefficient de permabilit de
Darcy encore appel coefficient de
permabilit.
En posant V =

Q
H
et I =
S
L
Figure 39 : exprience de Darcy

La relation se transforme en :
V = KI qui est l'expression la plus simple de la loi de DARCY.

Hydrologie - Hydrogologie Septembre 2008 Vronique MERRIEN-SOUKATCHOFF

51

V est la vitesse apparente moyenne, encore appele vitesse de Darcy ou vitesse hydraulique
ou vitesse de filtration. Le terme "vitesse de percolation" est aussi parfois utilis.
(en toute rigueur si H = H sortie H entre l'coulement s'effectuant dans le sens des charges
dcroissantes il faudrait crire V = K I ).
H 2 H1
si la diffrence de charge H 2 H 1 et la distance L deviennent trs petite,
L
dH 1
.
la limite la loi de Darcy peut tre crite sous la forme diffrentielle : V = K
dL
V = K

L'expression de la vitesse en fonction du coefficient de permabilit a t gnralise par

Schlichter2 (cit par Magnan) au cas dun coulement tridimensionnel dans un sol homogne
et isotrope, sous la forme :
r
r
uuuuur
v = - Ki = - K gradH
Le coefficient K dpend du terrain travers et du fluide qui traverse ce terrain. En se basant
sur les quations aux dimensions et par vrification exprimentale, on trouve que la constante
K varie en fonction inverse de la viscosit dynamique du fluide [ML-1T-1] ( 1 cPo :
centipoise pour l'eau 20C, 1 cPo =10-2 Po : poises; 1 cPo =10-3 Pa.s : Pascal.seconde = 103
N.s/m2).
Remarque : viscosit cinmatique =

= 10 -6 m 2 / s ou myriastokes [L2T-1]

La loi de Darcy peut donc s'exprimer sous la forme gnrale : V =

k
grad p + g grad z

V ou VD est une grandeur macroscopique de mme que , et p. On est amen caractriser

ces proprits sur un Volume Elmentaire Reprsentatif (V.E.R.) dont les dimensions
dpendent de la structure du matriau. La vitesse relle (ou vitesse intergranulaire, ou vitesse

d'aprs http://www.e-atlantide.com/thematiques/shadok/index.htm#dessins
2 Schlichter, C. S. 1905. Field Measurements of the Rate of Movement of Underground Waters. U.S. Geol.
Surv. Water Supply Paper 140. 122 pp.

Hydrologie - Hydrogologie Septembre 2008 Vronique MERRIEN-SOUKATCHOFF

52

de pore) du fluide Vr est suprieure la vitesse de Darcy VD.

Vr =

Q
ou SV est la surface des vides offerte l'coulement.
SV

Vr =

Q S VD V VD VD
=

=
nc : porosit cinmatique
SV
SV
VV
nc

Le coefficient k ou ki coefficient de permabilit intrinsque est reli au coefficient de

permabilit K par la relation :

K = ki

Le coefficient de permabilit intrinsque n'est pas fonction du liquide. Il est caractristique

du terrain travers. La permabilit intrinsque n'est dfinie qu' l'chelle macroscopique. Sa
dimension est celle d'une surface [L2]. On l'exprime souvent en :
DARCE : 10-12 m2
DARCY : 0,987 10-12 m2 ou MILIDARCY (10-3 DARCY)
Une permabilit intrinsque k de 1 Darcy correspond un coefficient de permabilit K de
10-5 m/s.
Ces units sont notamment employes dans le domaine ptrolier et les gisements ptroliers
ont des permabilits courantes variant de 1 quelques milliers de milidarcy.
Le coefficient de permabilit de Darcy K est fonction du fluide, mais dans notre cas nous ne
nous intresserons qu' l'eau. Il faut noter que ce coefficient est galement fonction de la
viscosit du fluide qui varie sensiblement avec la temprature ( est de 1,787 cPo pour l'eau
0C, 1,310 cPo 10C et 1,002 cPo 20C). En gnral on fait l'hypothse que la
temprature est constante, cependant des variations climatiques importantes pour les nappes
superficielles peuvent engendrer des variations non ngligeables de la permabilit.
Si l'on s'intresse la permabilit l'eau 20C on a donc K ( en m / s ) *107 = ki (en m 2 )
B.3.3 Exprience de Reynolds (pour mmoire)

L'coulement d'un fluide peut se produire de deux manires diffrentes selon les conditions
locales de vitesse. A faible vitesse les lignes de courants sont stables et ne se "mlangent"
pas. Dans cet coulement appel laminaire, les couches fluides glissent les unes sur les autres
et il n'y a pas de transfert de particules d'un filet fluide un autre. Lorsque la vitesse crot, les
filets fluides paraissent osciller et vibrer, puis ils perdent leur identit propre. Dans ce rgime
appel turbulent, les particules oscillent autour d'une trajectoire moyenne. Le passage d'un
rgime l'autre dpend de la valeur d'un paramtre adimensionnel, le nombre de Reynolds:
U : vitesse caractristique de l'coulement
D : est une des dimensions gomtrique
Re =

UD

: viscosit cinmatique du fluide =

(la viscosit cinmatique de l'eau 20C est de 10-6 m2/s)

Dans le cas d'une conduite circulaire le nombre de Reynolds critique est de 2000. Si
Re < 2000 le rgime est laminaire et si Re >> 2000 le rgime est turbulent (dans ce cas la

Hydrologie - Hydrogologie Septembre 2008 Vronique MERRIEN-SOUKATCHOFF

dimension est D est le diamtre hydraulique DH, DH =

53

4S
avec S : section mouille, P
P

primtre mouill).
Dans le cas des sols, on peut pratiquement considrer que l'on se trouve toujours en rgime
laminaire sauf aux abords d'ouvrages particuliers (puits de pompage par exemple) o les
vitesses peuvent devenir trs importantes.
B.3.4 Ecoulement dans les roches stratifies

La permabilit K dpend du matriau. Supposons un matriau anisotrope form par la

superposition de couches horizontales d'paisseur ei et de permabilit Ki ; l'coulement se
fait la vitesse V qui peut tre dcompose en VH + VV.
B.3.4.1 Permabilit horizontale :

He
q1

A travers chaque strate s'coule un dbit :

H
q i = Ki (e i 1)
L

Hs
e1

Le

dbit

total
est
donc
n
H H
Q = qi = Ki (e i 1)
=
(K i e i )
L L i=1
i =1
i =1
n

qi

ei

mais le dbit total peut galement s'crire :

H n
Q = KH
(e i 1)
L i =1
n

(K e )

en

Figure 40 : permabilit horizontale

quivalente
B.3.4.2 Permabilit verticale.

Il est ais d'tablir que : K v =

e i
e
i
Ki

donc K H =

i =1

e
i =1

Hydrologie - Hydrogologie Septembre 2008 Vronique MERRIEN-SOUKATCHOFF

54

Q
e1

H0
H1
Hi-1

ei
Hi

Le dbit vertical qui passe travers l'ensemble des couches:

(H - Hn )
Q = S Kv I = S Kv 0
est galement le dbit qui
ei
passe au travers de chaque couche. On peut galement crire
(H - Hi )
Q = S K i I i = S K i i-1
Hi-1 et Hi-1 tant la charge
ei
hydraulique respectivement au somment et la base de la
couche i.
n

Hn-1

Or H 0 Hn = (Hi -1 Hi )
i =1

en

Q e i

Hn

donc

i=1

S Kv

Q ei
=
soit K v =
i =1 S K i

Figure 41 : permabilit
verticale quivalente

e
i=1
n

ei

K
i=1

B.3.5 Mesures et estimation de la permabilit au laboratoire

Dans ce paragraphe, nous ne mentionnerons que les mthodes utilises en laboratoire. Il faut
garder l'esprit que ces mthodes ne permettent pas de mesurer correctement la permabilit
de l'ensemble des terrains. Pour estimer la permabilit d'un terrain dans son ensemble, des
mthodes in situ (notamment essai de pompage) sont utilises.
B.3.5.1 Problmes poss par l'chantillonnage

Pour estimer ou mesurer la permabilit au laboratoire, il est ncessaire de prlever un

chantillon de terrain. Cet chantillon, de petite taille, ne sera pas reprsentatif de l'ensemble
de l'aquifre :
-

Les caractristiques du terrain seront modifies du fait de l'chantillonnage ;

L'chantillon ne permettra pas de prendre en compte les variations de permabilit dues

aux failles ;

L'aquifre sera en gnral prlev l'affleurement (o le terrain est modifi par

l'altration). Pour constituer un chantillon caractristique, il faudrait faire des
prlvements diffrents niveaux de l'aquifre, ce qui est difficilement ralisable et serait
coteux ;

L'chantillon ne se trouvera pas dans les conditions de pression, de forces adjacentes et de

temprature qui taient primitivement les siennes et qui sont difficilement valuables.

Ces techniques de mesures ou d'estimation de la permabilit en laboratoire sont en fait plus

utilises par les mcaniciens des sols que par les hydrogologues (en effet si l'on travaille sur
des sols remanis, comme le sont les chantillons, l'ordre de grandeur de la permabilit
fourni peut tre acceptable).
B.3.5.2 Estimation de la permabilit

La permabilit au laboratoire peut tre estime partir de la granulomtrie (relation de

Hydrologie - Hydrogologie Septembre 2008 Vronique MERRIEN-SOUKATCHOFF

55

Hazen ou relation de Casagrande) dans le cas d'une roche meuble :

B.3.5.2.1 Relation de Hazen :

A partir d'expriences effectues avec des sables filtre, d'uniformit leve (Cu < 2;
D
Cu = 60 ) et peu compacts, Hazen a obtenu les quations empiriques suivantes :
D10
2
K (en cm/s) = C1 d 10

d10 : diamtre en de duquel il y a 10 % des grains ou diamtre efficace en cm,

C1 est un coefficient variant entre 100 et 150 s.cm-1.

B.3.5.2.2 Relation de Casagrande :

Pour des sols gros lments (> 1 mm) dont les grains sont supposs cubiques, on peut
exprimer la permabilit en fonction de l'indice des vides e :
K = 1,4 K 0.85 e 2
K0.85 est la permabilit pour e = 0.85. Il suffit donc de dterminer la permabilit
correspondant une valeur arbitraire de e et on obtient les valeurs de K correspondant
d'autres valeurs de e au moyen de l'quation.
Ces relations ne tiennent pas compte de la forme des grains. Elles ne doivent tre utilises que
pour les cas prcis pour lesquels elles ont t dfinies. Dans la pratique, elles sont
inutilisables pour les terrains naturels qui ont des structures diffrentes et plus complexes que
les sols tudis.
B.3.5.3 Permamtres

La permabilit peut tre galement mesure au moyen d'un permamtre sur un chantillon
de terrain.

B.3.5.3.1 Permamtre charge constante

L'chantillon prlev est ramen aux dimensions requises pour l'appareil de mesure. Il est
ensuite mis saturer, puis l'prouvette est traverse par un fluide de telle faon que la charge
au sommet de l'chantillon soit constante.
La dtermination de la permabilit se fait partir de la mesure du dbit d'coulement Q et du
H
gradient hydraulique
:
L
Q=KS

H
donc
L

K=

Q L

S H

Remarque : pour permettre de bloquer l'chantillon dans l'appareil, il est possible de fixer
chaque extrmit un matriau de trs forte permabilit et de trs faible paisseur. Nous
avons vu que la permabilit verticale d'un ensemble de couches de terrain s'exprimait par...

Hydrologie - Hydrogologie Septembre 2008 Vronique MERRIEN-SOUKATCHOFF

56
n

Kv =

e
i=1
n

ei

K
i=1

... si on a deux terrains dont l'un est de trs faible paisseur et de permabilit trs
e
e
importante : e i e et i .
Ki K

B.3.5.3.2 Permamtre charge variable

Dans ce type d'appareil la charge hydraulique applique au sommet de l'chantillon est
variable.
Pour une variation lmentaire de la charge dh, il est possible d'crire :

dh
: dbit d'coulement
dt

dh
S h(t )
= K
ou encore
s
dt
L
dh K S
dt
=
h
sL

h
S t
; ce qui
si t = 0 on a h = h0 ; ln = K
s L
h0
h
.S
devient ln = t en posant =
sL
h0

ou

encore h = h 0 e

Q = Q0 e

ou

comme

K S
Q = h
L

(loi de Maillet)

Figure 42 : Principe du
permamtre charge variable

B.3.6 Mesures in situ

Les mthodes de mesures in situ ne sont pas dveloppes dans ce polycopi, il faut cependant
rappeler que ce sont elles, en gnral, qui permettront d'valuer la permabilit des terrains en
place.
B.3.7 Ordre de grandeur de la permabilit

Pour fixer les ordres de grandeur, on rencontre frquemment les valeurs suivantes :
-

graviers, sables grossiers

K = 10-1 10-5 m/s

(Sables de Fontainebleau

K = 2 10-5 m/s)

sables fins

K = 10-5 10-6 m/s

silts

K = 10-6 10-8 m/s

marnes

K = 10-8 10-9 m/s

tourbe

K = 2 10-8

Hydrologie - Hydrogologie Septembre 2008 Vronique MERRIEN-SOUKATCHOFF

argiles

K <10-9 m/s

calcite

K = 10-11 m/s

granite (non fractur, non altr)

57

K = 10-11 m/s

B.3.8 Gnralisation en 3 dimensions

Jusque l nous avons suppos que les matriaux taient homognes et isotropes (mis part
dans le paragraphe B.3.4).
L'exprience de Darcy est ralise en dimension 1, mais lorsque nous avons crit :

V = grad p + g grad z , implicitement nous avons gnralis la loi 3 dimensions en

supposant que le terrain tait homogne et isotrope et que la loi tablie en dimension 1 tait
valable pour les 3 composante de l'espace i.e. que l'on pouvait crire

Vx = K

H
x

Vy = K

H
y

Vz = K

H
z

ou encore V = K grad H
Pour justifier cette gnralisation 3 dimensions, il est possible de raliser des calculs partir
des quations de Navier-Stokes qui montrent que les dplacements de fluide en milieu poreux
sont dus aux gradients de pressions et aux forces extrieures (ici la gravit).

Remarque :
On rappelle que les quations de Navier Stokes (conservation de la quantit de
mouvements) sont tablies pour l'coulement laminaire d'un fluide
incompressible

2u 2u 2u
p
+ 2 + 2 + 2 =
x
z
y
x

2v 2v 2v
p
+ 2 + 2 + 2 =
y
z
y
x

u
u
u
u
+u
+v
+ w
x
y
z
t

v
v
v
v
+ u + v + w
x
y
z
t

2w 2w 2w
p
+ 2 + 2 + 2 =
x
z
y
x

w
w
w
w
+u
+v
+ w
x
y
z
t

avec si on suit une particule dans son mouvement (mthode de Lagrange) qui
x
y
z
; v=
; w=
passe du point Pt0(x0,y0,z0) au point Pt(x,y,z) : u =
t
t
t
On sait qu'a priori un terrain n'est ni homogne, ni isotrope; il suffit de regarder une
succession de couches sdimentaires pour s'en convaincre. Dans une srie argilo-sableuse, la
permabilit horizontale sera suprieure la permabilit verticale. De mme dans des

Hydrologie - Hydrogologie Septembre 2008 Vronique MERRIEN-SOUKATCHOFF

58

alluvions, il existe en gnral des chenaux plus permables. Ces constatations conduisent
considrer la permabilit d'un point de vue mathmatique comme une proprit tensorielle.
=

On dfinit un tenseur de permabilit K que l'on admettra tre un tenseur du deuxime ordre
symtrique (c'est dire reprsent par une matrice 3*3 symtrique par rapport la diagonale
et ayant donc 6 coefficients indpendants).
Kxx
K = Kyx
K
zx
=

K xy
K yy
Kzy

K xz
Kxy = Kyx

K yz avec K xz = Kzy
Kzz
K yz = Kzy

La relation V = K grad H peut donc s'expliciter par :

Vx = K xx

1 p
H
H
H
1 p
1 p
Kxy
Kxz
= K xx
K xy
K xz
+ 1
g z
x
y
z
g x
g y

Vy = K yx

1 p
H
H
H
1 p
1 p
Kyy
Kyz
= Kyx
K yy
K yz
+ 1
g z
x
y
z
g x
g y

Vz = Kzx

1 p
H
H
H
1 p
1 p
Kzy
Kzz
= Kzx
K zy
Kzz
+ 1
g z
x
y
z
g x
g y

Il est possible de reprsenter ce tenseur dans un espace ayant comme axes les directions
principales du tenseur des permabilit, la matrice se rduit alors dans ce nouvel espace
(XYZ ou 123) 3 composantes diagonales:
=

K XYZ

K XX
= 0
0

0
KYY
0

0
0
KZZ

Mathmatiquement X, Y et Z sont les directions propres de la matrice et KXX, KYY, KZZ les
valeurs propres associes. Physiquement, X, Y et Z sont les directions pour lesquelles
l'coulement est parallle au gradient de charge.

Vx = K H

1 p
g x

Souvent le milieu est stratifi et on distingue

donc deux permabilits : une permabilit
1 p
horizontale (KH = Kxx = Kyy) et une Vy = K H g y
permabilit verticale (KV = Kzz) On a alors :
1 p
Vz = KV
+ 1
g z
Les paragraphes prcdents portent sur les coulements en milieux poreux "continus". Nous
reviendrons sur la possibilit d'coulement dans les milieux discontinus au B.11.

Hydrologie - Hydrogologie Septembre 2008 Vronique MERRIEN-SOUKATCHOFF

59

B.4 Les nappes

B.4.1 Conditions d'existence des nappes

Une nappe est une accumulation d'eau dans les pores d'un terrain permable.
B.4.1.1 Processus gnral

Par le jeu de la pesanteur, une partie de l'eau de pluie s'infiltre dans le sol, soit directement,
soit aprs circulation la surface de celui-ci. Nous tudierons plus en dtails au paragraphe
IV.2 l'alimentation des nappes.
Selon la permabilit des terrains rencontrs, elle descend une plus ou moins grande
profondeur. Cette circulation, approximativement verticale, est interrompue par la rencontre
d'un terrain de faible permabilit par rapport celle qui la prcde.
Sous des terrains permables, cette formation "impermable" reprsente en quelque sorte "le
fond du rcipient". L'eau s'y accumule en saturant l'ensemble des vides des terrains susjacents plus permables. Ainsi se constituent dans ces formations relativement permables
appeles aquifres (qui portent l'eau) des nappes.
Lorsque le "rcipient" est plein, il dborde vers l'extrieur ou en direction d'autres terrains
permables (cf. B.4.3).
B.4.1.2 Facteurs d'existence d'une nappe

L'existence d'une nappe est conditionne par la conjonction de trois facteurs :

- facteur lithologique : il doit exister une roche "rservoir" la fois poreuse et permable qui
constituera le terrain aquifre. Cette roche doit avoir sa base un mur impermable pour
soutenir la nappe ou plus simplement qui servira de fond "tanche" au rceptacle constitu
par le rservoir ;
- facteur d'alimentation : il faut que de l'eau puisse venir remplir les pores de la nappe ;
- facteur de structure : il est ncessaire d'avoir une structure favorable l'accumulation de
l'eau (une structure anticlinale est dfavorable).
Ces trois facteurs ont une influence variable selon les types de nappes rencontres. Les
conditions d'alimentation, lithologiques et structurales ne peuvent pas tre considres
isolment.

B.4.1.2.1 Cas d'une alimentation et d'une lithologie favorables :

Si l'eau parvient de faon abondante et continue une rgion constitue par un terrain trs
permable d'une part, au mur duquel se trouve un terrain trs impermable d'autre part
(conditions d'alimentation et lithologiques excellentes), une accumulation souterraine pourra
se former et persister.

B.4.1.2.2 Cas d'une lithologie et d'une structure favorables :

Si deux terrains, l'un trs permable et l'autre pratiquement impermable, sont disposs de
telle sorte que la limite de permabilit forme un creux (conditions lithologiques et
structurales excellentes), il suffira d'un faible apport d'eau, mme temporaire, pour qu'une
nappe, mme permanente, se forme.

60

Hydrologie - Hydrogologie Septembre 2008 Vronique MERRIEN-SOUKATCHOFF

B.4.1.2.3 Cas d'une alimentation et d'une structure favorables :

Si les apports d'eau sont abondants et continus dans une structure en creux (conditions
d'alimentation et de structure excellentes), il suffira d'une faible diffrence de permabilit
entre deux terrains pour qu'une partie de l'eau soit retenue et s'accumule au-dessus de la limite
de permabilit.
B.4.2 Alimentation

Les nappes peuvent avoir plusieurs modes d'alimentation :

-

par infiltration directe de la pluie lorsque l'aquifre est l'affleurement,

par drainage d'un cours d'eau de surface,

par dversement d'une nappe sus-jacente,

par drainance travers des terrains semi-permables lorsque deux nappes ne sont pas en
quilibre.

B.4.3 Exutoires

Si on reprend la vision schmatique du terrain permable rcipient qui se remplit par

alimentation verticale de la pluie, lorsque le rcipient est plein, l'eau s'coule de ce rcipient,
soit la surface topographique, formant ainsi une source ou une mergence, soit dans un
terrain sous-jacent
On appelle donc exutoires d'une nappe, les points privilgis o l'eau sort de la nappe. Les
exutoires sont de diffrents types.
B.4.3.1 Les sources

B.4.3.1.1 Classification des sources

On rencontre une source ou une mergence lorsque, l'affleurement du niveau aquifre, le
niveau pizomtrique est gal la cote du sol.
Si la sortie de l'eau est contracte, il s'agit d'une source ; le point de la sortie d'eau est appel
griffon de source. Si la sortie d'eau se fait sur une grande tendue, qu'il y a un suintement
diffus, l'exutoire est une mergence.
Les sources peuvent se classer selon leur position structurale :

Hydrologie - Hydrogologie Septembre 2008 Vronique MERRIEN-SOUKATCHOFF

61

B.4.3.1.1.1 Sources de dversement :

La totalit des filets liquides se meut en amont de la source, au-dessus du niveau de celle-ci.
Exemples :
S. de dversement

S. de dversement

a
Figure 43 : Nappe horizontale avec deux sources de dversement

S. de
dversement

b
Figure 44 : Nappe incline avec une source de dversement

62

Hydrologie - Hydrogologie Septembre 2008 Vronique MERRIEN-SOUKATCHOFF

B.4.3.1.1.2 Sources de dbordement :

Une partie au moins des filets liquides se meut en amont de la source, au-dessous du niveau
de celle-ci et l'eau se dplace dans le sens du pendage du terrain.
Exemples :

c
Figure 45 : Nappe de dbordement en tte de nappe captive

d
Figure 46 : Nappe de dbordement par faille

Hydrologie - Hydrogologie Septembre 2008 Vronique MERRIEN-SOUKATCHOFF

B.4.3.1.1.3 Sources de trop plein :

Mme dfinition que la source de dbordement, mais l'eau circule contre-pendage.
Exemples :

e
Figure 47 : Nappe de trop-plein en tte de nappe captive

f
Figure 48 : Double nappe de trop-plein

Trop-plein
Dversement

g
Figure 49 : Nappe de dversement et nappe de trop-plein jointes

63

64

Hydrologie - Hydrogologie Septembre 2008 Vronique MERRIEN-SOUKATCHOFF

B.4.3.1.2 Fonctionnement des sources

B.4.3.2 Les exutoires cachs

Les exutoires cachs sont par exemple le dversement souterrain dans une autre nappe,
l'alimentation du rseau de surface ou la drainance.
Remarques :
- Nous n'avons mentionn ici que les exutoires naturels ; il peut exister des exutoires
artificiels tels que forages pomps et drains.

- Nous avons dcrit, dans les derniers paragraphes, les points d'alimentation et les exutoires
d'une nappe. Dans un systme aquifre donn, il s'tablit un quilibre entre le volume d'eau
entrant, et celui qui en sort par les exutoires naturels ou artificiels.
B.4.4 Classification des nappes
B.4.4.1 Critres gologiques

B.4.4.1.1 Nappes de terrains sdimentaires stratifis

Les terrains sdimentaires stratifis sont des terrains de grande extension latrale qui forment
des rservoirs d'importance rgionale. Ces nappes peuvent tre classes en fonction de la
structure (synclinal, monoclinal).

B.4.4.1.2 Nappes de terrains sdimentaires mal (ou non) stratifis

Ce sont des nappes de terrains tertiaires ou quaternaires d'extension limite. Il s'agit des
nappes d'alluvions fluviatiles, d'boulis, de dpts oliens ou glaciaires.
Les alluvions forment souvent de petites units hydrogologiquement discontinues entre les
mandres d'une rivire.

Infiltration de la rivire vers la nappe

Alimentation de la rivire par la nappe

Figure 50

B.4.4.1.3 Nappes de terrains cristallins ou ruptifs

Il ne peut exister une nappe dans ce type de terrain que dans deux cas :
- terrain fissur (ce qui donne une porosit maximum de 4 %),
- terrain altr (un granite sain ne peut pas contenir une nappe).

Hydrologie - Hydrogologie Septembre 2008 Vronique MERRIEN-SOUKATCHOFF

65

B.4.4.2 Critres hydrodynamiques

On distingue deux types de nappe sous l'angle de la pizomtrie.

B.4.4.2.1 Nappe libre

On appelle nappe libre, une nappe dont la surface pizomtrique se confond avec la surface
(ou le toit) de la nappe.
B.4.4.2.1.1 Nappes perches ou non soutenues
La nappe s'coule dans ce cas l'air, c'est--dire par des sources. Toutes ces nappes sont
convexes et elles peuvent tre classes en fonction des exutoires :
-

nappe de dversement,

nappe de dbordement,

nappe de trop plein.

B.4.4.2.1.2 Nappes soutenues

La nappe s'coule dans une autre nappe souterraine, dans un lac, dans une rivire ou dans la
mer. Ces nappes peuvent tre soit convexes soit concaves.
Exemple :
Terrasse d'alluvions anciennes
Terrasse d'alluvions rcentes

Rivire

Nappe d'alluvions
anciennes soutenue
par la nappe des
alluvions rcentes

Nappe dans alluvions rcentes

soutenue par la rivire

Substratum
permabilit plus faible

Figure 51 : Nappes d'alluvions

B.4.4.2.2 Nappe captive

On appelle nappe captive, une nappe dont le toit est un niveau infrieur la surface
pizomtrique. Ceci suppose que la couche situe au toit de l'aquifre soit "impermable",
mais cette condition n'est pas suffisante. En effet, si le terrain permable est suffisamment
aliment la nappe est captive ; si cette alimentation est trop faible, elle reste libre.

Hydrologie - Hydrogologie Septembre 2008 Vronique MERRIEN-SOUKATCHOFF

66

(b)

(a)
Nappe captive

(I1)
Source

(I1)

(P)

SP

SP

(P)

Nappe libre

(I2)

SP : surface pizomtrique
I1, I2 : terrains faible permabilit
P : terrain forte permabilit

(I2)

Figure 52 : nappe captive (a) puis libre (b)

Hydrologie - Hydrogologie Septembre 2008 Vronique MERRIEN-SOUKATCHOFF

67

B.5 Etude de quelques coulement dans les nappes

B.5.1 Paramtres dont dpend l'coulement. Rseau d'coulement
B.5.1.1 Transmissivit

Si l'on considre un terrain d'une

permabilit donne K, le dbit passant
travers une section de ce terrain sera
fonction de la permabilit mais aussi de la
surface de la section traverse. On appelle
transmissivit le produit de la permabilit
par l'paisseur de la nappe.
T=Ke

L=1

(en m2/s)

Figure 53 : transmissivit

B.5.1.2 Coefficient d'emmagasinement

Lorsque qu'une nappe est en rgime transitoire c'est--dire quant au moins un des paramtres
varie en fonction du temps, la seule transmissivit ne suffit plus caractriser le milieu
aquifre. En effet, lorsque le niveau pizomtrique d'une nappe baisse, il y a dpart d'eau.
Pour caractriser ce phnomne, on utilise la notion de coefficient d'emmagasinement.
Le coefficient d'emmagasinement S est le
volume d'eau que l'on peut extraire d'une
tranche de 1 m2 de surface horizontale
pour une baisse de pizomtre de 1 m.
D'aprs cette dfinition, on constate que
S est sans dimension (m3/m2/m).

S = 1 m2
H

Dans le cas de nappes libres, le

coefficient d'emmagasinement reprsente
la porosit efficace. S est alors de l'ordre
de quelques %. Par contre, pour une
nappe captive, S dpend du coefficient de
compressibilit du fluide et du terrain. S
est alors beaucoup plus faible, environ
10-5 10-6.

Figure 54 : coefficient d'emmagasinement

Si e est l'paisseur de la nappe on dfinit le coefficient d'emmagasinement spcifique Ss par :

Ss =

S
= g ( + n )
e
dVT

: compressibilit du milieu poreux ( =

VT

d e

; VT volume total de terrain, d e

variation de la contrainte effective)

: compressibilit du fluide (4,8 10-10 m2/N ou Pa-1 pour l'eau)

68

Hydrologie - Hydrogologie Septembre 2008 Vronique MERRIEN-SOUKATCHOFF

Le coefficient d'emmagasinement spcifique reprsente le volume de fluide que l'on peut

extraire d'un volume de terrain donn, lors d'une baisse de pizomtre de 1 m. Sa dimension
est [L]-1.
Remarque : en premire approximation la compressibilit du milieu peut tre assimile
l'inverse du module de dformation volumique (Bulk modulus)

K=

E
3 (1 2 )

1 3 (1 2 )
=
K
E

Pour une roche trs peu poreuse comme un granite tel que n =1,5%, E = 50 GPa et = 0,3 :
SS = 3 10-7 m-1
B.5.1.3 Surface pizomtrique et isopizes

On appellera surface pizomtrique d'une nappe, la surface, relle ou fictive, sur laquelle la
pression des particules liquide est gale la pression atmosphrique. C'est encore le lieu de
l'ensemble des hauteurs pizomtriques de la surface de la nappe.
Une ligne isopize ou isopize, est le lieu des points ayant mme niveau pizomtrique. Une
ligne d'isocharge hydraulique est une quipotentielle. Dans ce qui suit, nous confondrons en
gnral isopizes, isocharges hydrauliques ou quipotentielles car les vitesses considres
sont faibles.
La connaissance des isocharges hydrauliques permettra de dterminer le sens des
coulements dans une nappe.
B.5.1.4 Lignes de courants

Nous rappelons que la ligne de courant est la ligne tangente au vecteur vitesse en chacun de
ces points l'instant considr.
Si le milieu est homogne et isotrope les lignes de courant seront perpendiculaires aux
isopizes. La connaissance de la pizomtrie permettra donc de tracer les lignes de courants.
Si le milieu est anisotrope les lignes de courant feront un angle diffrent de 90 avec les
isopizes (sauf si les isopizes sont orthogonales aux axes principaux de permabilit).
Le trac des isopizes et des lignes de courant est appel rseau d'coulement.
B.5.2 Effet des variations de la transmissivit

Supposons l'coulement d'une nappe qui ne reoit pas d'apports intermdiaires. Pour un mme
dbit, la transmissivit et le gradient hydraulique varient en sens inverse.

Hydrologie - Hydrogologie Septembre 2008 Vronique MERRIEN-SOUKATCHOFF

69

Isopizes

T1 >> T2

Figure 55 : variation de la transmissivit

B.5.3 Ecoulement entre deux tranches

Supposons une tranche (tranche de captage), creuse dans un massif permable renfermant
une nappe, jusqu'au substratum impermable horizontal. Un foss parallle (tranche
d'alimentation), situ une distance R, atteignant le fond impermable et rempli d'eau sur une
hauteur constante H1, assure la ralimentation.
Un dbit constant Q maintient une tranche d'eau de hauteur H2 dans la tranche de captage,
entranant un coulement permanent des eaux souterraines dans la tranche de terrain
considr.
B.5.3.1 Cas d'une nappe libre

Dsignons par h l'ordonne d'un point quelconque de la surface (galement surface

pizomtrique) situe la distance x de la tranche de captage. Soit l le chemin parcouru par
les filets liquides. Dans la section d'abscisse x, le dbit par mtre linaire peut s'crire :
Q=-Kh

dh
dl

soit

Q dl = - K h dh

En intgrant avec les conditions aux limites suivantes : pour l = 0, h = H2 et pour l = L,

h = H1 :

H 12 H 22
2L

Tranche
d'alimentation

mais il est difficile de calculer L.

Tranche de
captage

Q=K

Plan

70

Hydrologie - Hydrogologie Septembre 2008 Vronique MERRIEN-SOUKATCHOFF

y
A' A

Courbe
de Dupuit

h
Pente dh/dl

B' B

O
x

R
Coupe AB

Figure 56 : Ecoulement entre 2 tranches (nappe libre)

Formule approche de Dupuit : Pour simplifier, Dupuit remplace l'arc lmentaire dl par sa
projection horizontale dx, ce qui est une approximation valable si on ne considre que les
parties de la nappe dont la surface n'a qu'une faible inclinaison (dans les zones forte
inclinaison, comme au voisinage de la tranche de captage, de grosses erreurs sont ainsi
introduites). On a alors :

Q dx = - K h dh
qui par intgration nous donne :

H 12 H 22
Q=K
2 R
La surface de la nappe est donc approximativement parabolique.
B.5.3.2 Cas d'une nappe captive

Si nous reprenons le mme coulement dans le cas d'une nappe captive d'paisseur constante
e, en gardant les mmes notations, nous pouvons crire :
Q dx = - K e dh
y
Q

Courbe pizomtrique

dh
H

O
x

Figure 57 : Ecoulement entre 2 tranches (nappe captive)

soit en intgrant avec comme conditions aux limites : pour x = 0, h = H2 et pour x = R, h =

H1 :

Hydrologie - Hydrogologie Septembre 2008 Vronique MERRIEN-SOUKATCHOFF

Q=Ke

71

H1 H 2
R

La surface pizomtrique de la nappe est donc une droite.

B.5.4 Ecoulement radial circulaire en rgime permanent
B.5.4.1 Nappe captive

Examinons le schma suivant : un puits o l'on pompe un dbit

Q constant et un pizomtre (c'est--dire un forage o l'on
pratique une prise de pression statique), situ une distance r du
puits, o l'on observe un niveau pizomtrique h.
Considrons les hypothses suivantes:
-

l'aquifre a une paisseur e constante,

le substratum est suppos horizontal,

Le milieu infini (ou trs grand dans toutes les directions),

le puits est quip d'une crpine sur la totalit de l'aquifre,

le dbit Q de pompage est constant.

R
Surface du terrain

Puits
Hauteur pizomtrique initiale de la nappe
Hauteur pizomtrique aprs
un temps de pompage trs long

Surface isopize
Ligne de courant

Substratum impermable
Figure 58 : coulement radial circulaire
La vitesse la distance r est V = K

dh
dr

tube
pizomtrique

72

Hydrologie - Hydrogologie Septembre 2008 Vronique MERRIEN-SOUKATCHOFF

La surface offerte l'coulement est S = 2 r e

Le dbit est donc Q = VS = 2 eK
soit en intgrant : ln

dh
dr
K
d'o
= 2e dh (1)
dr
r
Q

r1 2Ke
(h 1 h 2 ) or T = Ke
=
r2
Q

Si h0 est la charge hydraulique initiale, on note s = h 0 - h (rabattement)

En supposant qu'il existe une distance RA (appel rayon d'action) partir de laquelle s 0,
l'quation (2) devient :
2,3 log10

R A 2T
=
s
r
Q

ou encore

s=

R
R
Q
0.366 Q
ln A =
log A (3)
2 T
T
r
r

Cette relation est galement appele formule de Dupuit

Remarque : cette notion de rayon d'action est contestable (et contest) car li la signification
de s 0 . Il vaudrait mieux le dfinir par rapport une valeur absolue, par exemple : s = 10-2.
Cette notion est galement fonction du moment o l'on considre que l'on a atteint le rgime
permanent. En effet le rabattement est fonction du temps de pompage et comme nous le
verrons plus loin ( B.7.3), l'expression du rabattement en rgime permanent peut aussi tre
dduite de l'expression plus gnrale du rabattement en rgime transitoire, quand les temps de
pompage sont trs longs.
Les courbes caractristiques [Q = f (s)] sont des droites.
B.5.4.2 Nappe libre

Nous supposerons que les conditions sont les mmes que prcdemment, mais dans ce cas, la
nappe est libre et a une hauteur initiale constante h0.
La vitesse la distance r est V = K

dh
dr

dh
, dl tant le trajet du filet liquide.
dl
Ici, nous faisons l'approximation de Dupuit en assimilant l'arc dl sa projection horizontale
dr. Ceci est admissible loin de l'axe du puits, mais ne l'est plus au voisinage.

Remarque : En toute rigueur nous devrions crire V = K

La surface offerte l'coulement est S = 2rh (si h est la hauteur de la nappe une distance r
du puits).
Le dbit est donc Q = V S = 2 r h K

dh
d'o
dr

dr 2 h K
=
dh
r
Q
dr 2 K
=
h dh (4)
Q
r
soit en intgrant entre deux distances r1 et r2 o les niveaux pizomtriques sont h1 et h2...

Hydrologie - Hydrogologie Septembre 2008 Vronique MERRIEN-SOUKATCHOFF

ln

73

r1
r
K 1
K

= 2 h12 h22 ou ln 1 = ( h12 h22 ) (5)

r2
Q 2
r2
Q

... ou encore, en supposant qu'il existe une distance RA (dite rayon d'action) partir de
laquelle le niveau pizomtrique est le mme que le niveau initial h0.
(h 02 - h 2 ) =

R
Q
ln A
K
r

(h 02 - h 2 ) = 0,732

R
Q
log10 A
K
r

or

h 02 - h 2 = (h 0 - h) (h 0 + h) = (h 0 - h) (2h 0 + h - h 0 )

et

s = h0 - h

donc

s(2h 0 - s) = 0.732

ou encore :

R
Q
log A
K
r

2h 0 s s 2 = 0.732

R
Q
log A (6)
K
r

Les courbes caractristiques [Q = f (s)] sont dans ce cas des paraboles.

Remarque : En supposant que le dbit Q est faible et que l'extension de la nappe est
importante, pour des pizomtres loigns du puits, il est possible de considrer que le terme
2h0 est nettement suprieur s (2h0 >> s).

Dans ce cas le terme : s(2h 0 - s) se simplifie en 2h 0 s et l'on retrouve... s =

R
0.366 Q
log A
Kh 0
r

... c'est--dire une formulation analogue celle d'une nappe captive.

D'aprs http://www.e-atlantide.com/thematiques/shadok/index.htm#dessins

Hydrologie - Hydrogologie Septembre 2008 Vronique MERRIEN-SOUKATCHOFF

74

B.6 Cas Gnral : quation fondamentale de l'hydrodynamique

B.6.1 Etablissement de l'quation de diffusivit

Les quelques exemples que nous venons de citer sont en fait des cas particuliers de rsolution
de l'quation fondamentale de l'hydrodynamique (ou de continuit ou de diffusivit).
Problme : on se propose d'tablir une relation qui permette de dcrire le mouvement d'une
particule d'eau en fonction de sa position, du temps, des caractristiques physiques du fluide
et de l'encaissant (le terrain).

On dispose de 3 types de relations :

* les quations d'tat isotherme du fluide :

= e
0

p p

: coefficient de compressibilit du fluide [M-1 L T2]

* L'quation de continuit ("rien ne se perd, rien de se cre")

: volume de frontire

M : masse du fluide
n : porosit

M
n
=

t t
uur ( n)
+ q = 0
div( VD ) +
t

Figure 59

r
r
* Les quations de la dynamique ( F = m )
Les quations de la dynamique se traduisent pour les fluides visqueux, dont les coefficients
de viscosit sont supposs constants, en quation de Navier Stokes.
A partir de ces 3 types d'quations on peut tablir l'quation de diffusivit

Hydrologie - Hydrogologie Septembre 2008 Vronique MERRIEN-SOUKATCHOFF

75

Equation de diffusivit :

div K grad (H ) = S s

H
+q
t

Ss : coefficient d'emmagasinement spcifique

S : coefficient d'emmagasinement (S : storage)

si la permabilit horizontale est

constante sur l'paisseur de la nappe : S =

div T grad (H ) = S

H
+Q
t

si T est isotrope l'quation devient

2H =

S H Q
+
T t T

toit

s
substratum

dz

Q > 0 dbit prlev par unit de surface de la nappe

Q < 0 dbit inject par unit de surface de la nappe
(donc en m/s) : par exemple pluie

qui est la forme de l'quation de

diffusivit que nous utiliserons dans la
suite.
2
En rgime permanent l'quation devient : H = 0

B.6.2 Proprits de l'quation de diffusivit

B.6.2.1 Unicit de la solution (Cf. Cours TCS11, Maths 1A et SG232 EDP))
En rgime permanent ou pour une nappe captive en rgime transitoire condition que S et T
soient constants et les CL, l'quation de diffusivit est linaire en h.

Dans ces cas, si D est un domaine d'intgration donn de l'quation de diffusivit, muni de
conditions aux limites et initiales, on dmontre, que si h vrifie ces conditions aux limites et
si h vrifie l'quation de diffusivit, h est, en gnral, la solution unique du problme (en
rgime permanent ou transitoire).
B.6.2.2 Principe de superposition

Dans les mmes conditions, dans un domaine D l'quation de diffusivit est linaire en h et en
q. Si (h1, q1) et (h2, q2) sont 2 solutions particulires de l'quation de diffusivit vrifiant des
conditions aux limites donnes, alors h1 + h2 est solution de la mme quation, avec les
dbits q1 + q 2 et les conditions aux limitent qui en rsultent

B.7 Solutions de l'quation de diffusivit lors de pompage dans un

puits
Un des ouvrages le plus frquent en hydraulique souterraine consiste en un puits atteignant un
aquifre. Ces puits permettent de prlever un dbit Q constant ou variable dans la nappe (ou
une partie de la nappe) ainsi atteinte.
Pour connatre les caractristiques hydrauliques des rservoirs il est frquent d'effectuer des
"pompages d'essais" qui consiste enregistrer les variations de la charge hydraulique (ou du
dbit) au puits et dans des pizomtres proches du puits en fonction du temps. Ces variations
sont interprtes selon diffrents modles thoriques de comportement. Dans ce qui suit nous
prsentons quelques-uns de ces modles.

Hydrologie - Hydrogologie Septembre 2008 Vronique MERRIEN-SOUKATCHOFF

76

B.7.1 Equation de diffusivit en coordonnes radiales

En coordonnes radiales l'quation de diffusivit s'crit :

1 H 2 H S H
=
r
+
r r r z 2
T t

B.7.2 Solution de Theis

Theis a propos une solution cette quation dans le cas d'un milieu infini, homogne et
isotrope de transmissivit constante dans le temps et dans l'espace (nappe captive ou nappe
libre peu rabattue), d'un forage captant la nappe sur toute son paisseur d'une injection ou
d'un pompage dbit constant, dans un puits infiniment petit :
s=

Q
4 T

W (u ) =

r2 S
ey
avec
u
=
dy
4 T t
4 T u y
Q

t : temps
r : distance radiale partir du centre du puits
s : rabattement de la nappe (diffrence entre le niveau pizomtrique initial de la nappe et son
niveau aprs un temps t de pompage, observe la distance r du puits)
Q : dbit de pompage
T : Transmissivit (produit de la permabilit par l'paisseur de la nappe)
S : coefficient d'emmagasinement
B.7.3 Solution de Jacob

Une solution approche de l'quation prcdente peut tre propose :

s=

2,25 T t
ln
(formule de Jacob)
4 T S r 2

Cette approximation reste valable pour des valeurs de u infrieures 10-1 ; c'est--dire, en
pratique, pour des temps de pompage long et des distances proche du puits de pompage.
En effet on peut crire que :
n

1
(
ey
y y2 y3
1) u n
u y dy u y 1 + 2 3! + 4! ...dy ln(u )
n (n!)
n =1

1 1
avec = lim n 1 + + + ... + ln n 0,5772156... constante d'Euler
n

2 3
pour des faibles valeurs de u on peut limiter le dveloppement au 2me ordre et crire :

ey
1
u y dy ln(u ) ln u 0,5772 ou encore

Q
2,25 T t
ln 0,5772 =
ln
soit en log10 :

2
4 T u

4 T S r

Hydrologie - Hydrogologie Septembre 2008 Vronique MERRIEN-SOUKATCHOFF

77

0,183 Q
2,25 T t
log10

2
T

S r

Remarques :

r2 S
correspondent des temps de pompage longs ou des
4 T t
points proches du puits (r faible)

1. Les faibles valeurs de u =

2. Le terme

2,25 T t
S

a les dimensions d'une longueur au carr. Si on pose

2,25 T t
0,366 Q
L(t )
log10
= L2 (t ) on peut crire le rabattement sous la forme : s

S
T
r
ce qui est proche de la formule de Dupuit crite au B.5.4.1 page 72 en posant :
L(t ) = R A .
On comprend aussi que cette notion de rayon d'action soit contestable, puisque ce rayon
dpend du temps de pompage. Cependant comme les nappes ne sont jamais infinies il existe
souvent une certaine distance laquelle il existe une ralimentation de la nappe.
B.7.4 Pompage d'essai

Un essai de pompage classique consiste mesurer (dans le puits de pompage ou dans un

pizomtre) les rabattements induits dans la nappe par un pompage de quelques heures
(ventuellement plus long) afin de dterminer les paramtres hydrauliques de l'aquifre : S et
T (les autres paramtres Q, t et r tant a priori mesurs). La solution de Theis est souvent
reprsente sous forme d'abaque dans des axes log-log, pour faciliter l'interprtation
graphique des essais de pompage.
Des essais plus spcifiques sont utiliss dans des terrains de faible permabilit.
Les essais sont interprts en comparant diffrents modles thoriques aux essais rels et les
proprits des terrains investigus sont obtenus par calage des courbes thoriques aux courbes
exprimentales. Les diffrents modles thoriques correspondent des gomodles diffrents.
B.7.4.1 Interprtation graphique de la mthode de Theis

L'interprtation graphique du pompage consiste ajuster la courbe exprimentale la courbe

thorique. En effet :
soit en log :

s=
u=

Q
4 T

r2 S
4 T t

W (u )

4 T

log[W (u )] = log s + log

Q
r2 S

log u = log t + log

4 T

La courbe exprimentale log(s), -log(t) peut donc se dduire de la courbe thorique

log[W (u )] log(u) 2 translations prs :

78

Hydrologie - Hydrogologie Septembre 2008 Vronique MERRIEN-SOUKATCHOFF

4 T
selon l'axe des x ou des s ou des W(u)
une translation de log
Q

r2 S
selon l'axe t ou u
une translation de log
4 T

4 T
permet d'en dduire la valeur de la
La valeur de la translation selon l'axe des x : log
Q
transmissivit.
r2 S
par dduction permettra d'en dduire la valeur du
La translation selon l'axe y : log
4 T
coefficient d'emmagasinement S.

Aprs superposition des courbes, on choisit donc un point arbitraire pour lequel on note les
coordonnes W(u), u de la courbe thorique et les coordonnes s, t de la courbe
Q
W (u ) et
exprimentale. L'introduction de ces valeurs dans les quations s =
4 T
r2 S
permet d'obtenir la valeur de la transmissivit et du coefficient
u=
4 T t
d'emmagasinement

Log10 HWHuLL
2
1
0
-1
-2
-3
-4

Log10 H1uL

Figure 60 : courbe thorique de Theis

B.7.4.2 Interprtation graphique de la mthode de Jacob

On peut tracer les rabattements en fonction du temps sur du papier semi-log. Les points
0,183 Q
. On peut donc en dduire la transmissivit.
s'alignent selon une droite de pente
T
S r2
permet d'en dduire le coefficient d'emmagasinement
L'abscisse l'origine t 0 =
2,25 T

Hydrologie - Hydrogologie Septembre 2008 Vronique MERRIEN-SOUKATCHOFF

(quand s =

79

2, 25 T t
2, 25 T t
= 1)
ln
=0
2
S r2
4 T S r

Remarque : les points correspondant des temps cours sont mal aligns car ils correspondent
r2 S
des valeurs de u =
> 10 1 . L'abscisse l'origine permet d'en dduire S.
4 T t
B.7.5 Remonte de la nappe

Si aprs pompage on observe la remonte du niveau de la nappe, il est possible d'interprter

les rabattements observs qualifis de rabattements rsiduels partir du principe de
superposition. En effet l'arrt du pompage peut tre considr comme la somme d'un
pompage dbit Q 'et d'une injection dbit -Q. le rabattement rsiduel s' s'exprime donc
sous la forme :

s' =

Q
4 T

[W (u ) W (u ' )] avec u =

r2 S
r2 S
t
et u ' =
=u
4 T t
4 T t'
t'

t : temps depuis le dbut du pompage

t' : temps depuis l'arrt du pompage : t = t p + t ' avec tp dure du pompage
Le rabattement rsiduel peut donc s'crire :

s' =

t
)
W (u ) W (u
t t p
4 T
Q

L'quation mise sous cette forme montre que le retour l'tat initial (rabattement nul) ne se
fera que pour un temps t trs grand par rapport au temps de pompage.
B.7.6 Mthode des paliers enchans

Il est possible d'interprter un pompage effectu avec une srie de paliers de dbits, sans
attendre entre chaque palier la remonte de la nappe.

"Je pompe donc je suis"

d'aprs
http://www.eatlantide.com/thematiques/shadok/index.htm#dessins

B.7.7 Quelques cas particuliers de nappe

Nous avons vu que la mthode de Theis permettait d'interprter le pompage dans une nappe
infinie, horizontale, non ralimente.
Il existe d'autres solutions thoriques dans des cas plus complexes.

80

Hydrologie - Hydrogologie Septembre 2008 Vronique MERRIEN-SOUKATCHOFF

B.7.7.1 Drainance
Il s'agit d'intgrer l'quation de diffusivit en prenant en compte les pontes.

B.7.7.1.1 Cas des aquifres superposs semi-permable : Schma de Hantush

La solution de Hantush, propose un modle thorique lorsque l'une des 2 pontes de l'aquifre
est constitu d'un horizon semi-permable.

B.7.7.1.2 Solution de Boulton

B.7.7.2 Ralimentation latrale

B.7.7.2.1 Mthode des puits images

La mthode des puits images est couramment utilise pour analyser des coulements situs
prs d'une limite impermable ou une limite charge impose. Cette mthode est base sur le
principe de superposition expos au B.6.2.2 page 75.

Hydrologie - Hydrogologie Septembre 2008 Vronique MERRIEN-SOUKATCHOFF

81

B.7.7.2.1.1 Principe de superposition

O'
X

Considrons 2 puits de centre O et O' pompant chacun avec un O

X
dbit Q et Q'. D'aprs le principe de superposition le rabattement
en un point M peut tre obtenu en additionnant les rabattements
que provoqueraient chacun de ces puits sparment.

r
r'

Si les conditions du milieu et du pompage permettent d'utiliser la

solution de Theis, le rabattement au point M s'crira :

s=

Q
4 T

W (u ) +

Q'
W (u ' )
4 T

M
Figure 61 : principe de
superposition

B.7.7.2.1.2 Limite potentiel impos

Supposons maintenant que dans le puits O' on injecte un

dbit Q (ou ce qui revient au mme, on pompe avec un
dbit -Q). Le rabattement en M s'crira :
Q
[W (u ) W (u ' )]
s=
4 T

O' Si on s'intresse aux points M situs gal distance de O et

de O' (c'est dire sur la mdiatrice OO'), pour ces points
r=r',
le
rabattement
s'crit
alors
:
r
Q
[W (u ) W (u )] = 0 ou exprim autrement, la
s=
r'
4 T
charge au niveau de ces points est constante, puisqu'il n'y a
M
pas d'volution du rabattement en fonction du temps. Nous
avons construit une solution au problme suivant : un puits
Figure 62 : limite charge
O' situ une distance d d'une limite rectiligne charge
impose
impos (ou potentiel constant). Cette solution est unique
(principe d'unicit) puisqu'elle vrifie les conditions aux
limites de l'quation de diffusivit ( condition qu'au forage
O le rayon r0 soit ngligeable devant la distance 2d)

B.7.7.2.1.3 Limite flux constant

Supposons que dans le puits O' on pompe avec le mme dbit Q qu'en O. Le rabattement au
niveau de la mdiatrice s'crira : s =

Q
4 T

[W (u ) + W (u )] =

Q
2 T

W (u ) . On peut

H
= 0 On a donc l la solution au problme d'un
x
forage situ une distance d d'une ligne rectiligne flux nul.

dmontrer que le long de cette ligne

B.7.7.2.1.4 Plusieurs limites

Par ce principe des images on peut reprsenter un problme plusieurs limites en multipliant
les puits images.

Hydrologie - Hydrogologie Septembre 2008 Vronique MERRIEN-SOUKATCHOFF

82

B.7.7.2.2 Application l'interprtation d'un pompage avec ralimentation latrale du puits

B.7.8 Essais Lugeon et Lefranc

Ces essais sont plutt utiliss en gotechnique.

B.8 Les essais spcifiques :

Dans ce paragraphe nous voquerons quelques tests plus spcifiques qui sont utiliss pour
dterminer des faibles permabilits ou des permabilits dans les milieux saturs.
Pulse test et slug test, essais en milieux non saturs (double anneaux).
Les slug et pulse se pratiquent gnralement dans une portion de forage isol par un systme
de double obturateur qui dlimite une chambre d'injection.
B.8.1 Essais par chocs hydrauliques : Pulse test

Le pulse test consiste appliquer une impulsion de pression trs brve en un point d'un
forage, entre deux obturateurs et observer en ce mme point la rponse en rgime
transitoire. L'interprtation varie en fonction de la forme de la cavit o lieu l'impulsion
(cylindre ou sphre). Cet essai permet d'estimer la transmissivit T et avec moins de prcision
le coefficient d'emmagasinement S.
Les pulse tests sont gnralement interprts en utilisant les courbes types tablies par
Cooper et al. (1967). Ces courbes type sont fonctions de 3 paramtres , et le temps t.
Le coefficient sans dimension, dfinit la forme de la courbe et est donn par :

rw2 S
avec
Cg

rw : rayon du forage

(m)

: coefficient sans dimension

(-)

S : coefficient d'emmagasinement

(-)

C : coefficient d'emmagasinement de la cavit de pompage ou injection

(m3/Pa)

: masse volumique du fluide

(Kg/m3)

g : gravit (9,81)

(m/s2)

Deux mthodes d'interprtation peuvent tre utilises :

-

l'ajustement de la courbe exprimental aux courbes types permet de dduire le coefficient

d'emmagasinement. Les valeurs de obtenues par ce type d'ajustement conduisent
souvent des valeurs irralistes du coefficient d'emmagasinement

Le calcul des valeurs de partir d'une estimation du coefficient d'emmagasinement bas

Hydrologie - Hydrogologie Septembre 2008 Vronique MERRIEN-SOUKATCHOFF

83

sur les proprits de la roche et du fluide (cf. B.5.1.2 page 67) et l'ajustement de la courbe
exprimental aux courbes types avec des formes prdfinies
La transmissivit T est calcule partir de l'ajustement de la courbe exprimental aux courbes
types par l'expression :

T=

C g
t

avec

: paramtre sans dimension

(-)

t : temps

(s)

Le paramtre C peut tre mesur pendant le pulse test en mesurant les variations de
pression lors de l'extraction d'un volume d'eau en supposant

La ralisation d'un pulse test, n'est pas possible si la permabilit est trop lev (il est alors
impossible d'appliquer la brve impulsion de pression, le fluide pntrant instantanment
dans la roche. Dans ce cas, il est possible d'effectuer un test d'injection dbit constant entre
obturateurs
B.8.2 Slug test (dbit variable partir d'une injection unique)

Le dispositif de ralisation d'un slug test est quivalent celui d'un pulse, mais convient des
transmissivits moyennes faibles. Pour ce test, on impose galement une pression donne,
mais aprs avoir fait varier la charge hydraulique dans le tubage raccord la zone de test, la
vanne entre la zone de test et le tubage reste ouverte, de sorte que, selon la diffrence entre la
roche et le tubage on observera un coulement d'eau du tubage dans la roche ou inversement.
Il est possible de calculer le dbit d'eau qui afflue ou sort, partir de la variation de pression
dans le tubage.

B.9 Mthodes analytiques de rsolution de l'quation de

diffusivit : fonctions harmoniques
B.10 Mthodes numriques de rsolution de l'quation de
diffusivit
La modlisation numrique des coulements souterrains permet de rsoudre l'quation de
diffusivit par approximation numrique.
Le problme consiste trouver la rpartition des charges hydrauliques dans un domaine dont
on connat la gomtrie, la permabilit, le coefficient d'emmagasinement (s'il s'agit d'un
problme transitoire), les conditions aux limites et initiales (pour un problme transitoire).
La reprsentativit de la modlisation dpendra de la plus ou moins bonne approximation de
la ralit et donc des hypothses plus ou moins restrictives qui seront faites. Parmi les
hypothses qui simplifient la rsolution numrique on peut citer :
-

milieu homogne et isotrope ;

gomtrie simplifie ;

conditions aux limites ;

eau incompressible ;

pas de dpendance par rapport au temps (rgime permanent)

84

Hydrologie - Hydrogologie Septembre 2008 Vronique MERRIEN-SOUKATCHOFF

Les modles, reprsenteront d'autant mieux la ralit qu'ils ne seront pas limits par ces
hypothses restrictives, mais ils seront de plus en plus complexes et lourds mettre en uvre,
ds lors que l'on s'affranchira de ces hypothses. De plus, plus ils prendront en compte la
complexit du milieu naturel, plus la quantit d'information ncessaire pour les "alimenter" et
les caler sera importante.
La prcision de la modlisation sera, entre autre, lie la finesse du maillage, qui se traduira
par une lourdeur informatique plus grande.
Dans ce qui suit nous considrerons 2 types de modles numriques : les modles diffrences
finies et les modles lments finis. Ces 2 types de modlisation imposent de dcouper
l'espace en lments en utilisant des points nodaux, situ soit aux extrmits des lments
dans le cas des lments finis ou des diffrences finies (mesh-centered nodes) soit au centre
des lments (block-centered nodes) dans le cas des diffrences finies centres.
Dans les 2 types de mthode, il s'agit de rsoudre un systme d'quations avec pour inconnues
la charge en un nombre fini de points.
Exemple : Un aquifre bord par une rivire. L'aquifre est aliment en surface par des
prcipitations. Horizontalement les seuls changes sont ceux avec la rivire.

Affleurement terrain
impermable

Lac

Rivire

Faille

Figure 63 : exemple de problme modliser

B.10.1 Equations rsoudre

2
En rgime permanent : quation de Laplace H = 0 ou div K grad (H ) = 0 , K
matricielle

Hydrologie - Hydrogologie Septembre 2008 Vronique MERRIEN-SOUKATCHOFF

En rgime transitoire div K grad (H ) =

85

S H Q
+
T t T

B.10.2 Conditions aux limites :

Pour pouvoir rsoudre les quations prcdentes, il est ncessaire de prciser les conditions
aux limites qui peuvent tre de diffrents types :
a) Charge connue le long de la frontire (condition de Dirichlet)
b) Flux (dbit) connu le long de la frontire (condition de Neumann)
c) combinaison de (a) et (b) (condition mixte)
B.10.3 Principe des diffrences finis :

L'espace est dcoup en une grille de points distants de x et .y (cf. Figure 64).

Figure 64 : Diffrences finies avec nuds centrs sur la grille (d'aprs Wang et
Anderson [33])

Si on connat la fonction charge hydraulique H en 2 points a et c, la valeur de H au point i

milieu de a et c peut tre approche par un dveloppement en srie de Taylor.

Figure 65 : dcoupage diffrences finies

2
H 1
2 H
H a = H i + ( xa xi )
x
x
+

(
)
a
i

2
x i 2
x i

Hydrologie - Hydrogologie Septembre 2008 Vronique MERRIEN-SOUKATCHOFF

86

B.10.3.1 En rgime permanent :

x 2

hi +1, j hi , j

hi , j hi 1, j

hi 1, j 2hi , j + hi +1, j

(x )2

2 h hi , j 1 2hi , j + hi , j +1

y 2
(y )2
2
dans le cas ou x = y H = 0 => hi 1, j + hi +1, j + hi , j 1 + hi , j +1 4hi , j = 0

Il existe plusieurs techniques numriques pour rsoudre le systme d'quation, en gnral des
techniques itratives (explicites) sont utilises.
On peut cite parmi les mthodes de rsolution :

Jacobi (la plus simple et la moins bonne)

Gauss Seidel (encore assez simple)

surelaxation

gradients

IADI

B.10.3.2 En rgime transitoire :

2H =

S H Q
+
T t T

Le temps est divis en pas de temps et entre le pas de temps n et le pas de temps n+1 on peut
n
n 1
h hi , j hi , j
hn + 1 hn

h
i, j
i, j
t

ou t
crire :
t
t
L'quation de diffusivit peut tre crite sous la forme :
n +1
n
n
S hi , j hi , j Qi , j
+
=

dans cette expression hin, +j 1 est

2
2
t
T
T
(x )
(y )
valu par les termes du pas de temps prcdent (mthode progressive).

hin+1, j 2hin, j + hin1, j

hin, j +1 2hin, j + hin, j 1

D'autres mthodes de rsolutions peuvent tre utilises

B.10.4 Principe des lments finis :

La mthode des lments finis consiste dcouper le domaine tudi en lments de forme
simple : triangles (le plus souvent), paralllogramme en 2D, ttradres, paralllpipdes en
3D. Les lments sont interconnects en des points appels nuds.
Le potentiel h dans chaque lment est dfinit par des fonctions d'interpolation partir des
valeurs aux nuds.
Principe : minimisation de l'nergie dans le domaine de rsolution de l'quation

Hydrologie - Hydrogologie Septembre 2008 Vronique MERRIEN-SOUKATCHOFF

Dans un lment : he ( x, y ) = a 0 + a1 x + a 2 y
Intrt de la mthode: maillage beaucoup plus souple qu'en diffrences finies
Inconvnients : plus lourd mettre en uvre d'un point de vue mathmatique

87

Hydrologie - Hydrogologie Septembre 2008 Vronique MERRIEN-SOUKATCHOFF

88

B.11 Ecoulement dans les milieux fracturs

Jusqu' prsent nous avons considr des coulements en milieux poreux "continus". Dans les
roches l'essentiel de l'coulement va s'effectuer dans les fractures. Les milieux fracturs
peuvent tre tudis par deux approches :
-

en prenant en compte les fractures une une;

en considrant un milieu continu quivalent.

B.11.1 Ecoulement dans une fracture

Dans un premier temps nous nous intresserons aux coulements dans une fracture en
supposant que la matrice c'est dire les blocs compris entre les fractures sont impermables.
Les vitesses d'coulement dans une fracture dpendent de son ouverture et de la rugosit de
ses pontes. La rugosit absolue correspond la hauteur moyenne des asprits dans la
fracture. La rugosit relative, sans dimension, est le rapport entre la hauteur des asprits et le
diamtre hydraulique moyen :
Rr =

DH

DH Diamtre hydraulique : DH = 4 RH = 4

RH rayon hydraulique

S : section mouille

P : primtre mouill

S
P

Remarques:
-

Pour une fracture infinie et plane d'ouverture e, on a D h 2e . En effet si L est la largeur

e
S
e L
e
de
la
fracture
si
<< 1
=2
D h = 4R h = 4 = 4
2e
L
P
2 (e + L )
e

+ 1
L
donc Rr =

DH

2e

B.11.1.1 Diffrents rgimes d'coulement

Les travaux de Louis ([27]) ont montr que l'coulement dans une fracture lmentaire
pouvait s'crire :

Hydrologie - Hydrogologie Septembre 2008 Vronique MERRIEN-SOUKATCHOFF

en rgime laminaire V = k f J f

en rgime turbulent V = k 'f J f

89

V : vitesse moyenne d'coulement dans la

fracture ([L/T] ou m/s) ;
kf : conductivit hydraulique de la fracture
([L/T] ou m/s) ;
kf' : conductivit turbulente de la fracture
([L/T] ou m/s) ;
Jf : projection orthogonale du gradient
hydraulique sur le plan de la fracture (sans Figure 66 : Ecoulement dans une fracture
dimension) ;
: degr de non-linarit sur le plan d'coulement (0,5 1. vaut 0,5 en rgime
turbulent compltement rugueux et 1 en rgime laminaire).
Louis a dfini, en fonction de ses exprimentations, 5 Rgimes d'coulement en fonction du
nombre de Reynolds (Re) et de la rugosit relative (Rr) de la fracture.
Re =

VD

: viscosit dynamique du fluide (cf. B.3.2 page 50) = 10 3 Pa s

: poids volumique de l'eau

Pour les 5 rgimes d'coulements Louis donne une expression de la vitesse V fonction de :
-

e : ouverture de la fracture

la viscosit

Remarque :
Le coefficient de perte de charge est dfini par la relation : J f =

V : vitesse moyenne d'coulement

1 V2

DH 2 g

Hydrologie - Hydrogologie Septembre 2008 Vronique MERRIEN-SOUKATCHOFF

90

Ecoulement
non parallle

Type 5

0,1

Type 4

0,033
Type 3

0,01

Ecoulement
parallle

Type 1
Type 2

0,001
1,00E+02

1,00E+03 2300

1,00E+04

1,00E+05

Figure 67 : diffrents rgimes d'coulements dfinis par Louis [27]

1. Type 1 : Ecoulement laminaire lisse (Poiseuille)

g e2
J f
V =
12

96
Re

2. Type 2 : Ecoulement turbulent lisse (Blasius)

4

1
7

5 4

g
e
2
g

Jf
V=

0,079

= 0,316 Re

3. Type 3 : Ecoulement turbulent rugueux (Nikuradse)

3,7
J f
V = 4 e g ln

R
r

DH
1

= 2 log
3,7

4. Type 4 : Ecoulement laminaire rugueux (Louis)

g e2
V =
1, 5
12 1 + 8,8 Rr

J f

k
96
1 + 8,8
=
Re
DH

5. Type 5 : Ecoulement turbulent trs rugueux (Louis)

1, 5

Re

Hydrologie - Hydrogologie Septembre 2008 Vronique MERRIEN-SOUKATCHOFF

91

DH
1

= 2 log
1,9

1,9
V = 4 e g ln J f
Rr

Dans le cas des coulements laminaires (type 1 et 4) la vitesse varie en fonction de l'ouverture
au carr, donc le dbit en fonction de l'ouverture au cube : on parle couramment de "loi
cubique"
Si la fracture n'est pas entirement ouverte (les 2 lvres se touchent par endroit) il faut
multiplier les expressions par le degr de sparation de la fracture F :
F=

Surface ouverte de la fracture

Surface totale de la fracture

L'expression de l'coulement laminaire de type 1 peut se dmonter partir des quations de

continuit et de Navier-Stokes
L'quation de continuit s'crit :
ur ( )
div ( V ) +
+ q = 0 (cf. page 74)
t
L'quation de Navier Stokes dans le cas d'un fluide incompressible :
ur
ur
ur uuuuur
ur
(V )
(V )
2
= f grad p + V en rgime permanent

= 0 donc l'quation se simplifie en

t
t
ur uuuuur
ur
f grad p + 2V = 0
z
Si l'coulement est unidirectionnel selon x
:
ur
V = (V ( x, z ),0,0)
y
Q
e
Si le fluide est incompressible et ici la
porosit dans la fracture est 1, l'quation
de continuit se ramne
ur
Vx
= 0 d'o V = (V ( z ),0,0)
x

x
Figure 68 : Ecoulement unidirectionnel
horizontal entre 2 plans parallles lisses

Hydrologie - Hydrogologie Septembre 2008 Vronique MERRIEN-SOUKATCHOFF

93

ur uuuuur
ur
La seule force agissant est la gravit oriente selon z donc f grad p + 2V = 0 est
quivalent
p
2Vx
+ 2 =0
x
z
p

=0
y
p
g +
=0
z
2Vx 1 dp
p dp
soit en intgrant :
=
=
la premire quation devient
dx
z 2
x dx
1 dp 2
z + C1 z + C2 Aux parois (pour z=0 et z=e) la vitesse est nulle donc
Vx =
2 dx

Comme

Vx =

1 dp
( z 2 + ez )
2 dx
e

1 dp
e3 dp
la vitesse
( z 2 + ez ) dz = b

2
12
dx
dx
0

Le flux travers une section b est q = b Vx dz = b

0
2

moyenne : V = b

e dp
12 dx

La relation q = b

e3 dp
est souvent appele "loi" cubique.
12 dx

Si l'coulement n'est pas horizontal, il faut prendre en compte la gravit dans la direction de
l'coulement et la relation s'crit en fonction de la charge hydraulique:
q = b

g e3 dh
g e3 dh
= b
12 dx
12 dx

[ML-1T-1] viscosit (dynamique) pour l'eau 20C 10-3 N.s/m2

[L2T-1] : viscosit cinmatique =

= 10 -6 m 2 / s ou myriastokes

B.11.1.2 Conductivit hydraulique, permabilit, transmissivit d'une fracture

La conductivit hydraulique d'une fracture est le rapport entre la vitesse et le gradient

hydraulique. Pour l'coulement de type 1 :
kf =

g e2
12

Bear et al. dfinissent la permabilit k par rapport la conductivit K par la

g
k . Cette permabilit k est une permabilit intrinsque. Ces auteurs
relation : K =

dfinissent donc la permabilit d'une fracture par k f =

e2
. Nous utiliserons le terme de
12

Hydrologie - Hydrogologie Septembre 2008 Vronique MERRIEN-SOUKATCHOFF

94

permabilit intrinsque pour cette expression

La transmissivit d'une fracture est le produit de la conductivit par l'ouverture de cette

dernire :
Tf =

g e3
12

Il existe une certaine confusion dans l'utilisation de ces termes et la transmissivit est souvent
appel par abus de langage conductivit hydraulique de la fracture.
B.11.1.3 Coefficient d'emmagasinement

Par analogie avec le milieu poreux il est galement possible de dfinir le coefficient
d'emmagasinement d'une fracture par la relation :

S = g 1 +eCw
Kn

: masse volumique (1000 kg/m3)

g :gravit

Kn : raideur normale du joint normal (en Pa/m)

e : ouverture

Cw compressibilit du fluide (4,4 10-10 m2/N ou Pa-1 pour l'eau)

Ce coefficient d'emmagasinement ne traduit que les variations dans la fracture.

Dans beaucoup de cas de massifs rocheux fracturs, l'coulement s'effectue principalement
dans les fractures, mais la matrice rocheuse joue un rle d'emmagasinement qu'il est possible
de prendre en compte, l'aide d'un coefficient d'emmagasinement suprieur celui qu'aurait
une fracture seule.
B.11.2 Ecoulement dans un rseau de fractures

Les coulements dans un rseau de fracture vont dpendre de la rugosit et de l'ouverture de

chacune des fractures mais galement de la rpartition gomtrique des fractures et des
intersections de fractures. Souvent on caractrisera la fracturation par un nombre fini de
famille de fractures partir de reprsentations strographiques et de considrations
gologiques plus qualitatives. Pour chaque famille on s'intressera la distribution des
orientations et des espacements de fracture. Cette distribution aura une influence importante
sur les coulements (par exemple s'il existe une seule famille de fracture d'orientation donne,
les coulements dans la direction orthogonale seront nuls)
Les modles conceptuels les plus simples considrent que l'coulement n'a lieu que dans les
fractures, mais l'influence de l'eau contenue dans les blocs dlimits par les fractures peut tre
non ngligeable surtout si l'on s'intresse aux phnomnes transitoires.
Dans ce qui suit nous prsentons quelques analyses effectues sur les rseaux de fractures,
mais la prsentation n'est pas exhaustive.
B.11.2.1 Analyse de la dimension de l'coulement

Hydrologie - Hydrogologie Septembre 2008 Vronique MERRIEN-SOUKATCHOFF

95

B.11.2.2 Modle d'coulement radial gnralis de Barker

B.11.3 Milieu continu quivalent

Certaines approches visent dfinir un milieu continu quivalent qui pourra tre trait par les
mthodes utilises dans les milieux continus. Pour pouvoir dfinir un milieu continu
quivalent, il faut pouvoir trouver un volume de terrain appel VER (volume Elmentaire
Reprsentatif) qui est dfini comme le volume minimal au-del duquel la conductivit
hydraulique moyenne reste constante quel que soit le volume. Il semble que pour les milieux
fracturs, il ne soit pas toujours possible de dfinir un tel VER

B.12 Transport de soluts en milieu poreux

B.12.1 Position du problme et mcanisme de transport

Le transport en milieu poreux est un phnomne difficile apprhender car il fait appel aux
domaines de la physique, de la chimie, de la biologie et de l'hydrodynamique. Il faudra en
effet prendre en compte :
-

la nature du produit transport (son caractre miscible ou non, sa solubilit) ;

la nature du sol (texture, structure, richesse en matire organique) ;

l'hydrodynamique du site (position et fluctuation de la nappe, zone sature, non sature) ;

l'activit biologique ;

l'activit anthropique (pompage, drainage, barrage,)

96

Hydrologie - Hydrogologie Septembre 2008 Vronique MERRIEN-SOUKATCHOFF

FLUX MASSIQUE
ENTRANT

Interaction
Convection-Dispersion
solut-matrice solide
PARAMETRES PHYSICO-CHIMIQUES
(pH, temprature, potentiel d'oxydo-rduction
teneur en eau, porosit, matire organique)

Solut

MO

Interaction
solut-autre solut

FLUX MASSIQUE
SORTANT

Bactries

Interactions
solut-biomasse

Figure 69 : Principaux phnomnes intervenant dans le transport de solut en milieu

poreux (d'aprs Benslimane [11])

Dans ce qui suit nous ne nous intresserons qu'au transport de produits "en solution" dans
l'eau par opposition aux coulements de 2 fluides non miscibles comme l'eau et l'huile ou
l'eau et l'air. Nous ne considrerons donc qu'une seule phase fluide et on dfinira la
concentration C d'une substance dans l'autre (en gnral concentration d'une substance dans
l'eau).
Les lments transports par l'eau peuvent tre :
-

minraux (radioactifs ou non) ;

organiques

Ces lments en solutions peuvent tre sous forme de sels solubles (ventuellement
partiellement ioniss), d'agrgats de molcules et/ou d'ions, de collodes.
La variation (dans l'espace et dans le temps) de la concentration d'un lment en solution sera
due aux principaux mcanismes suivants :
-

la convection (ou advection) ;

la diffusion ;

la dispersion cinmatique ;

les phnomnes d'change avec le milieu solide (phnomnes d'adsorption : rtention par
les particules solides et de dsorption) ;

Hydrologie - Hydrogologie Septembre 2008 Vronique MERRIEN-SOUKATCHOFF

97

les processus de dgradation physique, chimique ou biologique parmi les quels on peut
citer la dcroissance radioactive (que nous n'aborderons pas ici).

Les lments en solution peuvent tre conservatif (ou inerte ou non ractif ou traceur parfait)
s'ils ne subissent aucune modification et n'interagissent pas avec la phase solide du milieu lors
de leur transport. Dans le cas contraire les lments sont dit ractifs. Pratiquement aucun
lment n'est inerte, mais il existe des lments de ractivit ngligeable. Pour les lments
inertes la variation de concentration dans le temps et dans l'espace n'est due qu' l'advection,
la diffusion, la dispersion cinmatique et ventuellement des processus physiques tels que
filtration.
B.12.2 La convection (ou advection)

L'advection (parfois galement appele convection) correspond l'entranement du solut la

mme vitesse que le fluide qui se dplace. Pour l'valuer il est ncessaire de connatre la
vitesse relle du fluide, qui est diffrente de la vitesse de Darcy, car elle ne concerne que l'eau
"libre" (i.e. susceptible de s'couler sous l'action des gradients hydrauliques cf. B.2.7.1 page
45). La vitesse relle dpend donc de la porosit cinmatique c.
Si on suppose que l'advection est le seul phnomne de transport, le principe de conservation
de la masse permet d'effectuer un bilan des entres sorties et d'crire.
ur
C
c
= CU avec U vitesse de Darcy
t

( )

En effet (cf. Figure 70), considrons un volume lmentaire de milieu poreux (D) de frontire
extrieure ().La conservation de masse dans le volume D permet d'crire que lintgrale sur
() du flux massique de l'lment transport travers () est gal la variation de la masse
dans l'lment de volume (D) dans l'intervalle de temps unit.
r
"Le flux volumique de fluide traversant la surface () est l'intgrale de la vitesse de Darcy U
sur () ; ce flux volumique se transforme en flux massique de l'lment transport en
r
multipliant scalairement la vitesse de Darcy U par la concentration C" (De Marsily [8]). Le
premier membre de l'quation de conservation s'crit :
r r
CU .n d

ur r
r
n tant la normale sur () oriente vers l'extrieur de D, le flux entrant est U n

En transformant l'intgrale de surface en intgrale de volume par la formule d'Ostrogradsky :

r r
r
CU .n d = div CU .dv

( )

98

Hydrologie - Hydrogologie Septembre 2008 Vronique MERRIEN-SOUKATCHOFF

()

r
n

r
U

(D)

y
x

Figure 70 : Volume lmentaire reprsentatif (D) de milieu poreux (d'aprs Benslimane

[11])

"La masse de l'lment transport contenue dans D s'obtient en sommant les volumes
lmentaires de fluide c*dv contenus dans le milieu poreux, multiplis par la concentration
volumique C du fluide en lment considr" (De Marsily [8]) :

c C dv

"La variation de masse dans l'intervalle de temps unit s'obtient en drivant cette expression
par rapport au temps" (De Marsily [8]) :

( C dv ) =
D

C
dv
t

Soit en galant les membres :

r
C
div CU .dv = c
dv
D
D
t

( )

Comme D est arbitraire il est possible de supprimer les signes intgrale dans les deux
membres :
r
C
div CU .dv=c
dv
t
cqfd.

( )

B.12.3 La diffusion
B.12.3.1 La diffusion molculaire

La diffusion est lie l'agitation molculaire des particules

Fick a tabli que le flux massique de particules dans un fluide au repos est proportionnel au
gradient de concentration :
r
r
= D0 C (1re loi de Fick)
Le coefficient de proportionnalit D0, est le coefficient de diffusion molculaire. Il est
isotrope et peut s'exprimer pour les liquides par la relation de Stokes-Einstein :

Hydrologie - Hydrogologie Septembre 2008 Vronique MERRIEN-SOUKATCHOFF

D0 =

99

RT 1
N 6 r

D0 : coefficient de diffusion molculaire [L2T-1]

R : constante des gaz parfaits = 8,32 units SI [ML-1T-2K-1]
N : nombre d'Avogadro = 6,01.1023
T : temprature absolue (K)
: viscosit dynamique du fluide [ML-1T-1]
r : rayon moyen des particules qui diffusent [L]
Si le transport d'lments dans un fluide au repos est seulement du la diffusion de Fick (pas
d'advection), le principe de conservation de la matire et la loi de mouvement, permet
d'tablir comme prcdemment la relation :
r
C
( D0 C ) =
(2me loi de Fick)
t
Les valeurs du coefficient de diffusion des liquides sont comprises entre 1 10-9 2 10-9 m/s
25C (Robinson et Stockes, 1965). Elles ne varient pas beaucoup avec la concentration, mais
sont fortement dpendantes de la temprature, diminuant d'environ 50% 5C par rapport
25C.
B.12.3.2 La diffusion en milieu poreux

En milieu poreux, la diffusion s'opre dans la phase fluide (celle qui s'coule et celle sui est
immobile). Les grains solides arrtent ou ralentissent fortement le mouvement brownien des
particules. Pour un fluide statique en milieu poreux, le coefficient de diffusion effectif De
(Bear, 1972) est plus faible que D0.
De = D0 =

D0

: tortuosit du milieu
En pratique, le coefficient varie de 0,1 (argiles) 0,7 (sables).
B.12.4 La dispersion cinmatique

La dispersion cinmatique est lie l'htrognit du champ de vitesse microscopique (cf.

Figure 71). La dispersion dpend d'un trs grand nombre de paramtres classs en trois
rubriques :
-

paramtres caractrisant le milieu poreux :

courbe granulomtrique et en particulier le diamtre moyen des grains (dp) ;
la surface spcifique des grains ;
la densit apparente.

Ces paramtres contrlent :

la dimension des pores, donc le volume offert l'coulement et sa tortuosit ;
la quantit de solution lie la matrice solide.

paramtres caractrisant les solutions tudies :

Hydrologie - Hydrogologie Septembre 2008 Vronique MERRIEN-SOUKATCHOFF

100

variations de la masse volumique en fonction de la concentration en soluts ;

variations de la viscosit en fonction de la concentration et de la temprature ;
coefficient de diffusion molculaire dfini en solution libre.

Figure 71 : Facteurs lorigine de la dispersion longitudinale lchelle du pore, d'aprs

Fetter (1988)

Bear (1972) a tabli de manire thorique que la dispersion en milieu poreux tait analogue
la loi de Fick en considrant une distribution alatoire dans l'espace de canalicules
reprsentant les conduits travers les pores du milieu poreux :
r
r
= D C
Cette thorie est confirme par des expriences de laboratoire menes en laboratoire
l'Institut Franais du Ptrole (Pfannkuch, 1963).
D est un tenseur (suppos tre symtrique et du 2me ordre) dont les directions principales
sont la direction du vecteur vitesse et deux autres directions gnralement quelconques,
orthogonales la premire. Les coefficients du tenseur sont des fonctions de la vitesse
d'coulement.
Dans les directions principales (celle du vecteur vitesse et 2 autres orthogonales, le tenseur
s'crit :
DL

D= 0
0

DT
0

DT

DL : coefficient de dispersion longitudinale (dans le sens de l'coulement) [LT-1],

DT : coefficient de dispersion transversal (dans les deux autres directions) [LT-1].

Hydrologie - Hydrogologie Septembre 2008 Vronique MERRIEN-SOUKATCHOFF

101

Avec
DL = D0 + L u
DT = D0 + T u
L, T : coefficients de dispersion intrinsque ou dispersivits [L].
L : est de l'ordre de quelques centimtres en laboratoire, de 1 100 mtres sur le terrain.
T : est gnralement beaucoup plus petit.
Des traages in situ ont galement montr que le ratio L/T pour des milieux poreux varie
entre 0,05 et 0,3.
Afin destimer limportance relative de ladvection comparativement
dispersion/diffusion, il est courant d'utiliser le nombre adimensionnel de Peclet :
Pe =

la

U dp
D0

D0 : coefficient de diffusion molculaire du solut dans le solvant en m2/s (pour les liquides
10-9 2 10-9 m2/s)
U : vitesse du fluide
Quand Pe est faible (<<1) l'coulement est domin par la diffusion molculaire, quand Pe est
grand (>>1), c'est la dispersion cinmatique qui a l'impact le plus important, entre les deux se
situe une zone de transition.
100

1/Ped

10

1
1
0
10-2

2
10-1

3
10

4
102

103

5
104

105

106

Pe m

Figure 72 : Les diffrents rgimes de dispersion en milieu poreux satur (daprs

Pfannkuch, 1963)

Plus prcisment Pfannkuch (cf. Figure 72) dfinit 5 rgimes en utilisant un systme de
coordonnes adimensionnelles constitues par l'inverse du nombre de Pclet dynamique
U dp
U dp
en fonction du nombre de Pclet molculaire Pem =
Ped =
DL
D0
I

Diffusion molculaire pure

II

Superposition des effets de diffusion

III

Dispersion cinmatique prdominante

IV

Dispersion cinmatique pure

102
V

Hydrologie - Hydrogologie Septembre 2008 Vronique MERRIEN-SOUKATCHOFF

Dispersion cinmatique hors du domaine de validit de la loi de Darcy

B.12.5 Les phnomnes d'change avec le milieu solide

B.12.6 Les processus de dgradation

B.12.7 Expression de l'ensemble des termes du transport

En rsum L'quation ci-dessous donne l'expression du transport de solut en milieux poreux
(sous forme eulrienne)
q

1
C 1
C u i C
Dij
=

s (C C s ) liquide C + solide a S
R
t R xi x j R xi R c
c

Hydrologie - Hydrogologie Septembre 2008 Vronique MERRIEN-SOUKATCHOFF

103

Proprits physiques de l'eau : symboles, valeurs et units

Masse volumique

[M L-3]

Poids volumique

[M L-2 T-2]

Porosit

Sans unit

viscosit dynamique

[ML-1T-1]

Tenseur de
dispersion

D0

Coefficient de
diffusion
molculaire du
solut dans le
solvant

DL

coefficient de
dispersion
longitudinale

[LT-1],

DT

coefficient de
dispersion
transversal

[LT-1].

Indice des vides

Sans unit

Module d'Young

Pa [MLT-2]

gravit

[LT-2]

Coefficient de
permabilit de
Darcy

[LT-1]

Permabilit
intrinsque

Porosit

Sans unit

Pression

Pa [MLT-2]

VV
V

1 cPo :
centipoise pour
l'eau 20C, 1 cPo
=10-2 Po : poises;
1 cPo =10-3 Pa.s :
Pascal.seconde =
10-3 N.s/m2

e=

VV
VS

Page 41
Page 50+1

Page 41

Page 50

n=

VV
V

Page 41

Hydrologie - Hydrogologie Septembre 2008 Vronique MERRIEN-SOUKATCHOFF

104

Re

Nombre de
Reynolds

Sr

Saturation

Vitesse de Darcy

Volume Total

Re =

UD

Sr =

Vw
VV

Page 52

Page 41

m/s [LT-1]
[L3]

Page 41

Vitesse de Darcy

m/s [LT-1]

Page 50

VR

Vitesse relle,
intergranulaire

m/s [LT-1]

Vv

Volume des "vides"

[L3]

Page 41

Vv

Volume du solide

[L3]

Page 41

Vw

Volume d'eau

[L3]

Page 41

V ou VD

[L2T-1]

viscosit
cinmatique

= 10 -6 m 2 / s

Page 50+1

Masse volumique (masse par unit de volume)

M
[ML-3]
V

La masse volumique de l'eau varie avec la temprature et la pression; Schller [30] a propos
les rgressions empiriques suivantes :

= 999,9731

= 999,9731

(T 3,982)2
503570

T + 273
pour 0C T 42C
T + 62,26

(T 3,982)2 T + 273 350 T

466700

T + 67

pour 17C T 102C

365 T

Hydrologie - Hydrogologie Septembre 2008 Vronique MERRIEN-SOUKATCHOFF

105

Bibliographie Eaux Souterraines

[7] ANDERSON Mary P.,& WOESSNER William W., 1991, Applied groundwater
modelling (simulation of flow and advective transport), Academic Press
[8] BANTON O., BANGOY L., 1997, Hydrogologie, Multiscience environmentale des eaux
souterraines, Presses de l'Universit du Qubec, AUPELF
[9] BEAR, 1972, Dynamics of porous media, Elsevier
[10] BEAR J., TSANG C.-F., DE MARSILY G., 1993, Flow and Contaminant transport in
fractured Rock, Academic Press, Inc.
[11] BEN SLIMANE Frid, 1997, Mthodologie d'tude et de modlisation du transport de
polluants organiques en milieux poreux naturels : application l'tude de deux cas de
pollutions issues de l'industrie du bois, Thse INPL, LAEGO-EMN, Ecole Doctorale
PROMEN, soutenue le 12/09/97.
[12]

CASSAN, Aide mmoire d'hydraulique souterraine, Presses de l'ENPC

[13] CASTANY G., 1968, Prospection et exploitation des eaux souterraines, Dunod
Universit, Paris
[14]

CASTANY G., 1982, Principes et mthodes de l'hydrogologie, Dunod Universit

[15] CASTANY G., 1963, Trait pratique des eaux souterraines, reconnaissance et
exploitation, Dunod, Paris
[16] CASTANY G., MARGAT J., 1977, Dictionnaire franais d'hydrogologie, Editions
du BRGM, 249 pages.
[17] CHEREAU A., DEGALLIER R., DURAND A., GAILLARD P., DE MARSILY G.
Dtermination des paramtres hydrodynamiques des nappes souterraines partir de leurs
variations de niveau sous l'influence des pompage. Document BRGM 150, 1988.
[18] COOPER H.H. Jr., BREDEHOEFT J.D., and PAPADOPULOS I.S., 1967,"Response
of a finite-diameter well to an instantaneous charge of water: Water Resources Research,
v. 3, no. 1, p. 263-269.
[19]

DE MARSILY G., 1981, Hydrogologie Quantitative, Masson

[20] DE MARSILY G, 1986, Quantitative Hydrogeology: Groundwater Hydrology for

Engineers, Academic Press.
[21]

DE WIERT R. Hydrogeology.

[22]

DULHIEN, 1979, Porous Media : fluid transport & pore structure.

[23] ENGEBRISTEN S. E., SANFORD W. E.,1998, Groundwater in geologic processes,

Cambridge University Press
[24]

FILLIAT, 1981, La pratique des sols et des fondations, Edition du Moniteur

[25]

FREEZE R. A., CHERRY J. A., 1979, Groundwater, Prentice-Hall, Inc.

[26] HILLEL Daniel, 1988, L'eau et le Sol. Principes et processus physiques, collection
Pedasup 5, 2me dition.
[27]

LOUIS C., 1974, Introduction l'hydraulique des roches, Bull. BRGM deuxime

Hydrologie - Hydrogologie Septembre 2008 Vronique MERRIEN-SOUKATCHOFF

106

srie. Sect. III, n 4, 1974. pp.283-356

[28] MAGNAN J.P., 200??, L'eau dans les sols, C212, Revue des Techniques de
l'ingnieur, trait de construction, volume C22, 22 pages
[29]

SCHNEEBELI G., 1966, Hydraulique souterraine, Eyrolles

[30]

SCHOELLER H., 1962, Les eaux souterraines, Masson

[31]

STRACK O. D. L., 1988, Groundwater Mechanics, Prentice Hall, Englewood Cliffs

[32]

TODD D.K., 1959, Groundwater hydrology, John Wiley & Sons

[33] WANG H. F., ANDERSON M. P., 1982, Introduction to groundwater modeling,

Finite difference and finite element methods, Academic Press, San Diego

Quelques revues

La Houille Blanche

Les cahiers de l'ORSTOM

Water Ressource Research

Journal of Hydrology

Hydrology

Hydrogologie

Hydrologie - Hydrogologie Septembre 2008 Vronique MERRIEN-SOUKATCHOFF

107

Quelques sites web HYDRO & Environnement

http://waternet.com/

Welcome to WaterNet

http://water.usgs.gov/

USGS Water Resources Information

news:sci.geo.hydrology
http://www.usc.edu/dept/fccchr/ccvlib/

Newsgroup: sci.geo.hydrology
The World Wide Web Virtual Library: CrossConnection Control/Backflow Prevention

http://dir.yahoo.com/Science/earth_sciences/h
ydrology/

Yahoo! - Science:Earth Sciences:Hydrology

http://members.unine.ch/philippe.renard/hydr
ogen/hydrogen.html

Cours d'hydrogologie

http://www.ecologie.gouv.fr/rubrique.php3?id Ministre de l'cologie et du dveloppement _rubrique=4

Domaine de leau
http://www.cig.ensmp.fr/~thiry/courshydro/hydro00.htm
http://aih.brgm.fr/

Expliquer l'hydrogologie par Mdard Thiry

Comit Franais d'Hydrogologie
(Association Internationale des
Hydrogologues)

http://www.epa.gov/

US Environmental Protection Agency

http://www.cig.ensmp.fr/~hubert/glu/aglo.htm

Glossaire International d'Hydrologie

et
http://www.cig.ensmp.fr/~hubert/glu/HINDF
R.HTM
http://pagespersoorange.fr/stephane.guibert/hydrogeo/plan.htm
http://www.u-picardie.fr/beauchamp/

Notion d'hydrogologie et d'exploitation des

nappes
Plein de choses intressantes