ISSN 0458-5860

BULLETIN DE LIAISON DES LABORATOIRES DES PONTS ET CHAUSSÉES

Digues et barrages en terre

de faible hauteur

MINISTERE DE L'URBANISME, DU LOGEMENT ET DES TRANSPORTS

LABORATOIRE C E N T R A L DES PONTS E T CHAUSSÉES
58, boulevard Lefebvre- 75732 PARIS C E D E X 15
Tél. (1) 532-31-79 - Télex : LCPARI 200361 F

1984
Les articles rassemblés dans ce document sous la présentation de MM. Pilot et Cartier
ont été publiés dans trois numéros du Bulletin de liaison (130, 131 et 132), comme
précisé dans le sommaire détaillé ci-dessous. Ils ont conservé, en conséquence, leur pagi-
nation d'origine à l'intérieur de chacun de ces trois numéros du Bulletin de liaison.

Pages

Présentation générale : G. Pilot, G. Cartier 5

6
Bull. n° 130 • Calcul automatique des digues et barrages en terre
mars-avril J.-L. Bangratz, H. Josseaume, B. Lebrun
1984 • Exemples de problèmes spécifiques relatifs à des digues de faible hauteur 15
C. Mieussens

Bull. n° 131 • Les bassins de Maurepas - Courance et des Roussières 5

mai-juin Présentation : J.-P. Bogaert, J. Nallet
1984
• Constatations sur la digue de Maurepas-Courance 7
A. Delfaut, J. Jardin, R. Baldit
• Interprétation des mesures de pression interstitielle dans la digue de 23
Maurepas-Courance.
Anisotropie de perméabilité de l'ouvrage
H. Josseaume, C. Mahdavi, J. Jardin
• Constatations sur la digue des Roussières 29
J. Jardin, R. Baldit, A. Delfaut

Bull. n° 132 Construction et comportement d'une digue en terre à haute altitude 13

juil.-août M. Lagrange, J. Perrin, H. Josseaume
1984
Application de la méthode des éléments finis aux écoulements 21
en milieu poreux
P. Humbert

Abstracts 38
 Digues et barrages en terre
de faible hauteur

Présentation générale

Georges PILOT Gilles CARTIER

Chef de la division Géotechnique - Mécanique des sols 1 Chef de la section des ouvrages en terre
Laboratoire central des Ponts et Chaussées Laboratoire central des Ponts et Chaussées

Les collectivités locales sont assez fréquemment amenées à préparer des projets de
réalisation de bassins de retenue, que ce soit pour des travaux hydrauliques de
stockage d'eaux de pluie ou Secrètement de crues, pour la création de réserves
d'eau, ou encore dans le cadre daménagement de plans d'eau d'agrément. Les
dimensions généralement modestes des ouvrages en terre réalisés pour ces projets
interdisent la plupart du temps que Ton y consacre un budget d'étude très important.
Le comportement des digues s'apparente pourtant à celui que Ton connaît pour les
grands barrages, et le risque lié à leur stabilité reste élevé, notamment en zone
urbanisée.

Les laboratoires des Ponts et Chaussées ont souvent eu à se pencher sur ce dilemme
afin de tenter d'ajuster le besoin d'études de comportement et de stabilité aux
dimensions des ouvrages et, partant, au volume des crédits disponibles. Profitant
de Texpérience acquise au cours de ces nombreuses interventions, un effort de
synthèse a été récemment entrepris dans le cadre d'un programme d'études générales
des laboratoires des Ponts et Chaussées sous ï'autorité de M. Josseaume.

On a ainsi analysé le comportement de quelques ouvrages types à partir des

observations faites durant les différentes phases de travaux, et des mesures effectuées
pendant et après la construction. Les études théoriques ont, par ailleurs, porté sur
les conditions d'écoulement hydraulique dans les corps de digues après la mise en
eau.

On a, d'autre part, synthétisé Texpérience du réseau des LPC sous forme d'un
document de Recommandations qui sera bientôt mis sous presse et dont la sortie
sera annoncée dans un prochain numéro du Bulletin de liaison.

La série d'articles qui suit, et dont la publication s'étendra sur plusieurs numéros,
s'inscrit dans ce contexte. On trouvera ainsi la présentation des outils généraux
nécessaires à l'étude théorique des écoulements et à leur prise en compte dans les
analyses de stabilité, puis Texposé d'études de quelques digues illustrant les techni-
ques et les moyens mis en place pour suivre le comportement réel de ces ouvrages.
 Calcul automatique des digues
et barrages en terre
Jean-Louis BANGRATZ
Ingénieur
Laboratoire régional de l'Est parisien (Melun)

Henri JOSSEAUME
Attaché de recherche

Bernard LEBRUN
Technicien supérieur
Section des ouvrages en terre
Laboratoire central des Ponts et Chaussées

INTRODUCTION ET HISTORIQUE

L'étude d'un projet de digue en terre nécessite la vérification

de la stabilité de l'ouvrage aux principaux stades de son
histoire et notamment lorsque ses conditions de service sont
les plus défavorables, c'est-à-dire lorsque l'eau de la retenue
RÉSUMÉ atteint son niveau maximal. L'étude de stabilité correspon-
Le programme conversationnel DIGUE a pour dante comporte généralement l'étude de l'écoulement perma-
but de faciliter la mise en œuvre du programme nent à travers la digue et sa fondation et l'étude de stabilité
de calcul hydraulique NSAT et du programme en rupture circulaire du talus aval, cette dernière étant faite
de calcul de stabilité de pente PETAL, modifié
sous le nom STADIG, dans plusieurs configu- à partir des contraintes effectives et prenant en compte les
rations de digue hétérogène sur sol bicouche. pressions interstitielles obtenues à partir de l'étude d'écoule-
La digue est composée de quatre parties (talus ment.
amont, centre, talus aval, drainage aval) dont
chacune comporte plusieurs variantes. Leur Initialement, la procédure suivie dans les laboratoires des
assemblage permet de représenter un assez
grand nombre de configurations différentes Ponts et Chaussées pour l'étude de la stabilité en service
pouvant avoir pour éléments : des risbermes était la suivante : étude de l'écoulement par analogie électri-
amont et aval, un noyau avec ou sans drain- que (méthode du papier conducteur), relevé des charges
cheminée, un tapis drainant aval, un massif de
pied, un parafouille amont ou central, des puits hydrauliques en divers points du modèle analogique, calcul
de décharge à l'aval. des pressions interstitielles correspondantes et introduction
Les modifications du programme PETAL ont manuelle de celles-ci dans le programme de stabilité P E T A L
porté sur la définition des cercles de glissement ou de sa version originale N I X E S et T R O L L S .
potentiels et sur l'amélioration des performan-
ces dans la recherche des valeurs de la pression
interstitielle le long des cercles de glissement A partir de 1979, année où un projet de digue a nécessité
potentiels. l'étude d'écoulements complexes, ces derniers ont été calcu-
lés par la méthode des éléments finis au moyen des program-
MOTS CLÉS : 42 - Programme de calcul - mes F L O T puis N S A T [P. Humbert, 1984]*.
Barrage en terre - Stabilité des talus • Hydrau-
lique (fluide) - Surface de glissement - Pression
interstitielle -/Digue - NSAT - PETAL -
* HUMBERT P., Étude des écoulements par la méthode des éléments
STADIG.
finis. Programme F L O T - N S A T , Bull, liaison Labo. P. et Ch.,
1984, à paraître.

6
B u l l , l i a i s o n L a b o P. et C h . 130 - m a r s - a v r i l 1 9 8 4 - Réf. 2 8 8 5
 À

Afin de faciliter leur utilisation, trois programmes de — digue avec tranchée remblayée (dite encore : « clé
génération des modèles ont été successivement établis. d'ancrage ») et noyau, la digue étant équipée ou non
d'un drain-cheminée;
Deux de ces programmes, D I G U E 1 et D I G U E 2, — massif de pied aval prolongé sous la digue par un
associés au programme F L O T , permettent de mailler tapis drainant ;
automatiquement les modèles hydrauliques associés
aux configurations suivantes : digue homogène fondée — talus amont à risberme, laquelle peut représenter
sur un substratum imperméable et digue homogène à la limite un tapis amont étanche ;
fondée sur un sol formé de deux couches de caractéris- — talus aval à risberme ;
tiques hydrauliques différentes. — parafouille amont ou central;

Le troisième (et le plus récent), D I G U E 3, qui est — ligne de puits de décharge aval ;
associé au programme N S A T , permet le maillage des — en amont : configuration de « réservoir en
modèles précédents et de modèles beaucoup plus com- remblai » ou de « réservoir en déblai dans la première
plexes tels que ceux élaborés pour l'étude de digues couche de sol ».
homogènes ou à noyau, fondées sur un sol bicouche
partiellement interrompu par un écran parafouille ou
drainé par une ligne de puits de décharge incomplets.
SCHÉMA D E PRINCIPE :
L a substitution d'une méthode numérique d'étude des
L E S T R O I S ÉTAPES L O G I Q U E S
écoulements à la méthode analogique classique a per-
mis l'introduction automatique des valeurs de la pres-
sion interstitielle dans le programme de stabilité et, On a rappelé, en introduction, que le calcul d'une
de ce fait, a offert la possibilité d'en introduire un digue comporte deux phases : le calcul de l'écoulement
grand nombre. Il est alors devenu possible de définir dans la digue et dans sa fondation, et le calcul de sa
la distribution de pression interstitielle prise en com- stabilité vis-à-vis d'une rupture circulaire.
pte dans les calculs de stabilité par l'ensemble des
valeurs de ce paramètre, calculées aux nœuds du mail- Les programmes généraux correspondants, N S A T et
lage du modèle hydraulique. P E T A L , ont été conçus indépendamment l'un de l'au-
tre par des équipes différentes. E n conséquence, le
Toutefois, le programme P E T A L n'ayant pas été même modèle est décrit dans des langages différents
conçu pour traiter un champ de pression interstitielle propres à chacun des deux programmes.
défini par de très nombreuses valeurs, il s'est avéré
indispensable d'améliorer ses performances sur ce Prenons l'exemple de la géométrie de la digue. Pour
point. D'autre part, tant dans le but de diminuer le le programme N S A T , elle est décrite sous la forme
coût des calculs de stabilité que de faciliter l'exploita- d'un maillage d'éléments triangulaires affectés chacun
tion de leurs résultats, un nouveau mode de définition d'une perméabilité et d'un coefficient d'anisotropie.
des cercles de glissement potentiel a été adopté. O n a Pour le programme P E T A L , c'est le contour de la
ainsi constitué un programme dérivé de P E T A L , le digue et les intercouches de sols qui sont décrits par
programme S T A D I G qui est maintenant utilisé pour segments affectés chacun du numéro de sol qu'il sur-
les études de stabilité d'ouvrages hydrauliques (digues monte.
de retenue, digues de canaux, éventuellement talus de
bassin). L a mise au point de ce programme a été Afin de simplifier l'utilisation de ces programmes dans
faite pour l'essentiel, par P. Lebreton, ingénieur des le cas particulier du barrage ou de la digue en terre,
Travaux publics de l'État, au cours d'un stage au il était naturel de ne faire qu'une seule fois cette
Laboratoire central des Ponts et Chaussées. description.

Dans cet article, on présente le système de program- On distingue donc les trois étapes logiques schémati-
mes de calcul regroupant N S A T , D I G U E (*), S T A - sées à la figure 1 :
D I G ainsi que différents programmes de liaison,
actuellement utilisés pour les études courantes de — Étape I : description géométrique, hydraulique
digues en terre. et mécanique de la digue.
— Étape II : calcul de l'écoulement, notamment
sous la forme d'un réseau de pression
interstitielle.
LES OPTIONS D U P R O G R A M M E « DIGUE »
— Étape III : calcul de la stabilité à la rupture
circulaire.
De façon exhaustive, les options de D I G U E sont les
suivantes : Idéalement, l'intervention de l'utilisateur ne devrait se
faire qu'à l'étape I. E n pratique, comme le montre la
— sol de fondation mono ou bicouche ou sol de figure 2, les trois étapes logiques ont été par nécessité
fondation non pris en compte dans les calculs (cas décomposées en treize opérations effectuées au cours
d'une fondation rigide et imperméable) ; de quatre interventions conversationnelles, schémati-
sées par une console. L a première, c'est-à-dire la des-
cription de la configuration géométrique est prépondé-
(*) c'est-à-dire D I G U E 3 maintenant utilisé pour le mail- rante.
lage des configurations initialement traitées au moyen de
D I G U E 1 et D I G U E 2. On examinera dans ce qui suit le détail de ces étapes.

7
 ETAPE I :
DESCRIPTION D E L A DIGUE

L a description géométrique de la digue se fait en deux

temps :
— description qualitative de la configuration ;
— description quantitative de la géométrie.
PETAL
Le modèle de digue est divisé en quatre sous-modèles
comportant chacun plusieurs variantes (fig. 3) :
— A , talus amont (six variantes),
DESCRIPTION : — B, partie centrale (cinq variantes),
— géométrique
— hydraulique
— C , talus aval (six variantes),
— mécanique — D , drainage aval (quatre variantes).

Par une série de questions auxquelles on répond par

DIGUE oui ou par non, telles que : « L a digue a-t-elle un
noyau?, un parafouille amont? »... etc., le programme
sélectionne la combinaison A B C D parmi les i J k l

720 combinaisons possibles.

Géométrie (maillage)

Cercle
Conditions aux limites
NSAT
Pression
Interstitielle
PETAL critique Le deuxième temps de la description consiste à préci-
perméabilités
Cercles
ser les coordonnées des quelques points caractéristi-
ques de la figure que l'on vient de décrire qualitative-
ment. O n remarquera sur les figures 3 et 4 plusieurs
natures de points :

— les points noirs doivent être obligatoirement défi-

Fig. 1 — Schéma de p r i n c i p e d e la c h a î n e « D I G U E ».
nis par une abscisse et une ordonnée non nulles ;
— les points blancs sont facultatifs;

Fig. 2 — Les t r o i s étapes l o g i q u e s :

PRÉPARATION Dessin coté de la digue avec les approximations

DES C A L C U L S : I — D e s c r i p t i o n d e la d i g u e .
nécessaires pour entrer dans le cadre des sous-modèles
A1 a A 6 (talus amont) B1 à B5 (partie centrale) II — C a l c u l d e l'écoulement
C1 à C6 (talus aval) D1 à D 4 (drainage aval) III — C a l c u l de la s t a b i l i t é .

Elles s o n t d é c o m p o s é e s e n treize o p é r a t i o n s
regroupées, comprenant quatre interventions
conversationnelles (1, 2, 3 et 5) et deux
1/ Description de la
étapes a u t o m a t i q u e s (4 et 6).
T. ] configuration Mémorisation des
descriptions 1/ et 2/
11 Description géométrique

Fichiers :
« Super-éléments » et
if Conditions aux limites »
3/ Subdivision des super-
éléments

(Console graphique obligatoire)

Matrice de numérotation
des éléments
Coordonnées des nœuds
du maillage
4/ Titre
Affichage puis correction
5/ Données mécaniques
6/ Conditions aux limites
Surfaces à débit calculé
TX
\ eiPerméabilités 107 Couches ^
9/ Sorties graphiques 11/ Points
d'émergence
12/ Pied de talus
13/ a.ß.y

CERCLE

