PROPRIETES HYDRAULIQUES DES SOLS

1 - PERMEABILITE
1.1.- Définition de la charge hydraulique
Tous les sols sont plus ou moins perméables. L'eau filtre à travers l'argile comme à travers le
gravier mais ce phénomène sa manifeste avec des intensités différentes. On étudiera surtout
les écoulements qui correspondent à des régimes permanents. Dans de tels mouvements, les
particules fluides suivent des trajectoires invariables au cours du temps : les lignes de courant.
En supposant le liquide parfait (incompressible et non visqueux), on définit la charge hydraulique
2
V U
h à l'aide du théorème de Bernouilli : h= + + H ; Ce qui représente l'énergie d'une
2g γ w
2
V
particule fluide de poids unité. Mais dans le cas de filtration de l'eau par le sol la quantité
2g
est négligeable puisque les vitesses sont très faibles.

1.2.- Loi de Darcy

Les premières expériences sur la perméabilité ont été faite par Darcy en 1854 et il a remarqué
que le débit par unité de surface est proportionnel à la perte de charge et inversement
proportionnel à la hauteur : ce qui peut s'écrire :

P : pression de l'eau qui pénètre

Q : débit
Q P− P0 A : aire ou section de sol
=k.
A L
Po : pression de l'eau qui sort
l : longueur de filtration

Cette relation appelée loi de Darcy s'écrit :V = ki. Dans cette formule V est une vitesse fictive
Q
représentant la vitesse de décharge : , k représente le coefficient de perméabilité du sol. Les
A
valeurs sont très diverses suivant la nature du sol :
• gravier : de 10-1 à 102 cm/s.
• sable : de 10-3 à 10-1 cm/s.
• argile : de 10-11 à 10-7 cm/s.
Le coefficient i ( i= dh )est le gradient hydraulique, c'est-à-dire la perte de charge par unité de
dl
longueur :

1.3.- Mesure de la perméabilité en laboratoire :

Deux types d'appareils sont utilisés : le perméamètre à charge variable et le perméamètre à
charge constante.

1.3.1.- Perméamètre à charge variable

Le sol à étudier est placé dans un moule fermé à ses deux extrémités par deux couvercles. Le
couvercle supérieur permet l'évacuation de l'eau, le couvercle inférieur étant relié à un
réservoir et un tube fin.
 On commence par amener l'échantillon à saturation. Le moule
esst alors plein d'eau et relié au tube gradué. On constate
que le niveau baisse dans le tube et on mesure le temps T qui
s'écoule entre les niveaux h1 et h2, ces deux niveaux étant
repérés par rapport au niveau de l'ajutage de sortie. Si A
est la section de l'échantillon et a celle du tube gradué, le
débit q sera égal à : q = A.V

Or si à un instant t, l'eau dans le tube est au niveau h nous auront aussi : q = A.k.h (puisque le
L
gradient hydraulique est i= dh ). Le volume d'eau qui sort de l'échantillon est égal à la diminution
dl
de volume dans le tube : q.dt = -a.dh
A.k.h .dt =−a.dh
L
k.dt =− a .L. dh k.T =− a .L.log h2
A h A h1
Or la valeur du coefficient est donnée pour la valeur normalisée de la température de 20°C.
Celle-ci joue un grand rôle puisque la fiscocité de l'eau augmente dans des proportions
ηθ 1,79
importantes. On a la relation : k20 =kθ . avec ηθ =
η20 1+0,003368.θ +0,00022θ
2

⎧ θ en °C
⎨η en centipoise
⎩

Cet essai peut aussi être réalisé à l'aide de l'oedomètre.

1.3.2.- Perméabilité à charge constante

Cette deuxième méthode est utilisée pour des matériaux assez perméables (k > 10-3 cm/s). Dans
ce cas le tube gradué est déconnecté et on maintient un niveau constant dans le réservoir. Le
gradient hydraulique i est donc constant. La vitesse de décharge est déterminée en mesurant le
volume d'eau qui traverse un échantillon en un temps donné.

1.4.- Domaine de validité de la loi de Darcy

Les études de laboratoire ont permis de constater que la relation v = ki est bien vérifiée pour
tous les sols mais seulement à l'intérieur d'un certain intervalle de variation du gradient
hydraulique. Pour de très faibles vitesses de décharge, l'écoulement réel ne correspond pas à
l'écoulement théorique par suite des phénomènes colloïdaux. Pour les fortes vitesses, les
divergences sont dues aux forces d'inertie dans un mouvement non uniforme. Le nombre de
Reynolds Re correspond au rapport des forces d'inertie aux forces de viscosité : Re=V×d
η
Muskat a proposé de prendre comme limite de sécurité la vitesse critique correspondant au
nombre de Reynolds égal à 1, mais les expériences de Lewis et Barnes ont conduit à constater
que la loi était valable pour Re < 10. On peut donc dire que la loi de Darcy est une excellente
approximation pour les faibles nombres de Reynolds mais qu'elle devient de moins en moins
bonne lorsque V augmente.

