PR Erick Carlier Mcanique des sols Licence 3 Gnie Civil

Chapitre IV. Lois de Bernoulli et Darcy. Permamtres et permabilit

quivalente.
Ce chapitre comporte 6 exercices

I. Loi de BERNOULLI et charge (potentiel) hydraulique

1.1 Charge (potentiel) hydraulique

La charge hydraulique en un point est gale la somme de la pression en ce point (par exemple, lie la
hauteur deau au dessus de ce point) et de lnergie potentielle par unit de volume lie laltitude de ce point :

.g.h + .g.z

: masse volumique de leau et g : acclration due la pesanteur

Cest lnergie totale mcanique (par unit de volume) en ce point.
Leau se meut des points de haute nergie vers les points de basse nergie.
Lquation prcdente, pour tre complte doit inclure un troisime terme relatif lnergie cintique (par unit
de volume) qui dpend de la vitesse (1/2 ..V2). En divisant les termes par le produit constant .g, on exprime la
charge hydraulique en unit de longueur :

H = h + z + (V2/2g)

Mais, en rgle gnrale, les coulements souterrains sont trs lents et le terme relatif lnergie cintique peut
tre nglig. Considrons un tube rempli deau sous pression et inclin, avec un dbit deau Qx (m3/s) qui
circule ; la pression en un point du tube peut tre visualise en perant un trou et en y positionnant un petit tube
vertical ; la hauteur deau atteinte dans le petit tube verticale indique la pression en ce point exprime en unit de
longueur (m) :

La diffrence de charge est (Pe2 + Z2) (Pe1 + Z1) = H2 H1 = H)

On dfinit le gradient hydraulique comme la perte de charge (hydraulique) par unit de longueur :

Ix = H / x.
Remarquons, que sur une mme verticale, la charge hydraulique est constante en tout point.
Exemple : sur la verticale 1
PR Erick Carlier Mcanique des sols Licence 3 Gnie Civil

Un point la base du petit tube :altitude z=z1 et Pression = Pe1 (+ Patmos que lon prend gale 0 par
convention) soit H1=z1+Pe1
Un point au sommet du petit tube : P=pression atmosphrique soit 0 par convention ;altitude z=z1+Pe1
Soit H1=0+z1+Pe1

1.2 Application linterface eau douce-eau sale (zone littorale)

En zone littorale, il existe ce que lon appelle un biseau deau sale situ sous leau douce car leau de mer est
plus dense (lhuile surnage leau !). La frontire entre les deux types deau nest pas une simple ligne ; cest, en
fait, une zone de quelques mtres compose deau dite saumtre.

En considrant un tat dquilibre hydrostatique, les deux points a et b sont au mme potentiel. En prenant
comme origine des hauteurs la ligne ab, on peut crire :

0+ (Z+h).douce. g = 0+Z.sale. g
De la relation prcdente, il vient :

d
Z = h. 40 .h
s d

On remarque que pour une baisse dh de 1 m deau douce, la remonte dZ est de 40 m ! Une exploitation abusive
de forage en zone littorale conduit un envahissement par leau sale de ces derniers et une avance lintrieur
des terres du biseau.

II. Loi de DARCY

Nous allons la dfinir en partant de lanalogie entre phnomnes lectrique et hydraulique.