[MECA.j iGEOM.ilCERCLESl

10o\ 5 y m
STADI G

Résultat de l'analyse (PETAL modifié)

de stabilité

8
 B1

Ff
C1
A1 D1

D2
B2 C2
A2

D3
B3 C3
A3

D4
B4
C4
A4
Modèle D : drainage aval
D4 : pas de drainage aval

Zi
C5
A5
/ Z
3

B5

Modèle B : partie centrale

/ Zi B5 : pas de partie centrale
C6
A6
/ 2
3

Modèle A : talus amont

Modèle C : talus aval
A l , A4 : canal en déblai
A2, A5 : canal en remblai

Fig. 3 — Les quatre sous-modèles et leurs variantes. Chaque région de — les points marqués X doivent avoir obligatoire-
chaque sous-modèle est subdivisée par le programme en super-éléments
(groupements d'éléments de calcul) n o n représentés sur la figure.
ment une abscisse non nulle mais leur ordonnée est
facultative. Une coordonnée facultative, si elle est
donnée non nulle garde la valeur donnée, sinon elle
prend une valeur calculée par défaut en fonction de
données déjà connues du programme. C'est le cas de
l'ordonnée d'un point sur l'horizontale d'un point
amont déjà défini.

Cette description, qui est la plus longue des quatre

D4
étapes conversationnelles, est finalement assez simple
à faire. Elle a pour résultat (fig. 2) :
AVAL
— le fichier « S U P - E L E M » décrivant le super-
maillage pour N S A T ;
— le fichier « L I M I T E S » définissant les limites à
charge ou débit imposés ;
— le fichier « G E O M » définissant la géométrie de
Fig. 4 — Assemblage des quatre parties de la digue (A6, B5,
C6, D 4 est la plus simple combinaison possible). Si le parafouille la digue pour le programme S T A D I G .
amont ou le tapis drainant aval n'est pas déclaré lors de la
description qualitative, les super-éléments correspondants exis- U n exemple simple de constitution de maillage est
tent malgré tout mais il leur est affecté les caractéristiques
hydrauliques et mécaniques des milieux qui les contiennent. donné à la figure 4.

9
E T A P E II : Les limites à charge ou à débit imposés sont connues
dès l'étape I. Elles sont rappelées sur la console et
CALCUL D U CHAMP
l'on a, à ce stade, la possibilité de les contrôler et de
DE PRESSION INTERSTITIELLE
les modifier éventuellement. Par exemple, la limite
gauche du modèle, segment A B sur la figure 6, est
Cette étape logique se décompose en deux étapes considérée implicitement par le programme comme
une limite équipotentielle. Il est possible, à cette
conversationnelles marquées [T] et [3] sur la figure 2 étape II, de changer sa nature et de la déclarer comme
et une étape automatique marquée [TJ. Il s'agit de : limite imperméable.

— la subdivision des super-éléments créés par le pro-

gramme D I G U E en éléments plus petits, Charrie = h _ D E

4
y
. Charge = 0
— l'affichage de certaines données connues du pro-
gramme avec la possibilité de les modifier.

L a subdivision des super-éléments en éléments simples

c F ;
I
1
est faite au moyen du programme « M E S H ». Ce
Limite imperméable H
programme est utilisé par le service de Mécanique du
L C P C pour créer des maillages généraux de super-
éléments que l'on subdivise ensuite avec optimisation Fig. 6 — Conditions hydrauliques aux limites :
de la largeur de bande de la matrice de rigidité. O n — charge imposée : h sur ABCD et 0 sur FIJ et GH ;
— surface de suintement : DEF;
n'utilise ici que sa fonction de subdivision automati- — surfaces où l'on calcule le débit : débit total entré ABCD et puits de
que. décharge GH.

L a •figure 5 montre que, en raison de la propagation

de la subdivision entre éléments adjacents, tous les
côtés de super-éléments ne sont pas indépendants :
sur l'exemple simplifié de la figure 5 (une digue homo- D'une façon analogue, les perméabilités et les caracté-
gène sur fondation homogène et tapis drainant aval) ristiques mécaniques des sols ont des valeurs prédéter-
on dénombre sept super-côtés indépendants. minées dans le programme. Elles sont rappelées sur
la console, et il est bien sûr possible de les changer
A ce stade, le maillage est entièrement construit. Reste dans cette étape. Par exemple, tous les poids volumi-
à définir les limites à charge et à débit imposés, les ques sont de 20 k N / m par défaut; les perméabilités
3

perméabilités et les caractéristiques mécaniques des sont toutes isotropes, celle des drains, par exemple,
différents sols : poids volumique, cohésion et angle de est de 1 0 " m/s; la cohésion du corps de digue est de
2

frottement interne. 5 kPa, son angle de frottement interne est de 30°, etc.

A u terme de cette étape, les données mécaniques pour

S T A D I G sont prêtes dans le fichier « M E C A » et
les données du programme N S A T sous leur forme
standard sont prêtes dans le fichier « D O N - N S A T ».
Il n'y a plus qu'à déclencher l'exécution du pro-
gramme qui a pour résultat la génération du fichier
du champ de pression interstitielle « P.I. » et la repré-
sentation graphique de l'écoulement sous la forme
des réseaux de lignes équipotentielles et de lignes
équipressions (fig. 7).

(b)
ÉTAPE III : C A L C U L D E L A STABILITÉ

Les données géométriques, mécaniques et hydrauli-

ques sont maintenant préparées dans les formats du
programme P E T A L . A u lieu d'être dans un fichier
unique (programme P E T A L standard), elles sont ven-
tilées dans des fichiers distincts pour les raisons indi-
quées plus haut. Cela nécessitait un simple change-
1
/ / ' V Drain \ . 7
ment de numéro des unités d'entrée dans le pro-
2 gramme P E T A L et n'aurait pas justifié de rebaptiser
(c)
i ce programme.

3 6 Il restait à définir un jeu de cercles potentiels de

Fig. 5 — Subdivision des super-éléments en éléments simples. glissement dont le programme doit extraire les cercles
— Propagation de la subdivision par des super-éléments quadrilatères (s) critiques pour chaque couche de sol. Le mode de
et triangulaires (/>). définition des cercles potentiels a été changé par rap-
— Subdivision des super-côtés indépendants numéros 1 à 7 d'une
digue simple (c) : 2 et 5 en deux segments, 3 en trois segments; les port au programme P E T A L standard. C'est la princi-
autres ne sont pas subdivisés. pale raison d'être du nouveau programme S T A D I G .

10
 • C

L a deuxième raison est l'amélioration des perfor-

mances dans la recherche des valeurs de la pression
interstitielle à la base des tranches de calcul. Cette
amélioration a été rendue possible par une modifica-
tion mineure comme on va le voir maintenant.

M O D I F I C A T I O N S APPORTÉES
AU PROGRAMME PETAL
D
Fig. 8 — Définition d'un cercle de glissement potentiel.
Définition et sélection
des cercles de glissement potentiels
D'autre part, un cercle ne peut être considéré comme
un cercle de glissement potentiel que si son arc situé
Généralités
au-dessous de la surface du sol ne recoupe pas la
verticale en H, c'est-à-dire si l'ordonnée de son centre
Dans la plupart des programmes de calcul de stabilité
est supérieure à celle de H (fig. 10 b).
de talus en rupture circulaire, les cercles de glissement
potentiels sont définis par la position de leurs centres
(qui constituent les nœuds d'un quadrillage régulier), Sélection des cercles de glissement potentiels
par les valeurs maximales de leur rayon (distance des
centres au substratum rigide) et par le pas de variation L a méthode adoptée pour la sélection des cercles de
du rayon. Lorsque la configuration du talus est un glissement potentiels tient compte du découpage du
tant soit peu complexe, et lorsque le choix des centres talus en couches réelles ou fictives limitées intérieure-
et du pas de variation du rayon n'a pas été précédé ment par des plans horizontaux, et ayant pour limite
d'une étude géométrique approfondie (c'est-à-dire supérieure la surface du sol et/ou la limite inférieure
dans la majorité des cas), la sélection des cercles de de la couche sus-jacente. Chaque couche est repérée
glissement est assez arbitraire si bien que les résultats par un indice k dont la valeur augmente avec la
profondeur (fig. 9). O n considère que la limite supé-
de l'étude de stabilité ne sont fiables que s'ils ont été
obtenue partir du calcul de nombreux cercles. C'est
pourquoi un nouveau mode de définition des cercles
de glissement potentiels, permettant de mieux
« visualiser » ceux-ci et de limiter leur nombre, a été
introduit dans le programme S T A D I G .

Définition des cercles de glissement potentiels

k = 3
U n cercle de glissement potentiel est défini par ses
points d'émergence haut et bas, c'est-à-dire par ses
k = A
points d'intersection H et B avec la surface du sol et
par la valeur / de la flèche de son arc situé au-dessous Fig. 9 — Définition et repérage des couches prises en compte
de la surface du sol (fig. 8). pour la sélection des cercles de glissement potentiels.

11
rieure d'une couche, ou une partie de cette limite, fnoB étant la flèche de l'arc HOB et p une longueur
n'appartient à la couche que si elle coïncide avec la très petite, par exemple P = 0,01.
surface du sol. O n considère, d'autre part, qu'un
cercle appartient à une couche lorsque celui de ses — si B n'appartient pas à la couche k (fig. 11c) et
points D ou B, c'est-à-dire son point bas et/ou son s'il existe un cercle C appartenant à la couche
m a x

point d'émergence bas (fig. 8), situé le plus en amont, k — 1, qui ne soit pas tangent à la verticale passant
appartient à la couche. par H, C est le cercle ayant pour flèche :
m i n

/min ( * ) = / m » ( * " ! ) + Y
L a méthode de sélection des cercles de glissement
potentiels passant par H et par B (l'ordonnée de H y étant une longueur très petite, par exemple y=0,01.
étant supérieure à celle de B) consiste à définir les
faisceaux de cercles appartenant respectivement à cha- Si la dernière condition n'est pas satisfaite, il n'existe
que couche dont la limite inférieure se situe au- aucun cercle de glissement dans la couche k.
dessous de B, les caractéristiques de ces faisceaux
étant déterminées dans l'ordre des valeurs croissantes (a)
de k.

Le faisceau des cercles appartenant à la couche k est

limité par les cercles C et C ayant respectivement
m a x m i n

pour flèches / (k) et /

m a x (k). m i n

Le cercle C est le cercle tangent à la base de la

max

couche si l'ordonnée de son centre est supérieure à

(b)
celle du point H (fig. 10 a). Dans le cas contraire, C m a x

est le cercle admettant une tangente verticale en H

(fig. 10 b). Si, dans ce dernier cas, le point bas est dans
la couche k — 1 il s'agit en fait du cercle C (k — 1) et m a x

il n'y a pas de cercle C (k). m a x

(0

Fig. 11 — Définition du cercle Cmin appartenant à la couche k.

(a) B appartient à la couche k et se situe à l'amont du pied de talus.
(b) B appartient à la couche k et se situe à l'aval du pied de talus.
Fig. 10 — Définition du cercle Cmax passant par B et H (c) B n'appartient pas à la couche et il existe un cercle Cmax appartenant
et appartenant à la couche k. à la couche * — 1 .

L a méthode de définition du cercle C dépend de la min Une fois définis les cercles C et C tous les cercles
m a x m i n

situation du point B par rapport à la couche k et par du faisceau appartenant à la couche k sont déterminés
rapport au pied du talus : à partir de la valeur Af (k) du pas maximal de varia-
tion de la flèche dans la couche k. Ce sont les cercles
— si B appartient à la couche k et se situe à l'amont dont les flèches varient suivant une progression arith-
du pied de talus (fig. 11 a), C est le cercle ayant
m i n
métique de raison aussi voisine que possible de Af (k)
pour flèche : et au plus égale à cette valeur (fig. 12).
/ (k) = ct.BH
mln

a étant un coefficient réducteur assez petit pour qu'un L a sélection des cercles de glissement potentiels est
glissement le long de C puisse être assimilé à un
m i n
faite par le programme C E R C L E à partir des données
glissement plan, par exemple a = 0,01. suivantes :

— si B appartient à la couche k et se situe à l'aval — coordonnées des points d'émergence bas B(i) et
du point O représentant le pied de talus (fig. 11 b), des points d'émergence hauts H (j),
C est le cercle ayant pour flèche :
m i n
— coordonnées du pied de talus,
— valeurs Af (k) du pas maximal de variation de la
/min ( * ) = / I K » + P flèche,

12
 1

interstitielle au milieu de la base de chacune des tran-

ches du cercle de glissement potentiel étudié (c'est-à-
dire en des milliers de points). Il faut donc déterminer
à quel triangle du maillage appartient le milieu de la
base de chaque tranche considérée au cours du calcul.

Dans la version actuelle de P E T A L , la procédure

d'identification du triangle contenant un point M
donné, consiste :

— à déterminer le plus petit rectangle à côtés horizon-

Fig. 12 — Faisceau des cercles de glissement appartenant à la couche k.
taux et verticaux contenant tous les points où est
donnée la valeur de la pression interstitielle (rectangle
circonscrit au maillage),

— ordonnées y (k) des limites inférieures de couches, — à découper ce rectangle en douze zones rectangulai-
— valeurs des paramètres a, ß et y communs à toutes res égales (fig. 13) et à recenser les triangles apparte-
les couches. nant à chacune d'elles (un triangle appartient à une
zone lorsqu'au moins un de ses sommets y est situé),
Le programme calcule les caractéristiques (notam-
ment les coordonnées du centre et la valeur du rayon) — à déterminer dans quelle zone se trouve le point
des cercles de glissement potentiels passant par tous M et à étudier la position de ce point par rapport à
les couples de points obtenus en associant à chaque tous les triangles de cette zone.
point d'émergence bas B (i) tous les points d'émer-
gence hauts situés en amont de B (i) et à un niveau Cette procédure conduit à des temps de calcul accepta-
supérieur. bles dans le cas où le domaine maillé est relativement
compact, même si le champ des pressions interstitielles
est défini par une ou deux centaines de point. E n
Traitement du champ des pressions interstitielles revanche, le temps de calcul devient prohibitif lorsque
la procédure est appliquée à un cas où le champ des
pressions interstitielles est défini par un plus grand
Le champ des pressions interstitielles dans le talus est
nombre de points, et lorsque le domaine maillé ne
défini en un certain nombre de points. Ces points
recouvre qu'une faible proportion du rectangle circons-
sont reliés deux à deux par des segments de droite de
crit au maillage (fig. 14).
façon à former un maillage d'éléments triangulaires.
Lorsque le talus est le siège d'un écoulement perma-
nent et que celui-ci a été étudié par la méthode des C'est pourquoi la procédure d'identification du trian-
éléments finis, le maillage retenu est naturellement gle contenant un point M donné a été sensiblement
celui utilisé pour l'étude hydraulique. L a pression modifiée dans le programme S T A D I G . O n détermine
interstitielle en un point M du talus est donné par tout d'abord le rectangle utile c'est-à-dire le plus petit
l'expression : rectangle à côtés horizontaux et verticaux circonscrits
à l'ensemble des cercles de glissement potentiels
u = au +bu + cu
1 2 3
(fig. 14) et non plus à l'ensemble des points de pres-
sion interstitielle. Ce rectangle utile est ensuite divisé
u u , u étant les valeurs de la pression interstitielle
u 2 3
en zones rectangulaires égales dont le nombre au lieu
à chacun des sommets de l'élément triangulaire conte- d'être limité à 12, est approximativement égal à celui
nant le point M et a, b, c les coordonnées barycentri- des triangles du maillage. A ce stade de la procédure,
ques de M par rapport à cet élément. les opérations demeurent les mêmes que dans le pro-
gramme P E T A L :
L'étude de la stabilité du talus par une méthode de
tranches (méthode de Bishop, méthode des perturba- — avant tout calcul : établissement de la correspon-
tions) nécessite de connaître la valeur de la pression dance « zone -» triangles »;

• Point d'émergence haut

* Point d'émergence bas
Rectangle utile

Fig. 14 — Tracé du rectangle utile (plus petit rectangle

circonscrit aux cercles de glissement potentiels).

13
 Fig. 15 — Résultats complets
***CALCUL PAR METHODE JISHO° AVEC TEST RRT*«*
pour chaque cercle.

vvvvv. VVV'l VVVV'V

* * * * * * * * * * . . , * , . . . , * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * *
* * * * * * * * * * *"** * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * *
? O I N T
D ' EMERGENCE 2 AS \ 0 : 1
= OINT 0' £''£R3SVC£ HAUT NO : 1

COUCHE VAS 1 VAR 2

NO : X C E N T s. E YC ENTRE FLECHE RAYON 10<ENT POIOS J1AX FRRT FRRT
1 2 41 .24 ' 4 . 1 5 1 . 705 14.41 .594E+03 .113E+03 . 2 7 0 E + 01 1 .60 4 . 5 9
2 2 3 9.97 20.02 2.409 10.30 .343E+03 .264E+03 . 3 0 0 5 * 0 1 1.32 3.55
3 2 i 9 . 2 j 1 7 . i l 3.113 2. 95 .102E+04 .41oE*03 .142E»02 1 .91 3.26
4 2 33.72 1 5 . ' 7 2.516 7.97 . 1 1 3E-1-04 .57oE»03 .20SS*02 2.00 3.14
CALCUL N C 'i CONV.VAR . U t R CL £ C I 01 E S 0 U S
5 3 3 5. 4o 15.77 3.916 7.36 . 1 1 4E * 0 4 . sOOEH-03 .217E+02 2.01 3.10
0 3 33.36 14.31 4. 470 7.41 .120E+04 .732E+03 .271E+02 2.05 3.11
7 3 3 2.12 14.00 5.023 7.11 . 1 235E + 0 4 .370E+03 .325ET02 2.14 3.11
* * * * * * * * * * * * * * * * * * * * » * * ** * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * *
POINT 3 ' EMERGENCE 2 A S NO : 1
POIMT 5 ' EMER3ENCE HAUT NO: 2

COUCHE VAR 1 VAR 2

NO : X C E M T ti E y c ENTRE FLECHE RAYON •OIENT POIDS UMAX FRRT FRRT
3 2 41 .05 40.71 1 .773 30. 35 . 4 3 1 E + 04 .496E*C3 . 1 4 1 E •"•01 1.37 3.02
9 2 i d . 6 9 31.26 ' . 5 7 6 21.91 .454E+04 . 741E*03 . 680E*01 1 . 92 2.36

Fig. 16 — Récapitulatif des

TAHLEAUX C0^=«»ATIFS
cercles critiques par variante
mécanique.
........*......,,.*.
VARIANTE NO 2
023 VALEURS 32 FRRT " I N I H U I

.,.****,...,......*,.**,.**
PAR PAIRE DE POINTS D'EMERGENCE

* 1 2 : 3 : *..*.,***......*..*.*.*...*****
* 44 .00 4 0 . 0 0 : 3 7 . 0 0 :
EN COLONNE LE S POINTS »A3
* 1 0 . 0 0 1 0 . 0 0 : 1 1 . 0 0 :
t * * * * * * * * * * * * * * EN LIGNE LE S POINTS HAUTS
* * ** ** * * * * * *
3.12 : DANS CHAO UE CAS E : -FRRT
1 * 3.1C 2 . 5 3 :
-N'J.IEPO DU CERCLE
31 . 0 0 * >>>> :
. A A A A j 4 L E «IINTIUM PAR COLONNE/
1 4 .00« 5 1 ô : 23 :
S 0 > J

>>>> SI ON L'A PAR L I G N E /

2 « 2 . 6 6 2 . 6 9 : 2 . 9 3 :
24 . 0 0 * ----- AA A A . A A A A .
==== SI ON L'A PAR LIGNE ET PAR COLONNE
1 5 .00* 13 20 : 31 :

COUCHE FRRT.MIN CERCLE FRRT.MAX CERCLE NO

1 2.90 31 3.92 27
2 2 . 6 9 20 4 . 5 9 1
3 2.66 1 3 4 . 0 0 33

— en cours de calcul, pour chaque base de tranche : coefficient de sécurité F obtenues respectivement pour
• recherche de la zone rectangulaire qui la contient; tous les faisceaux de cercles de glissement potentiels
cette opération est immédiate puisqu'il s'agit de considérés dans le calcul. L a valeur minimale de F
calculer les coordonnées entières d'un point dans relative aux cercles du faisceau passant par les points
un maillage de rectangles égaux; B(i) et est lue (ainsi que le numéro du cercle
• recherche du triangle contenant le point, parmi correspondant) dans la case située à l'intersection de
les quelques triangles de la zone. Puisqu'il y a la colonne i et de la ligne j.
autant de rectangles que de triangles, chaque rectan-
gle intercepte en effet un très petit nombre de trian- Au-dessous de chaque tableau récapitulatif sont égale-
gles. ment indiquées les valeurs minimale et maximale de
F obtenues dans chaque couche ainsi que les numéros
Les essais effectués lors de la mise au point du pro- des cercles correspondants.
gramme S T A D I G ont montré que cette petite modifi-
cation du traitement des pressions interstitielles per-
mettait d'effectuer certains calculs de stabilité en qua-
tre fois moins de temps que celle mise en œuvre dans CONCLUSION
le programme P E T A L .

Le programme D I G U E a été réalisé dans le but

Présentation des résultats d'un calcul de stabilité d'assurer l'interconnexion d'un programme de calcul
d'écoulement hydraulique par la méthode des élé-
Outre le rappel des données du calcul, le listing des ments finis et d'un programme de calcul de stabilité
résultats comporte : de talus à la rupture par la « méthode des tranches ».

— un tableau général où sont reportés les paramètres O n dispose ainsi d'un moyen intégré pour l'étude
géométriques et mécaniques de chaque cercle de glisse- complète de la stabilité des digues ou des barrages en
ment potentiel (fig. 15); terre avec, pour l'utilisateur, le minimum de données
— un tableau récapitulatif par variante de calcul à fournir, sans redondance, pour les deux programmes
(fig. 16) où sont indiquées les valeurs minimales du de calcul.

14
 Digues
ef barrages en terre
07r/////////////fmmUi//, de faible haute

Exemples de problèmes spécifiques relatifs

à des digues de faible hauteur
Claude MIEUSSENS
Ingénieur ENSHG
Laboratoire régional de Toulouse

RÉSUMÉ

L'article décrit quelques exemples de problèmes

spécifiques rencontrés lors de l'étude ou du
suivi après la mise en service de cinq digues
situées dans la zone d'action du Laboratoire
régional des Ponts et Chaussées de Toulouse :
1. le cas d'une digue construite sur un terrain
trop perméable. L étanchement du barrage et
du sol support a été réalisé au moyen d'un
rideau injecté dont l'efficacité s'est avérée
insuffisante compte tenu des faibles débits de
fuite admissibles ;
2. celui d'un barrage de faible hauteur, mais
retenant un volume d'eau très important, pour
lequel l'évolution vers une rupture certaine du
parement aval a pu être mise en évidence grâce
à une instrumentation appropriée ;
3. les moyens mis en œuvre pour étudier la Pour les ouvrages de faible importance, l'ingénieur de
perméabilité du sol support d'un barrage, à bureau d'études est toujours à la recherche d'un équilibre
priori très défavorable d'après la géologie à
cause de la présence aléatoire d'horizons entre le niveau technique et le coût des études et l'impor-
perméables. Pour le corps de la digue, on a tance économique du projet.
mesuré l'influence du traitement à la chaux
sur la perméabilité du matériau compacté ;
Les cinq exemples présentés dans cet article illustrent les
4. la conception même d'un barrage écrêteur difficultés rencontrées à l'occasion des études ou de la
de crue afin qu'il remplisse au mieux la
fonction de stockage maximal. Les problèmes réalisation de digues de faible hauteur correspondant à des
d'érosion interne et de contournement sont ouvrages existants ou en projet pour lesquels, le plus sou-
plus particulièrement examinés ; vent, le Laboratoire régional des Ponts et Chaussées de
5. le cas particulier d'un plan d'eau situé à Toulouse est intervenu au niveau des études.
proximité d'une rivière, ce qui a nécessité
l'étude de la stabilité de l'écran étanche dans
l'hypothèse d'un fonctionnement inversé de la Les méthodes d'étude utilisées sont parfois simplistes, ou
digue (plan d'eau vide, rivière en crue). Ce au contraire disproportionnées avec l'ampleur des problè-
même exemple décrit aussi une simulation mes. L a première cause en est probablement que chaque
d'écoulement plan par analogie électrique.
étude nouvelle présente des caractères originaux devant
Enfin, à travers l'exposé de méthodes d'études lesquels un ingénieur même confirmé ne peut faire jouer
parfois simplistes, l'article pose le problème
de l'adéquation entre le niveau des études et utilement toute son expérience. L a principale difficulté est
l'économie du projet. toutefois d'ajuster ou de savoir exiger les moyens d'étude
les mieux adaptés potentiellement au projet.
MOTS CLÉS : Barrage en terre - Étanchéité -
Sous sol - Écran - Injection (mater.) - Perméabi-
lité - Rupture - Traitement des sols - Chaux • Les ouvrages étudiés sont volontairement décrits de façon
Érosion • Rivière - Simulation - Écoulement anonyme et uniquement du point de vue de l'ingénieur
(fluide) - Modèle analogique - /Plan d'eau -
Bassin écrêteur de crue - Digue. d'étude, afin de fixer l'attention sur les caractères originaux
des problèmes rencontrés en cours d'étude ou après la mise
en service.

15
Bull, liaison L a b o P . et C h . 130 • m a r s avril 1 9 8 4 Réf. 2 8 9 9
E X E M P L E D E P R O B L È M E D'ÉTANCHÉITÉ marne ou de conglomérat très compact et imperméa-
D U CORPS D EDIGUE E TD U S O L SUPPORT ble. O n observera sur la figure 1 la position de la
(barrage de C) nappe, quasiment horizontale et parallèle au substra-
tum en fond de talweg, ce qui met bien en évidence
l'homogénéité et la forte perméabilité des alluvions
Le barrage de C fait partie d'un aménagement touris- entre les forages F l et F4.
tique et résidentiel réalisé pour une commune. De
dimensions modestes, la digue de 6 m de haut au Malgré ces conditions défavorables, il était décidé de
maximum et 80 m de long retient un volume d'eau donner suite au projet; plusieurs solutions étaient
de 30000 m sur une surface de l'ordre de 2 hectares.
3
donc envisagées et proposées au maître d'oeuvre.

Les études préliminaires ont mis en évidence un triple

problème mettant en cause la faisabilité du projet : C A I L L CX X GRAVI ERS GROS SABLE S A B L E FIN
100

— le sol support est perméable, 90

— on ne trouve pas de matériaux imperméables pour 80

construire la digue à une distance compatible avec * \

I
l'économie du projet,
3
70
\
— l'alimentation pendant la période estivale peut être
quasiment nulle.
N
Une reconnaissance approfondie du site au moyen de
sondages a permis d'établir le profil géotechnique
(fig- 1).

Des essais de perméabilité de type Lefranc dans les

différents forages ont donné des valeurs de k de l'or- 200 100 50 20 10 5 2 1 0.4 0,2 0.1 50
dre de 2 . 1 0 ~ m / s pour les graves sableuses. L a
4 O mm U ^ ^
courbe granulométrique de la figure 2 précise la Diamètre équivalent (^m)

nature du matériau. L e substratum est constitué de Fig. 2 — Analyse granulométrique des alluvions et f o n d de talweg.

Fig. 1 — Profil en long géotechnique au droit de la digue de C.

16
 i

Constitution du corps de la digue de l'efficacité du tapis étanche peut être trouvé en

écrivant que le débit de fuite est réduit proportionnelle-
L'absence de matériaux imperméables pouvant consti- ment à la longueur du tapis étanche, soit :
tuer le corps de la digue a naturellement conduit
à proposer deux solutions permettant d'utiliser les
alluvions sablo-graveleuses du talweg :
Q : débit de fuite en l'absence de tapis étanche,
— la réalisation d'une étanchéité amont au moyen Q : débit de fuite avec le tapis étanche,
T

par exemple d'une feuille plastique (butyl, polyester, B : largeur de la digue à la base,
etc.); L : longueur du tapis étanche.
— l'exécution d'un écran étanche vertical avec un
coulis de ciment bentonite. O n a calculé ainsi qu'il fallait un tapis de 80 m de
long pour diviser les débits de fuite par 10.

Étanchement du sol support Dispositions constructives proposées

L a figure 1 montre que l'étanchement du sol support Pour des raisons économiques, il paraissait préférable
est nécessaire au droit de la digue mais aussi pour d'adopter une solution unique à la fois pour la digue
limiter les débits de fuite par contournement. Par et pour le sol support. O n pouvait donc imaginer
ailleurs la nature sablo-graveleuse des alluvions est deux solutions principales :
favorable à l'érosion interne pouvant à terme conduire
à la ruine de l'ouvrage. — un film plastique étanche unique à la fois pour la
cuvette et comme étanchéité amont de la digue,
Les débits de fuite sous-jacents à la digue ont été — un voile mince étanche, avec par exemple un coulis
calculés à partir de l'abaque de la figure 3. Pour ciment-bentonite, réalisé depuis la crête du barrage
l'ensemble de la digue, sans tenir compte des contour- jusqu'au substratum imperméable.
nements, on a calculé ainsi Q x 6 l/s. Cette valeur n'est
pas admissible car elle correspondrait, en été, à une Dans les deux cas, le corps de digue peut être exécuté
baisse du plan d'eau de l'ordre de 1 m/mois. avec les alluvions de la cuvette sans qu'il soit néces-
saire de prévoir un tapis drainant à l'aval.
L a figure 3 permet de vérifier que pour b/T=\,5 et
s/T=0,9 le débit de fuite est seulement divisé par 2, A u niveau du marché, la solution voile mince ciment-
une coupure partielle de la couche perméable est donc bentonite a été retenue avec les critères suivants :
sans intérêt. — épaisseur : e'=0,5 m,
1,6
— perméabilité : k' = 10~ m/s.
s

Compte tenu de ces valeurs, les débits de fuite à

travers la digue et le sol support sont divisés par 500,
soit 0,022 l/s au lieu de 11 l/s. Dans le marché, ce
faux de réduction des débits de fuite a été ramené à
100.

Dispositions constructives adoptées

et problèmes de fuites
L a solution finalement retenue a consisté à réaliser
un écran injecté, avec une seule file d'injection et un
espacement de 2,5 m entre points. Le coulis était
constitué d'un mélange ciment-bentonite.

Pour satisfaire les termes du marché il fallait donc

que l'injection réduise la perméabilité de 2 . 1 0 ~ m/s
4

à 2 . 1 0 " m/s sur une largeur de 2 m. Pour des allu-

7

0 0.1 0,2 0,3 0,4 0,5 0,6 0,7 0,8 0.9 1.0 vions, un tel résultat ne peut être obtenu qu'avec des
s/T
gels de silicate et à condition de réduire l'espacement
Fig. 3 — Abaque de calcul des débits de fuite sous une digue étanche
avec o u sans parafouille (d'après Polubarinova Kochina, Theory of ground des forages d'injection disposés sur deux lignes ou
water movement, 1962, Briston University Press). même trois lignes.

Dans l'ouvrage « Injection des sols », H . Cambefort *

Deux solutions ont donc été proposées : indique que, pour les alluvions, l'expérience montre
que l'injection « permet simplement de réaliser des
— une coupure étanche ancrée dans le substratum
écrans dont la perméabilité moyenne est comprise entre
imperméable et prolongée en rive droite jusqu'au
1 0 " et 1 0 " m/s », soit vingt-cinq fois plus que la
6 5

niveau du sondage F ;
valeur requise.
5

— un étanchement de bassin de la retenue au moyen

par exemple d'un film plastique. Cette solution n'a
pas été étudiée en détail, mais un ordre de grandeur * Eyrolles, 1964.

17
Lors de la mise en eau, des fuites sont rapidement
apparues en plusieurs points à l'aval. Ces différentes
Remblai
fuites ont été collectées, ce qui a permis d'évaluer les
débits à 1,5 1/s environ. C'est une valeur relativement Alluvions
~\ Rivière
faible qui pourrait être admissible si l'alimentation en
Surface de r u p t u r e étudiée
étiage était suffisante. Par ailleurs, compte tenu de
Q 1 ? ? 1

la nature du matériau, on peut toujours craindre le

développement d'une érosion interne dans la zone de
cheminement préférentiel de l'eau.
V
///////>^T?
_x N i v e a u p i é z o m é t r i q u e après r é a l i s a t i o n des forages
N a p p e c o n s i d é r é e d a n s les c a l c u l s
Une nouvelle ligne d'injection a donc été réalisée, en
quinconce avec la première, et un système permanent Fig. 4 — Profil en travers de la digue. Position de la surface de rupture
de contrôle des débits de fuite a été mis en place. O n observée et de la surface libre.
peut penser que ces injections complémentaires ont
eu peu d'influence sur la perméabilité moyenne de la
zone injectée; cependant, les premières observations
semblent mettre en évidence une diminution des fuites
L a figure 4 donne les résultats de calculs de stabilité
les plus importantes, mais les dispositions adoptées
qui ont permis de situer la surface de rupture les plus
sont trop récentes pour que l'on puisse en déduire
défavorables passant par les zones de mouvements
des conclusions plus précises.
observés. O n notera, sur cette même figure, la position
élevée de la surface libre malgré la présence du tapis
drainant.

E X E M P L E D E P R O B L È M E D E STABILITÉ Pour expliquer le glissement, on peut proposer l'ana-

D U PAREMENT AVAL D'UN BARRAGE lyse suivante : le phénomène de vidange rapide de la
(barrage de B) rivière, dû à la décrue, est probablement à l'origine
des mouvements.

Le barrage de B correspond à un aménagement touris- Avant la construction du barrage, les berges ont certai-
tique de conception assez particulière. Le plan d'eau, nement été soumises à des « vidanges rapides » sans
d'une superficie de l'ordre de 30 ha, n'a pu être réalisé que des désordres se soient déclarés. O n peut donc
sur un terrain pratiquement horizontal que grâce à la supposer que c'est la modification du régime hydrauli-
construction d'une digue périphérique de 3 m de haut que à l'aval du plan d'eau qui a déclenché la rupture.
en moyenne et de près de 2 km de circonférence. L'évolution régressive des fissures confirme bien cette
Les matériaux assez argileux et imperméables ont été hypothèse. Le « calage » des caractéristiques mécani-
prélevés dans le bassin même de la retenue. Des essais ques pour F = 1 sur la surface de rupture de la figure 4
de perméabilité en place avaient montré que le creuse- conduit à c' = 0 et cp' = 18°.
ment était possible sans risquer d'atteindre des hori-
zons plus perméables. Ce sont des valeurs très faibles pour ce type de maté-
riau, qui ne peuvent s'expliquer que par une évolution
Compte tenu de la nature des matériaux et de la vers des caractéristiques résiduelles à cause des gran-
faible hauteur de la digue il n'y avait a priori aucun des déformations.
problème majeur, en particulier, de stabilité.
L'hypothèse d'un phénomène de « renard » ne doit
Mais on le voit sur la figure 4, dans un souci d'aug- pas non plus être exclue. Une analogie électrique des
menter au maximum la surface du plan d'eau, la écoulements a permis d'estimer le gradient hydrau-
digue a été rapprochée le plus possible des berges lique à 1,5 environ au voisinage de la rivière, ce qui
d'une rivière qui contourne l'aménagement, l'alimenta- est très supérieur aux valeurs maximales conseillées
tion de la retenue se faisant par un affluent. dans les argiles, soit 0,6 à 0,7. Ainsi le phénomène de
renard peut expliquer l'amorce d'un glissement de
Après une forte crue de la rivière, des fissures ont pied qui serait à l'origine de la rupture progressive
commencé à apparaître sur les berges, puis elles ont observée.
progressé au cours du temps vers le talus aval de la
digue, tout en prenant une extension de plus en plus Cet exemple illustre bien l'importance des instabilités
grande en plan. locales, car il est tout à fait probable qu'un calcul de
stabilité générale sur la base de caractéristiques de
Afin de préciser l'ampleur et l'évolution des désordres, cisaillement mesurées en laboratoire aurait conduit à
une instrumentation a été mise en place. Elle compor- des valeurs tout à fait admissibles du coefficient de
tait (fig. 4) trois tubes piézométriques et un inclino- sécurité.
mètre ainsi que des nivelles qui permettaient un suivi
journalier des mouvements. Plusieurs systèmes de confortation ont été imaginés,
en particulier le drainage des berges au moyen d'épe-
L a figure 5 met bien en évidence l'existence d'une rons ainsi qu'un allongement du tapis drainant en
surface de rupture vers 7 m de profondeur et une pied de la digue. L'amélioration du coefficient de
évolution rapide en fonction du temps pouvant entraî- sécurité dans les zones ayant glissé n'est que de
ner la ruine de l'ouvrage. A u vu de ces résultats, dès 20 %. L a solution finalement retenue a consisté à
le 24 juillet 1979 ( l mesure), la décision de baisser
r e
éloigner la digue de 30 m environ de la crête de talus
le plan d'eau était prise. de la rivière.

18
 Fig. 5 — Relevés inclinométriques
MINISTERE DE L'ÉQUIPEMENT FORAGE I 1 en fonction du temps.
CETE DU SUD-OUEST
L A B O R A T O I R E DE TOULOUSE PROFIL INCLINOMÉTRIQUE

ee Description et légende Variations angulaires du tube (degrés)

géotechniques
• 2,00 1,00 2,00

Limon brun •M-

0,5- légèrement argileux
"''/ir
1,0- 1,00 'A''?-

1,5-
Argile limoneuse
brun clair assez
2,0- compacte
2,30 it**-
2,5- 5^

-r
3,0- Argile limoneuse
brun clair (molle et
3,5- saturée de 2,30 à
3,90 m)
4,0-

4,5-

5,0-
4,60
-r
Limon argileux et
vaseux gris bleu
5,5"
consistant

6,0- - •sa
6,5.-

7,0
7,10

7,5
Limon argileux brun •"-io !
et bleu petits rognons
8,0

8,5 • 8,50. -f Petites graves

argileusesjaune e_a
9,0 9,00 0/60 mm op e

9,5 Argile compacte

jaune et beige
10,0
10,20

10,5-

D a t e des mesures et 1. 24/07/79

c o t e d e la n a p p e : 2. 27/07/79
Date du zéro : 2 0 / 0 7 / 7 9 3. 31/07/79

E X E M P L E D E DIFFICULTÉS D'ÉTUDES Cette dernière couche correspond soit aux formations

D U SOL SUPPORT du Pontien, soit à celles du Donaù. Dans le cas du
Pontien, les dépôts sont constitués de sols argileux
comportant de gros éléments quelquefois décomposés.
Dans le cas du Donaû, qui est une formation ancienne
Le projet de barrage de C comporte une digue de
du Quaternaire, la granularité peut être localement
400 m de long et de 13 m de haut au maximum. L a
beaucoup plus grossière, et c'est généralement dans
surface du plan d'eau prévue est de 13 hectares.
ces niveaux que l'on rencontre de fortes venues d'eau.
Les données géologiques du site sont a priori défavora- A titre indicatif, pour illustrer le comportement hydro-
bles. E n effet, le terrain naturel est le siège de circula- géologique de ces sols, ce sont des matériaux que
l'on a pu drainer facilement au moyen de drains
tions d'eau importantes, généralement proches de la
subhorizontaux dans le cas de problèmes de stabilité
surface mais pouvant aussi intéresser des horizons
de pentes.
plus profonds. L a coupe géotechnique est la suivante :

— 1 à 3 m de tourbe, Une reconnaissance approfondie du sol support était

donc nécessaire. A cet effet, les opérations suivantes
— 1 à 2 m de graves argileuses (avec des éléments
ont été réalisées au droit de la future digue :
de dimensions supérieures à 400 mm),
— à partir de 3 à 5 m une forte épaisseur d'alternance — sept sondages à la pelle mécanique à 7 m de pro-
de graves argileuses et de sables argileux. fondeur. Outre l'identification des sols il a été possi-

19
ble, en cours de foration, d'observer d'éventuelles beaucoup trop fortes (31,5 % en moyenne, alors que
venues d'eau importantes ; la traficabilité n'est possible qu'en dessous de 26 %).
— quatre sondages carottés, dans lesquels des essais Seul un traitement entre 2 et 4 % de chaux permettrait
de type Lefranc ont permis de mesurer des coefficients la réalisation du chantier.
de perméabilité compris entre 1 0 et 10~ m/s;
- 8 7

Comme le montrent les résultats ci-après, le traite-

— un essai de pompage classique à partir d'un puits
ment à la chaux améliore très sensiblement les caracté-
de 10 m de profondeur et 1,5 m de diamètre, les
ristiques mécaniques des matériaux compactés (à
mesures de rabattement étant faites dans sept piézo-
95 % de l'OPN) :
mètres disposés sur le profil en travers correspondant
à la plus grande hauteur du barrage (fig. 6). Les - sol non traité c' = 2 0 k P a , 9 ' = 29°
piézomètres n'étaient crépines qu'en profondeur, au- - sol traité ( c ' = 160kPa, cp' = 33° (pic)
dessous de la couche perméable de surface. U n rabatte-
à2%deCaO ( c ' = 20kPa, (p' = 33°
ment moyen de 3 à 4 m a été maintenu pendant un
Pour le matériau traité, l'écart important entre la
mois au moyen d'un simple pompage quotidien. O n a
résistance de pic et le palier de grande déformation
pu ainsi mesurer k = 2.10" m/s et un rayon d'action
7

met bien en évidence la fragilité du matériau traité.

R = 10 m.

Par ailleurs, on a vérifié que le traitement à la chaux

Cet essai, un peu exceptionnel pour un petit barrage,
a pour effet d'augmenter sensiblement la perméabilité
a permis sur une assez grande échelle d'écarter le
du matériau :
risque de l'existence d'un horizon perméable continu
qui aurait pu être dangereux pour la stabilité de
matériau non traité k = 2.10" 8
m/s
l'ouvrage à cause des phénomènes d'érosion interne
25 j de cure :
et de risques de soulèvement en pied de barrage. Ce matériau traité fc = l , 3 . 1 0 m/s - 7

problème principal de perméabilité du sol support à 2 % de C a O 75 j de cure :

étant résolu, on pouvait donc concevoir une structure
fe = l , 2 . 10~ m/s 7

de digue relativement simple (fig. 7) :

27 j de cure :
— la purge des matériaux de surface, tourbes et gra- matériau traité k = 2.10" 7
m/s
ves argileuses perméables, est évidemment nécessaire, à 4 % de C a O 75 j de cure :
— la bêche, destinée à améliorer l'étanchéité au /c = 6 , 8 . 1 0 m/s
- 7

niveau des terrains de surface, permet aussi de vérifier

l'absence d'horizons drainants superficiels non détec-
tés lors de la reconnaissance,
PROBLÈMES P A R T I C U L I E R S POSÉS
— le tapis drainant en pied et la protection du pare-
P A R U N B A R R A G E ÉCRÊTEUR D E C R U E S
ment amont correspondent à des dispositions classi- (barrage de N)
ques pour une digue homogène,
— les pentes de talus proposées correspondent à un
coefficient de sécurité de 1,5 en vidange rapide et de U n barrage écrêteur de crues a un fonctionnement
1,6 pour le talus aval. hydraulique spécifique pour chaque site. Normale-
ment, cet aspect ne concerne pas le géotechnicien;
Le problème principal reste celui de la mise en œuvre toutefois une bonne définition des principes permettra
des matériaux qui en place ont des teneurs en eau généralement de mieux adapter le projet au site étudié.

(m) LÉGENDE
2
0 L j Piézomètre tube ouvert
10 20 (m)
Puits de 1,5 m de diamètre
-TN après décapage
TN -Niveau stabilisé de la nappe
Tourbes alluvions récentes après rabattement dans le puits

Argile limoneuse et Argile

graveleuse ocre jaune
limoneuse
Argile sableuse ou Fig. 6 — Profil en travers de la digue et
limoneuse ocre jaune implantation du puits et des piézomètres
avec gros éléments pour l'essai de pompage.

20
 À

Le fonctionnement idéal consiste à stocker le maxi-

mum d'eau tout en évacuant à l'aval un débit régulé
à la valeur maximale admissible par le chenal. L a
combinaison simple barrage en terre-vidange de fond-
évacuateur de crues est très éloignée de ce fonctionne-
ment idéal :

— si la vidange de fond ne comporte pas une vanne

régulée, le débit d'exhaure croît avec la charge, ce
qui ne permet pas d'utiliser au mieux les capacités
d'évacuation à l'aval et réduit d'autant le stockage
utile ;
— un évacuateur de crues largement dimensionné est
indispensable afin d'exclure tout risque de submersion
de la digue en terre. Lorsque l'évacuateur est constitué
d'un simple seuil déversant, ses dimensions et en parti-
culier la cote de la crête du seuil doivent être telles
que la crue maximale envisagée puisse s'écouler; mais,
dans ces conditions, la capacité totale du bassin ne
Fig. 8 — Profil en travers géotechnique au droit du barrage.
sera utilisée que pour la crue maximale, ce qui est
irrationnel. Dans ce cas aussi, une vanne régulée per-
met d'imaginer un fonctionnement plus satisfaisant
de l'évacuateur de crues.
Voie ferrée M u r
bajoyer
t Niveau max.
Selon l'importance du barrage et en particulier sa —k
s.
longueur, l'aménagement comprendra une digue et un
" * / / / , Remblai / / / /
ouvrage de décharge, ou bien un simple barrage Rembla^^^
Emplacement
mobile qui assurera la triple fonction barrage- ^.de Ja vanne J

décharge-régulation des débits et des hauteurs. \ </////>/

\ Ancien Ht; 5(\Radierv\X
Dans l'exemple du barrage écrêteur de crues de N
l'étroitesse du chenal à l'emplacement envisagé pour \ ^
. 10 m
V /
le barrage se prêtait bien à la solution « barrage \ Écran étanche /
mobile » dont quelques-uns des problèmes spécifiques
Q . A longueur \ ^-T-^-T-"/ ~ ~ " ~
vont être maintenant examinés.
2 m T T T T

1 hauteur ^^X'
Le ruisseau concerné a un régime de type méditerra- 2 m X^-' " >
néen. Il est le plus souvent à sec, mais la montée des
eaux peut être extrêmement brutale et atteindre pour
Fig. 9 — Profil en long schématique d'un barrage mobile.
la crue centennale un débit de 165 m /s. L'aménage-
3

ment du chenal à l'aval permet d'évacuer un débit

maximal de 21 m /s. Une autre contrainte dans le
3

projet est liée à la présence d'une ligne importante de

contournement et une érosion interne qui serait amor-
la S N C F , ce qui impose de façon stricte un niveau
cée à partir des berges de la rivière à l'aval du barrage.
d'eau à ne pas dépasser.
L'écoulement par contournement a été étudié de
Dans l'étude d'un tel ouvrage, on retrouve un certain manière simplifiée par analogie électrique en suppo-
nombre de problèmes spécifiques des digues qui sant au niveau du fond du lit actuel (fig. 9) un écoule-
demandent le même type d'investigation, en plus de ment horizontal, dont les conditions aux limites
celle plus adaptée à l'étude d'une fondation. seraient une charge de 6 m d'eau à l'amont et zéro à
l'aval. O n retrouve dans cette schématisation le pro-
L a figure 8 donne le profil géotechnique des sols au blème classique d'un écoulement autour d'un para-
droit de l'ouvrage. Il s'agit généralement de sables fouille mais avec un potentiel égal à la pression,
fins graveleux ou argileux. A priori, la perméabilité puisque la cote est constante.
du sol support est sans importance, de forts débits de
fuite étant sans conséquence sur le fonctionnement O n a montré ainsi que la longueur de l'écran latéral
d'un barrage écrêteur de crues. de la figure 9 est théoriquement insuffisante vis-à-vis
du gradient critique d'érosion interne pour ce type de
Mais, dans le cas décrit ici, un matériau sableux matériau. Le problème a été résolu non pas en allon-
pose surtout le problème de l'érosion interne. L'écran geant l'écran, mais en protégeant les berges à l'aval
étanche, tel qu'il est représenté sur la figure 9 n'a avec un masque qui, de toute façon, était nécessaire
donc été proposé que pour éviter ce risque. Cet écran pour assurer la stabilité vis-à-vis des forces hydrodyna-
pourrait être constitué d'un rideau continu en ciment- miques.
bentonite.
L a figure 10 donne un profil en travers schématique
O n notera sur la figure 9 que l'écran étanche déborde au droit du barrage mobile. O n notera la position à
de part et d'autre du barrage. Compte tenu de la l'amont de l'écran d'étanchéité. Cette disposition s'est
nature des matériaux on peut craindre en effet un avérée avantageuse du point de vue de la stabilité de
 NGF

Enfin, la proximité du remblai de la voie ferrée a

20 •
nécessité d'examiner deux problèmes, d'une part véri-
-*- Vanne fier l'insensibilité des matériaux aux variations de
teneur en eau et d'autre part dimensionner un masque
Masque étanche Events garantissant la stabilité en vidange rapide.
15
Radier béton

1 Béton poreux 'Filtre PROBLÈMES D E S B A R R A G E S P O U V A N T

Ecran étanche F O N C T I O N N E R DANS L E S D E U X SENS
10

Dans certains cas particuliers, un barrage peut être

; Calcaire
conçu pour fonctionner dans les deux sens. Par exem-
ple c'est le cas des digues de rétablissement des com-
Fig. 10 — Profil en travers schématique d u barrage mobile. munications de l'aménagement qui ont pour rôle de
maintenir un plan d'eau pour des raisons touristiques
lorsque la retenue principale est vide, mais qui peu-
Niveau d'eau maximal 21 N G F vent aussi exceptionnellement avoir un rôle inverse.

L a figure 12 donne un autre exemple d'aménagement

où le barrage, exceptionnellement, peut fonctionner
en sens inverse. Le projet comprend une voie sur
berge qui localement sert de barrage pour un petit
plan d'eau touristique d'un hectare et demi environ.
Le terrain naturel est constitué d'alluvions graveleuses
très perméables (k ~ 10~ m/s) reposant sur le substra-
2

tum compact assez imperméable dans la zone étudiée

(k = 2 à 6. 1 0 " m/s d'après des essais Lefranc dans
6

cinq forages).
Fig. 11 — Profil en travers type de la digue en terre.

Voie sur berge

l'ouvrage au glissement. E n effet, dans ce cas la pous-
Longueur 0 5 10m .61.70
sée hydrostatique ne s'exerce que sur la hauteur de la
Remblai compacté
vanne. Si, au contraire, l'écran étanche est placé sous
Masque étanche
le radier, la poussée hydrostatique s'exerce en plus
sur toute l'épaisseur du radier ce qui, dans le cas
étudié, multiplierait les efforts horizontaux par 1,8.

Les remblais situés de part et d'autre du barrage

mobile ont été étudiés de la même manière que pour
une digue en terre. L a figure 11 donne un profil en
travers de principe. Profil géotechnique d ' u n aménagement o u le barrage peut
fonctionner dans les deux sens.

Côté amont, un choix économique doit être effectué

entre un masque d'épaisseur importante ou une pente
de talus faible, qui a pour conséquence d'augmenter L'exploitation des graves jusqu'au niveau du substra-
la longueur des murs en retour. tum permet donc la réalisation du plan d'eau, à condi-
tion de réaliser un masque étanche entre le bassin et
Côté aval, compte tenu de l'étroitesse du chenal, le la rivière. U n film plastique étanche ou bien un mas-
niveau d'eau peut être important et varier rapidement que en argile compactée permettent de résoudre sim-
à la décrue, il est donc paradoxalement nécessaire plement le problème.
d'étudier la stabilité en « vidange rapide ». Le masque
en enrochements est nécessaire pour éviter les risques L a rivière bordant l'aménagement a un régime médi-
d'affouillements en pied de barrage. terranéen avec des crues extrêmement subites. O n

1 2m Fig. 13 — Protection d u film plastique étanche

contre les sous-pressions en cas de crue de la
rivière, le plan d'eau étant vide.

Sable
plastiqu Matériau < ì recharge
, (réutilisation i

Substratum

22
 i

peut donc se trouver dans la situation où le plan ramené à un problème plan dans lequel on suppose
d'eau est vide au moment d'une crue exceptionnelle l'épaisseur de la nappe constante, ce qui revient à
de la rivière. Pour éviter les désordres dus aux sous- négliger la courbure des lignes de courant selon la
pressions à l'arrière de l'écran étanche, on peut imagi- troisième dimension. O n peut ainsi simuler l'écoule-
ner soit un système de remplissage automatique de la ment par analogie électrique.
retenue pendant la montée des eaux, soit la mise en
place d'une surcharge qui assure la stabilité de l'écran. L a figure 14 donne un exemple de réseau d'écoule-
L a figure 13 donne le schéma d'une protection au ment qui permet d'estimer les débits de fuite compte
moyen d'une surcharge. O n peut vérifier la stabilité tenu des hypothèses simplificatrices faites.
du coin A B C en écrivant l'équilibre des forces; on
peut aussi utiliser des programmes de calcul en rup- Le débit total est donné par la relation :
ture non circulaire tel que P E T A L .

Enfin, il est nécessaire de préciser la longueur de En

la zone à étancher en fonction des débits de fuite
admissibles. Le problème tridimensionnel a été k : coefficient de perméabilité (10~ m/s),
2

H : épaisseur de la nappe (supposée constante,

H = 4,2 m),
N : nombre de tubes courant (ici N= 18),
E : nombre d'intervalles entre équipotentielles (ici
£ = 10),
m rapport géométrique constant de la largeur d'un
~ • tube de courant divisé par la distance entre deux
équipotentielles (en général on prend m/n = l).

O n calcule ainsi g ~ 3 2 0 1/s. Il ne s'agit bien sûr que

d'un ordre de grandeur qui dans le cas du projet
étudié est inadmissible. L a même méthode a permis
d'ajuster la longueur de l'écran à des débits de fuite
(approximative) beaucoup plus faibles.

23
 Digues
^^et barrages en terre
'-fWf//, de faible hauteur

Les bassins de Maurepas-Courance et des Roussières

Présentation

Jean-Pierre BOGAERT JeanNALLET*

I n g é n i e u r principal au syndical communautaire Directeur technique de l'Etablissement public
de l'aménagement de l'agglomération de la ville nouvelle de la ville nouvelle de Saint-Quentin-en-Yvelines
de Saint-Quentin-en-Yvelines (SCAAN)

La maîtrise du ruissellement des eaux de pluie demeure une contrainte pour le

développement de la ville nouvelle de Saint-Quentin-en-Yvelines. En effet, partant
d'un faible ruissellement en milieu rural, il faut contrôler et évacuer par des ouvrages
de régulation le fort ruissellement induit par l'urbanisation.

Les infrastructures d'assainissement étaient des plus succinctes car le cadre naturel
suffisait à lui seul à assurer les besoins. Tout était donc à faire pour permettre
l'installation et le développement de zones d'habitation, de commerce, d'activités,
de sports et de loisirs.

Le site se prête assez mal à l'élimination des eaux pluviales, puisque son éloignement
d'un cours d'eau à débit important (la Seine) représente une distance à vol d'oiseau
de 30 km.

D'autre part, sa position sur le plateau de Trappes ne favorise pas non plus de
grands écoulements. En effet, malgré la présence de nombreux thalwegs à son
pourtour, les possibilités sont faibles car ces thalwegs ne donnent naissance qu'à de
minuscules ruisseaux incapables d'absorber les débits importants provenant de
l'urbanisation.

Cette constatation est d'une importance capitale car elle entraîne des conséquences
lourdes tant sur le plan physique que sur le plan financier. Cela a donc conduit à
construire des réservoirs naturels destinés à écrêter les pointes pour ne restituer
ensuite qu'un faible débit compatible avec la capacité d'accueil des exutoires. Ainsi,
à terme, une trentaine de bassins seront construits afin de stocker et réguler environ
2 000 000 de mètres cubes d'eau. Ces réservoirs sont situés directement à l'aval des
zones urbanisées, soit sur le plateau, soit dans les vallées. Ils sont calculés pour
des orages de fréquence cinquantenaire ou centenaire et peuvent être soit des
« bassins secs » soit des « bassins en eau ».

Certains bassins ont été réalisés en « bassins secs » mais de plus en plus les
nouvelles retenues sont du type « bassins en eau ». Cela s'explique aisément car à
la fonction hydraulique viennent se superposer des fonctions de détente et de loisirs
dues à l'utilisation des plans d'eau.

* Actuellement au BCEOM.

5
Les sites dans lesquels sont implantés les ouvrages de Maurepas-Courance et des
Roussières nous a permis, dès le début, de retenir des bassins en eau.

Conçues par le même bureau d'étude, les deux digues sont de type homogène, en
limon, avec un dispositif de drainage identique : tapis, cheminée drainante et puits
de décharge.

Digues importantes, puisqu'elles se situent à la limite des grands barrages; elles

ont été équipées d'un système d'auscultation et de surveillance développé car leur
défaillance mettrait en jeu la sécurité publique.

Le Laboratoire régional de TOuest parisien a reçu pour mission de prévoir et de

suivre révolution du corps de digue. Cette mission a débuté avec la construction et
a été étendue, par la suite, à F exploitation.

Les articles de MM. Baldit, Delfaut, Jardin du LROP et MM. Josseaume et

Mahdavi du LCPC exposent les méthodes employées pour mener à bien ces études
et les résultats du suivi.

Ces constatations ont été complétées par des études théoriques, notamment par
Vanalogie électrique et par des méthodes numériques pour la détermination du
réseau d'écoulement à travers le corps de digue et de son support afin d'apprécier
ïanisotropie des sols.
 i,

Digues
et barrages en terre
7/7///////////////miUi//, de faible hauteur

Constatations sur la digue

de Maurepas-Courance
Abel DELFAUT
Assistant

Jean JARDIN
Ingénieur

Robert BALDIT
Technicien supérieur
Laboratoire régional de l'Ouest parisien

RÉSUMÉ

Le bassin de retenue de Maurepas-Courance

s'inscrit dans un programme d'aménagement
de la ville nouvelle de Saint-Quentin-en-Y vélines ;
il s'agit d'une retenue d'orages faisant office de
bassin tampon pour les eaux pluviales du
plateau dont l'urbanisation a accru les débits de L a retenue de Maurepas-Courance s'inscrit dans le plan
pointe du ru existant en fond de vallée. d'aménagement du réseau hydrographique du plateau de
La mission confiée au LROP a été orientée Trappes. L'urbanisation de la région, par la création de
plus spécialement vers le suivi du comporte- la ville nouvelle de Saint-Quentin-en-Yvelines, a entraîné
ment du corps de digue durant la phase de
construction et ensuite après la mise en service. progressivement l'imperméabilisation d'importantes surfa-
Pour cela, il a été nécessaire d'équiper deux ces conduisant ainsi à l'évacuation de quantités croissantes
profils en travers de la digue en cellules de d'eau de ruissellement vers les rus existants; or ceux-ci
pression interstitielle, de tassomètres et de
tubes piézométriques. Des mesures optiques sur n'ayant pas une capacité suffisante pour écrêter les pointes,
des repères bétonnés ont été effectuées afin il s'est avéré indispensable de créer des retenues d'eau le
d'apprécier les déplacements du corps de digue. long de leur tracé. Ces bassins tampons permettent l'accumu-
L'article rend compte des principaux résultats lation des eaux pendant les fortes pluies puis la restitution
obtenus lors de la construction de la digue, de lente dans le cours d'eau.
la mise en eau et de l'exploitation du bassin.
Pour essayer de comprendre le fonctionnement O n dénombre une quinzaine de retenues construites sur le
hydraulique réel de l'ouvrage, il a paru intéres- plateau de Trappes. Elles sont toutes de moyenne et petite
sant de rapprocher les constatations des
premières analyses théoriques faites dans le importance, les hauteurs des digues varient entre 3 et 15 m ;
cadre d'une recherche sur l'anisotropie des celles de Maurepas et du centre ville (digue C) sont les plus
limons utilisés en corps de digue. On verra que importantes par leur capacité et par la hauteur de la retenue.
l'on n'est pas parvenu à des conclusions nettes,
et combien il est difficile de tenir compte de
la complexité des phénomènes réels dans les Le bassin de Maurepas est caractérisé par une hauteur de
études théoriques à fortiori au stade d'un retenue normale en pied de digue (vidange de fond) de 10 m
projet.
et celle des plus hautes eaux (niveau du déversoir) de 14 m ;
MOTS CLÉS : Barrage en terre • Pression la capacité de stockage aux deux niveaux est respectivement
interstitielle - Perméabilité • Alluvions - Conso- de 213000 m et de 393000 m sur une superficie de plan
3 3

lidation - Tassement - Mesure • Déplacement - d'eau, à la retenue normale, de 4,9 ha.

Stabilité - Piézomètre - Auscultation - Sous sol -
Modèle analogique • Écoulement (fluide) •
Jsotropie - Limon - Drainage - /Digue - Bassin L a hauteur maximale de la digue est de 15,5 m (revanche
d'orage. de 1 m) avec une largeur en tête de 4 m et une longueur de
245 m.

B u l l , l i a i s o n L a b o P. et C h . - 131 • m a i - j u i n 1 9 8 4 - R é f . 2 7 9 4
CADRE D'INTERVENTION L'article a pour objet de décrire l'ouvrage et son
D U LABORATOIRE comportement pendant et après l'exécution, en s'atta-
chant particulièrement à comparer ce dernier avec les
prévisions.
Les missions confiées au Laboratoire régional de
l'Ouest parisien visaient :

— à définir le site géologique de l'ouvrage et plus

particulièrement celui de la future digue : il s'agissait
SITE D E L ' O U V R A G E
d'une étude préliminaire s'appuyant sur des investiga-
tions en nombre limité; cette étude géologique a été
toutefois complétée par quelques essais de laboratoire
L'ouvrage est implanté dans la vallée de la Courance,
(cisaillement et compressibilité);
thalweg orienté S E - N O qui entaille le plateau du
— à mettre en place et à suivre des appareils de Hurepoix au nord du village de Maurepas (fig. 1 et 2).
mesure : sondes de pression interstitielle, tassomètres, L a vallée dans l'axe de l'ouvrage accuse une asymétrie
piézomètres installés dans le corps de digue et dans assez marquée : sur le versant, du côté droit, exposé
le sol de fondation. au sud-ouest, la pente moyenne est de l'ordre de 20 %
alors que sur le versant du côté gauche, exposé au
Les études géotechniques complémentaires nécessaires nord-est, la pente n'est que de 12 % ; celle du fond de
à la conception de la digue et des ouvrages en béton thalweg peut être estimée à 2,5 % en moyenne.
armé ont été assurées par le Bureau d'ingénieurs-
conseils Coyne et Bellier. Cet organisme a d'ailleurs, Le versant le plus raide est boisé (bois de Maurepas)
par la suite, servi de conseil à l'entreprise adjudicataire tandis que le versant à pente plus douce est à vocation
(Viafrance) et au maître d'oeuvre ( D D A des Yvelines), agricole. A l'origine, il existait un étang dont le niveau
le laboratoire ayant obtenu parallèlement une presta- se situait à une cote voisine de l'émergence de la
tion directe avec le maître d'ouvrage délégué ( E P A nappe phréatique (nappe des sables de Fontainebleau)
de Saint-Quentin-en-Yvélines) et par la suite avec le et qui avait, jadis, desservi un moulin dont on pouvait
Syndicat Intercommunal ( S C A A N ) . encore apercevoir les ruines.

Fig. 1 —
Plan de situation.

Fig. 2 — Site de l'ouvrage.

8
L a présence de cette retenue a amplifié considérable- Le fond du thalweg est comblé par les alluvions du
ment le caractère marécageux de cette zone en ru. Leur structure est complexe, comme en témoigne
favorisant des dépôts de nature variée : vases, débris le profil de la vallée sur la figure 3.
végétaux...
Les premiers alluvionnements ont été déposés sur les
argiles à cérithes sous forme de lentilles provenant de
matériaux érodés à l'amont : on distingue dans la
SITE GÉOLOGIQUE
vallée, mélangés pèle-mêle, des sables fins plus ou
moins propres et organiques et des argiles sableuses
L a reconnaissance du site, si l'on regroupe toutes les (anciennes vases), le tout contenant des niveaux plus
phases d'intervention, y compris celle de l'entreprise, grossiers de débris siliceux (meulière) formant parfois
a nécessité l'exécution des sondages et essais suivants : de véritables lits de graviers inclus dans ces alluvions
fines. Ces derniers alluvionnements ont amené des
— dix sondages carottés profonds de 20 à 30 m dont niveaux tourbeux et des vases molles dans lesquelles
sept équipés de piézomètres; se trouvaient emprisonnés des débris végétaux.

— sept sondages à la tarière tous équipés de piézo-

mètres; Hydrogéologie
— douze sondages destructifs dont trois ont été
accompagnés d'essais pressiométriques; Le ru de la Courance sert d'exutoire, dans la partie
du projet, à une ligne de sources qui apparaît à la
— vingt-six sondages de pénétration dynamique; cote 111 N G F sur le versant droit, et à la cote
113 N G F sur le versant gauche. Cet écart du niveau
— quinze sondages au pressio-perméamètre;
provient certainement de la pente différente des ver-
— un essai de pompage en rive droite. sants, mais aussi de la présence d'horizons argileux
au sein des sables de Fontainebleau qui servent locale-
ment de guide support à la nappe aquifère la faisant
Lithostratigraphie déboucher à des endroits légèrement différents selon
ses ondulations.
O n trouve ici la série stratigraphique régionale com-
plète, à savoir : au sommet des flancs de la vallée, Les gradients hydrauliques caractérisant le rabatte-
les calcaires et meulière de Beauce qui déterminent ment de la nappe dans le fond du thalweg sont de
l'entablement du plateau à la cote voisine de l'ordre de 80/1000 sur les deux versants (en toute
170 N G F , puis les formations des sables de Fontaine- rigueur le gradient rive droite est légèrement plus fort
bleau (60 m d'épaisseur) dont l'affleurement est recou- que celui de la rive gauche, de 5/1000) et de 15/1000
vert d'un tapis d'éboulis argilo-sableux. Les sables en fond de vallée.
reposent sur un support marneux composé des argiles
à cérithes, des marnes à huîtres, d'une séquence cal- L a principale alimentation du ru provient du rabatte-
caire que l'on peut attribuer au calcaire de Brie, des ment de la nappe des sables de Fontainebleau, ce qui
marnes vertes et des marnes supragypseuses. maintient ainsi un débit permanent; le reste peut être

Prélèvements d'échantillons intacts

25m
Echelle horizontale
Niveau piézométrique

Fig. 3 — Profil géologique en travers de la vallée suivant l'axe longitudinal de la digue.

9
attribué aux eaux de drainage du plateau, aux eaux Le tableau I récapitule par familles de sol les perméabi-
de ruissellement le long des pentes et à celles de lités horizontales mesurées au pressio-perméamètre
l'artésianisme de la nappe des calcaires de Brie. Ce dans les sols rencontrés en fond de vallée.
dernier point mérite une attention toute particulière,
NOTA. — Un essai en laboratoire sur éprouvette œdométri-
car il conditionne pour une part le régime hydraulique
que (perméabilité verticale) d'un échantillon du complexe
du site de l'ouvrage : les piézomètres profonds (20 m
sablo-argileux (faciès d'argile sableuse) a donné une valeur
et plus) indiquent un niveau d'eau à la cote
de 0,7 x 10" m/s à la pression normale des terres (cote
7

109,20 N G F , tandis que les piézomètres courts (10 m) 108 NGF).

donnent un niveau proche du T N vers 108 N G F . Dans la couche alluvionnaire du complexe sablo-argileux,
on distingue des familles : les sables argileux de coefficient
Deux types d'essais pour la mesure de la perméabilité de perméabilité moyen fc = 1,6 x 10" m/s pour trente et
A
6

en place ont été réalisés : une valeurs et les passages graveleux ou sableux de coeffi-
cient de perméabilité moyen k =\,l x 10" m/s pour cinq
h
5

— des essais ponctuels au pressio-perméamètre valeurs. Le rapport du nombre de valeurs de chaque milieu
Ménard effectués dans quinze sondages ont permis reflète leur importance respective dans la couche.
de différencier les couches du point de vue de leur
perméabilité horizontale (fig. 4) ; A l'examen du tableau I, on peut émettre les conclu-
sions suivantes :
— un essai de pompage en rive droite a permis de
déterminer globalement la perméabilité des sables — les coefficients de perméabilité moyens des sables
roux légèrement argileux qui affleurent à l'émergence de Fontainebleau éboulés et en place sont dans un
de la nappe et qui ont servi d'appui à la digue. Cet rapport inférieur à 10;
essai a été réalisé par le bureau d'étude Burgeap.
— les coefficients de perméabilité moyens du com-
plexe sablo-argileux et des argiles à cérithes sont dans
un rapport voisin de 100. L'horizon des argiles à
io- 9
2 3 4 5 7
10- 8

2 3 4 5 7
10- 7

2 3 4 5 7
10" 6
,10-5k„
2 3 4 5 7 cérithes constitue hydrauliquement le support
Tourbes et « imperméable » ;
vases PP6 PP7
t — la grande dispersion des coefficients de perméabi-
lité du complexe sablo-argileux qui sert de fondation
Complexe * à la digue est significative de l'hétérogénéité du milieu.
argilo-sableux
L a présence de passages graveleux peut déterminer
-I des circulations préférentielles d'écoulements;

Argiles
a cénthes
10

15
m — les argiles à cérithes, base du sable de Fontaine-
bleau, sont deux cent cinquante fois plus imperméa-
bles que le sable de Fontainebleau argileux en place.

Comportement mécanique des terrains

Fig. 4 — Profils de perméabilité d u sol de f o n d a t i o n
d u corps de digue.
Les essais de comportement sur le site ont été peu
nombreux; ils concernent la couche alluvionnaire. Ils
L a nappe a été testée sur une hauteur de 5,50 m au ont été effectués à l'occasion de l'étude préliminaire,
moyen d'un forage équipé en puits : sept piézomètres dont l'objectif n'était pas de justifier la conception de
d'observations ont été mis en place à distance variable l'ouvrage en terre, mais d'entrevoir les dispositions
du point de pompage. à prendre pour assurer globalement la stabilité de
l'ouvrage.
L'interprétation des essais par la méthode de Theis a
conduit à un coefficient de perméabilité horizontal Les résultats des essais de cisaillement sont regroupés
compris entre 6,5 et 8 x 10" m/s; la valeur moyenne 6 dans le tableau II et les caractéristiques de compressi-
est donc voisine de 7 x 10" m/s. 6 bilité dans le tableau III.

TABLEAU I. — P e r m é a b i l i t é s h o r i z o n t a l e s ( p r e s s i o - p e r m é a m è t r e ) .

Nombre h h t
h maxi
Nature du sol h mini
h maxi
d'essais n

(m/s) h moyen

(m/s) (m/s) h
Sable fin en place A 4,5 x I O h mini
7
25 6xl0" 1,8 x IO" 5 -6 n

30
Argiles sableuses en place B 28 2,8x10" 7
IO" 5
2,1 xlO" 6
35
Sables fins éboulés C 12 9x 10" 7
2 x IO -5
6,1 x I O -6
22
Complexe sablo-argileux D 36 7xl0" 9
3,5 x IO" 5
1,8 x I O -6
5 000
Argiles à cérithes E 21 2xl(T 9
5xKT 8
2,4 x IO"
8
25

10
 TABLEAU II. Cisaillement Cisaillement
Cisaillement
w à la boîte au scissomètre
Paramètres Nature du sol Prélèvement à la boîte
(%) UU de chantier
CD
de cisaillement. (*) (kPa)

Tourbe CPS 1
limoneuse -1,00 m 105 10
CPS 2
-0,50 m 55 13
Tourbe CPS 1
-1,80 m 312 10

Sable CPS 1
vasard -2,00 m 30 = 34°
<f>HH
c„„ =10kPa
-3,00 m 30 = 31°
<f>uu
c = 10kPa
uu
Argile CPS 1
sableuse -4,00 m 19 = 30°
Complexe = 20kPa
sablo-argileux CPS 2
-1,50 m 20 = 32°
<f>uu
= 20kPa
-2,00 m 19 <f>HU = 34°
= 20kPa
Argile CPS 2
sableuse
avec graviers -3,00 20 <p' = 36° c' = 0

Argiles Sable 29 <p' = 35° c' = 0

sableuses argileux

(*) Essais UU effectués selon le mode opératoire. La détermination d'un angle de frottement indique un drainage au moment de l'essai
en raison de la nature sableuse des sols.

T A B L E A U III. — P a r a m è t r e s d e c o m p r e s s i b i l i t é .

Nature du sol Prélèvements

(kN/m ) 3
(kPa) c c

Argiles sableuses B C, 6,40 à 7,20 m 14,8 45 0,106

Complexe sablo-argileux D CPS 1 3,50 à 4,50 m 16,8 28 0,041
CPS 2 2,50 à 3,50 m 17,9 21 0,21
C 2 5,00 à 5,50 m 16,5 40 0,018 0,09
C 3 6,00 à 6,50 m 15,6 80 0,020 0,146
Argiles à cérithes E C 13,10 à 14 m
t 16,1 200 0,03
Marnes vertes C 22 m
t 15,1 130 0,145
25 m 15,3 160 0,140

C : indice de compression dans l'intervalle de contraintes a' à rs .

s v0 p

C : indice de compression au-delà de a .

c p

CONCEPTION D E L A DIGUE fondation et les matériaux filtrants, un limon ( 1) dont

l'indice de plasticité sera supérieur à 1 5 ( / > 15); p

Transversalement, la digue est de type homogène avec

tapis drainant aval et drain cheminée. Cependant, — au centre et à l'amont de la cheminée drainante,
l'insuffisance de limon argileux dans l'emprise du pro- un limon (2) de plasticité légèrement inférieure au
jet a nécessité un approfondissement de l'emprunt précédent (8 < I < 15); p

jusqu'au sable; cela explique la présentation sur les

coupes en travers de la digue de zones différenciées — à l'aval de la cheminée drainante, un sable
en fonction de la nature des sols (fig. 5). Cette évolu- limoneux (3) de faible plasticité (/ <8). p

tion de la nature du matériau étant connue avant le

début des travaux, la structure de la digue a été définie E n fait, les deux premières catégories de limons ont
de la manière suivante : été très difficiles à distinguer et l'on peut considérer
qu'il s'agit pratiquement du même sol, par rapproche-
— à l'amont, au contact du plan d'eau, et à la base ment des références proctor, qui appartiennent à une
du corps de digue de façon à recouvrir les sols de même famille de sol.

11
 NGF

120

Protection rocheuse
4
Niveau avant mise en eau

Marnes à huîtres

Calcaire de Brie

Marnes vertes

Fig. 5 — Structure d u corps de digue.

CAILLOUX GRAVIERS GROS S A B L E SABLE FIN

100
V)
-0)
90

E
80
O
tí)
CO
70 Pl.
c/>
•E 6 0
to

CO 50
Cl)
"D
tu A0

ro3 0
O)
w
C
tu
u 20
Po

1 0

• m m 200 loo 50 20 10 1 0,5 0,2 0,1 50>j 20>J

• nin 20 u 10)J 5 y 2» 1M 0,5 p 0,2p

Diamètres équivalents (Um) Diamètres équivalents (Um)

Fig. 6 — Courbes granulométriques des différents types de sols.

D i : faciès graveleux d u complexe alluvionnaire. 0 2 : faciès sableux d u complexe alluvionnaire.

Le corps de digue est caractérisé par une tranchée FIN

parafouille en pied de digue amont et par un drain

aval en matériaux filtrants choisis de manière à respec-
ter les conditions classiques des filtres vis-à-vis des
matériaux environnants (fig. 6). Ce filtre se retourne
à l'intérieur du corps de digue pour former une chemi-
née drainante inclinée à 1/1, dont la fonction est de
capter des infiltrations éventuelles pour les rabattre
vers le drain de pied. Le filtre horizontal est composé
en son centre d'un gravier propre F sur 0,40 m, 2

protégé sur ses deux faces par un gros sable F d'épais- t

seur 0,20 m ; les conditions de filtre de Terzaghi sont

respectées :
200 100 50 20 10 5 2 1 0.4 0.2 0.1 50« 20
D l 5 m a x
<5 et D l 5 m a x
^5 (fig. 7),
^ 8 5 min d 1 5 m i n
Fig. 7 — Fuseaux granulométriques des filtres.

D et d ouverture de la maille à 15 et 85 % de tamisât

du filtre et du matériau à filtrer.

Les matériaux F et F , reconstitués en centrale,

t 2 Le parement amont est protégé par une recharge en
jouent à la fois les rôles de drains et de filtre. petits enrochements de grès ( F ) disposés sur une 3

nappe de géotextile non tissé dont le rôle essentiel est

Les pentes générales de la digue ont été fixées à 3,23/1 d'éviter la pénétration du limon sous-jacent à travers
à l'amont et 2,32/1 à l'aval. les enrochements ou inversement.

12
Compte tenu de la complexité du site et de l'existence
possible de couches relativement perméables dans le
sol de fondation, il a été jugé prudent d'assurer la
dissipation d'éventuelles sous pressions par un massif Limon 1 Poids volumiques secs ( k N / m ) N b = 5 8
d'enrochements de même calibre que le précédent \ \ \ \\ Limon 2 Poids volumiques secs ( k N / m ) N b = 2 7
3

Sable 3 Poids volumiques secs ( k N / m ) N b = 69

3

ainsi qu'une ligne de puits filtrants en pied aval.

Les puits ont été exécutés par lançage sous pression

d'eau en diamètre de 250 m m et crépines en diamètre
de 150 m m ; l'espace annulaire entre crépine et forage
a été comblé en sable filtrant F , . Les puits sont dis-
tants de 5 m et ils sont ancrés de 2 m dans les marnes
à huîtres.

Mise en œuvre du corps de digue

Le calendrier des différentes phases des travaux est

représenté sur la figure 8. Les matériaux utilisés pour
la construction de la digue proviennent d'un emprunt
situé en rive gauche de la retenue. L a distribution des
déblais s'est faite en fonction de la nature des sols
extraits :

— C A T 1 et 2, ce sont des limons de pente dont

l'indice de plasticité varie entre 15 et 25, mis en place
à l'amont de la cheminée drainante, ils se classent
en A p (LPC) ou A - A ( R T R ) ,
2 3

— C A T 3, ce sont des sables limoneux ou argileux

dont l'indice de plasticité est inférieur à 10, mis en
place à l'aval de la cheminée drainante; ils se classent
en S A - S L (LPC) ou B - B ( R T R ) .
5 2

15,0 15,4 15,8 16,2 16,6 17,0 17,4 17,8 18,2 18,6 19,0 7 d sable kN/rri saole
15,3 15,7 16,1 16,5 16,9 17,3 17,7 18,1 18,5 18,9 19,3 T d limon k N / m 3
limon

J Travaux préparatoires
Enlèvement de la tourbe

0)P Vidange de fond

-I;
Pose des cellules
VTxTT w
OPN L i m o n
'

w
OPN S a b , e 3

2 Mise en œuvre du corps

de digue
ml 3 Mise en place des filtres
0>l et tapis de propreté

MM
^ Mise en place de la cheminée
drainante

Fig 8. —
Calendrier
des travaux.
COP
0>l

5 Ouvrage en maçonnerie
Évacuateur de crue
Tour de régulation

7 8 9 10 11 12 13 15 16 17 18 1 3 20 w (%)
£ Exécution des puits filtrants

Références OPN
Pose des tubes piézométriques

7 Mise en eau yd 18,9 18,3 18,1 -

1 Limon w(%) 15 15,5 15,5 -
yd 17,5 18,7 18,5 18,5
3 Limon
L a mise en œuvre des remblais s'est faite par couches sableux w(%) 15 12 14,5 14,5
de 0,30 m d'épaisseur compactées au moyen d'un
yd 16,4 16,6 16,3 16,5
rouleau à pieds dameurs.
3 Sable w(%) 15,5 14,5 15,5 15,8
L a mise en œuvre de la cheminée drainante a été
conduite par épaulements successifs des matériaux Fig. 9 — Contrôle de la mise en œuvre du corps de digue. Histogrammes
filtrants sur le talus du remblai amont. de poids volumiques et des teneurs en eau.

13
Le contrôle d'exécution a été à la charge de l'en- part la couche d'alluvions fines dont l'épaisseur est
treprise (autocontrôle). L'ensemble des résultats est de l'ordre de 3 mètres.
présenté sur la figure 9 sous forme d'histogrammes
par classe de matériaux. Les cellules implantées dans l'axe du remblai, donc les
plus chargées, enregistrent les tassements maximaux :
92 % des valeurs correspondent aux prescriptions du 30 cm au profil A et 23 cm au profil B.
CCTP:
E n outre, on constate sur ces courbes une vitesse de
Y ^0,95 y
d d o p n . tassement très rapide au cours de la construction de
l'ouvrage. Cela provient du fait que le coefficient de
consolidation c„ des alluvions est relativement élevé
SUIVI D E S M O U V E M E N T S D E L A D I G U E ( 3 . 1 0 ~ m /s) et que les tassements ont lieu surtout
6 2

ET DES PRESSIONS HYDRAULIQUES : dans cette couche; les tassements « secondaires » sont
L E DISPOSITIF D E SUIVI relativement faibles, de l'ordre de quelques centimè-
tres (le quart du tassement total pour les cellules
situées dans l'axe du remblai).
A u niveau des pièces du marché le Bureau d'études
avait prévu un système de surveillance par tubes piézo- L'étude du tassement primaire de la couche d'allu-
métriques en P V C crépines sur une hauteur de 1,50 m vions d'après les courbes oedométriques et la théorie
à la base et bouchonnés au-delà. de Terzaghi concluait à une amplitude de tassement
de 16 cm et un temps de tassement d'un mois après
L a disposition concentrée des tubes piézométriques la construction (fig. 13).
autour et à l'aval de la cheminée drainante avait
pour objectif essentiel la vérification du champ des O n s'attendait donc, d'après la théorie, au tassement
pressions interstitielles dans la zone du corps de digue suivant :
la plus critique du point de vue de la stabilité pour — alluvions fines 16 + 4 (tassement secondaire)
l'ouvrage en service. A la réalisation, ce système d'aus- = 20 cm,
cultation a été maintenu et complété par :
— couches argileuses sous-jacentes = 12 cm.
— des tassomètres classiques L P C de forme cylindri-
que, qui permettent la mesure des déplacements verti- Le tassement n'évolue plus du tout depuis la mise en
caux du sol sous l'effet de la surcharge engendrée par eau en juin 1977 pour deux raisons : la vitesse de
la construction de la digue; tassement s'était considérablement réduite (2 à 3 cm
par an de juin 1976 à juin 1977 pour les cellules centra-
— des capteurs des pressions interstitielles de type les) et le déjaugeage de la digue par remontée de la
P A C II à contrepression pneumatique, destinés à déce- nappe a pratiquement annulé le phénomène par la
ler les surpressions hydrostatiques ponctuelles dans suite.
une couche de terrain. Ces piézomètres ont l'avantage
d'avoir un temps de réponse très rapide (quasi immé- Dans ce cas particulier, malgré les hypothèses simplifi-
diat) et l'inconvénient de ne fonctionner qu'en milieu catrices du calcul de tassement, sol homogène et c v

parfaitement et constamment saturé. constant dans chaque couche, les résultats concordent
assez bien avec les observations.
Le dispositif d'auscultation a été groupé le long de
deux profils en travers, l'un à proximité du déversoir
et l'autre à proximité de l'ancien ru, dans l'axe de Capteurs de pressions interstitielles
l'ancienne vallée, à l'endroit où les hauteurs de rem-
blai sont maximales. Les coupes et plans des Les capteurs de pressions interstitielles ont obéi à
figures 10, 11 et 12 précisent l'implantation des cellu- trois types de sollicitations :
les et tubes piézométriques. 1. les fluctuations des nappes aquifères (Fontaine-
bleau et calcaire de Brie) ;
Pour des raisons pratiques, les cellules ont été implan-
2. la consolidation des couches;
tées après la mise en place des premières couches de
remblai, le fond de fouille n'étant pas accessible à 3. la variation du plan d'eau à l'amont.
pied donc a fortiori pour des machines de forage.
D'ailleurs, le curage de la tourbe et des vases molles
s'est opéré par terrassement à la pelle mécanique le Ainsi, quels que soient le niveau atteint et la nature
long de pistes préalablement construites. des matériaux entourant les sondes, le sens de varia-
tion des surpressions reste sensiblement le même, la
seule différence résidant dans l'importance prise par
Tassements l'une ou l'autre des trois sollicitations citées ci-dessus.
Les sondes situées dans les alluvions fines répondent
mieux à la variation du plan d'eau et au phénomène
Les deux profils d'étude ont été munis de cinq tassomè-
de consolidation, les sondes profondes à celles de
tres chacun, de manière à suivre les déplacements
l'artésianisme de la nappe du Brie.
verticaux. D'après l'analyse des courbes expérimen-
tales de tassement en fonction du temps, on constate
que l'amplitude de tassement du sol est restée modeste Les mesures de référence faites en juin 1975 ont été
vis-à-vis de la surcharge (270 kPa); cela est dû à la influencées par les rabattements de la nappe en fond
faible compressibilité des couches sous-jacentes, dont de fouille dus au terrassement des tourbes et vases et
la plupart se trouve en état de surconsolidation, à de l'action des lignes de pointes filtrantes destinées à

14
© Tube piézométrique

• Cellule de pression interstitielle

• Tassomètre Echelle

Fig. 10 — Vue en plan d u projet. Situation des cellules et tubes piézométriques de surveillance.

NGF 11 1 9 1 2 10
120 1

Niveau avant mjse en eau

Puits filtrant
de d é c h a r g e
«il

Marnes à huitres

Calcaire de Brie

Marnes vertes

(§) Cellules de pression interstitielle

M Tassomètre Fig. 11
Niveau d'eau dans le tube piézométrique
Pression interstitielle au niveau de la cellule
l avant mise en eau 15/5/1977
I 2 après mise en service à la côte 115 NGF (RN)
Fig. 11 et 12 — Profils en travers A et B d u corps de digue. Implantation des
cellules et tubes piézométriques de surveillance.
1 Limon plastique / p > 1 5
2 Limon peu plastique 8</p<15

NGF
3 Sable limoneux / p < 8
120 .
• Protection rocheuse

Puits filtrant
de d é c h a r g e

Marnes à huitres

Calcaire d e Brie

Marnes vertes

Fig. 1 2

15
 i Hauteur de
remblai il „
13m
, PROFIL A
/
*—PROFIL B

Cellule 5

o Cellule 4

* Cellule 3
(Cellule 2 détruite)

A Cellule 5

Cellule 4
o Cellule 3
* Cellule 2

J[J|A|S|O[N|D||J|F|M|A|M|J|J|A[S|O|N|D||J|F|M|A[M|J|J|A|S|O|N|D||J|F|M|A|M1J]J|A|S|O|N[D]

1975 1976 1977 1978

Fig. 13 — Évolution des tassements mesurés.

maintenir latéralement les talus, surtout en rive gau- culier sur les sondes 3, 4, 5 du profil A et 1, 3 du
che où l'affleurement était constitué par des sables profil B. Ces surpressions ont été insignifiantes dans
fins propres qu'il fallait consolider pour construire le la plupart des cellules, sauf à la cellule 3 du profil A
collecteur de vidange de fond de bassin. Dès que les où l'amplitude du pic est de 50 cm d'eau (fig. 14
premières couches de remblai eurent été montées et et 15).
que les systèmes de rabattement furent abandonnés, il
s'installa progressivement devant la digue une retenue Après une stabilisation des pressions interstitielles,
d'eau entre les cotes 105 et 108, ces eaux venant d'une on constate un fléchissement des courbes en avril et
part des terrains tourbeux situés en amont dans le juin 1976, surtout pour celles qui correspondent aux
bassin et d'autre part des versants (émergence de la cellules profondes et qui traduisent certainement l'inci-
nappe). dence des fluctuations saisonnières de la nappe du

Cela explique la remontée bru-

tale des pressions interstitielles Fig. 14 et 15 — Fluctuation des pressions interstitielles pendant les travaux et à la mise en eau
enregistrées par certaines sondes (cellules) selon le profil A o u B.

aux profils A et B; les sondes les

plus en amont et les moins pro-
fondes étant les plus sensibles. 115
L a sonde n° 3 mise à la base
Fluctuation du plan d'eau
du remblai a suivi sans retard 114
à l'amont de la digue
la montée du plan d'eau (fig. 14 113
et 15).
112

Les phénomènes de surpressions 111

dus à la consolidation des cou-

110
ches se sont confondus avec celui
de l'installation de la nappe. 109

Cependant, sur quelques sondes 108

fichées dans les alluvions ou à la
107 Rabattement dû au terrassement
base du limon, on note un pic du fond de fouille
dans la courbe d'évolution des 106 /
•J j j ¡ A l S | O | N | 0 ¡ J ¡ F l M Í A l M l j | j | 0 | s i O | N ) 0 ¡ j | F | M J A J M | J Í J | A Í S | O | N J D ¡ J J F | M | A ¡ M Í J ¡J|A¡S|Q|NJDF
pressions interstitielles, en parti- 1975 1976 1977 1978

16
 Surpression 3 NGF 110,51 pression Int.
interstitielle (cm) o 104,81 cellule
(Limon)

2 NGF 110,31
102,06
(Sables fins)
1
3004 NGF 110,06
102,31
(Sables et graviers)

P4 NGF 109,06
104,76
(Limon)
200_ P5 NGF 108,76
104.21
(Vases sableuses)

NGF 105,90
7 104,00
100. o (Sables et graviers)
NGF 107,10
103,60
(Sables fins)
NGF 106,01
104,66
(Limon)

J|J|À|S|O|N|D|J|F1M|A|M|J|J|O|S|O|N|D|J|F|M|A|M|J|J|A|S|O|N|D|J|F|M|A|M|J|J|A|S|O|N|D1 Fig. 14

1975 1976 1977 1978

Surpression
300. interstitielle (cm) P1 NGF 110,40 pression int.
104,10 cellule
280 (Sables et graviers)

260

240

220

P2 NGF 110,30
200- 102,35
(Argiles à Cérithes)
180

160

140

120

100.

80
NGF 107,85
60 A 99,86
(Argiles à Cérithes)
40

20

0 P4 NGF 109,35
-* 93,80
20 (Marnes à huitres)

40

60

80
J|J|A|S|O|N|D||J|F|M|A|M|J|J|A|S|O|N|D||J|F|M|A|M|J|J|A|S1O|N|D1J[F|M|A|M|J1J|A|S|O[N|D| Fig. 15

1975 1976 1977 1978

17
 NGF (m)
Fig. 16 — Fluctuation des
niveaux piézométriques dans les
tubes au profil A.

V a r i a t i o n d u plan d'eau
à la mise en service

109 .

o " ®

-©Sec

M J J A S O N D J F M A M J J A S O N D

Brie, mais également la proximité des puits aval de P V C 0 60 mm lanternes à la base sur 1,50 m avec
décharge construits à partir de janvier 1977 (sondes 3 bouchon de bentonite annulaire, ont permis de complé-
et 4 du profil B). ter les mesures de pressions interstitielles des sondes
en précisant dans les profils en travers la position de
O n se souviendra, à ce sujet, qu'une sécheresse d'im- la surface piézométrique à l'intérieur du corps de
portance centenaire a sévi de novembre 1975 à digue. Les résultats enregistrés avant et après la mise
septembre 1976 et que ses effets ont pu se faire sentir en service ont été reportés sur les graphes (fig. 17
avec un certain retard dans le temps. et 18) à leur cote N G F et non plus en surpression
comme pour les capteurs de pression interstitielle.
L a mise en eau du bassin à sa cote de retenue normale
décidée à la mi-juin 1977, a fortement agi sur les Les mêmes remarques sont à formuler que pour les
cellules disposées dans les alluvions et à l'amont de piézomètres : les tubes placés en avant de la cheminée
la cheminée drainante; celles-ci ont enregistré, pour ont bien répondu à l'élévation du plan d'eau, les plus
une élévation de plan d'eau de 7 m depuis la proches du plan d'eau étant évidemment les plus
cote 108 N G F , des surpressions allant jusqu'à 1,60 m sensibles (9); par contre, à l'arrière du plan d'eau, à
(cellule 3 profil A et cellule 1 profil B). cause de l'influence des drains, les réactions sont fai-
bles voire nulles. A u piézomètre 9, le niveau d'eau a
Certaines cellules ont peu ou pas réagi à la variation évolué de la même manière que les capteurs placés à
du plan d'eau : ce sont celles qui sont situées en l'avant de la digue : 1,60 m pour une élévation du
profondeur dans les couches argileuses (3-4 du plan d'eau de la retenue à 7 mètres.
profil B) ou dans les alluvions mais placées sous l'in-
fluence des drains et filtres (6-7 du profil A). Toutes ces données ont permis (fig. 11 et 12) d'estimer
le tracé de la surface libre (en effet avec le dispositif
d'observation utilisé, il n'est pas possible de détermi-
Piézomètres ouverts ner rigoureusement cette surface) en supposant le
milieu homogène et isotrope jusqu'au toit des argiles
Mis en place (fig. 16) juste avant la mise en eau, les à cérithes. L'influence du dispositif de drainage est
piézomètres ouverts, constitués par de simples tubes certaine.

18
 NGF (m)
Fig. 17 — Fluctuation des
niveaux piézométriques dans les

M | J | J | A | S | O | N | D | | J | F | M | A | M | J | J | A [ S | O | N | D |

1977 1978

Déplacements du corps de digue Analyse du comportement hydraulique de l'ouvrage.

pendant la mise en eau Comparaison des constatations avec les résultats
de modèles analogiques
Les repères disposés sur le remblai sont constitués de
plots en béton armé de section 0,50x0,50 m et de
1,50 m de profondeur (dont 1 m est ancré dans le Pour essayer de comprendre le fonctionnement
remblai). Dans ces bornes sont noyées des embases hydraulique réel de l'ouvrage, il a paru intéressant de
métalliques servant de repères topographiques; ceux- rapprocher les mesures et les observations des premiè-
ci permettent de suivre les tassements et les déplace- res analyses théoriques faites en 1979 dans le cadre
ments horizontaux du corps de digue ainsi que ceux d'une recherche sur l'anisotropie des limons utilisés
de la conduite sous remblai (repères 37-12-1). en corps de digue.

Les résultats sont donnés sur la figure 18, on constate Cependant, à cette époque, pour que cette étude com-
que les mouvements provoqués par la poussée du parative soit totalement satisfaisante, il manquait les
plan d'eau sont de l'ordre du centimètre et leur ampli- éléments suivants :
tude est maximale dans le tiers central du corps de
— des mesures de perméabilité dans le corps de digue,
digue (sur la vue en plan). Les mouvements verticaux
sont exclusivement des tassements, aucun soulèvement — des piézomètres dans le corps de digue et notam-
n'est à signaler; les déplacements horizontaux sont ment dans sa partie hydraulique.
orientés généralement dans l'axe de l'ancienne vallée
et vers l'aval, traduisant ainsi un resserrement et un De ce fait, nous avons été amenés à faire l'hypothèse
tassement de l'ouvrage dans l'axe de l'ancienne vallée. vraisemblable que le corps de digue possédait le même
Le long des ouvrages,, on peut observer quelques coefficient de perméabilité que le sol de fondation
mouvements anarchiques dus à un entraînement des (fc = 10~ m/s), moyennant quoi il était possible d'éta-
fc
6

terres lié à un mauvais compactage auprès des maçon- blir une comparaison entre les pressions interstitielles
neries. théoriques et mesurées et d'estimer la surface libre.

19
L'étude théorique hydraulique a été effectuée par la TABLEAU IV.
méthode de l'analogie électrique en utilisant successive- Comparaison des valeurs des potentiels théoriques
ment trois modèles. Les hypothèses et conditions aux et mesurées au d r o i t des cellules
des pressions interstitielles.
limites communes à ces trois modèles sont les
suivantes :
Cellules Modèle fc„/fc„=l Mesures in situ
— même perméabilité de la digue et de sa fondation, (potentiel en mètres) (potentiel en mètres)
celle-ci étant constituée par les alluvions (argile 1A 8,34 5,06
sableuse et argile sableuse plus graviers) ; 2A 9,10 5,31
3A 9,36 5,51
— le niveau imperméable est fixé au toit des argiles 4A 5,46 4,06
à cérithes (dans le modèle analogique une ligne de 5A 4,85 3,76
courant est imposée à cet endroit); 6A 2,40 1,01
7A 1,50 0,9
— une différence de potentiel est imposée entre le 8A 0 2,10
parement amont et l'exutoire du système drainant 1B 8,00 5,40
(tapis et puits aval). Elle correspond à la différence
de charge de 10 m entre le niveau de la retenue
normale (cote N G F 115) et l'exutoire (cote
105 N G F ) . moyenne de la couche, tendant à abaisser le réseau
de lignes de courant. Les différences s'atténuent lors-
que l'on se rapproche de l'exutoire, ce qui est normal.
Premier modèle (fig. 19)

E n utilisant la similitude des écoulements hydrauli-

Comme hypothèse particulière, nous avons supposé
ques et électriques, le débit de fuite par mètre linéaire
kJK=i.
capté par le tapis filtrant, en prenant
O n constate qu'il n'y a pas suintement le long de la k = lc = 1 0 " m/s (perméabilité moyenne), est :
h v
6

cheminée drainante dans ce cas, mais la ligne de

g = 3,2x 10~ 1/s par ml.
3

saturation se trouve plus haut que la ligne que nous

avons supposée réelle. Comme le montre le soit sur la longueur du barrage : Q =0,8 1/s.
tableau IV, où sont comparés aux niveaux des cellules
de pressions interstitielles les potentiels réels et les Les mesures faites sur seuil donnent un débit de 2 à
potentiels calculés, ces derniers étant beaucoup plus 3 1/s pour les fuites venant du tapis drainant. L a
forts, une différence de perméabilité doit exister entre valeur théorique est très proche de la réalité et l'écart
la digue et sa fondation. Des passages graveleux plus situé dans le bon sens confirme l'hypothèse du
perméables doivent augmenter la perméabilité contraste de perméabilité digue-sol de fondation.

20
 Fig. 20 — Deuxième modèle analogique.

t± = 25

Fig. 21 — Troisième modèle analogique.

Deuxième et troisième modèles (fig. 20 et 21) Le tableau V indique les valeurs de potentiels théori-
ques pour ces deux configurations et les valeurs des
Deux autres modèles ont été réalisés en milieu aniso- potentiels mesurés en place. Les valeurs théoriques et
trope dans un rapport de perméabilités horizontale et expérimentales sont en meilleure concordance qu'en
verticale de 25. Le premier modèle a été construit milieu isotrope, la meilleure adéquation étant obtenue
d'une façon homogène, le deuxième modèle a été dans le cas homogène avec k /k = 25 (k = 10~ m/s
h v h
6

construit en considérant l'influence d'un drain et k = 4 x 1 0

v m/s compte tenu de l'anisotropie).
- 8

« borgne » horizontal situé à la base de la digue, dans

le but de simuler la présence éventuelle d'une zone Le calcul du débit conduit à prévoir dans ce modèle
plus perméable dans la fondation. un débit de fuite sur la longueur de l'ouvrage de
0,22 1/s. On constate dans ce cas un écart important
Dans les deux cas, il apparaît un suintement le long avec la valeur mesurée, qui se situe entre 2 et 3 1/s.
de la cheminée drainante de 5 mètres de hauteur, avec
une ligne de saturation plus haute que dans le cas du E n conclusion de cette étude comparative, on constate
modèle homogène et isotrope. Le seul effet du drain que l'on ne peut faire coïncider totalement les observa-
« borgne » a été de modifier les équipotentielles de tions et les résultats théoriques quel que soit le modèle
part et d'autre du drain et de rabattre vers les limites analogique utilisé; toutefois avec la schématisation
du drain les lignes de courant venant de l'amont, le des modèles et l'imprécision des mesures le pourrait-
drain matérialisant une équipotentielle. on exactement ?

21
 TABLEAU V. Dans le domaine du comportement hydraulique de la
Comparaison des valeurs des potentiels t h é o r i q u e s
digue il existe, d'après les mesures faites sur les cellu-
et mesurées au d r o i t des cellules
des pressions i n t e r s t i t i e l l e s . les, une bonne corrélation des évolutions de pressions
interstitielles et du relèvement du plan d'eau : 1,60 m
Potentiel (m) entre la position initiale et finale pour les cellules les
plus proches de l'amont pour une élévation du plan
Cellules Modèle Modèle drain d'eau de 7 mètres à la mise en eau du bassin (cote 108
homogène borgne Mesures à 115 N G F ) , tandis qu'à l'aval les pressions sous
kjk. = 25 kj*. = 25 m situ
cheminée drainante ou sous le filtre ont très peu
1A 5,85 4,84 5,06 évolué (0,50 à 1,00 m environ), ce qui est normal. Les
2A 6,70 5,42 5,31 observations ne permettent cependant pas de mettre
3A 7,40 5,20 5,51 en évidence clairement le réseau d'écoulement réel et
4A 4,10 4,50 4,06
le comportement hydraulique des différentes couches,
5A 3,45 4,12 3,76
6A 1,52 3,35 1,01 même en s'appuyant sur des modèles analogiques
7A 1,32 2,08 0,9 qui prennent des hypothèses simplificatrices. Cette
8A 0,55 0,70 2,10 difficulté tient d'une part, au trop faible nombre des
points de mesures et, d'autre part et surtout, à la
1B 5,87 4,55 5,40
complexité du site du point de vue hydraulique en
particulier au niveau de la fondation. Il est donc
difficile dans ces conditions de porter un jugement
Stabilité des talus sur le rôle des différents éléments du système drainant.
C'est ainsi que celui de la cheminée drainante n'appa-
L a recherche des coefficients de stabilité du système raît pas clairement au niveau de la retenue normale.
digue-sol de fondation a été menée en rupture circu- Par contre, on peut émettre certains doutes sur la
laire. Tous les calculs ont été effectués en contraintes nécessité de puits aussi profonds à l'aval de l'ouvrage,
effectives avec des paramètres de cisaillement les mesures montrant en effet que les couches profon-
suivants : des sont très peu concernées par la mise en eau.
Cependant, s'il est aussi difficile d'analyser le compor-
tement hydraulique a posteriori, on conçoit que cela
le soit encore plus au stade du projet et que prévoir
Nature des couches
(°)
3
(kPa) (kN/m ) a priori un tel système drainant n'avait rien d'aberrant
Oet 5 en soi, surtout en se plaçant au niveau des plus hautes
Limon 1-2 32 19
(*)
eaux, les constatations ayant été faites à la retenue
Limon 3 0 36 21 normale.
Alluvions 0 34 2
Argiles a cerithes 5 25 19
Marnes a huitres 2 28 19 E n définitive le dispositif de surveillance prévu au
Calcaire de Brie 0 36 19 marché et complété par la suite a montré globalement
Marnes vertes 15 18 19
que le comportement de l'ouvrage était conforme à
ce que l'on pouvait attendre.
(*) 5 kPa pour la vidange rapide.
Malgré une certaine incertitude due à la schématisa-
tion des modèles analogiques sur la détermination des
Les valeurs des coefficients de stabilité minimaux réseaux d'écoulement à travers le corps de digue, il
sont : est probable que l'anisotropie du corps de digue soit
— parement amont en retenue normale : F = 1,81
en vidange rapide au plus hautes eaux : F = 1,06
— parement aval : F =1,41.
faible Q^<25^ et proche de l'isotropie compte tenu
de la méthode de mise en œuvre des matériaux et des
moyens de compactage (rouleau à pieds dameurs).

CONCLUSION

RÉFÉRENCES BIBLIOGRAPHIQUES
Les observations faites à l'occasion du suivi de la
digue de Maurepas-Courance ont mis en évidence
JOSSEAUME H., Digues en terre, Bull, liaison Labo. routiers
l'influence réciproque de l'ouvrage et du site. P. et Ch., hydraulique des sols, spécial N, avril 1970,
p. 199-214.
Les tassements dus aux poids des remblais (20 à
JOSSEAUME H., Le barrage de la Sorme. Étude et observa-
30 cm en moyenne) se sont produits essentiellement tions en cours de construction, Bull, liaison Labo. P. et
pendant la construction et leurs valeurs concordent Ch., 89, mai-juin 1977, p. 67-92.
assez bien avec les prévisions. NALLET J., Traitement des bassins de retenue d'eaux pluvia-
les dans la ville nouvelle de Saint-Quentin-en-Yvélines,
Les déplacements du corps de digue sous l'action de Les travaux publics, 1018, 3 trimestre 1977, p. 31-46.
e

la montée de l'eau sont de faibles amplitudes (de

POST G., LONDE P., Les barrages en terre compactée,
l'ordre du centimètre), montrant un léger affaissement Gauthier-Villars, Paris, 1953.
du remblai avec un mouvement général d'ensemble
vers le fond de l'ancienne vallée.

22
 Digues
et barrages en terre
/7/7////////////// 1 S de faible hauteur

Interprétation des mesures

de pression interstitielle dans la digue
de Maurepas-Courance
Anisotropie de perméabilité de l'ouvrage

Henri JOSSEAUME
Attaché de r e c h e r c h e

Cnahrokh MAHDAVI
Stagiaire
Section des ouvrages en terre
Laboratoire central des Ponts et Chaussées

Jean JARDIN
Ingénieur
Laboratoire régional de l'Ouest parisien

L a digue de Maurepas-Courance est une digue en terre de

15 m de hauteur maximale qui délimite un réservoir ayant
RÉSUMÉ une capacité de stockage de 390 000 m (hauteur maximale
3

de retenue 14 m). Ce réservoir, qui joue à la fois le rôle de

L'article rend compte de l'étude du comporte- bassin permanent et celui de bassin tampon, n'est rempli,
ment hydraulique en service d'une digue de rete-
nue homogène de 15 m de hauteur, construite en en période normale, qu'à environ la moitié de sa capacité
argile peu plastique compactée. Cette étude, (hauteur normale de retenue 10 m), le volume réservé per-
qui avait plus précisément pour objectif de mettant d'écrêter les crues du ru de la vallée de la Courance
caractériser l'anisotropie de perméabilité de
l'ouvrage, s'appuie : qui collecte une partie des eaux de ruissellement de la ville
nouvelle de Saint-Quentin-en-Yvélines.
- sur des mesures de pression interstitielle
dans l'ouvrage au moyen de l'équipement
piézométrique posé lors de sa construction et L'ouvrage, une digue homogène équipée d'un tapis drainant
d'un équipement piézométrique complémen- et d'un drain-cheminée, dont le projet a été établi par le
taire spécialement mis en place,
Bureau d'ingénieurs-conseils Coyne et Bellier, est décrit de
- sur des mesures de perméabilité en place
réalisées dans la digue et sa fondation, façon détaillée dans l'article de Delfaut et al. (1984) ainsi
- sur des études théoriques de l'écoulement que les études effectuées par le Laboratoire de l'Ouest
permanent à travers l'ouvrage, effectuées par parisien dans le cadre du suivi de son comportement pen-
la méthode des éléments finis en tenant compte
de diverses conditions d'anisotropie de perméa- dant et après sa construction. Outre des mesures de déplace-
bilité. ments, ce suivi comportait le contrôle des pressions intersti-
L'analyse des résultats expérimentaux et la tielles au moyen d'assez nombreux piézomètres répartis dans
confrontation des valeurs théoriques et expéri- deux profils transversaux, aussi Delfaut et al. ont-ils tenté
mentales de la pression interstitielle mettent en d'estimer l'anisotropie de perméabilité de la digue à partir
évidence un comportement quasi isotrope de la
digue vis-à-vis de la perméabilité. des pressions interstitielles mesurées. Cependant, la connais-
sance insuffisante de certaines données intervenant dans
MOTS CLÉS : Barrage en terre • Mesure - cette estimation n'a pas permis de conclure avec certitude,
Pression interstitielle - Interprétation - Piézo-
mètre - Perméabilité - Sous sol - Écoulement et une étude complémentaire s'est avérée nécessaire. Celle-ci
(fluide) - Éléments finis (méthode) - Isotropie - a été effectuée conjointement par le L R O P et la section des
Limon - Expérimentation - Calcul -/Digue - ouvrages en terre du L C P C dans le cadre des recherches
Bassin d'orage.
relatives aux digues de faible hauteur. L'article rend compte
de cette étude.

23
B u l l . Maison L a b o P. e t C h . • 1 3 1 - m a l - j u i n 1 9 8 4 - R é f . 2 9 1 1
P R I N C I P E D E L'ÉTUDE en principe, des perméabilités différentes, la zone 1
étant construite avec une argile un peu plus plastique
que celle constituant la zone 2 (fig. 1). Mais, comme
L'objectif de l'étude était de déterminer les valeurs du il s'est avéré très difficile de différencier ces matériaux
rapport d'anisotropie k /k (k h et k étant respective-
v h v lors de leur mise en œuvre, il est peu vraisemblable
ment les coefficients de perméabilité horizontale et que leurs comportements soient très différents.
verticale) des sols constituant le corps de digue et
sa fondation, à partir des mesures piézométriques Par ailleurs, les piézomètres mis en place lors des
effectuées dans ceux-ci au cours d'une période où travaux de construction, dans le cadre du suivi de
l'ouvrage est le siège d'un écoulement permanent, l'ouvrage, avaient principalement pour but de contrô-
c'est-à-dire en période normale, lorsque le niveau du ler les pressions interstitielles dans le sol de fondation
réservoir est stable depuis plusieurs mois. L a méthode et, de ce fait, un nombre limité d'entre eux a été posé
retenue pour l'interprétation de ces mesures consiste dans le corps de digue.
à comparer les valeurs expérimentales de la pression
interstitielle aux valeurs théoriques correspondantes Compte tenu de cette situation, il était nécessaire de
calculées en affectant successivement aux sols étudiés déterminer les perméabilités de la digue et de l'argile
différentes valeurs du rapport d'anisotropie. sableuse de fondation dans la zone des mesures piézo-
métriques ainsi que de compléter l'équipement du
L a mise en œuvre correcte de cette méthode corps de digue en piézomètres. Les mesures de perméa-
impliquait : bilité et la pose des piézomètres supplémentaires ont
été effectuées dans le profil de mesure A , c'est-à-dire
— d'une part, une connaissance suffisamment précise dans celui des deux profils de mesures instrumentés
des données nécessaires à l'étude théorique des écoule- lors de la construction, qui comportait le plus grand
ments et, en particulier, des contrastes de perméabilité nombre de piézomètres posés au-dessus du toit de
entre les différentes zones de l'ouvrage, notamment l'argile à cérithes. Les valeurs théoriques de la pres-
entre le corps de digue et sa fondation; sion interstitielle dans l'ouvrage en régime d'écoule-
— d'autre part, que la pression interstitielle puisse ment permanent ont été calculées par la méthode des
être mesurée en des points suffisamment nombreux et éléments finis au moyen du programme F L O T - N S A T
judicieusement choisis. [Humbert, 1984].

L'étude géotechnique du site [Delfaut et al, 1984]

donne une description précise du sol de fondation et
établit que l'argile à cérithes constitue un substratum PERMÉABILITÉS D E L A D I G U E
imperméable pour la couche d'argile sableuse sur ET D E SA F O N D A T I O N
laquelle est fondée la digue (fig. 1). Cette étude fournit
également des valeurs du coefficient de perméabilité Les essais complémentaires, ayant pour but de déter-
de l'argile sableuse mais celles-ci, très dispersées, miner les perméabilités respectives du corps de digue
devaient être précisées dans la zone des mesures piézo- et de la couche d'argile sableuse, étaient des essais de
métriques. D e plus, l'étude initiale ne donnait aucune perméabilité à charge constante effectués au pressio-
indication chiffrée sur la perméabilité du corps de perméamètre Ménard. Ils ont été réalisés dans deux
digue, constitué d'argile peu plastique compactée en sondages à la tarière hélicoïdale (diamètre 63 mm)
couches minces (0,30 m) au moyen d'un rouleau à pratiqués, l'un depuis la crête de la digue, l'autre
pieds dameurs. E n amont du drain cheminée, le corps depuis la surface du talus amont, un peu au-dessus
de digue est d'ailleurs subdivisé en deux zones ayant, du niveau atteint par l'eau dans le réservoir en période

T 119 N G F

Fig. 1 — Coupe transversale type de

l'ouvrage, et principales caractéristiques
d'identification de l'argile constitutive
du corps de digue.
Marne à huîtres

Caractéristiques d u corps de digue

Caractéristiques Caractéristiques optimales

Limites d'Atterberg
de compactage Proctor normal
Zone
w (%)
t /,(%) w (%) Y.ikN/m ) 3
w (%)
opl
ydmñx (kN/m ) 3

1 37<w <40 L 15</.<22 16<w<19 16,9<y <18,5i 15<w ,<15,5 18,l<y <18,9ÎMI

w = 38
L ï , = 18 w=17,8 Ï* =17,7 w ,= 15,3
op Y , ™ =18,4

2 26<w <28 t 10</,<12 ll,5<w<17 16,8 <y < 18,7

ä 12<w <15 opl 17,5< , <18,7
Yl mal

w = 27
t w=14,5 Yd = 17,8 w ,= 14
op Y*™* =18,3

24
 normale. Les essais ont dû être arrêtés à une profon- gène et isotrope. Il était alors possible de calculer les
deur de 17,50 m dans le premier sondage et à une valeurs k et fc des coefficients de perméabilité
h e q veq

profondeur de 4 m dans le second, cela en raison de horizontale et verticale équivalents et d'en déduire
.'instabilité des parois des sondages. le rapport d'anisotropie k Jk *. he O n a obtenu les
yeq

valeurs / c = 5 x 1 0 " m/s et fc = 8,5 x 10 m/s,

heq
6
veq
7

Les valeurs du coefficient de perméabilité mesurées qui correspondent à un rapport d'anisotropie de

au niveau des essais sont reportées sur la figure 2. Si l'ordre de 6.
l'on admet que chacune de ces valeurs est représenta-
tive de la perméabilité de la tranche de sol limitée par
les plans horizontaux situés respectivement à mi-
ÉQUIPEMENT COMPLÉMENTAIRE
distance du point de mesure considéré et des points de
E N PIÉZOMÈTRES
mesures immédiatement supérieur et immédiatement
inférieur, on obtient les profils de perméabilité repré-
sentés sur la figure 2.
L'équipement piézométrique mis en place dans les
profils de mesures A et B lors de la construction
1 0 " 6
10" Mm/s)
S
1 0 "7 1 Q - 6 -s10 1 Q - 4 fc(m/ )s
de la digue comporte quinze piézomètres classiques
I 1 —>- constitués de tubes de 60 mm de diamètre intérieur,
crépines sur une hauteur de 1,50 m, et de douze piézo-
2 mètres de type P A C [Levillain, 1976]. Comme indiqué
tenue normale o 4 précédemment, la plupart de ces appareils ont été
6 posés dans les couches de fondation et seuls cinq
•
8 piézomètres classiques permettaient d'effectuer des
- ^ ^ ^ 4
\ 10
mesures dans la digue.
12
' z (m) 14
•16 Six piézomètres supplémentaires ont été mis en place
18 dans la digue au niveau du profil de mesures A . Leur
Tz(m)
position dans ce profil a été fixée de façon à obtenir
la meilleure image possible de la distribution des pres-
Fig. 2 — Résultats des essais de perméabilité en place réalisés sions interstitielles dans la digue en période normale,
dans la digue et sa fondation. Profils de perméabilité.
c'est-à-dire pour une hauteur de retenue de 10 m. Ces
piézomètres sont des piézomètres classiques constitués
de tubes de 12 mm de diamètre intérieur, crépines sur
O n constate que, sur les onze essais effectués dans une hauteur de 0,50 m (ces caractéristiques ont été
le corps de digue, neuf fournissent des valeurs du adoptées afin d'obtenir un temps de réponse plus
coefficient de perméabilité comprises entre 1 0 " m/s 6
court que celui des piézomètres classiques déjà en
et 4 x l 0 ~ m / s (moyenne 2,4 x 10^ m/s), les deux
6 6
place et de réaliser des mesures mieux localisées).
autres conduisant à des valeurs qui se différencient Leur situation dans le profil A est représentée sur la
nettement des précédentes et les encadrent figure 3, ainsi que celle des piézomètres mis en place
(fe = 1,25 x 10 m/s etfe= 2 x 10~ m/s). O n constate
7 5
initialement dans ce profil et dans le profil B, au-
également que les essais ne mettent pas en évidence dessus du toit de l'argile à cérithes (les profils A et B
des différences significatives entre les perméabilités de étant pratiquement identiques aux plans de la géomé-
la zone amont (zone 1) et de la zone centrale (zone 2). trie et des propriétés géotechniques, on a tenu compte,
dans l'étude de l'anisotropie, des mesures effectuées
Les trois valeurs du coefficient de perméabilité mesu- dans les deux profils).
rées dans la moitié supérieure de la couche de fonda-
tion sont relativement groupées (10~ m/s</c 5

<4 x 1 0 m/s) et permettent une estimation correcte