1.5.- Perméabilité en terrain stratifié

Lorsqu'un terrain est constitué de plusieurs couches, il est évident que le coefficient de
perméabilité n'est pas le même pour chacune des couches et qu'il faudra distinguer les
perméabilités verticales et horizontale.

On cherche à déterminer le coefficient moyen de perméabilité horizontal kH, et le coefficient

vertical kV. Considérons d'abord un écoulement horizontal dont le gradient hydraulique est le
même à la traversée de chaque couche : V1 = k1.i, V2 = k2.i ; V3 = k3.i ;……….
Vn = kn.i

la vitesse moyenne sera alors : V = 1 .(V1H1+ V2 H 2 + V3H 3 +....+ Vn H n)=k H .i

H
V = i .(k1H1+ k 2 H 2 + k3H 3 +....+ k n H n)=k H .i
H
kH = 1 ∑ki.Hi
H i
Si l'écoulement est vertical, le principe de continuité exige que la vitesse de décharge soit la
même à la traversée de chaque couche : V = k1 i1 = k2 i2= k3 i3= k4 i4=…..= kn in= kV i
Or le gradient hydraulique i est égal à :
i = h = h1+ h2+ h3+...+ hn
H H
les hi représentant les pertes de charge.

=kv i =k1 1 =k 2 2 = h1 +h2 +...+hn

h h h
H H1 H 2 H1 + H 2 +...+ H n
k1 k2 kn
∑i hi ∑i hi
kV .i = d’où kV .i =
∑i Hkii ∑i Hkii
2 - HYDRAULIQUE SOUTERRAINE
2.1.- Définitions
Une nappe est une accumulation d'eau dans une roche fissurée ou un milieu poreux. Celle-ci est
retenue à sa partie inférieure par le "mur de la nappe", constitué d'un terrain imperméable.
 Si on enfonce des tubes capillaires
dans le sol, on pourra déterminer la
surface piézométrique, celle-ci étant
surmontée d'une zone appelée frange
capillaire due aux remontées
capillaires dans le milieu perméable.
Dans le cas de nappe captive, la
surface piézométrique se trouve au-
dessus du sol et on a des puits
artésiens.

2.2.- Généralisation de la loi de Darcy

En milieu homogène et isotrope, les vecteurs vitesses et gradient hydraulique sont colinéaires ;
or
V =k.i =−k. dh
dl
→
V =−k.grad H
→ → →
V =−grad (kH) = grad (−kH) = grad (φ)
Dans cette formule φ représente un potentiel des vitesses égal à -kH. Les équipotentielles sont
des surfaces d'égale charge hydraulique. Les composantes seront :
∂φ
Vx=
∂x
∂φ
Vy =
∂y
∂φ
Vz =
∂z
Or les conditions de continuité exigent que le liquide soit incompressible :

∂φ ∂φ ∂φ →
+ + =0 ⇒ divV ⇒ Δφ
∂x ∂y ∂z

2.3.- Ecoulement plans

Ce type d'écoulement est répandu et on choisira donc le plan vertical x o z comme plan
d'écoulement :
Vx = Vx(x,y)
Vy = 0
Vz = Vz(x,y)
La fonction potentiel φ est donc aussi une fonction de 2 variables ; on peut donc appliquer les
résultats de l'écoulement de l'électricité dans un conducteur qui sont analogues à ceux de
l'écoulement de l'eau dans le sol. On introduit alors la notion de potentiel complexe φ(x,z) +
iψ(x,z) de la variable complexe x + iz. Les courbes φ= constante sont des équipotentielles ; elles
admettent des trajectoires orthogonales ψ = constante qui sont les lignes de courant.

2.3.1.- Exemple : écoulement autour d'un rideau de palplanches

Le rideau est supposé de longueur infinie. Soit, dans une rivière, un rideau de palplanches fiché
dans une couche de sable de hauteur hs surmontant une argile. Les perméabilités étant très
différentes, on pourra supposer l'argile imperméable vis à vis du sable.
 Sur AB, la charge est constante et égale à :
U
hAB = AB + Z AB =hS + H1
γw

Sur CD, la charge est constante et égale à :

U
hCD = CD + ZCD =hS + H2
γw
La perte de charge globale est ∆ = H1 - H2

AB et CD sont des équipotentielles ; BEC et Ox sont des lignes imperméables donc des lignes de

courant pour lesquelles : ∂h =0 . (n : normale aux lignes). Ce type d'écoulement où la pression

∂n
n'est pas nulle à la surface du sol est un écoulement en charge.

