Hydraulique des milieux poreux

Ecoulement des eaux à travers un milieu poreux = déplacement de

particules d’eau, le long de trajectoires appelées lignes de flux ou
lignes de courant.

Ecoulement possible: différence de pression ou un

gradient hydraulique non nul entre deux points et si la
porosité de la roche est suffisamment bien connectée,

Q= f (plusieurs variables indépendantes x, y, z et le

temps) loi expérimentale de Darcy.
 h  h  h  h
 (K x ) + (K y ) + ( K z ) = S s
  x x y y z z  t
 Application de la loi de Darcy
Ecoulement en milieu poreux et en régime permanent : cas d’un Puits en nappe
libre

rabattement
cône d'affaissement

Rayon d'influence R

Q = (K i) A

A est l'aire d'alimentation à une position x, équivalente à l'aire d'un cylindre circulaire
de rayon x et de hauteur y.
H
H  y2 
A=2 π x y dy R Q Qln( x)r = 2 K  
 dx =  2Kydy
R
Q = 2 Kxy
dx r x  2 h
h
K ( H 2 − h 2 )
Q= formule de Dupuit
R
ln( )
r
En pratique on évalue le rayon d'influence par la formule empirique de Sichardt
R = 3000 ( H − h) K
 Ecoulement en milieu poreux en régime permanent cas d’un puits en
nappe captive

Q = (K i) A

A = 2xb
dy
Q = 2 Kxb
dx

Qln( x)r = 2 Kb  y h

H
R Q R H

r x
dx =  2Kbdy
h

2 K b ( H − h )
Q=
R
ln( )
r
 Application de la loi de Darcy

Écoulement en régime permanent d’une nappe captive

Débit unitaire
a h

 dL = − Kb  dh
'
Q q
q = '

l 0 h1

q = − Kb
dh q'
h = h1 −
'

dL a
Kb
Exercice : Un aquifère confiné présente les caractéristiques
suivantes : b = 33m (épaisseur), largeur de l’aquifère =
7km.
La conductivité hydraulique = 1,2m/j.
Deux forages d’observation sont creusés dans l’aquifère
suivant la direction d’écoulement de la nappe. Ces forages
sont distants de 1,2km.
Dans le premier forage la hauteur piézométrique h est de
97,5m.
Dans le deuxième forage h = 89m.
Calculez le débit total de la nappe
Calculer la hauteur piézométrique au point M situé à 0.3km
du forage f1.
 Écoulement en régime permanent d’une nappe libre

Pas de recharge ni de décharge

dh
q ' = − Kh
dL
Conditions aux limites L h2
dh q d = − K  hdh
'
- à P1 : x = 0
- à P2 : x = L
h = h1
h = h2
q ' = − Kh 0
x
h1
dx

K  h22 h12 
h2
 h2 
q x
' L
0 = −K   q =−
'

L 2
− 
2
 2  h1 Équation de Dupuit
 Écoulement en régime permanent d’une nappe libre
une nappe libre avec accrétion
aquifère homogène et isotrope h 2 = h12 −
(h
1
2
− h22 )
x+
W
(L − x )x
L K

le terme W représente la recharge par unité de surface

dh
q ' = − Kh
dx

dh K  dh 2
 ( h 2
− h 2
2) L
q ' = − Kh = −   = K 1 − W ( − x)
dx 2  dx  2L 2
Au niveau de l’axe du dôme q’x = 0

(h12 − h22 ) L
K − W ( − x) = 0 Avec x = d
2L 2

L K  h12 − h22 
d= −  
2 W  2L 

h 2
=h 2
−
(h 1
2
−h 2
2 ) d + W (L − d )d
max 1
L K
Un canal est creusé parallèlement par rapport à une rivière. La
distance qui sépare le canal et la rivière est de 1500 m. Le canal et la
rivière pénètre entièrement la nappe libre.
On donne les caractéristiques suivantes : K = 1.2 m/J
Infiltration 1.8m/an
Evaporation 1.3m/an
La charge au niveau de la rivière est de 31m et la charge au niveau du
canal est de 27m.
1- Calculez la position du dôme piézométrique existant entre la
rivière et le canal
2- Calculez l’altitude piézométrique maximale
m
Sol

surface eau

substratum
m
Sol

surface eau

substratum
avoir les données de terrain :
charge piézométrique et coordonnées
des points (x, y , z).

reporter les données sur une carte :

les noms des points et la charge
hydraulique

choisir l’équidistance de la carte

tracer les courbes selon la méthode

choisie
Méthode des triangles
Écoulement uniforme

• Q, V, sens et direction constants en tout point

• lignes équipotentielles parallèles

• profil piézométrique linéaire (aquifère à épaisseur constante et homogènes)

 Écoulement non uniforme Q, V, sens et direction variables

2 grands types d’aquifères : aquifère à nappe cylindrique et aquifère à nappe radiale

• aquifère à nappe cylindrique : lignes équipotentielles sont des droites parallèles à espacement variable

1-Resserrement des isopièzes dans le sens de l’écoulement un profil piézomètrique de type parabolique

2- Écartement dans le sens de l’écoulement le profil piézomètrique est de type hyperbolique

• aquifère à nappe radiale
formes coniques, bombées ou déprimées isopièzes présentent des courbures

L’orientation de la concavité par rapport au sens de l’écoulement différentie deux types d’aquifères

- aquifère à nappe convergente

Ouverture de la concavité vers l’aval

espacement croissant
nappe convergente à profil hyperbolique espacement est décroissant
nappe convergente à profil parabolique

40

35
30

25
Écoulement convergent et les courbes piézométriques fermées,
dépression piézomètrique
-aquifère à nappe divergente
concavité tournée vers l’amont

Espacement croissant  profil hyperbolique Espacement décroissant profil parabolique

25
30
35

40

Écoulement divergent et courbes piézométriques fermées dôme piézomètrique

Relation entre cours d’eau et nappe superficielle

drainage de la nappe par le cours d’eau

Écoulement vers le cours d’eau

lignes de courant convergent vers le cours d’eau.
Concavité vers l’aval de l’écoulement superficiel
 Relation entre cours d’eau et nappe superficielle

drainage du cours d’eau par la nappe

Le cours d’eau alimente la nappe

concavités tournées vers l’amont
lignes de courant sont divergentes par rapport à l’axe du cours d’eau
les 2 phénomènes associés suivant la période de l’année

les 2 phénomènes peuvent se passer sur un même système hydrogéologique

 Limites des nappes
lignes de courant parallèles à la limite géologique une limite imperméable

les lignes de courant sont perpendiculaires à la limite géologique, un flux entrant

 0m 20000m 40000m
1580000

Kaolack Birkilane
1560000
Djilor Saloum River

nd
1540000 sI la
Atla

um
la o
ntic

S Sokone
Oce
an

lon
1520000
Nioro du Rip

o
ob
Ba
Karang
300000 320000 340000 360000 380000 400000 420000 440000

Drainage network
Groundwater trench
Groundwater flow
Groundwater divide Water table contour