- 5

(tout au moins dans la zone des essais) du contraste ETUDE EXPERIMENTALE

de perméabilité entre la digue et sa fondation. E n D E L'ÉCOULEMENT.
effet, le coefficient de perméabilité moyen de cette DÉTERMINATION D E L A S U R F A C E LIBRE
dernière, égal à 2 x 1 0 " m/s, est à peu près dix fois
5

supérieur à celui du corps de digue.

L'écoulement à travers la digue et la couche d'argile
Les essais de perméabilité en place mesurant le coeffi- sableuse a été étudié à partir des mesures de pression
cient de perméabilité horizontal du sol ne fournissent, interstitielle effectuées six mois après la pose des piézo-
en principe, pas d'information sur l'anisotropie de mètres supplémentaires (les résultats des mesures sont
perméabilité. Cependant, celle-ci étant liée, pour une donnés sur la figure 3). L a surface du réservoir et les
grande part, à l'hétérogénéité du sol, une faible disper-
sion des mesures en place indique (si elles sont suffi-
* On rappelle que les valeurs de k et k relatives à un
samment nombreuses et si le sol n'est pas finement
hel< veq

ensemble de n couches homogènes et isotropes, d'épaisseurs

stratifié) une anisotropie modérée. L a relative homogé- respectives â et de coefficients de perméabilité respectifs k
t {

néité des résultats des mesures effectuées dans le corps (i variant de 1 à n), sont obtenues à partir des expressions :
de digue donne donc à penser que celui-ci est caracté-
risé par un rapport d'anisotropie assez faible. O n
a tenté d'apprécier sa valeur à partir du profil de
perméabilité relatif à la partie centrale de la digue, en *heq " et
faisant l'hypothèse que chaque tranche de sol prise
en compte pour l'établissement de ce profil était homo- 14

25
 -Profil A Piézomètre A A
Ai A 2 A 3
4
A 5 A 6 A, A 8 A t l 1 2
A» T 2 3 T 4
2 T 2 5 T 2 6 T 2 7 T 8 2

u(kPa) 77 82 57 44 45 15 22 43 0 45 22 0 15 0 25 6 8

Profil B
Piézomètre Bi B B B B B,
A Piézomètre 0 60 mm crépine sur 1,50 m
2 6 7 8

• Piézomètre tf> 12 mm crépine sur 0,50 m 62 80 7 51 0 52

u(kPa)
• Piézomètre P A C

Retenue normale

Fig. 3 — Situation des piézomètres dans les profils de

mesures A et B et valeurs des pressions interstitielles
mesurées en régime d'écoulement permanent.

pressions interstitielles n'ayant subi que de faibles représentant la digue et la couche de fondation (argile
fluctuations au cours de ces six mois (l'amplitude sableuse), les conditions aux limites imposées dans les
maximale de variation du niveau du réservoir était de différents calculs effectués sur ce modèle étant celles
0,30 m et la plupart des piézomètres indiquaient des qui auraient été prises en compte dans un calcul
variations de pression interstitielle inférieures à prévisionnel classique. Dans la seconde, les écoule-
2 kPa), le régime d'écoulement à l'époque des mesures ments respectifs dans la digue et la fondation ont été
pouvait être considéré comme permanent. étudiés à partir de modèles numériques distincts et
certaines conditions aux limites ont été fixées à partir
L a surface libre de l'écoulement a été déterminée à des valeurs expérimentales de la pression interstitielle.
partir des mesures piézométriques effectuées dans la
digue. Cette surface est située dans une zone limitée TABLEAU I. - V a l e u r s r e l a t i v e s d e s c o e f f i c i e n t s d e
à sa partie supérieure par les piézomètres les plus bas p e r m é a b i l i t é a t t r i b u é s a u c o r p s d e d i g u e e t à la f o n d a -
t i o n dans les é t u d e s d ' é c o u l e m e n t e f f e c t u é e s à p a r t i r d u
dans lesquels la pression est nulle ( T et T ) et, à
m o d è l e r e p r é s e n t é s u r la f i g u r e 5a.
2 3 2 5

sa partie inférieure, par les niveaux piézométriques

mesurés dans les piézomètres les plus élevés qui indi-
quent une pression non nulle ( T ) . Le point de pas-
2 8
Cas Corps de digue Fondation
sage de la surface libre au droit du piézomètre B 6
1 k„ = k„=l k„ = k„ = 10
correspondant, d'autre part, au niveau d'eau mesuré
dans ce piézomètre (puisque ce niveau d'eau a été 2 K=K=\ k =10,
H k„=l
relevé sur la hauteur crépinée), la surface libre a pu
3 fc„ = k„=l *. = *.= 15
être obtenue de façon assez précise, sauf dans sa
partie aval au voisinage du drain cheminée (fig. 4). 4 k = 25, k„=l
h k = k =
h v 250
L'allure de la surface libre suggère que l'anisotropie
du corps de digue, vis-à-vis de la perméabilité, est 5 k* = 25, fc„=l k„ = 250, fc„=10
assez faible.

Le modèle utilisé dans le cadre de la première procé-

dure est représenté sur la figure 5 a ainsi que les condi-
tions aux limites prises en compte dans les calculs. Les
valeurs des coefficients de perméabilité horizontale et
verticale successivement affectées à la digue et à sa
fondation sont consignées dans le tableau I. Le choix
de ces valeurs a été guidé par les considérations
suivantes :
Fig. 4 — Représentation de la surface libre déterminée à partir des mesures — les essais de perméabilité en place montrent que
effectuées dans les piézomètres supérieurs.
le rapport des perméabilités horizontales de la digue
et de la fondation est de l'ordre de 10;
— la forme de la surface libre indique que la digue
INTERPRÉTATION D E S M E S U R E S est isotrope ou faiblement anisotrope;
P A R L E C A L C U L . ÉTUDE D E L ' A N I S O T R O P I E — les résultats obtenus par Delfaut et al. ne permet-
D E PERMÉABILITÉ taient pas d'exclure l'éventualité de rapports d'aniso-
tropie de l'ordre de 25.

Deux procédures de détermination des valeurs théori- Dans aucun des cinq cas étudiés, la concordance entre
ques de la pression interstitielle ont été successivement les valeurs théoriques et expérimentales de la pression
employées. L a première a consisté à étudier l'écoulé- interstitielle n'est satisfaisante. D e ce fait, les résultats
ment dans l'ouvrage à partir d'un modèle numérique obtenus ne seront pas détaillés ici. O n se contentera

26
 À

(a) figure 5 b. L a distribution de la charge hydraulique a

été imposée le long de l'interface digue-fondation.
Puits de Cette distribution, définie par les valeurs de la charge
\ décharge mesurées le long de cet interface dans les
h =4 m I piézomètres A , B , A et A , a été supposée linéaire
3 9 4 6

h = 4 m dans les intervalles limités par ces piézomètres ainsi

qu'entre le pied amont de la digue et le piézomètre A 3

mperméable
et entre le piézomètre A et le pied du drain-cheminée.
6

L'étude théorique de l'écoulement à travers chacun

de ces modèles a été faite pour les valeurs 1, 9, 16,
25 du rapport d'anisotropie.
h = 14 m
Les valeurs de la pression interstitielle ainsi calculées
dans la couche de fondation sont pratiquement indé-
h = 9,7 m /i = 9,2 m / i = 8,4 m / ; = 5,5 m
pendantes du rapport d'anisotropie et, à l'emplace-
ment des piézomètres, sont très proches des valeurs
/7=14m \ I x"' * = 4 m
mesurées (l'écart varie de 1 à 7 kPa sauf au

II
-+-
/) = 4 m A, B Q piézomètre B où il atteint 22 kPa). Mais ces résultats
2