2.3.2. - Exemple d'écoulement à travers une digue en terre :

Soit une digue en terre reposant sur un sol imperméable. La retenue d'eau est H. AB est une
équipotentielle hAB = H et BC est la surface de suintement. AD est une ligne de courant : ∂h =0 .
∂n
Il s'agit ici d'un écoulement à surface libre.

2.4.- Ecoulement à 3 dimensions

On étudiera plus particulièrement, les écoulements de révolution (problème des puits). Ce sont
des problèmes que l'on rencontre à l'occasion des pompages, les deux applications en étant
l'alimentation en eau et le rabattement de nappes. Examinons le cas d'un massif perméable
reposant sur un massif imperméable, dans lequel on fore un puit de rayon r. La surface libre de la
nappe, présente une dépression en forme d'entonnoir. Ce phénomène a été étudié par Dupuit.

Le débit sur une surface de révolution (x,z) est :

q = V.z.2πx = k. dz .2.zπ
dx
R H

∫q.dxx =2πk ∫ zdz

r h
2 2
πk.(H −h )
q=
log R
r
Le rayon d'action R, correspond à la distance à partir de laquelle le pompage cesse de faire
sentir son action. Il a été déterminé par Sichard par la formule : R = 3000 (H - h)k
2πke.(H −h)
Dans le cas des nappes captives, on arrive au résultat suivant : q =
log R
r

3 - CONTRAINTES EFFECTIVES
Les efforts se transmettent dans le sol à travers les grains et l'eau interstitielle. Mais comme il
n'est pas possible d'étudier la répartition réelle des contraintes en fonction de la position et de
la taille des grains, on considère le milieu homogène du sol puisque les différentes phases qui
forment un sol (gaz + liquide + solide) sont régies par des lois différentes. On a donc envisagé
une loi unique et on introduit la notion de contrainte effective.

3.1.- Postulat de Terzaghi

P = P '+ u.( A − Ac )
P P ' u.( A − Ac ) P ' A
Divisons par l’aire A : = + = + u.(1 − c )
A A A A A P = Force normale totale
P’ = Force normale intergranulaire
A
Ac est très petit relativement à A, d’où : c 0 u = Pression interstitielle
A A = Aire totale
P P' Ac = Aire de contact entre les grains
= σ = Contrainte totale et = σ ' = Contrainte effective
A A
d’où σ ' = σ + u ou σ = σ '− u
C’est probablement l’équation la plus importante de la Mécanique des Sols.
On considère le milieu comme composé de 2 phases seulement : milieu sec (gaz + sol) ou saturé
(liquide + sol). Le postulat de Terzaghi consiste à formuler l'existence d'un nouveau tenseur de
contraintes, le tenseur des contraintes effective qui gouverne à lui seul le comportement du sol
et du squelette granulaire. Si σ et τ sont les composantes normale et tangentielle de la
contrainte totale en un point sur une facette quelconque, σ’ et τ’ les composantes de la
contrainte effective et u la pression du fluide (air ou eau) on a :
σ’ = σ + u
τ’ = τ

Dans le cas général d'un complexe de 3 phases, le calcul pratique

soulève de nombreuses difficultés. Parmi les tentatives les plus
intéressantes, on citera la formule de Bishop : σ’ = σ - ua + χ(ua – uw)
ua - uw représentant les pressions de l'air et de l'eau. χ est un
coefficient qui dépend du degré de saturation Sr.

3.2.- Applications
Nous allons appliquer cette théorie à des cas simples de sol saturés.
3.2.1.- Sol saturé, nappe en équilibre :

L'eau est en équilibre, la répartition de la pression est

donc hydrostatique. Sur une facette horizontale située
à une profondeur z dans le sol, on a :
σ = γw.D + γ.z
u = γw.(D + z)
alors
σ‘= σ -u = (γ - γw).z = γ‘.z

Donc on voit qu'au point de vue des contraintes effectives, le sol se comporte comme s'il n'y
avait pas d'eau à condition de remplacer par ' puisqu'il est saturé.

3.2.2.- Ecoulement verticaux descendants

Dans l'exemple ci-dessus, si l'eau est en mouvement, la répartition de la pression n'est plus
hydrostatique. La définition même de perte de charge, montre que l'écoulement transmet au sol
une contrainte. Si l'écoulement est descendant dans le sens des z croissants : u = γw.(D+z-h)
du =γ .(1− dh)=γ .(1−i) car =i = dh
dz w dz w dz
Donc entre la surface z = 0 et la cote z on a :
u = u0 = γw.z.(1 - i) avec u0 = γw.D
On a toujours : σ‘= γw.D + γ.z ⇒ σ‘= (γ‘ – i.γw).z

3.2.3.- Ecoulement verticaux ascendants