loi d'Ohm : si n est le nombre d'lectrons qui se dplacent par unit de volume et qui traversent une section s du
fil une vitesse moyenne u, la quantit d'lectricit est:
dQ = n .e. s.u.dt (e: charge de l'lectron),
c'est--dire le nombre de charges contenues dans le volume s.u.dt. L'intensit du courant est le nombre de
charges par unit de temps qui traversent la section s, soit

dQ/dt= i.
PR Erick Carlier Mcanique des sols Licence 3 Gnie Civil

i s'exprime par le rapport entre la diffrence de potentiel (ddp) entre a et b, lie aux nombres d'lectrons en a et
en b, et la rsistance du fil, elle-mme fonction de sa longueur, de sa section et de sa rsistivit dont l'inverse est
la conductivit ou conductance:
i =ddp/R = ddp.s/(.l)
s: section du fil lectrique
: rsistivit (1/: conductivit)
l: longueur du fil lectrique
R= .l/s (rsistance)

Electrique hydraulique
ddp= r.L.i/s soit
i =ddp. s/(L.r)

ddp h (m)

L, s L, s (m, m2)

i Q (m3/s)

1/ K (m/s)

Dans le dispositif hydraulique ci-dessus, qui est un permamtre charge constante, le mouvement de l'eau au
travers de l'chantillon de sol est cr par la diffrence des niveaux d'eau h, qui est l'analogue hydraulique de la
ddp, la section du fil a pour analogue hydraulique la section de l'chantillon, la longueur du fil a pour analogue la
longueur de l'chantillon. La conductivit a pour analogue une caractristique du sol qui exprime son aptitude
laisser passer l'eau: la permabilit. i a pour analogue un volume d'eau passant par unit de temps au travers de la
section de l'chantillon, c'est--dire un dbit Q.

Loi dOhm : = . =

Analogue hydraulique de la loi dOhm : = Q

PR Erick Carlier Mcanique des sols Licence 3 Gnie Civil

soit:

Loi de DARCY:
h/l
Q = - K . s .
Nous verrons la signification du signe un peu plus loin.
l: trajet des molcules d'eau dans le milieu poreux engendr par h, s surface traverse perpendiculairement par
les molcules d'eau
.
K a la dimension d'une vitesse (L.T-1).

Convention de signe :

Si le dbit est selon laxe des x croissant, alors Q>0, sinon Q<0. Selon le schma ci-dessus, laxe des x est
orient vers la droite, le terme h= h(x+x) - h(x) est ngatif et le dbit Q est positif. Il convient donc
dintroduire un signe la partie droite de lquation de Darcy. Cette formule avec le signe sera toujours
valable quelque soit lorientation des axes.

h/l est appel gradient hydraulique, il s'crit, en une dimension et en notation diffrentielle: dh/dx . En trois
dimension, il est reprsent par un vecteur:
grad (h) de coordonnes (dh/dx, dh/dy , dh/dz )

Le rapport Q/s = - K . dh/dx est appel vitesse de Darcy; c'est, en fait, un dbit par unit de surface, et elle
ne correspond pas la vitesse moyenne d'coulement dans le milieu (sauf dans le cas d'un tube vide ou d'un
cours d'eau).
La vitesse moyenne d'coulement dans un milieu s'exprime par:
u =vitesse de Darcy/n
o n est la porosit.

III. Permamtres charge constante.

Ce sont des appareils qui permettent de dterminer la permabilit d'un sol ou d'une roche.

volume d'eau
recueilli pendant
le temps de l'exprience
V/t = Q

section s

h. s)
K=Q.L/(-
Le permamtre de Darcy:
PR Erick Carlier Mcanique des sols Licence 3 Gnie Civil

section S

za=cte
l
ze
zs
0
En appliquant la loi de Bernouilli entre l'entre et la sortie de l'chantillon, on crit:

He = Ze + (Za - Ze) (Za-Ze reprsente la pression sur les points dentre de lchantillon exprime en unit de
longueur)

Hs = Zs +0 (le point de sortie est la pression atmosphrique donc gal 0 par convention)

donc He - Hs = Za - Zs

donc K = Q. l / (S.[Za-Zs])

IV. Permamtre charge variable.

Il est utilis pour la mesure de la permabilit de sol trs peu permable (argileux)
section s

h(t)

h(t+dt)

section S
Considrons un dbit positif, soit un axe des abscisses croissant vers la droite.
PR Erick Carlier Mcanique des sols Licence 3 Gnie Civil

Le gradient hydraulique est la diffrence de charge entre lentre et la sortie, soit h sur la figure ; il diminue avec
le temps (h(t+dt) < h(t) ce qui signifie que la variation du gradient hydraulique en fonction du temps
dh = h(t+dt)-h(t) est ngative. Le dbit Q peut sexprimer comme la variation instantane dun volume deau en
fonction du temps :

Q(t) = dV/dt= - s dh/dt (Q est positif et dh est ngatif)

Loi de Darcy: Q(t) = S . K . h(t)/l

=

En intgrant, on obtient:

= ( )

d'o:

s .l h
K= . ln 0
S .(t t 0 ) h

Exercice 1 :
En vue de la dtermination du coefficient de permabilit d'un chantillon, un essai au permamtre
charge constante a t effectu dans les conditions suivantes :
section de l'prouvette : 28,3 cm2
hauteur de l'prouvette : 10 cm
diffrence de charge hydrostatique : 50 cm
dure de l'essai : 10 s
quantit d'eau recueillie : 500 g
Calculer la permabilit (ms-1).

Exercice 2 :

Essai de permamtre charge variable :

section de l'prouvette : 23 cm2
hauteur de l'prouvette : 10 cm
section du tube de mise en charge : 1 cm2
charge hydrostatique initiale : ho = 60 cm
charge hydrostatique finale : h1 = 50 cm
dure de l'essai : 10 mn
Calculer la permabilit (ms-1).

Exercice 3:
Le filtre est rempli d'eau. Calculer le temps pour que le "caf soit pass" (filtre vide).
PR Erick Carlier Mcanique des sols Licence 3 Gnie Civil

filtre

10cm
90

2cm
caf (k = 0,001m /s )
4cm

V. Permabilit quivalente.

5.1 coulement perpendiculaire linterface.

On cherche dterminer la permabilit verticale quivalente au milieu stratifi.

Q
l
e1 I
I
R1
K1 e2

K2
R2 R
en
Kn

Rn

n
ddp= I . Ri= R.I
i =1

Ri a pour analogue hydraulique:

1 ei
.
K i l .L

R a pour analogue hydraulique:

1 ei
.
Kv l.L
Or :
PR Erick Carlier Mcanique des sols Licence 3 Gnie Civil

R = Ri

donc:
ei
1
.
e = K
i i
Kv l.L l .L

donc

Kv=
e i
e
K i

Kv est la permabilit verticale du milieu poreux quivalent au milieu stratifi.

5.2 coulement parallle linterface

On cherche dterminer la permabilit horizontale quivalente au milieu stratifi.

l
e1
I1
K1 e2
I I
K2
Q
en
In
Kn

ddp = R1 . I1 = R2 .I2 = .... = Rn. In = R . I ( R: rsistance quivalente)

I = I1 + I2 + .....+ In = R/R1. I + R/R2. I + ....+ R/Rn. I donc :

1 n 1
=
R i =1 Ri

1/Ri a pour analogue hydraulique

e .l
Ki . i
L
1/R a pour analogue hydraulique

l. ei
Kh.
L
Soit :

Kh=
K .e
i i

e i

Kh est la permabilit horizontale du milieu poreux quivalent au milieu stratifi.

PR Erick Carlier Mcanique des sols Licence 3 Gnie Civil

Exercice 4 :
h1-h2 = 1,30 m
Le dbit par unit de longueur tait, avant le colmatage, de 1,82 x 10-5 m3/s/m, la permabilit de la zone
colmate est de 1,4 x 10-5 cm/s.
Calculer le dbit par km de canal qui transite entre les deux canaux.

10m

h1
h2 4m colmatage
(paisseur = 4cm)

Exercice 6 :
Soit un versant de valle avec deux terrasses alluviales embotes : Donner l'expression du dbit unitaire
alimentant la rivire en fonction de L1, L2, K1, K2, h1 et h2. (nappes libres) .K1 > K2

h1 K1
hs K2
h2

L1
L2
PR Erick Carlier Mcanique des sols Licence 3 Gnie Civil

VI. Zone non sature

La zone non sature est un milieu gnralement triphasique : solide (roche, grains) + eau + air.

6.1 Remonte capillaire.

A cause des tensions de surface, leau peut slever au dessus de la surface de la nappe. Les remontes seront
dautant plus importantes que le diamtre moyen des vides sera petit. En effet, lquilibre entre force de gravit
et force de tension met en vidence ce fait. (cf ci-dessous). La pression de leau lintrieur dune remonte
capillaire est ngative (infrieure la pression atmosphrique prise comme rfrentiel gale 0).
PR Erick Carlier Mcanique des sols Licence 3 Gnie Civil

Le mnisque dans le petit tube est incurv vers le bas, cela signifie donc que la pression dun point situ juste au
dessous du mnisque est infrieure la pression atmosphrique ; elle est donc, par convention, ngative.
En prenant un point juste sous le mnisque du petit tube, donc une hauteur hc et un autre point situ sur la
surface de leau du rcipient et utilisant le thorme de Bernoulli qui exprime la conservation de lnergie
mcanique, on peut crire :
0+0 =hc+P
P tant la pression du point juste au dessous du mnisque. On en dduit donc que :
P=-hc
Si on exprime cette pression en pascal :
=

La figure ci-dessus montre que, dans un sol, les remontes capillaires, sur une hauteur totale h1 comprennent une
partie sature sur une hauteur h2.

Pression : u = w .h si partie sature de la remonte.

u = w .h.S si partie non sature avec S : indice de saturation (rapport entre volume de leau et
volume des vides)

La figure ci-dessous prsente les diffrents domaines du sol la nappe : Pe : pression de leau
PR Erick Carlier Mcanique des sols Licence 3 Gnie Civil

6.2 Exemples dexprimentation en zone non sature

6.2.1 Colonne horizontale

Tube de gauche : on utilise la loi de conservation de lnergie mcanique (Bernoulli) entre un point situ sur la
surface de leau du rcipient (altitude -20 cm) et un autre point situ lentre du sol (altitude 0) :

-20+0 =0+ P => P = Hp,1 = -20 cm

La charge hydraulique (potentiel) est donc :

H1 = 0 + (-20) = -20 cm
De mme, pour le tube de droite de sortie :

-30+0 =0+P =>P= Hp,2= -30 cm

H2 = 0 + (-30) = -30 cm

La pression moyenne H dans la colonne est -25 cm.

p

Si la permabilit K =1.3 cm/hr par exemple, alors le flux (dbit par unit de surface) est, par application de la
loi de Darcy:

6.2.2 Colonne verticale

Considrons une colonne verticale avec les mmes valeurs de pression aux limites :
Le niveau deau dans le rcipient de gauche est 20 cm sous lentre du sol donc une altitude de 30-20 =10 cm
La pression lentre du tube est :

+10+0 =30+P =>P=Hp2= - 20 cm

PR Erick Carlier Mcanique des sols Licence 3 Gnie Civil

La char ge hydraulique lentre est donc :

H2=30-20 =10 cm
Par le mme raisonnement, on obtient la charge la sortie :

0+P =-30+0 => P = Hp1 =-30 cm

H1 = -30 cm
On remarque que laxe des ordonnes est dirig vers le haut alors que le dbit dcoulement est vers le bas, selon
la convention, il sera donc ngatif.

Par rapport la colonne horizontale, on remarque que le flux dans la colonne vertical est plus important ; ceci
est du lexistence, dans ce cas, dun Potentiel gravitaire.

6.2.3 Colonne verticale bis

Considrons le cas particulier o laltitude du niveau de leau dans le rcipient connect au tube dentre est le
mme que celui de la sortie du sol, soit 0.

On remarque que dans ce cas, la pression lentre et la sortie du sol est la mme (-30 cm) :

Entre : 0+0 = Hp2+30 => Hp2 = -30 cm

Sortie : 0+Hp1 = -30 +0 => Hp1 = -30 cm

PR Erick Carlier Mcanique des sols Licence 3 Gnie Civil

Le gradient de pression est donc nul. La pression est uniforme dans le sol ; ceci implique que la teneur en eau est
constante. Ce dispositif particulier est donc trs utile pour dterminer la conductivit hydraulique en fonction de
la pression en zone non sature. Le mouvement de leau est donc seulement du au gradient gravitaire

6.3 deux tensiomtres proches

densit du mercure : 13,5

PR Erick Carlier Mcanique des sols Licence 3 Gnie Civil

Examinons le dispositif de gauche. Leau contenue dans le tube, traverse la paroi poreuse situe la base (point
2) une altitude de -100 cm, compte partir de la surface du sol. Leau est donc aspire par le sol (phnomne
de succion). Ce dpart deau provoque une monte du mercure dans le tube dont le niveau se stabilise 35 cm
au dessus du niveau du mercure contenu dans le rcipient. Le niveau du mercure dans le petit tube sest stabilis
45 cm au-dessus du sol.
Considrons un point linterface eau-mercure et un point la surface du mercure contenu dans le rcipient ; on
peut crire :

10 + 0 = 45 + %& '(

On en dduit que la pression linterface est :

= 35

Au point 2, le potentiel total (potentiel gravitaire+pression) scrit :

100 + *+,- = *-

Comme on est en quilibre statique, le potentiel total dun point de linterface est gal au potentiel total au point
2:
100 + *+,- = (10 + 35) 35

En divisant par la masse volumique de leau et par lacclarration due la pesanteur g, on exprime les potentiel
en unit de longueur (hauteur deau) et on en dduit la valeur de la pression en 2 :

*+,- = 145 35 = -327 cm

Avec d, la densit du mercure (13,5)

Le potentiel total au point 2 est donc :

*- = 100 327 = 427

En tenant le mme raisonnement, on trouve que le potentiel total pour le tensiomtre de droite au point 1 est :

* = 465

Comme H > H , c lcoulement se fait vers le bas.

2 1
La pression moyenne est :

Si K = 0.01 cm/hr, alors:

-331