Fig. 5 — Schémas des modèles utilisés pour l'étude des écoulements par
ne sont guère significatifs, car le fait de fixer les
la méthode des éléments finis. pressions interstitielles le long du toit de la couche de
a) la digue et la couche de fondation sont représentées par un modèle fondation caractérisée par une épaisseur très faible,
unique.
b) la digue et le sol de fondation sont représentés par des modèles revient pratiquement à imposer le champ de pression
distincts. interstitielle dans cette couche.

E n revanche, la distribution théorique des pressions

interstitielles dans la digue varie avec le rapport d'ani-
de signaler que la moins mauvaise concordance entre sotropie et la comparaison des valeurs mesurées et
les mesures et le calcul est obtenue dans les cas 1 et des valeurs calculées à l'emplacement des piézomètres
3 (digue isotrope sur fondation isotrope). (tableau II) montrent que la meilleure concordance
est obtenue entre celles-ci dans le cas où la digue est
Les écarts entre les distributions théorique et réelle supposée isotrope. L a concordance entre mesures et
de la pression interstitielle doivent être imputés princi- calcul peut alors être considérée comme satisfaisante
palement à un manque de représentativité du modèle pour les valeurs correspondant à sept piézomètres sur
au niveau de la fondation, notamment au fait que dix. Il apparaît également que, dans la zone où elle a
cette dernière, considérée comme homogène, est, en pu être définie avec une précision suffisante, c'est-à-
fait, hétérogène et, éventuellement, à une mauvaise dire sur la majeure partie de sa longueur, la surface
appréciation des conditions imposées le long des limi- libre observée est très proche de la surface libre calcu-
tes verticales du modèle. Aussi, pour éviter que les lée dans l'hypothèse d'un ouvrage isotrope (fig. 6 a)
erreurs liées à la modélisation de la fondation ne se et ne présente guère de similitudes avec la surface
répercutent sur les valeurs théoriques de la pression libre théorique correspondant à un rapport d'anisotro-
interstitielle dans la digue, a-t-on étudié séparément pie de 9 (fig. 6 b). Il s'ensuit que le comportement
l'écoulement dans la digue et dans la couche de fonda- hydraulique global de l'argile constitutive du corps de
tion à partir des deux modèles représentés sur la digue est celui d'un matériau isotrope ou faiblement