Dans ce cas, c'est le phénomène inverse qui se produit, le courant agit en sens contraire du
poids. Le terrain est en partie porté par l'écoulement et sa force portante décroît. La pression
est plus forte en bas.
z = z0 = 0 u0 = γw.D du =γ .(1+ dh)=γ .(1+i)
dz w dz w
0
u0 – u = ∫γ w.(1+i).dz
z

u = γ w.(1+i).z +γ w.D ⇒ σ'=σ −u =(γ'−i.γ w).z

 γ'
Lorsque le gradient i atteint la valeur , les contraintes effectives s'annulent. Le squelette
γw
solide flotte au milieu de l'écoulement. Il peut en résulter des accidents graves pour les
constructions : digue, fond de fouille. Pour les sables (n = 40 % ; γs = 2,65 g/cm3) le gradient i
est voisin de 1. Le sable devient boulant (phénomène de Renard). Pour éviter ces phénomènes on
utilisera des filtres judicieusement choisis.

4 – CAPILLARITE de L’EAU
4.1.- Loi de Jurin
Lorsqu'on plonge l'extrémité d'un tube fin dans un réservoir rempli d'eau, on constate que l'eau
monte à l'intérieur du tube jusqu'à une certaine hauteur h. C'est l'ascension capillaire qui
s'explique par l'attraction des molécules liquides par celles constituant la paroi. On peut
considérer que le phénomène est dû à une tension superficielle qui s'exerce sur le ménisque le
long de la ligne de contact avec le tube. Ecrivons les conditions d'équilibre :

2
d
πd.T.cosα = π. .h.γ w
4
h=π. 4T cosα
d.γ w
Or la tension superficielle de l'eau a fait l'objet de
nombreuses mesures. Sa valeur est voisine de 74,7
dyn/cm.
0,3×cosα ⎧h, d en cm
Alors h= avec ⎨ α en rd
d ⎩
Dans le sol, le phénomène est identique sans être aussi simple car les vides forment un réseau
complexe ; mais tous les vides ne sont pas occupés et même si le sol est saturé, il n'en est pas de
même avec la frange capillaire. La hauteur totale d'ascension dépend donc de l'indice des vides,
⎛ 0,1 < c < 0,5 cm 3
et de la granulométrie. On a sensiblement : h= c ⎜⎜
e.d10 ⎝ h, d10 cm

4.2.- Moule kh :
Si les sols fins ont une ascension capillaire importante ils ont aussi une faible perméabilité.
D'après Terzaghi, la valeur de d10 étant 20 µm, on conçoit que l'ascension capillaire est maximale
au bout de 24 h. On pourra donc déterminer plus facilement la vitesse d'ascension capillaire qui
dépend du produit kH ; d'où le nom de l'appareil utilisé : moule kH. Un échantillon cylindrique
dont on connaît la porosité est préalablement desséché et pesé. L'eau se propage dans le sol
sous l'action de la succion capillaire h en suivant la loi de Darcy. A l'instant t l'eau est au niveau
z. La particule dans le ménisque à l'altitude z est donc à une pression u = -h.γw par rapport à la
pression atmosphérique puisqu'il y a succion. D'après le théorème de Bernouilli, sa charge est
alors : h'= u + z or u = -h.γw ⇒ h’ = z - h
γw
A la base du moule, les particules sont à une charge nulle. La perte de charge est donc égale à la
charge initiale moins la charge finale soit : 0 – (z – h) = h – z ⇒ i = h −1 ⇒ q = A.k.( h −1)
z z
Le débit capillaire q fait monter l'eau à une hauteur dz telle que
q.dt = A.n.dz (n : porosité) ⇒ q.dt = A.k.( h −1).dt = A.n.dz
z
k .dt = z.dz ⇒ log( 1 )− z = k .t
n h− z 1− z h n.h
h
2
z
Si z <<< h on fait un développement limité du Log ce qui donne : 2 )= k .t ⇒ z = 2.k.h .t . Le
n .h n
z.h
volume d'eau absorbée V à l'instant t est donc : V = n.A.z. En effectuant différentes
pesées, on peut déterminer les volumes absorbés. On trace les résultats en fonction du t en
remarquant qu'au début la loi de Darcy ne s'applique pas car i est trop grand. On détermine
graphiquement V0 :
V −V0 = z = A. 2n.k.h.t
kH = z =(V −V0 ) . 1
2

A 2n.t
On peut classer les sols en 3 catégories :
• - sols à fortes remontées capillaires
kH > 1
• - sols à remontées moyennes
0,1 < kH < 1
• - sols à faibles remontées capillaires
kH < 0,1

Remarque
On a supposé que l'échantillon initialement sec était saturé après la remontée capillaire. Ce n'est
pas toujours vrai. Il faut alors remplacer n dans les formules par une fonction de wf - wi (wf et
wi étant les teneurs en eau initiales et finales).