TABLEAU II — C o m p a r a i s o n d e s v a l e u r s d e la p r e s s i o n i n t e r s t i t i e l l e m e s u r é e s d a n s la d i g u e
et des valeurs théoriques correspondantes calculées pour différentes valeurs d u rapport
d ' a n i s o t r o p i e kh/kv, à p a r t i r d u m o d è l e r e p r é s e n t é s u r la f i g u r e 5b.

Pression Pression interstitielle calculée (kPa)

Piézomètre interstitielle
mesurée (kPa) k /K=l
h kjk. = 9 kjk = 16
v Mt„ = 25

A„ 0 0 0 0 0,5

A,3 22 5 11 13,5 13,5

B 6
7 3,5 9 12,5 12,7

B 8
0 0 1 3 7,5

T 2 3
0 0 7,5 12 14,5

T 2 4
15 21,5 24,5 26 26,5

T 2 5
0 7,5 17,5 20,5 22,5

T 2 6
25 37 40 41 41,5

T 2 7
6 24 29 31 32

T 2 8
8 10 17,5 20 22

27
 (a)
CONCLUSION

L'interprétation des mesures piézométriques et des

mesures de perméabilité en place effectuées dans la
digue de Maurepas-Courance et dans sa fondation, a
permis de préciser le comportement hydraulique de
l'ouvrage.

L a dispersion relativement faible des mesures de

perméabilité réalisées dans l'argile peu plastique consti-
tuant le corps de digue, l'allure de la surface libre
de l'écoulement déterminée expérimentalement et la
comparaison des valeurs mesurées de la pression inter-
stitielle aux valeurs théoriques correspondantes calcu-
lées en tenant compte de différentes valeurs du rap-
Fig. 6 — Représentation de réseaux d'écoulement théoriques dans le port d'anisotropie de perméabilité, mettent en évi-
corps de digue et comparaison des surfaces libres théoriques et observées. dence le comportement quasi isotrope de l'ouvrage.
a) cas ou kh/kv = l E n revanche, l'anisotropie de la couche de fondation
n'a pu être précisée.
b) cas ou kh/kv=9
Ces réseaux o n t été établis à partir d u modèle de la figure 5 b.
R É F É R E N C E S BIBLIOGRAPHIQUES

DELFAUT A., JARDIN J . , B A L D I T R., Les bassins de Maurepas-

anisotrope. Ce comportement peut sembler assez sur- Courance et des Roussières. Constatations sur la digue
prenant en raison du fait que l'argile ne peut être de Maurepas-Courance, Bull, liaison Labo. P. et Ch.,
considérée comme homogène (en particulier, sa plasti- 131, mars-avril 1984.
H U M B E R T P., Application de la méthode des éléments finis
cité varie dans d'assez larges limites) mais est cohérent
aux écoulements en milieu poreux, Bull, liaison Labo. P.
avec la dispersion assez faible des perméabilités mesu-
et Ch., 1984, à paraître.
rées dans la digue, ainsi qu'avec le mode de compac- L E V I L L A I N J.-P., Mesures pneumatiques. La sonde piézomé-
tage de l'argile (compactage en couches minces au trique LPC, Bull, liaison Labo. P. et Ch., 76, mars-avril
rouleau à pieds dameurs). 1975, p. 19-24.

28
 Digues
et barrages en terre
r/////////////////
~7L de faible hauteur >

Constatations sur la digue des Roussières

Jean JARDIN
Ingénieur

Robert BALDIT
Technicien supérieur

Abel DELFAUT
Assistant
Laboratoire régional de l'Ouest parisien

L a digue des Roussières, située sur le territoire de la ville

nouvelle de Saint-Quentin-en- Y vélines, permet de réguler
RÉSUMÉ l'évacuation d'une partie des eaux pluviales du centre ville
par la création d'un bassin tampon dont la capacité est de
Située sur le territoire de la ville nouvelle de 400000 m aux plus hautes eaux ( P H E à 150,40 N G F ) .
3

Saint-Quentin-en-Yvelines (département des

Yvelines) la digue des Roussières est un barrage
de hauteur moyenne (environ 15 m) qui ferme L'ouvrage a 238 m de long, 17 m de haut au maximum,
150 m (en pied) et 65 m (en crête) de large, la voie
3
une retenue d'eau pluviale de 400 000 m .
C'est une digue homogène en limon avec un
dispositif de drainage (cheminée et tapis drai- routière E 2 étant liée à la digue proprement dite (fig. 1).
nant). Bien que de dimensions relativement A u niveau normal (143 N G F ) , le plan d'eau présente une
modestes, l'ouvrage qui intéresse la sécurité surface de 3 ha et la hauteur d'eau est d'environ 7 m en
publique a été équipé d'un important dispositif
d'auscultation. Des constatations sur son pied de digue. L'ouvrage est de type homogène en limon
comportement hydraulique et mécanique ont avec des pentes amont de 3/1 et aval de 2/1. Il comporte
été effectuées lors d'un remplissage aux plus un dispositif de drainage constitué d'une cheminée centrale,
hautes eaux, réalisé en fin de construction et
durant deux années après la mise en service de trois tranchées et d'une ligne de puits filtrants au pied
au niveau permanent. Plus particulièrement, aval. Le parement amont est protégé par des enrochements.
après mise en place des piézomètres complé- Le bureau d'étude Coyne et Bellier a été chargé des études
mentaires, on a étudié au niveau permanent
l'anisotropie hydraulique du limon constituant préalables et de la conception de l'ouvrage. Sous sa direc-
le corps de digue. Après un rappel des études tion, le L R O P a réalisé la reconnaissance géotechnique.
effectuées au stade du projet, l'article rend
compte des résultats de ces constatations. Le
comportement quasi isotrope (k /k «s 1) du Bien qu'elle ne relevât pas du contrôle du comité technique
n v

limon constitue le résultat le plus important permanent des grands barrages, la digue, située directement
de ces observations. Les auteurs pensent que à l'amont d'une petite agglomération, a été l'objet d'un soin
ce fait est lié pour l'essentiel au mode de mise
en œuvre du matériau. particulier du point de vue de la sécurité (application de la
circulaire ministérielle 70.15 du 14 août 1970 (fascicule du
MOTS CLÉS : Barrage en terre - Limon - ministère de l'Équipement et du Logement n° 94/70
Drainage - Auscultation - Isotropié • Perméa- texte 755 (70. 94) concernant l'inspection et la surveillance
bilité - Alluvion • Sable - Sous sol - Stabilité - des barrages intéressant la sécurité publique). Construite
Drainage • Mesure • Déplacement • Modèle
analogique • Écoulement (fluide) • Pression entre 1976 et 1979, elle a été équipée d'un système de
interstitielle-/Digue - Bassin d'orage. contrôle comprenant des repères topographiques et des pié-
zomètres. Soumise en fin de construction à l'épreuve d'une

B u l l , l i a i s o n L a b o P. e t C h . - 1 3 1 - m a i - j u i n 1 9 8 4 - R é f . 2 8 9 1 29
Cote 3
Niveau de la nappe du sable de Fontainebleau
NGF Plus hautes eaux Enrochement ~j avant travaux de la digue (avril-juin 1974)
149,50 du 30 avril 1979 / 1 ^ Cheminée drainante
9
150 au 4 mai 1979 / n

ll\
Niveau piézomètre enregistré le 6/4/74 avant Tranchées drainantes Puits filtrants
la mise en eau

Fig. 1. — Profil en travers de la digue des Roussières.

première mise en eau à la cote maximale, elle a ensuite Géologie du site

fait l'objet d'un suivi pendant deux ans après sa mise
en service. O n a saisi l'opportunité de ces opérations,
d'une part, pour effectuer des constatations sur le D u point de vue géologique (fig. 2) le projet est situé
comportement hydraulique de l'ouvrage lors de la en bordure du plateau de Trappes, constitué, sous
première mise en eau et, d'autre part, après mise en une couche superficielle de limon, par la formation
place des piézomètres supplémentaires, pour étudier, des argiles à meulière et calcaires de Beauce reposant
au niveau permanent, la question de l'anisotropie sur les sables de Fontainebleau dont l'épaisseur est
hydraulique du limon constituant le corps de digue. de l'ordre de 50 m. Le toit des sables forme une dune
L a présentation des résultats de ce travail constitue au relief bien marqué, orientée N O - S E , que la petite
le sujet principal de cet article. A u préalable, on Bièvre entaille en créant un versant abrupt au nord
présentera un ensemble d'informations concernant et un versant à pente douce au sud. Dans la partie
l'ouvrage : centrale de l'ouvrage, le sol est constitué par des
— les informations directement nécessaires à la com- alluvions reposant directement sur les sables de Fontai-
préhension du sujet, nebleau. Les sols sont baignés par la nappe des sables
— les informations présentant un intérêt du point de de Fontainebleau, qui est en charge sous les alluvions.
vue plus général de l'étude des digues.

30
 1

Fig. 3. — Analyse granulométrique d u limon de corps de digue. Les constatations et études hydrauliques ayant été
faites dans la partie centrale de l'ouvrage, on indique
ci-après les caractéristiques géotechniques des trois
catégories de sols concernées.

Limon constitutif du corps de digue

Extrait sur le plateau à proximité du site, il s'agit
d'une argile peu plastique (classification L P C : Ap)
d'une granularité particulièrement homogène à
l'échelle de l'emprunt (cf. l'étroitesse du fuseau
— fig. 3) et dont les caractéristiques d'identification
sont les suivantes :
2 1 0,4 3,2 0,10,08 5 0 i , m 20,.ml0, m
J 5 m 7 .< m li.mQ.5ym 0,2,,
— Teneur en eau : 16 %<w<25 % , moyenne=22 %
^ Diamètres équivalents (f, m l .
— Limite de liquidité : 30 % < w < 3 8 % L

— Indice de plasticité : 1 0 < / < 1 7 p

— Optimum proctor normal : y = 17,5 k N / m , d

3

Nombre de w = 1.6,5%
mesures
70.
L'utilisation de ce matériau posait un problème en
Teneur en eau de raison de ses fortes teneurs en eau. Grâce à un assèche-
compactage après ment lors du terrassement (aération par pulvimixer),
aération à l'emprunt
N = 294 la dispersion des teneurs en eau, et donc des poids
w % moyen = 1 7 , 9 volumiques secs, a pu être limitée (fig. 3) lors de
1 —
| Écart-type = 1,9
la mise en œuvre, sans pouvoir toutefois satisfaire
complètement les règles de l'art en la matière (45 %
des valeurs du poids volumique sec sont inférieures à
50
95 % de l'OPN) et cela malgré une énergie de compac-
Teneur en eau in situ tage correcte (couches de 0,30 m d'épaisseur maxi-
(mai 1974) male compactées à l'aide d'un rouleau à pieds
w % moyen = 21,9
40
Écart-type = 1,67
dameurs type P D 1 et P D 2 dans la classification
RTR*.
T
Histogramme
30
des teneurs
en eau d u limon. Couche constituant la fondation de la digue

D'une épaisseur maximale de l'ordre de 4 m, cette

couche est constituée d'un mélange (colluvions) d'allu-
vions et de sols en place. L a largeur du fuseau granulo-
métrique (fig. 4) traduit cette hétérogénéité bien qu'en
10 -. moyenne la partie supérieure (2 m) de la couche soit
plus fine (elle tend vers un sol type Ap) et la partie
inférieure, plutôt constituée d'un sable argileux renfer-
mant des petits blocs et grains de meulière, plus gros-
10 12 14 16 18 20 22 24 26 sière (sol type SA). L'ensemble est peu plastique
Teneur en eau (%)

Nombre
de mesures Fréquence^

Histogramme
des densités sèches
du limon.

Références OPN

7 OPN
d
7d 17,5 17,6 17,8

.95 % w 16,5 16,7 16

OPN U-'
1
[s*
20 10 5 0.4 0.2 0.10.08 5 0 p m 20^m10pm 5(im 2u

_ Diamètres équivalents 1^ m)

^1 Fig. 4. — Analyse granulométrique des alluvions.

I L, Td ( k N / m )4

* RTR : Recommandation pour les terrassements routiers, SETRA-

LCPC, 1976.

31
(J max = 14, w
P = 42), peu compressible et relative-
i r a a i

ment consistant, les essais au scissomètre et au pressio- Nombre Première médiane Deuxième médiane
de mesures
mètre ayant donné les résultats moyens suivants :

c„ = 35 kPa, /?, = 380 kPa, £ M = 2 500 kPa.

Sable de Fontainebleau 10" 6

6.10" 6
10' 5
5.10 5
10"
k (m/s)
h
Nombre total de mesures : 68
Il s'agit d'un sable fin (classification L P C Sm—SL)
compact (/?,>2000 kPa), dont la partie supérieure
(les dix premiers mètres approximativement) est en Fig. 6. — Histogramme des perméabilités
du sable de Fontainebleau.
moyenne plus fine que la partie inférieure (fig. 5).

GROS SABLE SABLE FIN LliVlON

4J -
i •I I Í . ! I
z • T V *

tt/h" inférieure, caractérisée par la valeur médiane

4 - k = 5.10"
h m/s. U n essai de pompage, réalisé dans
5

V-
y ~f *
*
- '• i : i - la partie inférieure, a également donné la valeur de
i Partie s u p é r i e u r e
-A d u sable d e 5.10 m/s. Quatre essais au perméamètre de labora-
5

-t-h-— '•1 Fnntainphlpai i

. .i -
toire ont donné les valeurs indiquées dans le tableau I
dbsi i it-cuiui i
pour la perméabilité verticale.
"{
Sable d e
Fontaineblea
U PC :' S m -
—'—l~~
r--
S ,—*

classification
1
1
y
— 1 PC • S m
Fondation
.[-i 4 , -. - i —
-4-
- . __
! •
S 30
'-' - T--T- Seuls deux essais Lefranc ont été effectués dans cette
couche.
—_
—• r-

r~--
Limon du corps de digue
.t
• 2 1 0,4 0,2 0,1 0 , 0 8 5 0 > . m 20/jmlO/jm 5/. m 2/.m
L a perméabilité a été déterminée lors d'essais oedomé-
m m

D i a m è t r e s é q u i v a l e n t s (¡1 m )
triques réalisés sur des éprouvettes de limons compac-
Fig. 5. — Analyse granulométrique d u sable de Fontainebleau. tés à l'énergie Proctor à des teneurs en eau proches de
cet optimum ( O P N , O P N + 2, O P N + 4) et soumises à
des contraintes verticales croissantes (100 à 400 kPa).
O n obtient / c „ < 1 0 ~ m/s. 8

Connaissance de la perméabilité au stade de l'étude du Le tableau I fait clairement apparaître les insuffi-
projet sances de la connaissance de la perméabilité au stade
de l'étude, à la fois du point de vue du contraste entre
le corps de digue et les sols en place sous l'ouvrage
Sable de Fontainebleau
et également du point de vue de 1' anisotropie (les
Mesurée par soixante-huit essais Lefranc, la perméabi- quelques valeurs de kjk obtenues dans le sable en v

lité du sable de Fontainebleau confirme l'existence comparant les essais au perméamètre sur échantillons
de deux niveaux dans la couche (fig. 6) : la partie intacts aux essais Lefranc effectués à proximité des
supérieure (environ les dix premiers mètres) caractéri- niveaux de prélèvement des échantillons ne sont pas
sée par la valeur médianefe = 5 . 1 0 m/s, et la partie A
- 6 significatives).

TABLEAU I

Connaissance de la perméabilité au stade de l'étude du projet.

k (m/s)
h k (m/s) v kJK

Pressions
verticales | OPN OPN + 2 OPN + 4
d'essai (kPa) |
Néant ?
Limon (corps de digue) 100 ; 10" 8
2. H T 1 0
4.10" u

(aucun essai)
9 1 0 11
200 2.10" LUT 6.10"
10 1 1 11
400 10.10" 5. H T 3.10"

Fondation <10" 8 Í 7
i
er 5
1 ensemble: 5.10"
e 6 7 -6 22, 4,
Sable de Fontainebleau 2 ensemble : 6.10 3,5.10" < <2.10
5 1,5, 11
5.10 ~ essai pompage

32
Niveau s t a t i q u e d e la
Hypothèses ont la même perméabilité.
• Régime permanent au niveau des plus hautes • L'ensemble des sols est isotrope.
nappe d u sable de
eaux. • Le bassin ne communique pas avec la nappe
Fontainebleau • Les limons du corps de digue et les alluvions du sable de Fontainebleau.
Ligne de
Plus hautes eaux
/ saturation Niveau s t a t i q u e de la
NGF 150,40
nappe d u sable de
Fontainebleau
Puits f i l t r a n t

Fig. 7. — Réseau d'écoulement construit par analogie électrique.

STABILITÉ DE L'OUVRAGE Étude hydraulique

L'étude des conditions hydrauliques de la vidange
L'ouvrage a été conçu par le bureau Coyne et Bellier
rapide a nécessité la construction préalable du réseau
sur la base de son expérience. Le maître d'oeuvre a
d'écoulement en régime permanent avec un niveau du
demandé au L R O P de vérifier sa stabilité. O n pré-
bassin correspondant à celui des plus hautes eaux
sente ci-après les résultats de cette étude de stabilité
(150,40 N G F ) . Construit par analogie électrique en
et, en particulier, l'étude hydraulique effectuée dans
supposant que les sols étaient isotropes et que le corps
ce cadre et qui a servi de référence aux constatations
de digue et la fondation avaient la même perméabilité,
et études ultérieures.
ce réseau est présenté sur la figure 7. O n constate :

— l'efficacité apparente du système drainant. O n

Stabilité à court terme peut remarquer toutefois que la nécessité de prolonger
le tapis drainant par des tranchées ancrées au toit du
Avec une cohésion c de 35 kPa pour la couche de
u sable de Fontainebleau (surtout pour ce qui concerne
fondation, le coefficient de sécurité vis-à-vis du poin- la tranchée amont et la tranchée centrale) n'apparaît
çonnement déterminé par la méthode de Mandel- pas clairement;
Salençon était de 2,7. E n tenant compte d'une cohé-
— le rôle de drain joué par le sable de Fontainebleau;
sion c„ du limon du corps de digue égale à 80 à
100 kPa et déterminée par des essais triaxiaux U U — le fait que la cheminée drainante n'est pas
sur matériau compacté, le coefficient de sécurité mini- concernée par le réseau hydraulique.
mal des talus de l'ouvrage était de 1,5.

Stabilité à long terme CONSTATATIONS FAITES AU COURS

DE LA PREMIÈRE MISE EN EAU
Celle du parement amont a été étudiée dans le cas
d'une vidange rapide totale du bassin à partir des Le dispositif d'auscultation de l'ouvrage (fig. 8)
plus hautes eaux. comprend :

Avec les caractéristiques mécaniques suivantes : — des repères topographiques : constitués d'un embout
- limon : cp' = 32°, c' = 5 kPa, métallique arrondi, ces repères sont fixés soit sur des
- fondation : <p' = 35°, c'=0 bornes en béton de 0,50 x 0,50 m ancrées de 0,50 m
dans le remblai (on compte trente-six repères de ce
on obtient : F>2,10.
type), soit sur les parties en béton de l'ouvrage (douze
E n fonctionnement permanent, avec les mêmes carac- repères de ce type). Quatre piliers en béton érigés sur
téristiques, le talus aval présente un coefficient de un massif de béton de 1 x 1 x 1 m, lui-même solidarisé
sécurité : F > 1,50. au sol par des fers à béton de 25 mm, constituent

33
 B A S S I N

0 20 m
I I I I p

• Bornes topographiques © P o i n t sur c o n s t r u c t i o n # E n r e g i s t r e u r d u d é b i t des e a u x d e d r a i n a g e

• P o i n t en galerie A Dispositif Vinchon • P i é z o m è t r e de p e t i t d i a m è t r e d a n s d i g u e o u f o n d a t i o n
A Pilier de référence t o p o g r a p h i q u e O Enregistreur du niveau du plan d'eau © P i é z o m è t r e d e gros d i a m è t r e

Fig. 8. — Plan d'auscultation de l'ouvrage.

le quadrilatère de référence. Effectué par l'Institut —des appareils Vinchon : scellés sur la clé de voûte de
géographique national, le suivi topographique a com- la galerie sous remblai, ils permettent de mesurer
porté des mesures fines (précision : ± 2 mm en plani- l'ouverture et la femeture des joints;
metrie et ± 1 mm en altimétrie) et des mesures couran-
— une échelle limnimétrique fixée le long de la tour
tes (précision : +1 cm en planimetrie et + 2 mm en
et un limnigraphe à flotteur, permettant de suivre
altimétrie);
l'évolution du niveau du plan d'eau du bassin ;
— des piézomètres : chaque profil (nord et sud) com- — un seuil situé dans le collecteur des eaux de drai-
porte dix-neuf piézomètres, dont quinze prévus dans nage, qui permet d'en mesurer le débit et son évolu-
le cadre de l'auscultation et quatre (n° 18, 19, 20, 21) tion dans le temps.
mis en place au titre de l'action de recherche sur
l'anisotropie. Certains (piézomètres de type 2 et 3, S'inspirant de la circulaire ministérielle n° 70.15 du
fig. 17) sont crépines sur 0,50 m (limon-fondation) ou
14 août 1970, l'essai de mise en eau, réalisé entre le
1 m (sable) d'autres (piézomètres de type 1, fig. 17)
6 avril et le 10 juin 1979, a consisté en un suivi du
sur une hauteur plus grande (fig. 17). Ils présentent
comportement de l'ouvrage pendant la montée et la
un diamètre de 56 mm pour ceux placés dans le sable
descente progressive du plan d'eau.
de Fontainebleau et de 12 à 19 mm pour ceux placés
dans la fondation et le corps de digue (piézomètres Le programme de cette auscultation est présenté sur
type Casagrande) ; la figure 9.

34
Les principaux résultats de cette expérience sont les Cote NGF (m)
du plan d'eau
suivants : 150
149 •
— Le temps de réponse court des piézomètres montre 148 •

que les fluctuations du plan d'eau se répercutent rapi- 147 •

146
dement à la nappe du sable de Fontainebleau et
145
au corps de digue. A titre d'exemple, la figure 10
144
représente, durant la montée du plan d'eau, les varia- Niveau — 143 •
tions des piézomètres 12 et 11532 situés respective- normal + Appareils
142
ment dans l'ouvrage et en rive du bassin. Les Vinchon
141 •
figures 11 et 12 illustrent également ce phénomène en 6 avril 1979

•y
Bassin 140
présentant les variations de la ligne « surface libre » vide 139,5 10 20 30 40 50 I I 6 n
Temps (j)
reliant les niveaux d'eau relevés dans les deux profils I I I I I I I I I
de mesures. Cette réponse hydraulique rapide s'expli- Mesures piézométriques
que en grande partie par la présence de la nappe des
Examen visuel journalier
sables de Fontainebleau qui baigne en permanence le
fond et les rives du bassin. Les niveaux n'étaient Fig. 9. — Programme de l'essai de mise en eau aux plus hautes eaux.
toutefois pas complètement stabilisés à la fin de la
montée du plan d'eau, et on a constaté une superposi-
tion des écoulements lors de la descente (niveaux de Montée du plan
25 % maximum plus élevés qu'à la montée pour un d'eau et des
piézomètres (m
même niveau du plan d'eau).

— Le système drainant fonctionne correctement, bien

qu'on ait constaté des niveaux d'eau anormaux dans
certains piézomètres de contrôle, en particulier en aval
du système drainant (piézomètres 4, 7 et 11 — fig. 1).
L'analyse montre que ces niveaux n'ont aucun rapport
avec le plan d'eau car ils ne réagissent pas à ses
fluctuations. Il s'agit de zones saturées (niveaux
« perchés ») dont la présence est liée à la pluviométrie
et à l'irrégularité de la perméabilité. Bien qu'elle ne
soit pas complément stabilisée, on constate (fig. 12)
qu'aux plus hautes eaux la « surface libre » ne
recoupe pas la cheminée drainante. Le débit des eaux
de drainage n'a pu être mesuré correctement en raison
d'un apport d'eau parasite en provenance de l'évacua- Temps (j)
teur de crue.
Fig. 1 0 . - Variations des niveaux de deux piézomètres
lors de la montée d u plan d'eau.
— L'examen visuel quotidien a permis de s'assurer
du bon comportement de la digue et de ses abords
au cours de l'essai, mises à part quelques défaillances
mineures (érosion rapide du sol par les vagues dans
la partie du talus non protégé par les enrochements, bassin est resté vide à la cote 139,50 ou plein à une
petites fuites diverses au droit des ouvrages en béton, cote d'environ 145 pendant plusieurs mois), ce qui
etc.). crée vraisemblablement des déformations supplémen-
taires. Par rapport au « tassement » du remblai, l'am-
plitude des déplacements supplémentaires « d'origine
— Déplacements :
• le déplacement vertical de la voûte de la galerie hydraulique » reste cependant faible. E n effet pour la
sous remblai est resté faible (au maximum environ période de mesures (environ deux années) on peut
2 mm); estimer que le déplacement vertical ou horizontal
cumulé est de l'ordre de 1 cm au maximum, alors que
• la mise en eau provoque un déplacement de la le « tassement » maximum atteint 5 cm.
zone amont de la digue, c'est-à-dire du remblai situé
approximativement sous la crête et le parement amont O n a recherché l'origine des déformations autres que
(fig. 13 et 14). Ce déplacement est dirigé à la fois vers celles provoquées par la mise en eau. Bien que la
le bas et vers l'amont. Il s'additionne à un déformation verticale soit prépondérante, ce qui
« tassement » en cours qui se manifeste seul dans la conduit à utiliser le terme de « tassement », les mesu-
partie aval de l'ouvrage, là où le remblai n'est pas res mettent également en évidence des déformations
concerné par la mise en eau (remblai routier). Comme horizontales en particulier en bordure de talus. Le
le montre la figure 15, l'effet de la mise en eau, qui phénomène n'est pas imputable à un tassement rési-
correspond en pratique à une montée temporaire du duel de la couche d'alluvions car il se produit égale-
plan d'eau de la cote N G F 141 à la cote N G F 149,50 ment sur les rives, là où le remblai repose directement
avant sa mise au niveau normal à la cote N G F 143, sur les sols en place (limon, argile à meulière, sables).
diminue rapidement avec le temps mais est encore Il s'agit donc pour l'essentiel d'un tassement du rem-
sensible deux ans après. Il faut remarquer toutefois blai sous son propre poids. Sachant que la mesure
que pendant cette période les variations du plan d'eau topographique initiale (4 avril 1979) a été faite dix-
autour de la cote N G F 143 ont été importantes (le sept mois après la fin de la construction (qui a duré

35
36
 Échelle des profils ? 15m
Échelle des déplacements ? 1B m m

,. 0 5 mm (sur b o r n e n ° 14)
Echelle des déplacements i l
(sauf sur borne n° 14) Déplacement mesuré au niveau
des plus hautes eaux

Fig. 14. — Déplacements dans les profils nord et sud durant la mise en eau.

elle-même deux ans) et que les tassements mesurés Anisotropie du limon constituant la digue
décroissent d'une façon plus ou moins linéaire en
fonction du logarithme du temps, on a estimé le
Principe de l'étude
tassement total depuis la fin de la construction, c'est-
à-dire sur une période de 3,5 années. Fonction de la Il consiste à tenter de retrouver les caractéristiques
hauteur du remblai (17 m environ au maximum), la de l'écoulement réel (pressions interstitielles, surface
valeur maximale serait de l'ordre de 15 cm. Ce phéno- libre) par calcul à l'aide d'un modèle dans lequel on
mène est lié à la teneur en eau de mise en œuvre en peut faire varier les paramètres et en particulier le
moyenne un peu trop élevée du matériau. A proximité rapport k /k„ du corps de digue. Utilisant le modèle
h

des talus, ce phénomène s'accompagne vraisemblable- théorique construit par analogie électrique et établi
ment de déformations provoquées par la mobilisation lors de la vérification de la stabilité de l'ouvrage, une
du cisaillement. E n conclusion, on peut dire que le première tentative d'analyse de l'anisotropie avait été
comportement mécanique traduit beaucoup plus faite lors de la première mise en eau en supposant un
l'aspect remblai de l'ouvrage que son aspect digue régime permanent établi aux plus hautes eaux. L'étude
ou ouvrage hydraulique. n'avait pu aboutir, principalement en raison du nom-

37
 Perméabilité
K • Piézomètre type 1
• H
• Piézomètre type 2
20 19 15 17 ' 14 • 10
- • — • — • • — • • • Piézomètre type 3
® 16 g©
2

s *•
: 35
6

, Trace du
plan d'eau
(NGF 143)

• Piézomètre
• Perméabilité
Profil sud © Dlagraphle
/
10 m
1 Fig. 17. — Coupe de l'expérimentation complémentaire.

Fig. 16. — Plan de situation de l'expérimentation complémentaire.

obtenue par essai en laboratoire. L a perméabilité hori-

zontale k de la partie supérieure de la couche d'allu-
h

bre trop restreint des points de mesure de la pression vions (colluvions fines sur la fig. 18) est plus forte
interstitielle dans le limon, les possibilités de l'analogie ( 5 . 1 0 " m/s)
6
que celle du limon susjacent
électrique étant par ailleurs limitées. L'étude a été (5.10 7
m/s) et comparable à la valeur mesurée par
reprise à la retenue normale ( N G F 143) après qu'une essai Lefranc dans la partie supérieure du sable de
expérimentation complémentaire ait été faite en place. Fontainebleau. Cela est en accord avec les variations
de texture des matériaux mises en évidence sur le
Expérimentation complémentaire graphique de la radioactivité naturelle (fig. 19) (celle-
ci est en effet en relation directe avec le pourcentage
L'expérience de la mise en eau ayant montré que les d'argile). Cm peut donc penser que la partie inférieure
caractéristiques de l'écoulement étaient quasi identi- de la couche d'alluvions (colluvions grossières sur la
ques dans les deux profils de mesure, seul le profil figure 19) est encore plus perméable, et on conçoit
nord a été utilisé par la suite. L'expérimentation com- que du point de vue de la perméabilité, l'ensemble de
plémentaire a comporté une étude de la perméabilité, la couche de fondation, au moins localement, puisse
la pose de piézomètres complémentaires et l'exécution s'apparenter plus au sable sous-jacent qu'au limon du
de sondages avec diagraphies (fig. 16 et 17) de radio- corps de digue.
activité naturelle, gamma-gamma (variation de la
densité) et neutron-neutron (variation de la teneur Piézomètres complémentaires
en eau).
Il s'agit de piézomètres de petit diamètre, crépines sur
Étude de la perméabilité 0,5 m et numérotés de 1 à 7 sur la figure 17. Ce
complément porte à 20 le nombre des piézomètres
Elle a consisté à effectuer, à partir de la risberme, un placés en amont du dispositif de drainage.
forage avec essais de perméabilité tous les mètres au
pressio-perméamètre Ménard, destiné à fournir un Sondages avec diagraphies
profil des perméabilités sur les trois couches de sol
concernées (limon, alluvions, sable de Fontainebleau). Deux sondages ont eu pour but la détermination de
Des difficultés d'exécution (éboulement des parois) la position de la surface libre, et un troisième a permis
n'ont pas permis d'atteindre la base du forage mais de compléter la reconnaissance géologique et hydrauli-
l'interprétation a pu être menée à bien grâce à l'appui que (voir précédemment).
des résultats d'un sondage diagraphie situé à proxi-
mité immédiate ( a l m) du forage (fig. 18 et 19).
ÉLABORATION DU PROFIL D'ÉTUDE
L a perméabilité horizontale du limon du corps de
digue varie entre l , 4 . 1 0 ~ m / s et 2 , 4 . 1 0 ~ m / s ce
9 6 A PARTIR DES RÉSULTATS EXPÉRIMENTAUX
qui peut surprendre si l'on songe à la constance de
sa composition granulométrique. E n fait, on peut Pour être sûr que les relevés piézométriques correspon-
penser que les écarts sont liés à la variation du volume dent à un niveau de retenue constant à la cote
du gaz occlus, elle-même traduite par la variation de N G F 143, il était nécessaire de s'assurer que celui-ci
la teneur en eau. Il semble que, plus la teneur en eau (enregistré en continu au moyen d'un limnigraphe à
est faible, donc plus il y a de gaz occlus, plus la flotteur) était effectivement stabilisé à cette cote sur
perméabilité est faible. une durée suffisante antérieurement à la mesure. Ce
n'est qu'à la mesure du 25 mai 1982, c'est-à-dire après
L a perméabilité horizontale moyenne ou perméabilité un an et six mois environ de suivi de l'ouvrage au
équivalente, est égale à 5. 10~ m/s, valeur cinquante 7
niveau théoriquement permanent (période pendant
fois plus grande que la plus grande valeur de k v laquelle douze séries de relevés piézométriques ont été

38
 1

RAN. WIDCO
Log de synthèse

faites), qu'on a constaté que le niveau du plan d'eau fondation (couche de fondation proprement dite et
était resté constant (à +10 cm près) depuis environ sable de Fontainebleau sous-jacent) et d'autre part le
trois mois et que cette mesure pouvait être utilisée degré d'anisotropie de chaque couche. Compte tenu
dans l'interprétation. des résultats rappelés ci-dessus dans le paragraphe sur
l'étude de la perméabilité, cinq profils d'étude ont été
E n ce qui concerne la perméabilité, il s'agissait de retenus. Sur la figure 20, on a représenté d'une façon
fixer, dans le profil d'étude, d'une part les rapports schématique les caractéristiques hydrauliques des cinq
de perméabilité entre le corps de digue et le sol de profils.

//////

O © O ©
© © © Rapport d'anisotropie k /k h t

Repère
de la zone A B C D
du profil

Corps de digue Couche 10 premiers Reste de la couche

Nature du sol
en limon de fondation mètres de sable de sable

Définition de la zone en fonction

de sa perméabilité verticale k
v
'/////,
Fig. 20. — Caractéristiques hydrauliques
Perméabilité k v 1 10 100 des cinq profils théoriques.

39
 Fig. 21 Niveau p e r m a n e n t Niveau piézo 8 N Niveau piézo 2 N
Conditions aux limites. N G F 143 136 N G F 135 N G F

200 m 100 m

.2
3
Niveau p u i t s .
134 N G F

-Imperni éable-

Construction du modèle théorique

L a détermination des pressions interstitielles ou des fils théoriques aux pressions interstitielles mesurées le
charges aux points correspondant à la prise de pres- 25 juin 1982. Les piézomètres repérés par le signe G
sion des piézomètres a été faite en utilisant le pro- sont crépines sur une faible hauteur (0,5 et 1 m), et
gramme F L O T . C'est un programme de calcul par la les piézomètres repérés par le signe 0 sont crépines
méthode des éléments finis des problèmes d'écoule- sur une hauteur plus importante. Les pressions sont
ment plan ou à symétrie de révolution en régime en kPa.
permanent avec détermination de la surface libre.
Les conditions aux limites sont représentées sur la De façon à rendre la comparaison plus lisible, on
figure 21. a déterminé en pourcentage, pour chaque pression

p calculée \ et
Résultats interstitielle calculée, son rapport ^ son

mesurée /
Comparaison des pressions interstitielles
p mesurée — p calculée \
écart avec la pression
Le tableau II et la figure 22 permettent de comparer p mesurée
les pressions interstitielles calculées dans les cinq pro- mesurée.

TABLEAU II
Comparaison des pressions interstitielles calculées et mesurées.

Piézomètres

Pression 1 2 3 4 5
Numéro Type mesurée
(kPa)

20 0 53,9 25,36 26,35 28,71 32,32 32,61

19 0 51,6 29,42 26,55 30,01 37,82 37,00

18 0 37,7 37,59 25,59 31,26 49,03 47,47

16 0 36,3 36,60 25,77 31,20 47,76 46,19

15 0 34,7 33,67 24,18 29,44 44,83 43,65 •

17 0 83,2 83,05 59,59 67,47 94,72 94,01

2 e 0,7 0 9,02 10,27 4,85 9,21

3 e 11,7 3,59 13,14 15,15 10,29 15,13

4 9 26,0 13,61 16,88 19,83 22,02 24,18

5 9 26,7 23,10 20,52 24,36 33,27 32,92

6 9 43,1 38,52 23,79 29,92 49,20 48,52

14 9 35,7 32,83 22,45 27,52 42,94 42,39

13 0 36,0 3422 22,64 27,69 43,54 43,13

12 9 43,1 31,88 21,95 25,93 39,35 39,10

7 9 18,0 10,64 8,49 9,79 14,21 14,12

10 9 27,3 28,44 22,15 23,90 32,63 32,54

9 9 51,8 54,51 42,61 47,23 58,21 58,10

40
 1

1 e r
Cas Hypothèses de perméabilité

• Tous les sols sont isotropes. Limon et alluvions *A = 5 . 1 0 " 7

khlkv
• Les alluvions ont la même perméabilité que le ^=5.10-'
limon. Sable de Fontainebleau sur */, = 5 . 1 0 - 6

khlkv kh sable (10 m)/*/, limon = 10

10 m *v=5.10" 6

Sable de Fontainebleau reste */, = 5 . 1 0 " 5

*» = 5.10" W*v=1 kh sableMA limon = 100

de la form. 6

2 Cas
e Hypothèses de perméabilité

• Anisotropic de 50 dans le limon et les alluvions. Limon et alluvions */, = 5.10~ 7

khlkv = 50
• Le sable de Fontainebleau est isotrope. *„=1.10" 8

• Le limon et les alluvions ont la même perméabilité. Sable de Fontainebleau sur *A = 5.10" 6
kt, sable (10 m)/*/, limon = 10
10 premiers mètres *,=5.10- 6
=1
Sable de Fontainebleau reste */,= 5 . 1 0 '5
kh sable/Ar/j limon = 100
kh/kv=1
de la form. *v=5.10~ 5

3 Cas e Hypothèses de perméabilité

Anisotropic de 50 dans le limon et les alluvions et Limon et alluvions */, = 5.10" 7

khi k* = 50
de 9 dans le sable de Fontainebleau. Av=1 1 0 " 8

> Les alluvions ont la même perméabilité que le Sable de Fontainebleau sur */, = 5.10" 6

khlkv = 9 kh sable (10 m)lkh limon = 10

limon. 10 m *v=5/9.10" 6

Sable de Fontainebleau reste */,= 5 . 1 0 - 6

kh/K =9 kh sable reste /kh limon = 100

de la form. *»=5/9.10 -5

Fig. 22.— Comparaison des pressions interstitielles mesurées et calculées pour chaque modèle théorique.

41
 4 e
Cas Hypothèses de perméabilité

• Tous les sols sont isotropes. Limon *A = 5.10 7

• Les alluvions ont la même perméabilité que les 10 *„=5.1CT 7

premiers mètres du sable de Fontainebleau. Alluvions et sable de Fontai- A/i = 5 . 1 0 - 6

=1 kh alluvions sable/fa limon = 10

nebleau sur 10 m *v=5.10~ 6

Sable de Fontainebleau reste k = &. 10~

n
5

ki,/kv = 1 kfi sable/An limon = 100

de la form. *»=5.1CT 5

5 e
Cas Hypothèses de perméabilité

i Anisotropic de 10 dans le limon Limon */, = 5 . 1 0 " 7

= 10
*„=5.10- 8

• Alluvions et sable de Fontainebleau Isotropes.

Alluvions et sable de Fontai- */, = 5 . 1 0 " 6

nebleau sur 10 m *„=5.10^ 6

*/,/*„ = 1 kh alluvions sable/fa limon = 10
Sable de Fontainebleau reste * = 5.10" 5

de la form.
n

*v=5.10~ 5 khlkv=-\ kh sableMn limon = 100

Fig. 22. — Comparaison des pressions interstitielles mesurées et calculées pour chaque modèle théorique.

On obtient les valeurs moyennes suivantes

TABLEAU III (tableau III) dans chacun des cinq cas étudiés.

Piézomètres Numéro du profil Bien que la comparaison des valeurs moyennes du

pris en rapport valeur calculée/valeur mesurée puisse être
compte 1 2 3 4 5 jugée peu représentative en raison d'un phénomène de
compensation entre valeurs supérieures et inférieures
Valeur moyenne
o+e 80,5 67 78,5 106 107,5 à 100 %, elle fait toutefois apparaître nettement le
du rapport (%)
modèle 2 comme très éloigné de la réalité et des quatre
0+ 8 20,5 34,2 25 21,0 20,5 autres modèles. Si l'on raisonne sur les écarts, où le
Valeur moyenne
de l'écart (%) phénomène de compensation ne joue pas, un modèle
e 25 32,5 25 16,5 17
supplémentaire est éliminé, il s'agit du modèle n° 3.
Les numéros 1, 4 et 5 semblent très proches les uns
R E M A R Q U E : les piézomètres 1 et 2 n'ont pas été pris en
des autres. Si, enfin, on ne s'appuie que sur les
compte en raison de leur très faible pénétration dans piézomètres G les plus précis, seuls les profils 4 et 5
l'écoulement. Ils ont été utilisés comme repère de la resteraient en compétition pour la meilleure représen-
surface libre au stade ultérieur de l'analyse. tation de l'écoulement réel.

42
Pour approfondir l'analyse on ajoutera aux points de matériau. Dans le cas de la digue de Maurepas mise
comparaison des pressions, ci-dessus, ceux particuliers en œuvre dans les mêmes conditions, on a obtenu
formant la surface libre, la surface libre réelle ayant un résultat identique même avec un matériau à la
été repérée en quatre points particuliers. granularité plus étendue (cf. article de M M .
Josseaume, Mahdavi et Jardin dans ce bulletin de
liaison, p. 23 à 28).
Comparaison des surfaces libres

Sur la figure 23 sont représentées les cinq surfaces

libres calculées et la position des points où la surface CONCLUSION
libre réelle a été repérée, ce qui correspond à des
piézomètres (1 et 2) ou à des sondages avec diagraphie
neutron-neutron (a et b) qui pénètrent très peu sous
la surface libre (de l'ordre de 1 m au maximum). Cet article conduit à la réflexion suivante
lorsqu'on songe à l'importance, pour la conception
O n constate que c'est la surface libre du modèle 4 de l'ouvrage, des informations apportées par les mesu-
qui, des deux modèles 4 et 5, se rapproche le plus res que l'on n'a pu faire qu'une fois Fouvrage construit,
des points réels, cela étant d'autant plus net que la on conçoit qu'il soit toujours difficile de prévoir, au
comparaison se fait dans la zone (entre parement stade du projet, le comportement d'un barrage en
amont et risberme) où les positions relatives des surfa- terre, et cela en dépit d'une reconnaissance géotechni-
ces libres théoriques sont les plus distinctes. Ce résul- que approfondie.
tat est confirmé si, tenant compte du fait que dans
ce modèle 4 les équipotentielles sont proches de la C'est ainsi qu'il était difficile de prévoir que le tasse-
verticale, on considère que les niveaux des piézomètres ment propre d'un remblai en limon de 17 m de hau-
sont ceux de la surface libre. Ce raisonnement étant teur, construit lentement (deux ans) et dans les règles
également valable pour le modèle 1, dont d'ailleurs la de l'art (quoique avec une teneur en eau moyenne un
surface libre est également proche des points réels 1 peu trop élevée malgré les efforts faits pour assécher le
et a, on conclut que la surface libre réelle se situe en matériau au terrassement) n'aurait pas cessé 3,5
moyenne dans un « fuseau » approximativement années après la fin de la construction, et atteindrait
limité par les surfaces libres des modèles 1 et 4. Cette à cette date environ 15 cm, ou que l'amplitude des
vraisemblable ondulation de la surface libre réelle est déplacements provoqués par la mise en eau aux P H E
à notre avis une conséquence de l'hétérogénéité de et orientés à la fois vers le bas et vers l'amont serait
nature et donc de perméabilité (cf. texte précédent sur aussi faible (de l'ordre du centimètre).
ces deux points) de la couche de fondation dans le
sens vertical mais également, et surtout, dans le sens Mais c'est surtout dans le domaine hydraulique que
longitudinal c'est-à-dire suivant le profil en travers de les prévisions étaient les plus délicates à faire. Par
l'ouvrage. exemple déterminer la perméabilité moyenne d'une
couche aussi hétérogène que la couche de fondation
E n conclusion, la prise en compte de l'ensemble des était aussi malaisé que de cerner le contraste de
observations sur le réseau d'écoulement réel montre perméabilité entre le futur corps de digue et cette
qu'il est proche de celui correspondant aux modèles 1 même couche. Les constatations montrent que ce rap-
et 4. O n peut donc dire que le limon du corps de port se situe entre 1 et 10. Cependant, c'est bien dans
digue a un comportement hydraulique isotrope ou l'appréciation du comportement hydraulique du corps
très proche de l'isotropie. de digue lui-même que résidait la plus grande diffi-
culté. Sur la base des valeurs de la perméabilité mesu-
A notre avis, ce résultat est essentiellement lié au rée sur échantillons de limon compacté, toutes inférieu-
mode de mise en œuvre, le facteur principal étant la res à 1 0 " m/s, on ne pouvait guère imaginer que la
8

faible épaisseur des couches successives (0,30 m au perméabilité globale réelle du corps de digue soit
maximum) liée au compactage à l'aide d'un rouleau beaucoup plus forte telle que les fluctuations du plan
à pieds dameurs, technique qui présente l'intérêt, sous d'eau se transmettent ainsi rapidement à l'ouvrage.
réserve d'une teneur en eau favorable (proche de D'ailleurs, on note aussi de la part du sable de Fontai-
l'optimum), d'assurer une grande homogénéisation du nebleau un temps de réponse court dans les rives du

43
bassin. Ce comportement moyen n'est cependant pas Après avoir présenté les principaux résultats des cons-
le fait d'un remblai hydrauliquement homogène, mal- tatations, il est intéressant, pour conclure, de revenir
gré la granularité particulièrement uniforme et cons- sur la réflexion initiale en se posant la question
tante du limon qui le constitue : le profil de mesures suivante : est-il possible de tirer des enseignements de
effectuées tous les mètres montre que k varie entre portée générale des constatations faites dans ce cas
1 0 m / s et 1 0 m / s (fe équivalent = 5 . 1 0
_ 9 _ 6
m/s).
_ 7
particulier, de telle sorte qu'ils puissent être pris en
compte au stade de la conception des ouvrages de ce
O n aurait pu penser qu'une telle plage de variation type ? A notre avis, l'incidence des conditions de mise
de k refléterait le caractère éminemment anisotrope en œuvre sur le comportement hydraulique du corps
du corps de digue avec la présence de niveaux d'écoule- de digue est le fait de caractère général le plus impor-
ment préférentiel plus ou moins horizontaux, justi- tant qu'apportent ces constatations (il est confirmé
fiant ainsi la prudence du concepteur qui a prévu une dans le cas de la digue de Maurepas) : en mettant en
tranchée drainante oblique pour pallier ce risque. Les œuvre un matériau du type limon par couche de faible
constatations montrent que de telles circulations épaisseur (0,30 m maximum) à Vaide d'un rouleau à
n'existent pas et qu'on est en présence d'un réseau pieds dameurs dans des conditions de compactage
hydraulique continu et homogène correspondant à un proches de l'Optimum Proctor, on obtient à travers
comportement globalement isotrope du corps de l'ouvrage un réseau d'écoulement continu et homo-
digue, résultat qui sur l'ensemble des constatations gène correspondant à un comportement hydraulique
est le plus intéressant et le plus difficilement prévisible. global isotrope du corps de digue.

44
 Digues
et barrages en terre
7L de faible hauteur

Construction et comportement
d'une digue en terre à haute altitude
Michel LAGRANGE
Technicien

Jacques PERRIN*
Ingénieur
Laboratoire régional de Lyon

Henri JOSSEAUME
Attaché de recherche
Section des ouvrages en terre
Laboratoire central des Ponts et Chaussées

Vue générale de la retenue.- En été.

RÉSUMÉ

L'aménagement d'un bassin de stockage en vue

de l'amélioration de l'alimentation en eau d'une
station de sports d'hiver a nécessité la construc-
tion d'une digue de faible hauteur (9,50 m) à
haute altitude.
Après avoir présenté les principales caractéris-
tiques géotechniques du site qui ont conduit à
la construction d'une digue homogène et après
avoir décrit cette dernière, les auteurs rendent
compte des problèmes qui se sont posés
pendant et après les travaux • problèmes liés
à la mise en place d'un géotextile entre le corps
de digue et les enrochements de protection,
fuites importantes lors de la mise en eau du
bassin qui ont nécessité l'injection de la fon- En hiver.
dation rocheuse.
L'article rend également compte des mesures
de pression interstitielle effectuées dans la Afin d'accumuler les eaux nécessaires à l'amélioration de la
digue et de leur interprétation. distribution en eau potable des stations de sports d'hiver
MOTS CLÉS : Barrage en terre - Construction de L a Plagne et de ses satellites, une retenue collinaire en
(exécution] • Comportement - Montagne • terre, haute de 9,50 m, a été édifiée au cours de l'été de
Non tissé - Protection - Roche - Filtre - Gravier - 1978 et 1979 dans la dépression des Blanchets, à 2 340 m
Sable - Étanchéité - Injection (mater.) - Coulis-
Bentonite - Ciment - Mesure • Pression intersti- d'altitude (fig. 1). Pendant l'été, il est prévu de stocker
tielle - Interprétation - /Digue - Étanchement - 150000 m qui seront utilisés durant la saison hivernale.
3

Enrochement.

Au Laboratoire de Lyon à la date de la rédaction de l'article.

13
B u l l . Maison L a b o P. et C h . - 1 3 2 - j u M . - a o û t 1 9 8 4 - R é f . 2 8 9 5
L'analyse hydrologique, le choix du site et le dimen-
sionnement des ouvrages annexes ont été réalisés par 0,10 CO Terre végétale

la Direction départementale de l'Équipement de la

Savoie. L a conception et le suivi de construction du Sols organiques
très humifies
corps de digue ont été assurés par le Laboratoire
0,60
de Lyon, dont les interventions sont décrites dans
l'article.
Argile glaciaire bleue
graveleuse compacte

DONNÉES GÉOLOGIQUES
o
ET GÉOTECHNIQUES
û

Grave sableuse propre

L a dépression des Blanchets est consécutive à un ' <$

0

Au droit du verrou 2 m
surcreusement glaciaire peu à peu remblayé par les (2 m en cuvette)
apports torrentiels et glaciaires, puis par les alluvions
Substratum quartzitique
récentes superficielles organiques. Le substratum quar-
tzitique en constitue le fond géologique et le verrou
(fig- 2). Fig. 2 — Coupe géologique représentative.

Les sondages de reconnaissance sommaires (tarière à

main, mini-terrassement, prélèvement à la boîte car-
rée) et les observations visuelles ont permis : Corps de digue, éboulis à fines (0-80)

— de constater que la cuvette présentait une — perméabilité verticale K„ m o y 2 . 1 0 " m/s

8

imperméabilisation naturelle très importante du fait — poids volumique sec :

du placage continu d'argile glaciaire bleue, graveleuse y : 20,5 k N / m
d m o y
3

et compacte; sur la fraction 0-20

— teneur en eau :
— d'orienter en conséquence les extractions de maté-
riaux imperméables hors du site ennoyé, d'autant que w m o y : 6,5 7o(w OPN = 9 %)]
les graves sableuses sous-jacentes, de forte perméabi- — caractéristiques intrinsèques à la boîte de
lité, pouvaient provoquer des pertes éventuelles Casagrande sur la fraction 0/2,5 m m : <p' = 3 9 ° et
importantes ; c ; = 0,12.10 kPa.
o y
2

— de mettre en évidence les grandes possibilités gra-

Filtres et drains (0-100)
nulométriques qu'offrent les éboulis quartzitiques
dans différents couloirs et cônes à proximité du site. — perméabilité verticale K^, m o y : 3,5.10~ m/s 6

Ce qui conduit à envisager une digue homogène ne

— poids volumique sec : y : 20,8 k N / m 3

nécessitant pas d'apports très onéreux de matériaux

dm o y

en provenance de la vallée (tableau I et fig. 3); — teneur en eau : w m o y : 4,9%

— de vérifier la bonne qualité hydraulique et mécani- — caractéristiques intrinsèques à la boîte de

que de ces sols dans les conditions suivantes de mise Casagrande sur la fraction 0/2,5 mm : <p' : 40° et moy

en œuvre.

14
 100 100 Fig. 3
Granulante des sols
90 90 constitutifs de la digue.

80 80

1 70
1 70

S 6°
Ê
2 50 2 50

"O
03 •D

§ 30 § 30
3 D
O O

" 20 1
20

10 10

• m m 200 100 50 20 10 5 2 1 0,4 0,2 0,1 0 , 0 8 5 0 ( j m • mm 200 100 50

Diamètres équivalents (/jm)

Corps de digue, grave limoneuse avec beaucoup de sable. Enrochements 100-500.

TABLEAU I CONCEPTION

Origine du matériau h ES
Les possibilités en matériaux et l'imperméabilisation
Corps de digue Éboulis à « fines » suffisante (estimée par une reconnaissance visuelle)
(0-80) du col de Forcle et des
22 5 19 du socle quartzitique ont conduit le laboratoire à
Quillis (distance de proposer finalement une digue homogène fondée au
transport «s 1,7 km)
rocher, soigneusement décapé, sans injection a priori,
Éboulis « gravelo-sableux » avec un système de drainage interne par tapis mono-
Filtres et couche et une protection générale par enrochements.
drains des Quillis (distance de (') C) 30
(0-100) transport »0,5 km) Cette protection était rendue nécessaire par les condi-
tions atmosphériques difficiles liées à l'altitude (eau,
Enrochements Éboulis « crus » neige, gel) et par l'amélioration de la protection du
(100-500) du Roc du Diable talus aval pour le cas où des filets d'eau y resurgi-
(distance de raient. Les pentes finales amont et aval retenues ont
transport « 1,3 km)
été de arctg 1/2 après une vérification des stabilités à
pleine eau et en vidange rapide, dans les conditions
(') indéterminé
géométriques de construction indiquées sur les
figures 4 et 5.

C o t e f i n a l e d e service

a. Profil P.2 Fig. 4 — Constitution de la digue

et repérage des capteurs de mesure de
pression interstitielle.

i i F i l t r e d e t r a n s i t i o n et d r a i n aval
F" / /\ Enrochements
É b o u l i s q u a r t z i t i q u e s à fines
/ / / / , Enrochements
Nappe textile (Sodoca A S . 4 2 0 )
Matériaux fins i d e n t i q u e s R o c h e r d'assise
à corps de digue P1 : p i é z o m è t r e o u v e r t
Non-tissé (crépine sur t o u t e la hauteur)
'. : ; . Substratum quartzitique Pac 01 et P a c 11 : S o n d e d e pression
Appareillages
interstitielle à c o n î r e p r e s s i o n p n e u -
d e mesures
matique
L P C 21 : S o n d e d e p r e s s i o n interstitielle
à contrepression hydraulique

b. Profil P.3

15
CONSTRUCTION D E L A DIGUE

Généralités

L a digue a été construite directement sur le substra-

tum quartzitique après un décapage soigné à la
main (*) (fig. 6).

Lors des travaux préparatoires de décapage, une

source, apparue sous le talus amont, à l'Est (près du
profil P2), a été bétonnée pour permettre la continua-
tion des travaux (mise en place des conduites et des
remblais).

Les contrôles de compactage ont été exécutés par le

laboratoire avec une vérification à l'aide du gamma-
densimètre D R 30 pour les sols fins (filtres et corps
de digue) (tableau II). Fig. 6 — Corps de digue en construction.

Le remplacement du filtre par le non-tissé (Sodoca

type A S 420) n'a pu être mis en œuvre pour les raisons Problèmes liés à l'utilisation d'un géotextile
qui seront évoquées dans le paragraphe suivant.
E n 1978, il a été décidé, à la suite d'une proposition
Lors de la première mise en eau (septembre 1980), il
du maître d'ouvrage et pour des raisons de « facilité »
est apparu une fuite importante à l'aval de l'éperon
de mise en œuvre, de remplacer la couche de filtres
rive gauche de la digue ; le traitement de ce problème
de transition par un géotextile entre enrochements et
est décrit dans le paragraphe « colmatage des fuites ».
corps de digue.

Le maître d'ouvrage a proposé l'utilisation du non-

(') Sauf très localement, en pied de falaise rive droite sur environ
2
50 m , compte tenu de la finesse des sols reconnuslors de la mise tissé de marque « SODOCA » type A S 420
en place des conduites de vidange et d'exploitation. (370 g/m ) (matériau comparable au Bidim U 44).
2

16
T A B L E A U II Épaisseur
Type Nombre ^moy
des Yd moy
de de O-D
couches (kN/m ) 3

(m) compacteur passes (%)

Corps de digue (0-80) Pieds

(Classe de terrassement : B 5) 0,40 dameurs 20 20,6 6,6
(PD1)

Filtres (0-100) Pieds

(Classe de terrassement : D 3) 0,40 dameurs 20 20,8 5
(PD1)

Enrochements (100-500) Vibrant

(Classe de terrassement : D4) 0,70 tracté 10
(V3)

A la fin des travaux réalisés lors de l'été 1978, les Dans tous les cas, les enrochements furent compactés
parties de la digue construites furent protégées de à l'aide d'un rouleau vibrant tracté type V 3.
l'hiver par le non-tissé, lui-même recouvert de polyes-
ter (immobilisée par enrochements et bastaings). Le sol support était lui-même constitué par une cou-
che de matériaux fins identiques à ceux du corps
A la reprise du chantier en 1979, on a constaté que de digue et présentait les mêmes caractéristiques de
le non-tissé présentait des déchirures et des perfora- compactage.
tions à l'endroit où les enrochements avaient été mis
en œuvre. A la suite des résultats obtenus, la solution
« non-tissé » a été abandonnée pour ne pas risquer
Après ces constatations sommaires, il fut décidé de l'entraînement des fines et les graves désordres dus
réaliser des planches d'essai pour définir une méthode aux déchirements prévisibles (quel que soit le mode
de construction qui n'altère en aucun cas le non-tissé. de mise en œuvre). L a méthode initialement envisagée
L a figure 7 résume ces différentes planches. du filtre gravelo-sableux a été reprise.

Principe de mise en œuvre

0 Perforations et déchirures
du non-tissé sur toute la
(Figuré du corps de digue) surface testée.
Environ 50 t r o u s / m
2

Enrochements
compactés (perforations et déchirures)
N.B - La chute d'un bloc
(p 200 d'une hauteur de 1,50 m
provoque la perforation du
Sol support compacté Sodoca

Poinçonnement du Sodoca par

Enrochements quelques blocs d'enrochements
compactés et principalement par les maté-
40 crr
riaux filtrés.
Environ 10 t r o u s / m
2

20 cm . , . w >-> ,r~_ ^ O ^ O & ' ^ ^ ' - O f a O C p o ^ Filtres non compacté (perforations et déchirures)
• Sodoca légèrement tendu
A W W Sol support compacté

©
Les enrochements arrivent
Enrochement* à traverser les quartzites
compactés à fines et à poinçonner le
géotextile (fig. 8).

. Quartzites à fines Environ 5 t r o u s / m

2

non compactées (perforations et déchirures)

T ^ ^ ^ ^ ^ ^ Sodoca sans tension
Sol support compacté

©
Le compactage des enroche-
Enrochements ments entraîne quelques blocs
compactés (100/200) (se trouvant à
l'intérieur des quartzites à fines}
qui perforent le non-tissé.
Quartzites à fines Environ 2 t r o u s / m
2

compactées et triées (une déchirure et une perforation).

Sodoca plissé
' Sol support compacté Fig. 7 — Recherche d'un mode de mise en
œuvre du géotextile.

17
 MESURES D EPRESSION INTERSTITIELLE
E T INTERPRÉTATION

Équipement de la digue en piézomètres

Des mesures piézométriques ont été effectuées dans

deux profils transversaux de la digue situés dans la
zone où celle-ci atteint sa hauteur maximale. Ces
mesures avaient pour but, d'une part, de contrôler les
conditions de fonctionnement de l'ouvrage (notam-
ment l'efficacité du tapis drainant), d'autre part, d'étu-
dier l'anisotropie de perméabilité de la digue, dans le
cadre des recherches sur les digues de faible hauteur,
effectuées par les laboratoires des Ponts et Chaussées.

Les profils de mesures [profils 2 et 3 espacés de 10 m

Fig. 8 — Perforations d u géotextile à la mise en œuvre.
(fig. 5)] ont été équipés d'appareillages de plusieurs
types : cinq piézomètres pneumatiques P A C et cinq
piézomètres L P C à double tubulure ont été mis en
place dans le profil 2 (fig. 4 a) alors que trois piézomè-
Colmatage des fuites tres P A C et quatre tubes de 26 mm de diamètre
intérieur ont été posés dans le profil 3 (fig. 4 b).
Le niveau d'eau dans la retenue étant à la cote
Les piézomètres L P C n'ayant pas fonctionné de façon
2346,40 m ( e « 6 , 4 0 m d'eau), des fuites d'un débit
satisfaisante, seules les mesures ponctuelles effectuées
important ont été observées (15 1/s mesuré en (C)
au moyen des piézomètres P A C ont été considérées
dans le thalweg, à l'aval de la digue).
comme fiables. D e même, en raison de l'étanchéité
Le suivi « à l'oreille » a permis au maître d'œuvre et défectueuse du bouchon d'argile réalisé à la partie
au laboratoire de localiser la fuite sous l'enrochement supérieure du forage recevant le tube crépine PI, les
en rive gauche (secteur A ) . relevés de niveau d'eau effectués dans celui-ci n'ont
pas été retenus.
U n e reconnaissance sommaire, à la pelle mécanique
effectuée dans le secteur A jusqu'à proximité du socle
Étude des conditions d'écoulement
rocheux (sans l'atteindre pour des raisons de sécurité),
n'a pas mis en évidence l'entraînement de fines. Les
essais à la fluorescéine n'ont pas permis de localiser Les mesures piézométriques, effectuées à partir de la
une entrée d'eau éventuelle à l'amont. Le débit impor- première mise en eau, montrent que les conditions
tant et ponctuel de la fuite et le fait que le rocher d'écoulement sont très différentes dans les profils 2 et
décapé présentait une fracturation très superficielle 3, bien que ceux-ci soient voisins. Dans le profil 3
importante ont conduit le laboratoire à proposer une (fig. 10 et 11), le comportement hydraulique de la
imperméabilisation par injection au coulis de bento- digue ne présente pas d'anomalies : sauf en période
nite ciment (c/e de 0,3 à 1) de la surface du rocher de vidange du réservoir, la charge hydraulique est
dans l'axe de la digue (fig. 5 et 9). constamment inférieure à la charge dans la retenue
et décroît d'amont en aval pour devenir très faible au
Les fuites ont été réduites de 15 à 3 1/s (sous 7,60 m voisinage de l'extrémité amont du tapis drainant. E n
d'eau dans la retenue). Les fuites résiduelles provien- revanche, dans le profil 2 (fig. 12), la charge hydrauli-
nent du secteur rive droite. L a source rencontrée lors que mesurée au niveau de certains piézomètres est
du décapage initial alimente probablement pour partie supérieure à la charge dans la retenue : en particulier,
cette fuite résiduelle. à la fin de la première mise en eau, seul le niveau
piézométrique correspondant au piézomètre 3 se situe
nettement au-dessous de la surface du plan d'eau.

Les valeurs élevées de la pression interstitielle mesu-

rées dans le profil 2 ne peuvent, a priori, s'expliquer
Tube d'injection (50/60) dans le
coulis de gaine (vingt.sept forages
que si l'on se place dans l'une des hypothèses
distants de 1 m dans l'axe de la digue) suivantes :

— soit, de fortes pressions interstitielles se sont déve-

Corps
de digue loppées dans la digue pendant sa construction et leur
vitesse de dissipation a été suffisamment faible pour
que, à la fin de la première mise en eau, la charge
a) injection du contact corps de digue
rocher (à partir d'une manchette d'extrémité)

Substratum
quartzitique
b) injection directe au rocher dans N.B. La durée des travaux d'injection a été de trois semaines (y
l'avant-trou {<p 48 mm)
compris amenée et repli) pour un coût en 1982 de 335 kF TTC)
3
(longueur de forage «225 m, volume de coulis « 3 3 m ). Le coût
Fig. 9 — Principe d'imperméabilisation par injection. de la construction de la digue a été de 1800 kF TTC (1980).

18
 i

Charge hydraulique Im) Tubes crépines : K, L

Piézomètres P A C : J, M , N hydraulique dans certaines parties du profil soit supé-
Charge dans la retenue
rieure à la charge dans la retenue ;

— soit, l'écoulement à travers le profil est alimenté,

non seulement par la retenue, mais aussi par une
source indépendante de la retenue, située sous la base
de la digue.

Les pressions interstitielles mesurées en fin de cons-

truction dans chaque piézomètre du profil 2 (elles
sont comprises entre 16 et 50 kPa) étant inférieures à
celles mesurées en fin de première mise en eau (c'est-à-
dire plus de dix mois après la fin de la construction)
ou du même ordre, la première hypothèse ne peut
être retenue. E n revanche, la seconde est cohérente
1980 1981 1982 avec les observations faites au début des travaux. Il
Fig. 10 — Variations de la charge hydraulique dans le profil 3 est, en effet, extrêmement probable que la source,
entre le début de la première mise en eau (juillet 1980) et la fin de 1982. découverte à l'emplacement de la conduite de vidange
et qui est vraisemblablement alimentée par des eaux
provenant du massif constituant la rive droite, per-
turbe les conditions hydrauliques dans les profils voi-
sins, notamment dans le profil 2 situé à proximité
immédiate de cette conduite (à environ 4 m).

Influence des injections sur la pression interstitielle

Bien qu'elles aient été peu nombreuses et qu'elles
aient intéressé essentiellement le profil 3, les mesures
effectuées après la fin des travaux d'injection mettent
en évidence l'effet bénéfique de ces derniers sur le
comportement hydraulique de l'ouvrage à l'extérieur
de la zone injectée (comprise entre les profils 4 et 7,
fig. 5). Les mesures réalisées dans le profil 3 montrent,
en effet, que les pressions interstitielles et les niveaux
d'eau dans les tubes se sont abaissés pendant les
deux mois suivant les travaux d'injection, alors que le
niveau d'eau s'élevait dans la retenue (fig. 10). Cela
implique, soit que le coulis d'injection ait colmaté des
fissures s'étendant jusque sous le profil 3, soit que,
avant traitement, l'écoulement dans le profil ait été
b. l'écoulement théorique a été étudié alimenté partiellement par la partie de digue située à
dans l'hypothèse d'un rapport d'anisotropie de 9. sa gauche et que les injections aient eu pour effet
Fig. 11 — Comparaison des écoulements réel et théoriques dans le profil d'annuler ou de réduire cette alimentation latérale.
3 à la date du 22 novembre 1982. Les réseaux théoriques ont été
déterminés en supposant que la fondation de la digue est imperméable.

Piézomètres P A C : A , B, C, D, E , F
Anisotropie de perméabilité du corps de digue
A
° o-° A o
Charge hydraulique (m) L a méthode la plus fiable de détermination du rapport
k
d'anisotropie de perméabilité — (k et k étant respec-
h v

K
tivement les coefficients de perméabilité horizontale
et verticale) du sol constitutif d'un ouvrage hydrauli-
que en service consiste :
— à mesurer les valeurs de la pression interstitielle
en divers points de l'ouvrage lorsque celui-ci est le
siège d'un écoulement permanent;
— à calculer en ces points les valeurs de la pression
interstitielle en régime permanent pour diverses
valeurs de k ;
h

k
—i—i—i 1—i—i—i 1 1 1—l—r-T—1 1 1—I—V-i—r 1 — > -
— à retenir la valeur de — assurant la meilleure
1

J A S O N D J J A S O N D J A S 0 N
1980 1981 1982 K
Fig. 12 — Variations de la charge hydraulique dans le profil 2 entre le concordance entre les valeurs théoriques et expérimen-
début de la première mise en eau (juillet 1980) et la fin de 1982.
tales de la pression interstitielle.

19
L'application de cette méthode à la digue de L a Pia- Ce résultat, qu'il convient de retenir avec une grande
gne se heurte à plusieurs difficultés. D'une part, prudence compte tenu du peu de mesures significati-
aucune série de mesures de pression interstitielle ne ves effectuées, est sensiblement différent de ceux pré-
correspond à un régime d'écoulement permanent, d'au- sentés par Shérard et al. (1963)** qui recommandent
tre part, les conditions aux limites de l'écoulement d'adopter un rapport d'anisotropie de 25 pour l'étude
sont très mal définies dans le profil 2 (existence d'une des projets de barrages en terre construits au moyen
source dans la fondation à proximité du profil) et de matériaux graveleux extraits d'emprunts homogè-
ne sont pas connues avec certitude dans le profil 3 nes.
(possibilité d'alimentation latérale ou par des fissures
de la fondation). Dans ces conditions, on ne pouvait
procéder qu'à une estimation assez peu précise du CONCLUSION
rapport d'anisotropie.
L'exemple de la digue de L a Plagne fait clairement
Cette estimation a été faite à partir des mesures effec- ressortir la nécessité d'injecter le massif rocheux consti-
tuées dans le profil 3 à la fin du mois de novembre tuant la fondation d'une digue fermant un bassin de
1981, c'est-à-dire plus de deux mois après l'achève- stockage, même si l'ouvrage est de faible importance
ment des travaux d'injection, alors que le niveau de et si une reconnaissance succincte semble indiquer que
la retenue n'avait pas subi de fluctuation sensible le massif est peu fissuré.
depuis un mois.
L'interprétation des mesures piézométriques effec-
A cette époque, l'injection de la fondation avait amé- tuées dans un profil de la digue semble montrer que
lioré l'étanchéité de cette dernière au droit du profil 3, celle-ci est isotrope ou faiblement anisotrope vis-à-vis
et/ou réduit l'alimentation latérale de l'écoulement de la perméabilité, bien qu'elle soit construite à partir
dans ce profil. Les résultats des mesures sont compa- d'un sol graveleux; ce résultat est toutefois à considé-
rés sur la figure 11 aux résultats de l'étude théorique rer avec prudence compte tenu du faible nombre de
de l'écoulement permanent effectuée en supposant que mesures significatives effectuées.
le sol de fondation est imperméable et en affectant
successivement les valeurs 1 et 9 au rapport d'anisotro-
pie. O n constate que si les résultats expérimentaux ne
* Des réseaux d'écoulement théoriques (non représentés ici) établis
concordent pas parfaitement avec les résultats théori- pour des rapports d'anisotropie supérieurs à 9 montrent que les
ques obtenus dans l'hypothèse d'isotropie, les écarts pressions interstitielles correspondantes sont toutes inférieures à
restent relativement limités. E n revanche, les résultats
celles obtenues pour — = 9 et, par suite, diffèrent davantage des
expérimentaux et théoriques divergent fortement lors-
K
que ces derniers correspondent à un rapport d'aniso- valeurs mesurées.
tropie de 9. Cela indique que le comportement hydrau-
* * S H E R A R D J . - L . , W O O D W A R D R . J . , GIZIENSKI S. F . , C L E V E N G E R
lique de la digue s'apparente à celui d'un ouvrage W. A., Earth and earth-rock dams, John Wiley and sons, Ed.
isotrope ou faiblement anisotrope *. (1963).

20
 .IL

Digues
et barrages en terre
T^A^f/A de faible hauteur

Application
de la méthode des éléments finis
aux écoulements en milieu poreux
Pierre HUMBERT
Chef de la section modèles n u m é r i q u e s
Laboratoire central des Ponts et Chaussées

Si les premières applications de la méthode des éléments

finis ont été orientées vers la résolution des équations d'équi-
RÉSUMÉ libre en élasticité ou en élastoplasticité, il est très vite apparu
que cette méthode pouvait aussi être utilisée pour la résolu-
Après un rappel des équations de base régissant
les écoulements en milieu poreux, une techni- tion de nombreux problèmes physiques régis par des systè-
que de résolution par la méthode des éléments mes d'équations aux dérivées partielles.
finis de l'équation de Poisson non linéaire en
régime transitoire est présentée. E n particulier, les problèmes de répartition de la charge
Les adaptations à apporter à cet algorithme hydraulique dans les écoulements en milieu poreux qui, tout
général pour résoudre certains problèmes
particuliers sont exposées et illustrées par des comme les problèmes de répartition de température ou de
exemples dans le cas des écoulements en potentiel électrique, se ramènent à la résolution d'une équa-
milieu poreux saturé avec surface libre, des tion du type A(u)=f, peuvent être traités aisément par la
écoulements en milieu poreux non saturé méthode des éléments finis. Après avoir posé les équations
et dans le cas des calculs de nappes aquifères.
du problème, on abordera la technique générale de résolu-
MOTS CLÉS : 42 • Éléments finis (méthode) - tion et on verra sa particularisation pour trois importants
Écoulement (fluide) - Calcul - Nappe aquifère • types de problèmes : les écoulements en milieu poreux non
Sol - Saturation - /Milieu poreux. saturés, les écoulements en milieu poreux saturés avec sur-
face libre et les calculs de nappes aquifères.

21
B u l l , l i a i s o n L a b o P. et C h . 132 - j u i l . - a o ù t 1 9 8 4 Réf. 2 9 2 0
I - ÉQUATIONS DES ÉCOULEMENTS
EN MILIEU POREUX Le terme p n g ^ P y — P + - J est le coefficient d'em-
s

magasinement spécifique C .
1.1. Cas du milieu poreux saturé
s

L'équation générale de l'écoulement en milieu poreux

L'équation générale qui régit les écoulements en saturé est donc :
milieu poreux saturé est obtenue à partir des trois
équations suivantes : dh
di\(K grad h) = C s Yq (6)
— équation de continuité du fluide : dt

div(p (p«) + p î = 0 (1) I I . 2. Cas du milieu poreux non saturé

et
Si S est le coefficient de saturation, on remplace n
— équation de continuité du solide :
par, nS dans les équations (1) et (3) et k par
div [p (1 -n)s f?J+ | [(1 -n) p ]=0
s (2) k .k dans (3)(k est la perméabilité relative comprise
T r

dt entre 0 et 1). Le coefficient S est relié à la teneur en

— loi de Darcy : eau 9 et à la porosité n par 0 = n S.

Une démarche analogue à celle du cas du milieu

U-nU = s [grad p + p g grad z] (3) poreux saturé conduit à l'équation :

avec : 9 _ ÔQ\Ôh
div (k K grad h) --
- 7Tiir.
r
Cs+ +q (7)
p : masse volumique du fluide,
n dty J dt
n : porosité, avec : \|/ = h — z.
q : débit volumique,
L a perméabilité relative k et la teneur en eau 0 sont
r

0 : vitesse apparente du fluide, des fonctions de O n détermine expérimentalement

U s : vitesse réelle moyenne du solide, la courbe K(G) (fig. 1) et la courbe v|/(0) ou courbe
p : masse volumique du solide, de succion (fig. 2). Cette dernière courbe présente
s

p : pression,
z : cote,
g : accélération de la pesanteur,
k : perméabilité intrinsèque,
p : viscosité dynamique.

Si l'on porte dans (1) la valeur de U donnée par (3) Fig. 1

et si l'on fait intervenir la notion de dérivée particu- Courbe de perméabilité.

— = 1- Û,. grad / ) ,on

dt dt )
obtient :

div p - (grad/j + p g grad z)

dp p dn P"dp
= n — + — — +P 1 (4)
dt l-n dt
Sans entrer dans le détail des calculs, il apparaît
que, moyennant certaines hypothèses simplificatrices,
l'équation (4) peut se mettre sous la forme :

dh
d i v ( X grad h) = pngl p>-p +
et + q (5)
s
Fig. 2
Courbe de succion.

avec :

K= p g - : tenseur de perméabilité,
M*
P r : compressibilité du fluide,
P s : compressibilité des grains solides,
a : compressibilité du milieu poreux.

22
 À

souvent des phénomènes d'hystérésis. Si on les néglige, conduction électrique,

on peut déterminerfe (\|/)et 0(\|i) qui sont alors intro-
r électrostatique,
duits dans l'équation (7). torsion cylindrique, etc.

E n régime permanent, l'équation (7) devient : L'équation (11), valable sur un domaine Q (fig. 3),
est associée à trois types de conditions aux limites sur
le contour :
div (k (v|/) K grad h) = q
r (8)
— condition de potentiel imposé ou condition de
Dirichlet
O n retrouve le cas du milieu poreux saturé en faisant
h = h(P) sur S u

— condition de flux imposé ou condition de Neuman

d'où : div (K grad h) = q qui peut s'écrire : -k'Jp- «, = / ; (P) sur S , 2

dxj
— condition mixte, dite condition d'échange ou de
dxj ( " ôxt ' ^ Cauchy

-k ^-n
ii
i = E(P).[h-h ]snrS ,
e 3

Si, de plus, il s'agit d'un milieu isotrope de OXj

perméabilité k, on obtient l'équation de Poisson :
avec :
-àh=f (9) h(P) : potentiel imposé au point P,
E n régime transitoire, si l'on pose n{ : composantes du vecteur normal au contour,
fi (P) : flux imposé au point P,
C(v|/) = E (P) : coefficient d'échange au point P,
he : potentiel extérieur.
l'équation (7) s'écrit :

e«0 ôh
div (kr K grad h) = C + C(v|/)
s
+q
~ôt

ôh
Le coefficient de — comporte un terme dépendant

des compressibilités du fluide, des grains solides et du

milieu poreux, et un autre dépendant de la présence
d'une zone non saturée, qui est souvent prépondérant Fig. 3
par rapport au premier.

E n faisant k (\|/) = 1 et pour un milieu isotrope, on

r

obtient une équation de la forme :

ôh
-Ah=f+c (10)
~dt

Cette équation n'est autre que l'équation de la chaleur

en régime transitoire. L a démarche suivie pour arriver à une formulation
« éléments finis » est une démarche en deux étapes :

1 — O n ramène le problème de la résolution d'un

II - RÉSOLUTION PAR LA MÉTHODE système différentiel à celui de la recherche du mini-
DES ÉLÉMENTS FINIS mum d'une fonctionnelle (principe variationnel).

2 — O n ramène ensuite le problème de la détermina-

L a méthode décrite ici est une méthode de résolution tion d'une fonction inconnue minimisant cette fonc-
de l'équation de la chaleur : tionnelle à celui du calcul d'un nombre fini d'incon-
nues, par décomposition du domaine Q en un nombre
ôh fini d'éléments sur chacun desquels on impose le type
- d i v ( K grad h)=f+c (H) de variation de la fonction inconnue.
~ô~t

qui pourra s'appliquer aux problèmes suivants : O n traitera d'abord le cas du régime permanent en
milieu isotrope puis anisotrope. O n abordera ensuite
— transfert de chaleur, le régime transitoire et la prise en compte de la non-
— écoulements en milieu poreux, linéarité éventuelle.

23
II. 1. Cas du régime permanent
<">(«)= j" Z ( ^ ) d î î - 2 J ufàn 2
(15)

Principe variationnel

Fonctionnelle associée à un système différentiel

— condition — =f (P) sur S (Neuman);
t

E n régime permanent, l'équation —Ah=f est de la dn

forme plus générale : Su dty ,
on a — = — = / , ; il en resuite :
Au=f (12) dn dn

où A est un opérateur défini, positif et symétrique

agissant sur un espace de fonctions D définies sur Q. A

• Conditions aux limites homogènes - 2 w/dn-2 M / Ì dS (16)

Si on se limite aux fonctions de D nulles sur le A

contour, on démontre que, si ( 1 2 ) a une solution, elle

est unique et minimise la fonctionnelle ôu
— condition \-E (u—u„) = 0 sur S (Cauchy);
F(u) = (Au,u)-2(u,f) (13)
ôn
où (u, v) est un produit scalaire. O n prend en général ôu
elle est de la forme h E u=f 2 (P) sur S.
(u, i;) = u.vdCi. Réciproquement, si u rend F(u) dn
du
minimum, il vérifie ( 1 2 ) .
Comme \-Eu= — +E v|/ = / 2 sur S, on a :
dn du
• Conditions aux limites quelconques dites « non
homogènes » - 2 f M,ifàfì .
n i \ Sx, Jn
O n montre que si \|/ est une fonction de D qui vérifie A

les conditions aux limites, résoudre ( 1 2 ) revient à Eu 2

dS-2 uf dS
2 (17)
rendre minimum

(A u, « ) - 2 («, / ) + (u, A*)- (*|r, A u)

Généralisation à Téquation — div{K grad u)—f
E n général, cette fonctionnelle peut s'écrire :
Dans le cas d'un tenseur non isotrope
a>(u) = (Au,u)-(2u,f)+ f N(u)dS (14)

Application à Téquation — Au =f K
~ [ k 1 2
k y 22

D'après ce qui précède, la fonctionnelle O(w) est de

la forme : on remplace l'opérateur laplacien —À par
* ( « ) = ( - A « , « ) - 2 ( u , / ) + («, -Av|/)-(v|/, - A M )

Sxj \ Sx t )
E n intégrant par partie

(-A«,«)«- uu, dQ,

it (u, — A\|/) et -Au), O n prend alors pour fonctionnelle :

1
$ («) s'écrit : J(u) =
n Sxt Sxj

<D(u) = M u,idn-2
; i ufdCl
fu dCì- N (u) dS.
Jn Jn

- f u^dS- f „*+dS + f
Js <5n J s 5n JJs 5 Sn Le dernier terme dépend des conditions aux limites
(flux imposé ou condition d'échange) :
Selon le type de conditions aux limites et sachant
que i|f vérifie ces conditions, on obtient :

— condition u = û~(P) sur S (Diriehlet),

| N(u) dS= j
Js 2
uf dS-
x | ( ^ -uf2 ) dS

comme u = v|/ = u s u r S, on peut écrire flux échange

24
Formulation matricielle Si l'on pose :

Formulation locale de; dri dÇ

Après décomposition du domaine à étudier en sous-

domaines élémentaires appelés éléments (fig. 4), on
'FVE= f / N | j\àZ,àr\ dÇ (20)
approche sur ces éléments le potentiel inconnu M par
un développement polynomial :

u=Nl/ £
G (N) dS
,
FSE= f E

avec :
N : matrice des fonctions d'interpolation de on peut écrire :
l'inconnue,
U E : vecteurs des valeurs nodales du potentiel. JE ( « ) = \ 'U E KE U E - ( ' F V E + 'F ) SE U E

Pour exprimer que J (u) est minimum, calculons SJ : E E

6J = >U K
E E E 5t/ -CF £ K £ +'F S £ ) 6U E

E n écrivant que 8J est nul pour tout bU admissible

E E

satisfaisant les conditions aux limites et en transpo-

sant l'expression obtenue, il vient :
Fig. 4 K U —F + F, (21)
E E VE
- SE
Décomposition d u domaine.

avec :
KE : matrice de rigidité élémentaire,
F VE : vecteur « forces de volume » élémentaire,
F SE : vecteur « forces de surface » élémentaire.

L'équation (21) exprime la relation très simple entre

le vecteur des inconnues nodales U et le vecteur des E

forces nodales F + F . VE SE

L a géométrie de l'élément est définie par les coordon-

nées des nœuds de l'élément contenues dans le vecteur Formulation globale
X et par les fonctions de forme F qui sont les fonc-
E

tions d'interpolation des coordonnées. Si r\, Ç sont O n définit les potentiels nodaux U de l'élément E

des coordonnées locales, on a : numéro k à partir du vecteur U des potentiels G

Z=[x,v,z] = F ( c ; , q , Q Z nodaux relatifs à l'ensemble des éléments par une

£

relation de la forme :
Lorsque F et N sont identiques, l'élément est dit
U = P" U
E G (22)
« isoparamétrique ».
où P représente une matrice « d'extraction » formée
k

Si Jf et Jf c sont les dérivées des fonctions

L de 0 et de 1. De la même manière, on a :
d'interpolation N par rapport aux coordonnées globa-
les et locales, on a : F = k
E PF k
G

J ( u ) étant en fait l'énergie potentielle élémentaire,

E

écrivons que la variation de l'énergie des « forces »

J désignant la matrice jacobienne. internes de l'ensemble des éléments est égale à la
variation de l'énergie des « forces » externes F pour G

tout 8U admissible :
G

F, TI S 'UÌ K k
E SU* = £ «F* ÒU = 'F k
E G W G

D'après (22) :
L'élément de volume dCl s'écrit d£î = | J | d£, dr\ dÇ, où
| J | est le déterminant de J. X 'U X* ÒU = X 'U 'P
k
E
k
E G
k
P k
bU G

k k
L a fonctionnelle s'écrit alors au niveau élémentaire : = 'I/ (£'F*KÏ/*)5I/
0 G

k
x
Jf c KJf d£ic E 1 U R

O n a donc :

•U K 6U G G G = 'F ÒU G G

- f fNd£l U - E E ! G(N)dS U E E (19)

JaE Js E pour tout 8U admissible G

25
et par conséquent :

KoU G =F G (23)

avec :

IK^X'P* K P
k k

Fig. 5
E

L'équation matricielle (23) représente un système

linéaire dont les coefficients sont les termes de la
matrice de « rigidité » globale KQ, les inconnues les
composantes du vecteur des potentiels nodaux U et G

le second membre les composantes du vecteur des

« forces » ou débits nodaux F . G Les dérivées de N par rapport aux coordonnées loca-
les s'écrivent :

Application au triangle à interpolation linéaire T 3 r-N,u = r i o- i i

|_N,r|J l_0 1 - 1 _
O n choisit pour coordonnées locales les coordonnées
barycentriques \, r\, Ç (fig. 5) :

^aire(MJK) 6 [ ( U ]
L a matrice jacobienne a pour expression

aire (IJK) r = T IKX

yiK "I
aire ( M K I ) L JK yjK J
e[0, 1]
X

11 =
aire (IJK) 3VCC X / j f =
Xj ^K.'

Ç
= ^ Ï Ï K T 6 [
' 0 1 ]
Son déterminant a pour valeur | J | = 2 S où S est l'aire
aire (IJK) de l'élément IJK. O n a d'autre part :
T|, Ç sont reliés par la relation £ + n. + Ç = 1.
yjK
O n choisit pour fonction d'interpolation N et pour 2S X
IK J
fonctions de forme F les coordonnées barycentriques
elles-mêmes : O n en déduit l'expression de la matrice élémentaire
K d'après (20) :
F=N=K,n,a
E

'JTL'J-'KT 1
Jf \ L J\ dei dp

1 0~
yjK X
JK "ir*, o i r ^ -y iri IK o - i dSl
0 1
4 S 2

-1 -1
_ - y ¡ K X
1K JL o k
J L-*« X
IK J Lo 1
-i J

Soit

: k y yKi
x JK
: k y y,j
x J K

+ x2
JK ' +k Xj x¡j
'. + k x y JK xKI y K

KyKiyu (24)
K = - symétrie
E
E

4S k yli
x + ky ll x

+ ky X XjJ
K I

_ symétrie symétrie • k yj + k
x
2
y
x
ÌJ_

où e représente l'épaisseur de l'élément.

Cette expression exprime directement les coefficients analytiquement et nécessite très souvent le recours à
de la matrice élémentaire K en fonction des coordon- E
une technique d'intégration numérique. L a figure 6
nées des nœuds de l'élément et des coefficients du représente quelques types d'éléments utilisés dans les
tenseur de perméabilité ou de conductibilité. Le calcul systèmes PPR et R O S A L I E et dans le programme
de l'intégrale sur l'élément n'est pas toujours possible CESAR.

26
 Fig. 6. — Types d'éléments
utilisés dans les PPR,
ROSALIE et CÉSAR.

II. 2. Cas du régime transitoire Après assemblage, on obtient l'équation matricielle

globale :
Principe variationnel et équation matricielle
KQ V 0 + C G UG = FG (25)

E n utilisant les mêmes procédés que pour le cas du Procédé d'intégration dans le temps
régime permanent, on montre que résoudre :
Il existe de nombreuses méthodes d'intégration dans
du le temps, tant implicites qu'explicites. Parmi les métho-
- d i v ( X grad u)=f+g des implicites, on a choisi celle des différences finies
dt dans le temps et celle de Galerkin.

équivaut à rendre minimale la fonctionnelle : Méthode des différences finies dans le temps

1 A partir des valeurs de U aux temps t _ , et t , on f {

y (M, t) = k« * * dfì calcule les valeurs de U et Ù au milieu de l'intervalle

2J dxi dxj de temps :

[ N(M
)dS

où, comme pour le régime permanent, le dernier terme

dépend des conditions aux limites (voir (18)). At\_ UjtQ-U&.J
Ul f,_ 1 +

At
OU
A t fixé, si l'on pose Ù = — - , on obtient l'équation
E E n reportant ces valeurs dans l'équation (25), on
dt obtient la formule de récurrence de Crank-Nicholson :
matricielle élémentaire :

KE UE + C U = F K C
- + —
E E E

2 At
où K , U et F ont la même signification que dans
E E E

le cas du régime permanent, et où apparaît en plus la + (26)

matrice C « d'emmagasinement », qui a pour
E

expression :
Méthode de Galerkin

C E =- g ' N N dn £ C'est une méthode d'éléments finis dans le temps qui

nécessite le choix de fonctions d'interpolation du

27
 temps dans l'intervalle [ t - i , t j . O n prendra des fonc- f O n en déduit
tions d'interpolation linéaires définies par :

où X =
N tt =X

1
F 51/,. = 0,
Les valeurs interpolées des vecteurs U et F dans l'inter- >-i • f r

6 * * 3
valle de temps sont données par : quel que soit 5(7;.

Û(t) = N _ U(t _ )+N U(t )

u l i 1 ti l
E n regroupant au premier membre les termes en l/j,
(27) on obtient au niveau global après assemblage la for-
F(t) = N , _ F ( t _ ) + i V F ( 0 | 1 ( 1 ( |
mule de récurrence de Galerkin :

L'application du principe variationnel revient à écrire L3 2 Al J

que sur tout intervalle t j , on a :

f
l 5

^-i+2F,
5<D («, t) = 0 (28)

pour tout 5u(t) admissible de à t.;

Technique de calcul

Les deux méthodes décrites ci-dessus se ramènent

Or on a :
toutes deux à une formule de récurrence de la forme :

5<D (u, t) = 'U K E E dU + 'Û C 8U -'F

E E E E E oU ,
E

Û E étant calculé à partir de (27) : Les matrices K et K sont fonctions de At. L'algori-
t 2

thme de calcul est décrit sur la figure 7. O n constate

BU(t) BN . . t BN,. que à chaque modification de At, il faut recalculer ou
pour le moins réassembler les matrices K et K et t 2

triangulariser K Les deux méthodes donnent de très

v

Posons :

L / ( r _ ) = t/,-_ et
i 1 1 U(td=U : t

ù = -u _
E { i 1 + u ^ = l
^i
État initial
Bt At V (to)

D'autre part, 8U est interpolé par E

Calcul de K t

5t7 = ( l - > . ) 8 [ 7 _ + X,S[/

£ ( 1 i

Calcul du second membre

Parmi tous les 3U admissibles, choisissons ceux pour

E

lesquels :
9

§[/,._! =0 => dU = X6Ui E

Non Oui
O n peut alors écrire :

Jt,' ¡ - 1
SO (u, t) =

At [(1 — À.) 'U _ i 1 + 'k 'C7J K X SfJi dX

+ At C X SfJj àX
o At
1
• At I [(1 -X) Fi.^X F J X oUi dX At 1+1 = Af,
9
Jo

1 Non Oui
Or, X(l-X)dX= - et X dX= - 2

o 6 o 3 Fig. 7. — Algorithme de calcul en régime transitoire (modèle linéaire).

28
bons résultats pour les petits pas de temps et sont — la prise en compte de caractéristiques variables
inconditionnellement stables. Toutefois, pour des en chaque nœud rend obligatoire l'utilisation de la
valeurs croissantes des pas de temps, la méthode de technique d'intégration numérique pour le calcul des
Galerkin oscille moins que celle de Crank-Nicholson. matrices élémentaires K et C , même si les éléments
E E

C'est pour cette raison que nous l'avons retenue pour sont géométriquement « droits ».
les groupes 03 et 04 du système R O S A L I E et pour
les programmes B I D I F et N S A T du système PPR. L'utilisation de méthodes explicites, pour lesquelles il
n'est pas nécessaire de faire une résolution à chaque
pas de temps, supprimerait le premier inconvénient,
mais ces méthodes ne sont pas inconditionnellement
Prise en compte de la non-linéarité de comportement stables et nécessitent l'emploi de très nombreux pas
de temps.
Présentation du problème
Technique de calcul
Que l'on soit en régime permanent ou transitoire, les
méthodes présentées aux paragraphes II. 1 et II. 2 L'algorithme de calcul est décrit sur la figure 8. Par
de ce chapitre supposent que les caractéristiques des rapport au cas linéaire, on constate la présence d'une
matériaux ne varient pas, en particulier en fonction boucle supplémentaire à l'intérieur de la boucle sur
de la valeur du potentiel (charge, température) et les pas de temps. D'autre part, la matrice est triangula-
éventuellement du temps. S'il n'en est pas ainsi, les risée à chaque itération de chaque pas de temps.
matrices K et C dépendent de U et de t, et il est alors
nécessaire de faire un calcul itératif (à chaque pas Le processus itératif cesse lorsqu'une certaine norme
de temps en régime transitoire) en réactualisant ces de la différence entre la solution à l'itération j et la
matrices, jusqu'à stabilisation. solution à l'itération j—l est inférieure à e donné.

L'équation à résoudre est de la forme :

III - APPLICATIONS
K, (u,t)U(t )=K (u,
t 2 1)^^ + 6(10.

L a technique d'actualisation (méthode dite « du point Les algorithmes décrits précédemment permettent de
fixe ») est coûteuse en temps de calcul, surtout pour traiter les problèmes d'écoulement en milieu poreux
les modèles tridimensionnels, et ce pour deux raisons : saturé sans surface libre, et les problèmes d'écoule-
— K étant modifiée à chaque itération, il est néces-
t
ment en milieu poreux non saturé sans surface de
saire de refaire sa triangularisation à chaque fois : suintement. Pour pouvoir prendre en compte la pré-
sence d'une surface de suintement, d'une surface libre,
ou pour les calculs de nappes aquifères, les algorith-
mes généraux nécessitent quelques adaptations qui
sont décrites dans les paragraphes III. 2 à III. 4 de
cette partie.
État initial
U(t ) 0

III. 1. Milieu poreux saturé sans surface libre

\ ' =0 I
U n exemple classique est le problème de l'écoulement
i = i+l
autour d'une palplanche. Dans ce type de calcul, les
chargements se réduisent en général à des flux nuls à
; =0
travers certaines surfaces; ce type de condition aux
limites est satisfait implicitement sur un bord libre.
Les conditions aux limites sur la charge hydraulique
ne sont en général des conditions de charge constante
j=J+l sur des parties de contour, ou de charge égale à la
cote (condition de pression nulle) sur des surfaces de
Calcul de ^ (l/j,) suintement connues. Les figures 9 et 10 représentent
le maillage et les réseaux d'équipotentielles et d'équi-
Calcul du second membre pressions obtenus pour un calcul d'écoulement autour
K I O ^ - I + GCI-I) d'une galerie circulaire. L'alimentation est assurée à
la partie supérieure du modèle par une condition de
Triangularisation charge constante. L a surface de la galerie constitue
1 la surface de suintement ou de pression nulle.

Ce calcul peut-être effectué avec les programmes

Oui Non P P R - B I D I F ou P P R - N S A T . Dans ce dernier cas, le
contour de la galerie est déclaré en surface de suinte-
ment, ce qui évite à l'utilisateur de rentrer point par
Nouveau pas Nouvelle point la condition h=y, cette condition étant automa-
de temps itération tiquement mise en place à la deuxième itération. L e
milieu restant saturé, le processus itératif s'arrête
Fig. 8. — Algorithme de calcul en régime transitoire (modèle n o n linéaire). après deux itérations.

29
 Fig. 9. — Maillage pour un calcul d'écoulement autour d'une galerie. Fig. 10. — Équipotentielles ( ) et équipressions ( ).

III. 2. Milieu poreux saturé avec surface libre

Formulation du problème en régime permanent

Considérons le cas du barrage représenté sur la

figure 11. Dans le domaine Q limité par les surfaces
S , S„ S , S et S , l'écoulement est régi par :
h r h f

Fig. 11
- d i v ( X grâd h)=f (29)

où / correspond à un débit volumique en général nul.

Les conditions aux limites sont de deux types : sur la charge et une condition sur le flux : S, est
simultanément une surface à pression nulle et une
— celles qui s'appliquent sur des limites connues : surface de courant.

h = h(P) sur S (charge imposée)

h

h = h' (P) sur S' (charge imposée)

h (30) Méthodes de résolution
ôh
— k — n , = / i (P) sur S
iJ
(flux imposé).f L'idée générale consiste à se fixer au départ de façon
! dxj arbitraire la position de S„ à lui imposer la condition
de flux nul (cette condition est implicite sur un bord
— celles qui s'appliquent sur la surface libre S, et la libre) et à ajuster sa position à chaque itération en
surface de suintement S„ qui ne sont pas des données cherchant à satisfaire la condition sur la charge
du problème : h = z(P). Cela s'obtient en déplaçant les points de S,
de manière à ce que leur cote soit égale à leur charge
/ h = z (P) sur S, et S r (pression nulle) à l'itération précédente.

—k iJ
— n, = 0 sur S, (flux nul). ( 3 1 )

Méthode de suppression d'éléments

Toute la difficulté consiste à résoudre (29) sur un Cette méthode, qui n'est jamais employée seule,
domaine qui n'est pas connu au départ, tout en satis- consiste à obtenir une première approximation de la
faisant les conditions aux limites (30) et (31). O n voit surface libre en supprimant à chaque itération les
que sur la surface libre, on a à la fois une condition éléments pour lesquels l'ensemble des sommets ont

30
 Fig. 1 2 . - Méthode de suppression d'éléments. Fig. 14. — Maillage initial.

- T

Fig. 13. — Ajustement final. Fig. 15. — Maillage déformé.

une charge inférieure à la cote. L a surface libre obte-

nue n'est pas régulière, surtout si le maillage est gros-
sier (fig. 12). Cette première phase est en général
suivie d'une phase d'ajustement consistant à déplacer
les nœuds du maillage en les amenant à une cote égale
à leur charge (fig. 13). Les points dont la cote est
amenée à diminuer peuvent rejoindre leurs voisins
immédiatement inférieurs et il faut à ce moment là
supprimer les éléments « plats » ainsi générés. Cette
technique est utilisée dans le groupe 06 du système
Fig. 16. — Détermination de R par extrapolation.
R O S A L I E pour les recherches de surfaces libres en
domaine bidimensionnel ou axisymétrique.

Méthode de déformation d'un ensemble d'éléments position des points de ruissellement après déplacement
des autres points, par extrapolation linéaire ou parabo-
Contrairement à ce qui se passe pour la méthode lique à partir des deux ou trois derniers points de S,
précédente, tous les éléments restent actifs lors du (fig. 16).
processus itératif. Le principe qui régit le déplacement
des nœuds reste le même, mais on déplace non seule- Exemple de calculs avec le programme FLO T
ment les points de la surface libre, mais aussi un
certain nombre de nœuds situés en dessous d'eux, en Le programme F L O T du système PPR utilise la tech-
les répartissant régulièrement sur ce que l'on peut nique de déformation d'un ensemble d'éléments. L a
appeler des « droites mobiles ». L a limite inférieure figure 17 représente le résultat d'un calcul d'écoule-
de la zone déformable doit être choisie de façon à ne ment à travers un barrage homogène.
pas perturber les conditions aux limites sur la charge.
Les figures 14 et 15 représentent le maillage d'un Limitations des méthodes précédentes
barrage avant et après déformation. Cette technique
est utilisée dans le programme F L O T du système Dans la mise en œuvre pratique des méthodes décrites
PPR. ci-dessus, on peut se heurter à deux types de difficultés
en milieu non homogène :
Traitement des points de ruissellement
— celles qui sont dues à une trop grande hétérogé-
Les points de ruissellement éventuels sont situés à néité des perméabilités ;
l'intersection de la surface libre et des surfaces de — celles qui sont dues à l'obligation de respecter la
ruissellement. Comme ils appartiennent à ces derniè- géométrie des interfaces.
res, on leur applique la condition h=z, ou condition
de pression nulle. Comme ils appartiennent aussi à la U n trop grand contraste de perméabilité dans le cas
surface libre S,, il est évident que le traitement appli- d'un noyau peu perméable (fig. 18) obligera les élé-
qué aux autres points de S, (faire z = h) ne les fera ments situés juste en aval du noyau à se déformer
pas bouger. Pour éviter ce blocage, on recalcule la beaucoup, car la surface libre subit une chute rapide

31
 Fig. 1 7 . — Écoulements à travers un barrage homogène (FLOT).

Cas du régime transitoire

Les techniques de déformation de maillage sont d'em-

ploi délicat en régime transitoire, car il faut pouvoir
traiter des mouvements importants de la surface libre,
tant en descente qu'en remontée. Ces difficultés étant
supposées résolues, on peut utiliser des formules de
récurrence du type Galerkin (28), où l'on itère sur la
géométrie à l'intérieur de chaque pas de temps; de
plus, on doit à chaque pas de temps prendre en
kx k2 k 3
compte dans le chargement le flux qui a traversé la
Fig. 18. — Utilisation d'un point de ruissellement fictif.
surface pendant ce pas de temps.

Tout cela conduit à des programmes très complexes,

et il est préférable d'utiliser une méthode de maillage
fixe et de « rigidité » variable, la zone située au-dessus
dans le milieu perméable. Cela a souvent pour consé-
de la surface libre étant traitée comme une zone non
quence d'entraîner des perturbations dans la partie
saturée.
de la surface libre située dans le noyau. O n supprime
cet inconvénient en calculant la position du point R*
comme celle d'un point de ruissellement, par extrapola-
tion à partir des deux ou trois points qui le précèdent. I I I . 3. Milieu poreux non saturé
O n utilise cette technique dans le groupe 06 de
k Bases de la méthode de calcul

ROSALIE pour des rapports de perméabilité — O n a vu en I que la non-saturation du milieu poreux

pouvait être traduite par une non-linéarité du tenseur
supérieurs à 100 environ. ôfo
L'autre type de difficulté correspond surtout au cas de perméabilité et du coefficient du terme en —
où la surface libre traverse un interface horizontal ôt
ou sub-horizontal. L a technique de déplacement des (équation (7)). U n algorithme permettant la prise en
points sur des lignes sub-verticales ne peut s'appliquer compte de ces non-linéarités est décrit en II. L a don-
aux points R et R (fig. 19), sous peine de déformer
t 2
née sous forme de fonction de v|/=h—z de la perméabi-
la géométrie des interfaces. Une technique possible en lité relative k et de la teneur en eau 0 permet à
r

cours de test consiste à appliquer la méthode classique chaque itération d'actualiser les coefficients du tenseur
aux points autres que R et R et à calculer la position
t 2
de perméabilité K et le coefficient d'emmagasinement
des points R et R par extrapolation à partir des
t 2

points de la surface qui les précèdent et qui les suivent.

Traitement des surfaces de suintement

Le long d'une surface de suintement, les conditions

aux limites changent de nature au passage du point de
ruissellement, l'intersection de la surface de suintement
et de l'équipression correspondant à la pression nulle.
Au-dessus du point de ruissellement, on a une condi-
tion de flux nul, en dessous, une condition de pression
nulle se traduisant par h=z. Lors des itérations per-
mettant l'actualisation des coefficients, les conditions
aux limites sur les surfaces de suintement sont modi-
fiées selon la méthode suivante :

Première itération

O n impose la condition de flux nul sur tous les points

Fig. 19. — Traversée de limites de zones subhorizontales. de la surface de suintement.

32
Passage de l'itération i à l'itération i+1

• Cas des points à condition de flux nul à

l'itération / :
— si h>z, on imposera h=z;
— si h_z, on maintiendra la condition de flux nul.

• Cas des points à condition h=zk l'itération i :

— si le flux est rentrant, on imposera la condition de
flux nul.
— si le flux est sortant, on maintiendra la condition
Fig. 20. — Infiltrations dans un massif.
h=z.

Pour éviter certaines instabilités, on ne modifiera la

condition aux limites que sur un seul point de chaque
surface de suintement à chaque itération. Cette techni-
que est employée dans le programme N S A T du sys-
tème P P R .

Exemples de calcul en régime permanent

Le premier exemple (fig. 20) concerne un problème

d'infiltration dans un massif à partir d'une alimenta-
tion sur une partie de la surface [Vauclin et al., 1979].

Les fonctions retenues pour représenter les variations

de k et de 0 en fonction de v|/ = /i—z sont de la
r

forme :
lsii|/ 0
=

k
' = X
P—— si xtr<0

( e si vh^o
s

( Y+ M 8

Dans l'exemple traité, on a pris :

a = 3.10 6

p= 4
- 1,5 - 1. -0,5 0
y = 4.10 4

Fig. 2 1 . — Variations de kr et de 8.
5 = 2,9
6 = 0,3
S

L'allure des courbes correspon-

dantes est donnée sur la
figure 21 pour diverses valeurs
de p. L a courbe en pointillé
représente le coefficient d'emma-
ÔQ
gasinement — .
t3v|/
Ce type de courbe est bien
adapté à des matériaux sableux
ou argileux et a été introduit
dans le programme N S A T du sys-
tème P P R .

Fig. 22. — Vecteurs vitesses,

équipotentielles et équipressions.

33
Compte tenu des importantes variations de perméabi- Le troisième exemple (fig. 25) concerne un barrage
lité que l'on peut avoir au-dessus de la zone saturée, avec drain horizontal. Contrairement à ce qui se passe-
la finesse du maillage dans cette partie du modèle rait avec le programme F L O T , le maillage ne nécessite
doit être d'autant plus grande que la perméabilité aucune adaptation particulière.
relative décroît rapidement pour des valeurs négatives
de la pression. Avec des éléments de type T3 (triangles Exemple de calcul en régime transitoire
à interpolation linéaire, dont l'usage est recommandé
dans le programme N S A T ) , il faut éviter que la L'exemple représenté sur les figures 26 à 29 représente
perméabilité relative soit divisée par un facteur supé- un calcul de barrage homogène en vidange rapide.
rieur à 3 sur la hauteur d'un élément. L a figure 22 Les valeurs des paramètres utilisés sont les suivantes :
représente les vecteurs vitesse, les équipotentielles
et les équipressions pour un flux entrant de ,fc =fc,= 2.10- 'm/s
x

4 , 1 . 1 0 " m/s et une perméabilité saturée de

5

9 , 7 . 1 0 " m/s.
5
perméabilité a=2.10 7

P=4
Le deuxième exemple (fig. 23 et 24) reprend l'exemple
du barrage homogène déjà traité avec le programme e =o,o6
s

F L O T (fig. 17). Comme il s'agit d'un calcul en régime

permanent, seuls les coefficients a et P interviennent teneur en eau \ y = 3 980
dans les résultats du calcul. Comme pour l'exemple 8 =1,8
précédent, y, 8 et 8 ne servent qu'à déterminer à titre
S

indicatif la teneur en eau en tout point.

Les charges sont initialisées par un calcul en régime
O n a pris a = 3 . 1 0 et P = 4,5. Une étude paramétri-
6 permanent (fig. 27). O n ramène ensuite la hauteur
que a montré que pour des valeurs de P comprises d'eau amont de 5 à 1 m en 8 pas de temps de
entre 3 et 5, les résultats variaient peu. L a comparai- 3 heures. Cette hauteur est maintenue constante jus-
son des figures 17 et 23 montrent une très bonne qu'à (=10 jours. Les figures 28 et 29 représentent les
concordance des deux méthodes de calcul. résultats pour t = 30 heures et t = 10 jours.

Fig. 23. — Maillage d'un barrage homogène (NSAT).

Fig. 24. — Écoulements à travers un barrage homogène (NSAT).

Fig. 25. — Barrage avec drain horizontal (NSAT).

34
 JL

Fig. 27. — Régime permanent.

Fig. 28. — Vidange rapide ( f = 3 0 h).

Fig. 29. — Vidange rapide ( f = 1 0 j ) .

III. 4. Calcul des nappes aquifères

Formulation du problème div (T grad h) = C — + q

s (32)
ôt
L a détermination des écoulements dans une nappe
où T représente le tenseur de transmissivité. Si e (x, y)
nécessite théoriquement la résolution d'un problème
est l'épaisseur de la nappe au point M , le tenseur de
tridimensionnel. Toutefois, si l'on fait l'hypothèse sim-
transmissivité est relié au tenseur de perméabilité par
plificatrice que la charge hydraulique h est constante
la relation simple :
sur une verticale (hypothèse de Dupuit), on est ramené
à résoudre en domaine bidimensionnel l'équation : f=e(x,y)K

35
L'épaisseur e dépendant directement de la charge h, L a deuxième difficulté est d'ordre théorique : dans la
il s'agit à nouveau, comme pour les écoulements en formulation exposée ci-dessus, on suppose connu le
milieu poreux non saturé, d'un cas particulier de non- tenseur de perméabilité K en tout point du modèle.
linéarité, qui sera traité par la méthode générale Autrement dit, on considère le modèle calé au départ,
décrite en IL et on l'utilise directement pour déterminer l'influence
de diverses perturbations (infiltration, pompage, modi-
Technique de résolution en régime permanent fication locale de perméabilité pour représenter un
ouvrage, etc.).
Pour pouvoir déterminer l'épaisseur de l'aquifère e, il
faut connaître en tout point les cotes du substratum E n fait, les reconnaissances sur le terrain fournissent
imperméable et du toit de la formation poreuse, z surtout des hauteurs d'eau et rarement des perméabili-
inf

tés.
et z -sup

A chaque itération du processus d'actualisation de la L a première partie de l'étude numérique consistera

transmissivité, l'épaisseur e est calculée de la façon donc à caler le modèle, de façon à retrouver au mieux
suivante selon la valeur de la charge h : les hauteurs d'eau mesurées avant perturbation. Ceci
sera obtenu par un ajustement, souvent délicat, des
si / i ^ z s u p , e= z s u p -z i n f (nappe captive), perméabilités. Il s'agit de résoudre un problème du
si z inf <h<z , sup e = h-z in{ (nappe libre), type « problème inverse », à savoir : quelle est la
répartition de perméabilités qui correspond à la piézo-
si h^z , e = 0 (domaine hors d'eau).
métrie observée. L a solution d'un tel problème n'est
inf

pas unique. E n effet, on peut multiplier les perméabili-

Cet algorithme est programmé dans le groupe 03 du tés par un même coefficient à l'intérieur d'un tube de
système R O S A L I E . L a convergence est en général courant sans changer la répartition des charges. Dans
obtenue pour un nombre d'itérations inférieur à 10. les méthodes numériques de résolution de ce problème
inverse, on rajoute des contraintes d'homogénéité qui
Difficultés de mise en œuvre permettent de lever cette difficulté.

L a première difficulté est d'ordre pratique : il s'agit,

une fois le maillage réalisé, de fournir pour chaque Exemple de calcul
nœud du modèle la cote du substratum et l'épaisseur
maximum de l'aquifère ( z —z ). sup inf
Il s'agit d'étudier l'influence sur une nappe de la
construction d'une rocade dont certains tronçons sont
L'utilisation d'un pré-processeur d'interpolation (pro- construits en tranchée, entre deux parois fichées dans
gramme T E R P O L ) permet de générer ces deux le substratum imperméable. Ces tronçons constituent
tableaux à partir de valeurs connues en un nombre donc un écran pratiquement étanche. Le maillage
réduit de points différents des nœuds du maillage. utilisé est représenté sur la figure 30. Il comporte des

Fig. 30. — Impact hydrogéologique d'une rocade (maillage). Fig. 3 1 . — Calage du modèle.

36
 CONCLUSION

Après avoir présenté un algorithme général

pour la résolution de l'équation de Poisson
non linéaire en régime transitoire, on a étu-
dié quelques modifications qui permettent de
traiter les cas particuliers que représentent
les écoulements en milieu poreux saturé avec
surface libre, les écoulements en milieu
poreux non saturé et les calculs de nappes
'Il aquifères. Les techniques numériques utili-
/ sées sont celles qui ont été programmées dans
les « groupes » 03 et 04 de R O S A L I E et
dans les programmes F L O T et N S A T du
système PPR. Les divers exemples de calculs
présentés montrent l'étendue des champs
d'application de ces programmes, en particu-
lier pour les problèmes d'écoulements dans
les digues ou les barrages.

Fig. 3 2 . — Impact de l'ouvrage.

RÉFÉRENCES BILIOGRAPHIQUES

GUELLEC P., ROSALIE, Système de calcul des massifs et des

éléments quadratiques à 6 et 8 nœuds et des éléments structures. Notice de présentation et d'utilisation, Rapp.
de continuité (*) sur le tracé de l'ouvrage. Pour L C P C , 1976.
simuler la coupure correspondant à la partie étanche ISAACS L. T., Adjustement of Phreatic Line in Seepage Ana-
du tracé, il suffit de rendre ces éléments de continuité lysis by Finite Element Method, Research Report
inactifs. N ° CE2, Université de Queensland, 1979.
NEWMAN S. P., Saturated-Unsaturated Seepage by Finite Ele-
Les conditions aux limites sont du type charge impo- ments, J. Hydraulics Div., A S C E , 1973.
sée sur A B (rivière) et B C (équipotentielle) et du type VAUCLIN M . , KHANJI D . , VACHAUD G., Experimental and
flux nul sur C D et D A . L a rivière est bordée sur une Numerical Study of a Transient, Two-Dimensional,
grande partie de sa longueur par une bande mince Unsaturated-Saturated Water Table Recharge Problem,
très peu perméable. L'alimentation est assurée par Water Resources Research, vol. 15, 5, 1979.
une infiltration sur toute la surface. L a figure 31 repré- W A N G S. Y. et al, Finite Elements in Water Resources, C. R.
sente les équipotentielles du modèle calé. L a zone 3 Conf. Inter, sur l'utilisation des elements finis en hydro-
e

située à l'extrême gauche est hors d'eau. L a figure 32 géologie, Université du Mississipi, 1980.
représente l'effet de la construction de l'ouvrage, ZIENKIEWICZ O. C , The Finite Element Method in Enginee-
rendu étanche sur les parties en trait fort. ring Science, Me Graw-Hill, Londres, 1971.

37
The automatic design calculation of earth dams and dykes Interpretation of measurements of pore pressure in the
J.L. BANGRATZ - H. JOSSEAUME - B. LEBRUN Maurepas-Courance dyke. Anisotropy of permeability of the
structure
The purpose of the DIGUE conversational programme is to facilitate the H. JOSSEAUME - C. MAHDAVI - J. JARDIN
application of the NSAT hydraulic design calculation programme and the
PETAL programme for calculating slope stabilization (modified under the The authors give an account of the study of the hydraulic behaviour in
name STADIG) in several configurations of heterogeneous dykes built on service of a homogeneous impounding dam 15 metres high, constructed
two-layer soil systems. in compacted clay of low plasticity. This study, whose main purpose
The dyke is in four parts (upstream slope, centre, downstream slope, was to characterize the anisotropy of permeability of the structure, was
downstream drainage), and there are several variants of each part. Fitting based on:
these variants together makes it possible to represent quite a large number - Measurements of pore pressure in the structure, using piezometric
of different configurations which may have upstream and downstream instruments set in place during its construction, and additional
berms, a core with or without vertical chimney drains, a downstream piezometric instruments specially installed.
drainage blanket, a toe fill, an upstream or central cut-off wall, and - Measurements of permeability in situ, made in the dyke and its
downstream relief wells. foundation.
Modifications to the PETAL programme relate to the definition of the - Theoretical studies of the permanent flow through the structure, using
potential slipping circles and the improvement of performance characteris- the method of. finite elements and allowing for various conditions of
tics in seeking values of the pore pressure along the potential slipping anisotropy of permeability.
circles. Analysis of the experimental results and comparison of the theoretical and
experimental values of pore pressure reveal an almost isotropic behaviour
Bull, liaison Labo. P. et Ch.. n°130, mars-avr. ¡984, p. 5-14.
of the dyke in respect of permeability.
Bull, liaison Labo. P. et Ch., n°13l, mai-juin 1984, p. 23-28.

Examples of specific problems relating to low dykes

Observations on the Roussières dyke
C. M1EUSSENS
J. JARDIN - R. BALDIT - A. DELFAUT

The author relates some examples of specific problems encountered

during the design or the follow-up after putting into service of five dykes Located on the territory of the new town of Saint-Quentin-en-Yvelines (in
located in the area covered by the Laboratoire Régional des Ponts et the Yvelines département) the Sourrières dyke is a dam of medium height
8
Chaussées in Toulouse. (about 15 metres) holding back a 400,000 m rainwater reservoir. It is a
1. The case of a dyke built on excessively permeable soil. The dyke and homogeneous silt dyke with a drainage system (drain shaft and drainage
underlying soil waterprofing was realized by means of a grout curtain blanket). Though of relatively modest dimensions, the structure, which
which proved inadequate in view of the low admissible leakage flow. concerns public safety, was equipped with important inspection
2. A low dam holding back a very considerable quantity of water, whose facilities. Observations of its hydraulic and mechanical behaviour were
downstream face was shown by appropriate instrumentation to be made duringfillingof the reservoir to the highest level on completion of
tending towards certain failure. its construction, and for two years at the permanent level after the dyke
3. Investigative methods employed to study the permeability of the soil had been put into service. In particular, after additional piezometers had
underlying a dam, which was a priori geologically very unfavourable been installed a study was made, at the permanent level, of the hydraulic
because of the hazardous presence of permeable horizons. For the body anisotropy of the silt constituting the body of the dyke. After reviewing
of the dyke, the influence of lime treatment on the permeability of the the studies carried out at the planning stage, the authors report the results
compacted material was measured. of these observations. The almost isotropic behaviour (k /kv = 1) of
n

4. The design of flood control dam so that it could best ensure maximum the silt is the most important result of these observations. The authors
storage. The problems of internal erosion and by-passing are more believe that this fact is mainly linked with the way the material was
particularly examined. applied.
5. The special case of an expanse of water located near a river,
necessitating a study of the stability of the sealing curtain should the Bull, liaison Labo. P. et Ch., n°131, mai-juin 1984, p. 29-44.
dyke function in reverse (pond empty, river in flood). This same
example also includes a description of a simulation of plane flow by
electrical analogy.
In conclusion, in giving an account of sometimes simplistic methods of The construction and behaviour of an earth dyke at high
study, the author raises the problem of adequation between the level of altitude
studies and the economy of the project. M. LAGRANGE - J. PERR1N - J. JOSSEAUME

Bull, liaison Labo. P. et Ch., n°130, mars-avr. 1984, p. 15-24.

The provision of a storage reservoir to improve the water supply of a
winter sports resort necessitated the construction of a low dyke (9.50 m) at
high altitude.
After presenting the principal geotechnic characteristics of the site, which
led to the construction of a homogeneous dyke, and a description of the
Observations on the Maurepas-Courance dyke dyke itself, the authors give an account of the problems which arose
A. DELFAUT - J. JARDIN - R. BALDIT during and after the work; problems connected with the laying of a
geotextile between the body of the dyke and the protective enrockment,
and considerable leakages when the reservoir was filled, necessitating the
The Maurepas-Courance storage basin is part of the development injection of the rock foundation.
programme of the new town of Saint-Quentin-en-Yvelines. It is a storm The article also reports measurements of pore pressures in the dyke, and
water reservoir serving as a buffer basin for the rainwater of the plateau, their interpretation.
whose urbanisation has increased the peak flow of flood water at the
bottom of the valley. Bull, liaison Labo. P. et Ch., n° 132, juil.-août 1984, p. 13.
The task assigned to the LROP was more particularly to keep track of the
behaviour of the body of the dyke during its construction, and
subsequently after it had been put into service. For this purpose it was
necessary to fit two transverse profiles of the dyke with pore pressure Application of the method offiniteelements toflowin a porous
cells, settlement recorders and piezometric tubes. Optical measurements medium
were made on concreted reference points so as to assess movements of the
body of the dyke. P. HUMBERT
The authors report the principal results obtained during the construction
of the dyke, during pondage, and during the operation of the basin. After a reminder of the basic equations governing flow in a porous
To attempt to understand the real hydraulic functioning of the structure, it medium, a technique of solution, using the method of finite elements, of
appeared worth-while to compare the observations with the initial Poisson's non-linear equation under transitory conditions is presented.
theoretical analyses carried out in the context of research on the The adaptations to be made to this general algorithm in order to solve
anisotropy of the silts used in the body of the dyke. Clear conclusions certain special problems are explained and illustrated by examples in the
were not reached, and it is difficult to take account of the complexity of case offlowin unsaturated porous media and in the case of water-table
real phenomena in theoretical studies, the more so in the planning stage. calculations.
Bull, liaison Labo. P. et Ch., n°131, mai-juin 1984, p. 5-22. Bull, liaison Labo. P. et Ch., n" 132, juil.-août 1984, p. 21.