Remerciements
Remerciements
Remerciements
Tout d'abord, nous tenons à exprimer notre gratitude envers Dieu pour nous avoir accordé la
santé, l'inspiration et la patience nécessaires pour mener à bien ce travail.
Nous aimerions également exprimer notre profonde reconnaissance envers tous ceux qui ont
contribué, de près ou de loin, à l'élaboration de ce rapport. Votre soutien et vos efforts sont
hautement appréciés, et nous tenons à vous témoigner toute notre gratitude et notre respect.
Nous espérons que notre reconnaissance sincère est transmise à tous ceux qui ont contribué à
notre succès. Nous leur sommes profondément reconnaissants pour leur aide précieuse.
Page | 1
Résumé
Le calcul de la stabilité des barrages aux séismes s'est beaucoup développé ces dernières
années. L’amélioration des techniques de mesure a montré, lors des séismes, que les niveaux
d’excitation sismiques sont souvent supérieurs à ceux retenus lors du dimensionnement. Ceci a
conduit à l’obligation de considérer l’accroissement du niveau d’excitation lors du
dimensionnement de nouveaux ouvrages et lors des nouvelles évaluations de tenue au séisme
des structures existantes.
Cette étude repose sur une méthodologie rigoureuse qui combine des analyses numériques et
des simulations d'éléments finis. Les logiciels géotechniques tels que Géo-Slope et Plaxis sont
utilisés pour modéliser les propriétés du sol, simuler la construction des barrages, et analyser la
stabilité des pentes. Dans le cadre de l'étude sismique, des logiciels spécialisés tels que
ABAQUS est employé pour reproduire les effets dynamiques.
Comprendre le comportement des barrages en remblai et des barrages poids dans des
conditions variées, y compris en réponse aux séismes, revêt une importance cruciale pour
garantir la sécurité, la durabilité, et l'efficacité de ces structures essentielles à l'infrastructure
hydraulique mondiale. Les résultats de cette étude fourniront des informations précieuses pour
l'ingénierie géotechnique et contribueront au développement de meilleures pratiques de
conception et de construction des barrages.
Ce rapport est structuré de manière à présenter en détail les différentes étapes de l'étude, les
paramètres utilisés, les modèles de simulation, ainsi que les résultats obtenus à partir des
analyses effectuées. Chaque section offre un aperçu approfondi des aspects spécifiques étudiés,
offrant ainsi une compréhension holistique du comportement des barrages en remblai et des
barrages poids.
Page | 2
Table de matières
Remerciements..............................................................................................................................1
Résumé..........................................................................................................................................2
Introduction générale.....................................................................................................................7
Chapitre Ⅰ : Bibliographie...........................................................................................................9
I. Introduction...........................................................................................................................9
II. Généralités.............................................................................................................................9
1. Définition......................................................................................................................9
2. Rôle :.............................................................................................................................9
4. Ouvrages hydrauliques................................................................................................11
1. Cadre géologique.........................................................................................................29
Page | 3
1.2.1 Stratigraphie...............................................................................................................33
2. Le relief.......................................................................................................................37
3. Topographie.................................................................................................................37
4. Réseau hydrographique...............................................................................................38
1. Introduction.................................................................................................................44
1. Cas de charges.............................................................................................................49
I. Présentation de l’ouvrage....................................................................................................74
1. Situation géographique................................................................................................74
Page | 4
Page | 5
Liste des figures
Figure 1 : Le rôle du barrage.......................................................................................................12
Figure 2 : Évacuateur des crues...................................................................................................14
Figure 3 : Vidanges de fonds.......................................................................................................15
Figure 4 : Plan du barrage en Terre à noyau argileux..................................................................16
Figure 5 : Barrage en terre compactée, homogène......................................................................17
Figure 6 : Barrage à noyau..........................................................................................................18
Figure 7 : Barrage à masque........................................................................................................18
Figure 8 : Barrage poids en béton...............................................................................................19
Figure 9 : Barrage voûte en béton...............................................................................................19
Figure 10 : Barrage poids à contreforts.......................................................................................20
Figure 11 : Vallée encaissée avec des flancs presque verticaux..................................................21
Figure 12 : Vallée étroite en V.....................................................................................................21
Figure 13 : Vallée étroite en U.....................................................................................................22
Figure 14 : Vallée large................................................................................................................22
Figure 15 : Site du barrage Kharroub dans la wilaya de Tanger.................................................29
Figure 16 : Coupe type de la digue principale du barrage...........................................................30
Figure 17 : Coupe type de la digue de col...................................................................................30
Figure 18 : Coupe type de l'évacuateur de crues.........................................................................31
Figure 19 : Coupe longitudinale de la galerie de vidange de fond..............................................31
Figure 20 : Coupe type de prises d'eau potable...........................................................................32
Figure 21: Le bassin Tahaddart au nord du Maroc dans un extrait simplifié de la carte
géologique du Maroc à 1/1000.000 « Trait noir : limites du bassin ».........................................32
Figure 22 : Log lithostratigraphique de la zone d’étude réalisé par le département...................37
Figure 23 : Relevé géologique de la tranchée réalisée en rive gauche........................................38
Figure 24 : Diagramme de répartition de la fracturation dans la région.....................................39
Figure 25 : Carte des pentes et failles du bassin de Tahaddart (Réf : science Lib Editions
Mersenne)....................................................................................................................................40
Figure 2626 : La courbe H-S-V du barrage kharroub.................................................................41
Figure 2727 : Carte géologique du bassin versant.....................................................................42
Figure 2728 : Les cours d'eau (classe1) du bassin de Tahaddart. (Réf : Science Lib Editions
Mersenne)....................................................................................................................................43
Figure 2829 : Evolution des apports du site du barrage KHARROUB.....................................44
Page | 6
Figure 2930 : Hydrogrammes de crues du site KHARROUB....................................................44
Figure 3031 : Carte topographique de la région Ayacha, site du barrage....................................46
Figure 3132 : image représentative de la structure du sol argileux (Réf : prise de photo sur
chantier).......................................................................................................................................46
Figure 3233 : Les propriétés physiques des matériaux dans les différentes zones du barrage.. 49
Figure 3334 : Fuseaux granulométrique des matériaux de construction.....................................50
Figure 3435 : Courbe granulométrique du noyau et du filtre, et critères du filtre......................50
Figure 3536 : Calcul de stabilité, cercles de glissement et description......................................54
Figure 37 : Ensemble des forces appliquées sur une tranche (Méthode ordinaire)....................55
Figure 38 : Ensemble des forces appliquées sur une tranche (Méthode de Bishop simplifiée). .56
Figure 39 : les fonctions typiques de f(x)...................................................................................58
Figure 40 : Menus disponibles sur Slope/W de Géostudio.........................................................62
Figure 41 : Coupe géologique du site du barrage Kharroub.......................................................65
Figure 42 : Résultats du calcul de stabilité du parement amont en régime permanent avec
séisme (Digue principale) :.........................................................................................................66
Figure 43 : Résultats du calcul de stabilité du parement aval en régime permanent avec séisme
.....................................................................................................................................................67
Figure 44 : Génération de maillage dans le barrage en remblai..................................................68
Figure 45 : Déplacements totaux du versant en présence de la nappe........................................68
Figure 46 : L’effet du poids propre du barrage en remblai..........................................................69
Figure 47 : dgdrghtujyj................................................................................................................69
Figure 48 : variation du coefficient de sécurité en fonction du poids volumique.......................70
Figure 49 : Variation du coefficient de sécurité en fonction de la cohésion................................71
Figure 50 : Variation du coefficient de sécurité en fonction du Module de Young.....................71
Figure 51 : Variation do coefficient de sécurité en fonction du coefficient Poisson...................72
Figure 52 : Variation du coefficient de sécurité en fonction de l’angle de frottement Interne....73
Figure 53 : Variation du coefficient de sécurité en fonction de la hauteur du talus....................73
Figure 54 : Variation du coefficient de sécurité en fonction du niveau de la nappe....................74
Figure 55 : Situation géographique du barrage de l’Oued Fodda...............................................77
Figure 56 : Profil type du barrage de l’Oued Fodda...................................................................79
Figure 57 : Propriétés des matériaux du barrage-poids de l’Oued Fodda...................................79
Page | 7
Liste des tableaux
Tableau 1 : : Données des précipitations et températures mensuelles de 2014...........................46
Tableau 2 : prédimensionnement des digues..............................................................................58
Tableau 3 : Les paramètres géotechniques des différents matériaux utilisés..............................63
Tableau 4 : Les résultats de calcul de stabilité du parement en amont sans et avec séisme de 0,14g
.....................................................................................................................................................64
Tableau 5 : Les résultats de calcul de stabilité du parement en aval sans et avec séisme de 0,14g65
Page | 8
Introduction générale
L'eau, élément essentiel à la vie, devient de plus en plus rare, faisant de la sécurisation de l'eau
potable une priorité au Maroc. Face à la saturation des eaux souterraines, au climat semi-aride
et aux défis de la sécheresse, le pays a adopté la stratégie de mobilisation des eaux de surface
depuis les années 30, principalement à travers la construction de barrages.
Les barrages ont une grande importance pour le développement rurale d’un pays, mais après
construction ces barrages sont exposés au risque de rupture qui est certes assez rare mais quand
elle arrive, les conséquences peuvent être catastrophiques et occasionner des dommages
considérables sur les hommes et les biens.
Le présent rapport de projet propose une étude approfondie sur le comportement de deux types
de barrages, à savoir les barrages en remblai et les barrages poids, en mettant particulièrement
l'accent sur leur stabilité des pentes et leur réaction vis-à-vis des séismes. Cette étude est menée
à l'aide de divers logiciels géotechniques et d'éléments finis, offrant ainsi une analyse complète
et détaillée des aspects fondamentaux de ces structures cruciales dans le domaine de l'ingénierie
hydraulique et géotechnique.
Page | 9
Données Utilisations
Chapitre Ⅰ : Bibliographie
Concepts généraux
Objectif
Généralités sur les barrages et leurs fonctionnalités
Page | 10
Chapitre Ⅰ : Bibliographie
I. Généralités
1. Définition
Les barrages sont, par définition, des ouvrages naturels ou artificiels construits sur le lit d’un
cours d’eau afin de barrer le passage de l’eau sur toute la largeur d’une section de la vallée,
créant ainsi une cuvette artificielle géologiquement étanche.
2. Rôle :
Par la construction de barrages, l’homme influence de manière prépondérante l’écoulement
naturel des eaux de ruissellement. trois raisons principales peuvent justifier cette intervention :
La création d’une retenue : Selon le volume utile de la retenue et le débit des apports, on
distinguera les accumulations journalières, hebdomadaires, saisonnières ou inter-saisonnières.
La régulation des apports : Dans la plupart des régions du monde, les précipitations sont
concentrées sur des périodes courtes. Ces apports sont souvent très irréguliers d’une année à
l’autre. Les besoins en eau sont répartis de manière beaucoup plus homogène sur l’année. Il
s’ensuit donc une succession de périodes de pénurie et d’excès que le seul moyen de compenser
est la réalisation d’une retenue.
La surélévation du plan d’eau d’une rivière : La mise en place d’un barrage en travers d’un
cours d’eau a pour effet de surélever le plan d’eau à l’amont. Cet effet est bien entendu utilisé
pour la production hydroélectrique, mais également pour gérer la dérivation des eaux d’une
rivière vers une prise d’eau puis un canal d’amenée pour l’irrigation ou l’alimentation en eau.
Page | 11
Figure 1 : Le rôle du barrage
L’étanchéité propre du barrage qui est liée à sa constitution, et donc peut aisément être
maîtrisée et contrôlée ;
L’étanchéité du bassin de retenue qui dépend de la géologie du site (nature des terrains, état
de fracturation, réseaux karstiques éventuels, ...), elle est établie lors du choix du site et peut
éventuellement être traitée sur des zones de faible étendue ;
L’étanchéité de la liaison barrage - sol de fondation ; la zone de pied de barrage est celle des
plus fortes pressions, et le terrain peut être en partie remanié lors de la construction de
l'ouvrage, cette zone doit donc être contrôlée et traitée soigneusement.
b. Stabilité :
Les barrages sont soumis aux efforts liés à l'action de l'eau :
On doit prendre en compte le poids propre du barrage et les actions de liaison du sol de
fondation. On devra vérifier :
Page | 12
La stabilité d'ensemble de l'aménagement (barrage et massif de fondation) qui dépend
des qualités du massif de fondation ;
La stabilité propre du barrage sous l'ensemble des actions extérieures ;
La stabilité interne du barrage sous les sollicitations.
4. Ouvrages hydrauliques
a. Ouvrage de sécurité et d'entretien :
Il s'agit de prises d'eau commandées par des vannes et protégées par grilles permettant
d'évacuer rapidement tout ou partie de la retenue. Le problème est d'éviter les dégradations,
dues à la puissance de l'écoulement, aussi bien pour l'ouvrage que pour les rives en aval de
l'ouvrage.
Évacuateurs de crue :
Il s'agit de trop plein permettant de limiter le niveau d'eau dans le barrage lorsque le débit
d'alimentation devient trop important. L'évacuation se fait par le haut du barrage soit le long
d'un canal à l'air libre où une partie de la puissance sera consommée par les remous, soit par-
dessus le barrage où une partie de l'énergie est absorbée par frottement et pulvérisation dans
l'air.
Vidanges de fonds :
Il s'agit de prises d'eau situées au point bas du barrage permettant la vidange pour l'entretien.
Les conduites traversent le plus souvent le barrage. Pour éviter la cavitation, qui entraînerait
Page | 13
une érosion importante des conduites, on utilise des formes très progressives et des parois les
plus lisses possibles (blindages métalliques). En sortie, on cherche à réduire la vitesse de
l'écoulement par élargissement des conduites et éventuellement projection dans l'air.
b. Ouvrages hydroélectriques :
Les équipements hydroélectriques peuvent se décomposer en 3 parties :
Les prises d'eau doivent être convenablement protégées (criblage, dessablage) pour éviter la
dégradation des équipements.
c. Ouvrages de circulation :
Les barrages sont équipés de galeries permettant les accès aux équipements hydrauliques et
électriques et permettant l'auscultation et le contrôle du barrage.
Pour permettre la navigation fluviale, des écluses ou canaux de dérivation peuvent être
aménagés (exceptionnellement on peut utiliser des ascenseurs à bateaux pour de fortes
dénivellations).
Pour favoriser la circulation des poissons migrateurs (saumons), on équipe les barrages de
passes à poissons (canaux à faible vitesse d'écoulement).
Page | 14
Il est très fréquent que les barrages soient utilisés comme voie de circulation en crête, on
équipe donc le sommet d'une chaussée et des dispositifs usuels de sécurité.
6. Types de barrages :
Le principe général d'un barrage est de barrer une vallée, plus ou moins encaissée, pour
accumuler l'eau. Cette accumulation peut avoir plusieurs objectifs. Leurs caractéristiques
géométriques et le principe de leur conception permettent aussi de regrouper les barrages en
grandes familles. On distingue donc, de façon schématique, divers types de barrages : Barrage
en béton :
Suivant leur fonction, sachant qu'un même ouvrage sert le plus souvent à plusieurs
fonctions :
Page | 15
Barrages réservoirs, de grande hauteur ;
Barrages au fil de l'eau, sur le cours d'une rivière.
Barrages poids, pour lesquels, comme leur nom l'indique, c'est la masse de l'ouvrage qui
s'oppose à la pression de l'eau ;
Barrages voutes, pour lesquels la pression de l'eau est transmise par l'ouvrage sur les
parois latérales.
Homogène,
À noyau,
À masque.
Barrage homogène
Un barrage en terre est dit homogène lorsqu’il est constitué d’un même matériau à dominante
argileuse, relativement imperméable. Selon les ouvrages, la pente des talus sera plus ou moins
forte, en fonction notamment des caractéristiques du matériau employé.
Page | 16
Figure 5 : Barrage en terre compactée, homogène
Barrage à noyau
Dans un barrage à noyau, les fonctions de résistance et d’étanchéité sont en quelque sorte
séparées. La résistance est assurée par les recharges placées sur les flancs de l’ouvrage, et
l’imperméabilité par le noyau central. Le noyau au centre de l’ouvrage va être constitué de la
terre la plus imperméable possible. Il sera tenu de part et d’autre par des recharges composées,
selon les cas, de terre plus perméable, d’alluvions ou d’enrochements.
Barrage à masque
Il peut aussi exister des sites où aucune terre n’est disponible, mais seulement des
enrochements. Ceux-ci sont alors employés pour réaliser le corps du barrage, tandis que
l’étanchéité est assurée par un masque de béton, ciment ou béton bitumineux posé sur l’ouvrage
lui-même, côté amont.
Page | 17
Figure 7 : Barrage à masque
Barrage poids
Comme son nom l’indique, ce type de barrage oppose son poids à l’eau pour la retenir. En
fonction des propriétés de résistance du matériau, la forme triangulaire à l’aval de l’ouvrage
s’est peu à peu imposée.
Barrage voûte
Page | 18
l’arc central sur les rives, permettant de construire des ouvrages moins volumineux, à
performance égale. En revanche, les fondations, sur lesquelles se reporte une grande partie de
l’effort, doivent posséder des caractéristiques mécaniques élevées afin de supporter celui-ci.
Les autres formes de barrages voûtes sont des variantes des deux premiers types. Le barrage à
contrefort est ainsi un barrage poids allégé de l’intérieur. Le barrage à voûtes multiples, quant à
lui, repose sur des appuis intermédiaires, lorsqu’une trop grande portée sépare les deux rives
d’un cours d’eau.
Page | 19
minimisant les impacts causés par l’ouvrage et le chantier. Les critères principaux à prendre en
compte sont :
a. La forme de la vallée :
La géométrie de la vallée permet d’exclure d’emblée certains types de barrages :
Vallée étroite en V :
Le barrage voute peut être envisagé si la géologie et l’intégration des ouvrages annexes le
permettent D’autres critères topographiques doivent être vérifiés :
L’élancement ;
Les courbes de niveau du rocher d’appui doivent être parallèles à l’axe de la vallée,
voire légèrement convergentes ;
La géologie doit être adéquate (résistance, fracturation, pendages).
Page | 20
Figure 12 : Vallée étroite en V
Vallée étroite en U :
L’élancement
Les courbes de niveau du rocher d’appui doivent être parallèles à l’axe de vallée, voire
légèrement convergentes
La géologie doit être adéquate (résistance, fracturation, pendages)
Vallée large :
L’élancement du barrage est très important, éliminant le barrage-voûte. Tous les autres types
de barrages peuvent être envisagés pour autant que les autres critères de choix soient satisfaits.
Page | 21
Figure 14 : Vallée large
b. La séismicité :
Les types de barrages les plus résistants aux sollicitations dynamiques sont :
Les barrages à contreforts ne supportent que peu les sollicitations transversales, par la forme
même de la structure.
c. La géologie :
On peut classer la géologie d’un site dans les catégories suivantes :
Page | 22
Un rocher de moyenne qualité permet tous les types de barrages sauf les barrages –voûtes pour
lesquels les contraintes à la fondation sont trop importantes. La tenue du rocher doit être
étudiée minutieusement pour vérifier que les déformations de fondation seront supportées par le
barrage. Un comportement homogène de la fondation est très important.
La déformabilité du rocher est trop grande pour y fonder une structure rigide comme un barrage
en béton. On préférera sur ce type de terrain un barrage en remblai, à noyau ou éventuellement
à masque amont.
Sol meuble :
Seuls les barrages en remblai à noyau peuvent supporter la déformabilité de tels sols. Une
attention particulière sera portée à la continuité de l’étanchéité de la retenue dans la fondation
sous le barrage.
Page | 23
est directement dépendant de cette problématique. Les barrages en béton peuvent supporter
sans dommages excessifs un éventuel dépassement des débits de dimensionnement et donc un
déversement par-dessus le couronnement. En revanche le déversement par-dessus le
couronnement des barrages en remblai serait catastrophique et pourrait conduire à la ruine de
l’ouvrage et à des dégâts considérables à l’aval.
Page | 24
Chapitre Ⅱ : Présentation d’un barrage en remblai
Concepts généraux
Objectif
Collecte des données géographiques, hydrologiques, et topographiques
Données Utilisations
Les barrages en remblai, utilisés comme ouvrage de protection contre les crues le long des
rivières ou comme barrages d’accumulation, comptent parmi les ouvrages de génie les plus
anciens dans l’histoire de l’homme. En effet, on retrouve les traces de tels barrages vieux de
plus de 2000 ans Inde, au Sri Lanka, au Yémen et en Chine.
Bien que des milliers de barrages en remblai aient été construits au cours des siècles avant
l’époque de l’industrialisation, les méthodes analytiques pour les problèmes de stabilité sont
très récentes.
II.Présentation du projet :
1. Destinations principales du barrage :
Le site du barrage KHARROUB se situe sur l’oued KHARROUB dans la Wilaya de Tanger à
environ 22 Km à l’Est de la ville d’ASILAH, et à environ 9 Km à L’est du village ARBA
AYACHA.
Page | 26
- Maître d’œuvre : Direction Des Aménagements Hydrauliques
- Province : Tanger
c. Barrage
- Type : Digue homogène à noyau central argileux
- Pentes du parement amont : 4.8 H/1V -Pentes du parement aval : 4.5 H/1V
d. Dérivation provisoire
- Fonctions : Dérivation provisoire en phase de dérivation, vidange en phase finale,
- Description : Galerie
- Type : circulaire
- Diamètre : 6.5 m
Page | 27
- Cote aval : 41 NGM
e. Evacuateurs de crue
- Type : Tulipe en rive gauche (Rs=8.5 m) ;
f. Vidange de fond :
- Type : Galerie de 6.5 m de diamètre ;
g. Digue de col
- Type : Digue homogène à noyau central argileux
- Hauteur maximale : 30 m
Page | 28
- Volume de remblai : 815 000 m³
3. Situation géographique :
Le site du barrage KHARROUB se situe sur l’oued KHARROUB dans la Wilaya de Tanger à
environ 22 Km à l’Est de la ville d’Asilah, et à environ 9 Km à L’est du village ARBA
AYACHA.
Page | 29
Une protection du talus amont, d’une pente de 4,8H/1V, par rip-rap et transition ;
Une protection du talus aval, d’une pente de 4,5H/1V, par des galets
b. Digue de col
Le col, situé en rive gauche à la cote 68 NGM environ, est fermé par une digue homogène à
pseudo-noyau central argileux avec des recharges en limons et tout venant limoneux d’une
longueur en crête de 600 m. La digue de col a une hauteur maximale sur fondation de 30 m
avec un volume de 815.000 m3 de remblai.
c. L’évacuateur de crues
Etant donné l’absence d’une fondation rocheuse en rive, la solution d’un évacuateur à seuil
libre latéral est écartée. On prévoit donc un évacuateur de crues de type tulipe implanté en rive
gauche est constitué de l’amont vers l’aval de :
Un seuil circulaire de 8,5 m de rayon calé à la cote 85 NGM d’une longueur déversant
de 33 mètres,
Un puits de 7,5 m de diamètre et une hauteur de 30 m dont le fond est calé à la cote 54,4
NGM,
Page | 30
Une galerie souterraine de 8,5 m de diamètre et de 480 m de longueur qui débouche à
l’aval à la cote 41 NGM.
L’évacuateur de crues est dimensionné pour faire transiter un débit de 270 m3/s au passage de
la crue de projet de fréquence centennale.
d. Vidange de fond :
La galerie de dérivation provisoire sera transformée en fin de chantier en vidange de fond
permettant d’évacuer un débit de 270 m3/s sous la cote de la retenue normale. Elle est
constituée d’un pertuis équipé d’une vanne de garde de type wagon à brimballes manœuvrable
à l'aide d'un servomoteur à simple effet disposé à la chambre de commande, et d’une vanne de
réglage de type segment manœuvrable à l'aide d'un servomoteur à simple effet disposé dans la
chambre de visite à la cote 57,50 NGM.
Le puits de vannage est d’une conception classique et regroupe la vidange et les prises AEP.
Les prises alimentent un collecteur vertical qui vient rejoindre une conduite placée en calotte
de la galerie de vidange de fond.
Page | 31
Figure 20 : Coupe type de prises d'eau potable
Figure 21: Le bassin Tahaddart au nord du Maroc dans un extrait simplifié de la carte géologique du Maroc à
1/1000.000 « Trait noir : limites du bassin »
Les zones internes situées à l’intérieur de l’arc structural du Rif entre Sebta et Jebha et
dans les Bokoyas. Elles se composent des sébtides, ghomarides et dorsale calcaire.
Le domaine des nappes des flyschs jalonne l’arc rifain à l’ouest et se compose des
nappes Tisirène, Beni Ider, Chouamat – Melloussa et Numidienne.
Page | 32
Le domaine externe, où se situe la zone du projet, est formé d’épaisses séries méso-
cénozoïques et envahit par certains lambeaux des nappes de flyschs.
Le Rif interne se compose d’un ensemble d’unités allochtones alpines, empilées les unes
sur les autres. Du Nord -Est au Sud -Ouest, les affleurements sont successivement comme suit :
Les Sebtides :
Les Ghomarides :
Elles sont formées d’une série terrigène à sédiments fins schisto-gréseux et grossiers
conglomératiques, de roches détritiques immatures et de flyschs carbonatés épais du
Paléozoïque inferieur et moyen. Le Paléozoïque supérieur est composé d’une série turbiditique
épaisse. Les Ghomarides se terminent par des séries triasiques détritiques et grossières et des
dolomies, reposant en discordance sur le Paléozoïque plissé.
La Chaîne calcaire :
Elle est prédominée par des dépôts carbonatés et dolomitiques, séparés par des accidents
transverses. Ces dépôts carbonatés constituent la source potentielle des matériaux de
construction en agrégats et enrochements dans la région.
Les unités Tisirène, Beni Ider, Chouamat – Melloussa et la nappe Numidienne sont formées de
séries sédimentaires détritiques.
La Nappe de Tisirène :
Page | 33
Un flysch du Crétacé inferieur formé d’une alternance répétée de bancs de grès
verdâtres, finement micacés et arkosiques et d’argiles marneuses ou pélitiques. Des
niveaux de grès et quelques lits de micro-brèches calcaires se rencontrent aussi.
La Nappe de Béni Ider :
La Nappe Numidienne :
Elle est constituée principalement par un flysch gréseux grossier d’âge Oligocène -
Miocène inférieur et se termine par des argiles versicolores à nodules. Le flysch oligomiocène
est partout composé de gros bancs de grés et de fins lits argilitiques. Les dépôts issus des
flyschs gréseux intéressent la zone du site.
c) Le Domaine externe :
Il comprend l’unité intrarifaine de Ketama - Tanger, les unités mésorifaines, les unités
pré-rifaines et nappes rifaines et l’avant-pays plissé de la chaine rifaine.
La disposition actuelle très complexe des unités du Rif externe résulte de phases
tectoniques qui se sont succédé depuis le paroxysme miocène.
Page | 34
L’unité de Ketama est constituée surtout de terrains mésozoïques allant du Trias au
Crétacé moyen. Ces unités constituent l’autochtone relatif à l’ensemble des nappes de flysch
précédemment décrites.
Unités mésorifaines :
Elles correspondent à des formations complexes où les séries sont caractérisées par des
faciès à prédominance marneuse dès le Crétacé. Ces formations constituent un complexe
tectono-sédimentaire d’âge Miocène supérieur. Elles sont considérées comme un méga-
olistostrome dont la matrice essentiellement marneuse est d’âge Tortonien supérieur à
Messinien.
Contexte structural
Ainsi cette phase majeure dont le paroxysme est daté du Miocène supérieur se traduit
par :
Le charriage des zones internes sur la zone des flyschs et le domaine externe ;
Page | 35
La mise en place de la nappe prérifaine ;
Dans leur ensemble, les événements miocènes amènent l’empilement des écailles et des
nappes dans les zones internes avec remobilisation du socle. Ces structures sont à leur tour
scellées par le Miocène supérieur « post – nappes » composé d’un ensemble de terrains eux-
mêmes déposés en régime compressif.
1.2.1 Stratigraphie
Les terrains qui affleurent dans la zone s’étalent du crétacé au quaternaire. Le Crétacé est
constitué de pélites grisâtres à rares intercalations de calcaire gréseux comme indiqué sur le log
lithostratigraphique suivant :
Page | 36
Figure 22 : Log lithostratigraphique de la zone d’étude réalisé par le département
Une tranchée a été réalisée dans la rive gauche et a montré que le calcaire gréseux a un
pendage moyen de 80°.
Page | 37
Figure 23 : Relevé géologique de la tranchée réalisée en rive gauche
Dans la rive gauche, l’étude géologique réalisée par une campagne de sondage atteignant
une profondeur de 60 m a révélé trois types de pélites d’âge Crétacé données comme suit :
Il convient de signaler que la rive gauche de l’oued Kharroub abritera les ouvrages
annexes tels que l’évacuateur de crues et la dérivation provisoire qui sera équipée en vidange
de fond. Ces ouvrages souterrains seront réalisés dans la série pélitique sur presque leur
longueur totale à l’exception de leurs sorties du côté aval et quelques zones très localisées qui
traverseront les formations du calcaire gréseux intercalées avec des pélites.
1.2.2 Fracturation :
Les travaux antérieurs montrent que la région de l’oued Kharroub a été affectée par 4
familles de fractures :
Page | 38
La 1ère famille F1 : c’est la famille la plus dominante, elle montre un plongement fort
vers un azimut compris entre 285 et 295° ;
La 2ème famille F2 : il s’agit d’une famille dont le plongement est fort à subvertical,
montrant un azimut variant entre 5 et 30° ;
La 4ème famille F4 : il s’agit d’une famille dont le plongement est fort à subvertical et
montrant un azimut variant entre 145 et 155°.
2. Le relief
Les reliefs présentent au niveau du bassin de tahhddart ont des niveaux d’altitudes plus hautes.
Ils sont essentiellement au SE du bassin, tel que (1065m Jebel Ez Zaouia, 1020m Sehan Daria),
ainsi que les basses altitudes de 50 a228 m.
Page | 39
Le bassin caractérisé vers l’est par des pentes plus raides dans la partie amont (orientale et sud
orientale), et par des pentes faibles à l’aval (occidentale et nord occidentale). Néo moins les
zones à pentes raides sont orientées selon les directions privilèges ; les directions des
principales failles en particulier les directions NW_SE du bassin, dans la bordure orientale du
bassin et le long de la bordure sud de la nappe numidienne.
Figure 25 : Carte des pentes et failles du bassin de Tahaddart (Réf : science Lib Editions Mersenne)
3. Topographie
L’étude topographique de l’environnement du barrage s’est basée sur plusieurs documents
topographiques dont :
Ces cartes nous ont aidées aussi pour l’étude topographique afin de caler des structures de
dérivation provisoire.
Page | 40
Les caractéristiques géomorphologiques obtenues à partir de l’analyse des levés topographiques
au 1/50000ème et 1/1000ème réalisés au près du site sont les suivantes :
4. Réseau hydrographique
Le bassin de l’oued EL HACHEF totalise un bassin versant de 600 Km². Il est composé de
deux affluents principaux, oued EL KHOBZ qui draine un bassin versant de 242 Km² au niveau
du barrage de 9 avril 47 occupant la moitié Nord et l’oued KHARROUB qui draine un bassin
versant de 314 Km², les 44 Km² restant sont en zone de plaine à l’aval de la confluence. Les
oueds se rencontrent à environ 5 Km à l’aval du barrage de 9 avril 47.
Page | 41
Le site du barrage objet de notre projet est situé sur l’oued KHARROUB, contrôlant un bassin
versant de 187 Km².
Le site projeté contrôle un bassin versant de 191 Km2 qui culmine à 1065 NGM. L’oued prend
naissance à 16.5 Km environ en amont de la zone du projet, avec une dénivelée de 1025 m
entre le début du réseau hydrographique et la zone du site du barrage (exutoire). L’oued
Kharroub en traversant ce bassin versant atteint une pente de 2 % en moyenne.
La majeure partie du bassin versant est située dans le domaine externe, notamment l’unité de
Tanger externe, formée de marne et marno-calcaire attribués au Crétacé supérieur. Vers les
bordures du bassin, les nappes de flyschs sont représentées surtout par la nappe du grès
Numidien d’âge Oligocène-Aquitanien et un seul affleurement de la nappe de Béni Ider, d’âge
Eocène-Oligocène dans l’extrémité Nord de ce bassin versant. Les terrains de couverture se
résument à des faciès dunaires du Villafranchien, des éboulis de pente et des alluvions du
Quaternaire.
Page | 42
a. Hydrologie
Les résultats de l’étude hydrologique des séries d’apport mensuel et annuel, entre 1945 et 2002,
de la station de JBEL LAHBIB et du barrage 9 avril 1947, ont permis l’extrapolation des séries
d’apport du barrage KHARROUB. Les résultats obtenus sont les suivants :
Le graphe figure suivante illustre l’évolution des apports au niveau du barrage KHARROUB.
Figure 27 : Les cours d'eau (classe1) du bassin de Tahaddart. (Réf : Science Lib Editions Mersenne)
Page | 43
Figure 28 : Evolution des apports du site du barrage KHARROUB
b. Hydrologie de surface :
L'objectif d'un barrage étant de stocker un certain volume d'eau pour diverses utilisations, il
faut donc se préoccuper, d'une part, des conditions de remplissage de la réserve envisagée et,
d'autre part, de veiller à ce que l'ouvrage ne soit pas un obstacle au passage des crues qui
Page | 44
risqueraient de le submerger. Pour cela, il faut étudier les apports mensuels, saisonniers et
annuels et les débits instantanés pour définir les crues maximales pouvant entrer dans la
retenue. Si les données sont suffisantes, on utilise des méthodes statistiques directes. Si cela
n'est pas le cas, il faut étendre l'information par corrélation avec d'autres données disponibles
comme la pluviométrie ou le régime de bassins similaires. Toutes ces études demandent de
grandes précautions pour couvrir les incertitudes des données.
Pour étudier le comportement de ce bassin versant, on doit faire une étude hydrologique
concernant les caractéristiques physiques et morphologiques, ainsi une étude climatique,
concernant les précipitations et le mode du climat de la région par le bilan hydrique.
Pour une étude hydrologique, il est nécessaire de déterminer certains paramètres physiques, tels
que :
La pente : ce paramètre influence sur le débit de la retenue d’eau vers l’exutoire. Ci-dessous
une assimilation de topographie du bassin versant.
Ainsi, lorsque le temps de concentration est atteint, ceci veut dire que toutes les régions du
bassin participent au débit.
Pour le cas du bassin versant d'oued Kharroub, la topographie n'est pas complexe, dont la
retenue calculée par le bureau d'étude est très importante, ceci confirme que le temps de séjour
des gouttes de précipitation est très faible.
Page | 45
Ci-dessous une carte topographique simplifié, montre une basse altitude des reliefs constituent
ce bassin.
(Réf : http://fr-fr.topographic-map.com)
Pour le cas du bassin versant Kharroub, le sol est majoritairement argileux, alors le
ruissellement est très important, donc une grande retenue d’eau superficielle, ainsi l’infiltration
est très faible à cause de la granulométrie des argiles ( inf à 2µm).
Figure 31 : image représentative de la structure du sol argileux (Réf : prise de photo sur chantier)
Bilan hydraulique : Les apports annuels du bassin versant Kharroub Selon les équations
suivantes, on peut calculer les apports annuels de ce bassin :
A=ETR+R+Ie et Ie=Qex+Dr+F
A : apport pluviométrique ;
ETR : évapotranspiration ;
Page | 46
R : Ruissellement ;
Ie : l’infiltration efficace ;
Dr : drainage souterrain.
(Réf : accuweather.com)
Les matériaux de construction des barrages en remblais, sont réalisés généralement par des
matériaux naturels. Principalement par des composés situés au voisinage du cours d’eaux, dont
les matériaux grossiers s’installent toujours au lit de l’oued par la force de gravité et paramètre
de la pente, cependant les autres matériaux tels que les sables fins sont fabriqués dans la STM
de l’entreprise pilote du chantier, alors que les géo ressources de ces matériaux sont liés
toujours par le fait de l’étude des géo-matériaux par des spécialistes en géotechniques.
En outre, les barrages en remblais ont une forte emprise au sol qui leur permet d’être construits
sur des fondations de qualité médiocre. Ils résistent à la pression de l’eau grâce à leur poids.
Trois catégories principales de barrages en remblais existent. Les barrages en remblais sont des
barrages constitués d’un matériau meuble qui peut être fin (terre) ou grossier (enrochements).
De forme trapézoïdale, les barrages en remblais sont des barrages à forte emprise au sol. Grâce
à leur forte emprise au sol, ils peuvent reposer sur des fondations de qualité moyenne, là où un
barrage en béton ne pourrait pas être construit.
Page | 47
Les barrages en remblais résistent à la pression de l’eau grâce à leur poids, une caractéristique
qu’il partage avec les barrages poids.
La pente des barrages en remblais est généralement faible, mais elle peut évoluer en fonction
du matériau utilisé.
Il existe trois types de barrages en remblais, qui diffèrent en fonction des matériaux utilisés et
de la méthode utilisée pour assurer l’étanchéité du barrage.
Le barrage homogène est construit avec un seul matériau étanche l’argile ou le limon. Il s’agit
du type de barrage le plus ancien.
Qualité :
Non organique ;
Non altérable ;
Extraction transport et lise en place possible ;
Compactage possible ;
Résistance au cisaillement ɸ et cohésion c élevées.
Economie.
Les exigences sur les propriétés physiques des matériaux dans les différentes zones du barrage
sont présentées sur la figure ci-après :
Page | 48
Figure 32 : Les propriétés physiques des matériaux dans les différentes zones du barrage.
Pour garantir une faible perméabilité, la teneur en argile doit être en minimum 5% en cas des
limons argileux. Souvent des matériaux morainiques sont utilisés à la place des limons argileux.
Ces matériaux sont utilisables au cas où ils se composent de plus de 15% de limons et plus de
2% d’argiles.
b. Corps d’appui :
Comme matériaux des enrochements non cohésifs de provenance des :
- Formations alluvionnaires
- Moraines
- Carrières
Sont utilisés. Une granulométrie étendue des limons aux graviers est favorable. Le diamètre
maximum est limité par l’épaisseur des couches de compactage :
Page | 49
Les fuseaux granulométriques du matériau du noyau A (cohésifs) et du matériau des corps
d’appui B (non-cohésifs) sont représentés sur la figure 30.
c. Filtres
Le rôle du filtre est d’empêcher :
Un filtre est toujours inévitable comme zone de transition entre des matériaux cohésifs et
non cohésifs. Pour éviter la migration des grains fins vers la zone grossière, des critères de
construction du filtre sont pris en considération. Les plus couramment utilisés ont été définis
par Terzaghi-peck comme suit :
Page | 50
Page | 51
Chapitre Ⅲ : Etude de stabilité du barrage kharroub
Concepts généraux
Objectif
Sollicitations et approches de calcules
Données Utilisations
Page | 52
I. Etude de stabilité d’un barrage en remblai
1. Cas de charges
On distingue entre :
Pendant la construction :
Mise en charge de la digue pendant les crues (la hauteur finale de la digue n’est pas encore
atteinte). (CS)
Fin de la construction :
Tremblement de terre à lac vide : les pressions interstitielles dans le noyau sont maximales au
moment où les remblais viennent d’être achevés. (CE)
Page | 53
Selon DIN 4084 les coefficients de sécurité doivent être observés :
Le facteur de sécurité est défini dans la plupart des méthodes de tranches par le rapport :
Forces stabilisantes
FS=
Forces Motrices
Les valeurs des matériaux à la fin de construction doivent être utilisées en considérant la
variation des essais (+- écart type)
2. Principe d’analyse :
On calcule la sécurité au glissement pour diverses surfaces de rupture qui traversent plusieurs
zones de la digue et éventuellement le sous-sol. Pour aboutir à un résultat à savoir calculer le
facteur de sécurité qui assure la stabilité requise quelques hypothèses sont à respecter :
T =C+(n−ρ)tg ϕ
Avec :
On contrôle plusieurs surfaces de rupture orientées soit vers l’amont, soit vers l’aval qui
traversent
Page | 54
3. Méthodes d’analyse :
Pour les massifs hétérogènes comme les digues zonées on utilise les méthodes analytiques bien
connues de la mécanique des sols notamment celle d’équilibre limites.
Avec :
h : hauteur
b : épaisseur de la tranche (𝑏𝑏 = 𝑙𝑙. 𝑐os(α))
W : poids total de la tranche
𝑁 : Force normale à la base de la tranche
𝑇 : Force de cisaillement à la base de la tranche
XL, XR : Force horizontale inter-tranches
Page | 55
EL, ER : Force verticale inter-tranches
α : angle que fait la base de la tranche avec l’horizontale
R : rayon du cercle de rupture de centre O
L : longueur du plan de glissement de la tranche
X : bras de levier du poids des terres
i. La méthode de Fellenius :
C’est la plus simple pour l’analyse de stabilité des talus. Fellenius suppose que le volume de
glissement est délimité par la surface de glissement et la topographie du talus est subdivisée en
n tranches. Chaque tranche est considérée comme un solide indéformable, en équilibre sur la
ligne de glissement. Considérons un talus recoupant un certain nombre de couches de sols de
caractéristiques différentes i. La stabilité est étudiée en considérant le problème 2D, c'est-à-dire
en analysant l'équilibre d'une masse de sol d'épaisseur unité dans le sens perpendiculaire à la
figure.
La force Figure 36 : Ensemble des forces appliquées sur une tranche (Méthode ordinaire)
Page | 56
susceptibles de mobiliser un frottement. N est donc la composante normale du poids
immergé c.à.d. le poids calculé à partir de la densité saturée diminuée de 1 appelée N’, R le
rayon du cercle de glissement.
Cette méthode calcule le facteur de sécurité par la vérification de l’équilibre global des
moments tout en négligeant les forces entre tranches.
Hypothèses :
Figure 37 : Ensemble des forces appliquées sur une tranche (Méthode de Bishop simplifiée)
Cette méthode
calcule le facteur de
sécurité par la vérification de l’équilibre global des moments ainsi que l’équilibre global des
forces verticales pour chaque tranche (i). Cependant, elle néglige l’équilibre horizontal des
forces.
Page | 57
En conséquence,
1
∑ [ c . lcos α + (W −u . lcos α ) tan φ ]
' '
m∞
F=
Wsin α
Avec : Où,
m∞=cos α ¿ )
Il est remarqué que le facteur de sécurité (Fm) se trouve dans les deux parties de l’équation ci-
dessus, par conséquent, la résolution doit être passée par la technique itérative. Le point de
départ préféré pour cette itération ; est celui du coefficient de sécurité estimé par la méthode
ordinaire (Fellenius) puis le résultat est introduit dans le calcul par Bishop ce qui donne une
nouvelle valeur plus fine.
La méthode de Bishop étant plus rigoureuse que Fellenius nécessite plus de calcul (3 à 4
itérations). Les coefficients de sécurité obtenus sont en général un peu plus élevés. Le plus
souvent pour éviter d’alourdir exagérément la masse de calcul, on détermine le cercle critique
par la méthode de Fellenius et on vérifie que la valeur du facteur de sécurité calculé par la
méthode de Bishop est supérieure à la valeur du coefficient de sécurité calculé par la méthode
de Fellenius. Si ce n’est pas le cas, il faut reprendre la recherche du cercle critique par la
méthode de Bishop.
Avec :
Page | 58
La figure ci-dessous montre les fonctions typiques de f(x). Dans la première itération, les
forces de cisaillement verticales sont égales à 0. Les coefficients de sécurité sont obtenus par la
détermination de deux séries de coefficients de sécurité pour diverses valeurs de λ. La première
série correspond à l’équilibre des moments et la deuxième à l’équilibre des forces. Ces
coefficients de sécurité sont mis en graphique en fonction de x par rapport à chaque λ. Le point
d’intersection satisfait les deux conditions d’équilibre. Cette méthode est précise, mais
demande une certaine pratique pour donner le bon f(x).
A l’aide des logiciels, les surfaces de glissement critiques peuvent être trouvées assez
facilement.
Page | 59
glissement des surfaces de rupture critique. En pratique, on accepte des déformations et
tassements limités dues aux instabilités pendant des tremblements de terre.
Matériaux de remblai
Element étanche Pente amont Pente aval
(Corps d’appui)
Page | 60
Une fois la géométrie du modèle du référence définie, le matériau géologique est mis en place.
Les caractéristiques géotechniques spécifiques aux calculs sont regroupées dans le tableau
suivant :
Avec :
- Les actions considérées sont le poids du massif P et le chargement induit par les secousses
sismiques (effort pseudo – statiques KH et ± KV) on base sur la réglementation parasismique
marocain.
c. Procédure de calcul :
Le but de la série des calculs effectués est de déterminer les valeurs du coefficient de sécurité
au glissement du talus par l’approche pseudo-statique par la méthode " Phi-c réduction " à
l’aide du code PLAXIS).
Page | 61
- La seconde combinaison dans lequel l’effet de la nappe, dont le niveau est supposé
atteindre le terrain naturel (massif entièrement saturé), est pris en compte.
- La troisième combinaison dans lequel seul l’effet du séisme, avec une accélération
horizontale de référence égale à 0.10g (g=10 m/s² étant l’accélération de la pesanteur),
est pris en compte.
- La quatrième combinaison dans lequel les effets combinés de la nappe (saturation
Complète) et du séisme (accélération horizontale de référence égale à 0.10g) sont pris
en compte.
L’état initial est caractérisé par des contraintes de type géostatique, le calcul du modèle de
référence considérées se fait définie en deus phase :
L’exploitation détaillée des résultats des calculs sont représentés sur les figures au-dessous.
La figure 45 représente le maillage déformé du talus, les déplacements totaux (figure 46),
et les contraintes de cisaillement (figure 47).
Page | 62
Dans le présent travail, seul SLOPE/W parmi tous ces programmes sera utilisé. Le
programme SLOPE/W est l'un des programmes intégrés dans la fenêtre du Géo Studio et qui
est basé sur la théorie de l'équilibre limite dans le calcul du coefficient de sécurité des pentes.
L'élaboration de ce programme permet d'analyser facilement à la fois des problèmes simples et
complexes de stabilité de pente en utilisant une variété de méthodes pour calculer le coefficient
de sécurité. Ce programme est de type graphique et compatible avec les versions de Microsoft
Windows XP. L'interface du logiciel SLOPE/W : Les menus disponibles et la fonction de
chacun sont les suivants :
- File : Ouvre et enregistre les fichiers, les importations des photos et imprime le dessin.
- Edit : permet de copier le dessin dans la Presse-papiers du Windows.
- Set : définir les paramètres de la grille, le zoom, les axes et de la page.
- View (aperçu) : Options d'affichage des contrôles, des informations des sols et des
points, et nous permet aussi d'afficher les forces inter-tranche sur un schéma corporel de
chaque tranche et le polygone force.
- Keyin (menu) : Permet d'introduire automatiquement à l'aide de tableaux les
paramètres géométriques de la pente (coordonnées et couches), les caractéristiques
mécaniques des différentes couches constituant le talus, les conditions complémentaires
nécessaires pour le calcul, le niveau de la nappe, la sélection des différentes méthodes
de calcul, les surcharges, l'effort sismique et les confortements.
- Draw (dessin) : permet d'introduire manuellement les paramètres géométriques de la
pente (coordonnées et couches), des conditions complémentaires nécessaires pour le
calcul, le niveau de la nappe, les surcharges et les confortements.
- Sketch (croquis) : permet de définir les objets graphiques, d'étiqueter, d'améliorer et de
clarifier les résultats du problème.
- Modify (modifier) : Permet de déplacer, de modifier ou de supprimer les objets
graphiques, les textes et les photos.
Page | 63
des parements amont et aval du barrage principal. Les calculs de stabilité ont été effectués en
considérant des coupes perpendiculaires aux axes des digues, correspondant aux hauteurs
maximales en fond de vallée, tout en se basant sur les données géomécaniques et en utilisant les
méthodes des équilibres limites. Dans ce qui suit les résultats seront donnés par la méthode de
Morgenstern-Price. La méthode de rupture choisie est la rupture circulaire.
a. Paramètres géotechniques :
Les matériaux de fondation du site de la digue principale Kharroub sont dominés de pélites
altérée à rares blocs. La direction globale de la stratification des couches du substratum est rive-
rive, avec un pendage moyen à fort généralement vers l’amont.
Le fond de vallée est tapissé d’alluvions grossières à éléments de grés numidien au-
dessus du substratum pélitique.
La zone de la digue de col située au large de la rive gauche, montre un substratum
pélitique remanié en surface.
Plusieurs essais ont été réalisés pour la reconnaissance des caractéristiques
géotechniques des matériaux de la fondation. En effet, des sondages de
reconnaissances ont été réalisés le long de l’axe de la digue
Densité : une moyenne de 2.45 (t/m3), avec un minimum de 2.22 (t/m3) et un maximum
de 2.64 (t/m3).
Angle de frottement résiduel : une moyenne de 22.33° avec un minimum de 19° et un
maximum de 28°.
Cohésion résiduelle : une moyenne de 11.33 kPa avec un minimum de 8 kPa et un
maximum de 16°.
Dans la suite des calculs et compte tenu de l’aspect hétérogène de la fondation, on procèdera à
une étude de sensibilité portant sur l’angle de frottement et la cohésion de la fondation. Les
valeurs suivantes sont retenues pour cette étude :
Page | 64
Quant aux matériaux de construction nous retenons les valeurs suivantes :
b. Calcul de stabilité :
Profil examiné :
Page | 65
Figure 40 : Coupe géologique du site du barrage Kharroub.
c. Résultats :
Page | 66
Parement Aval (RN=85 NGM)
Cas de charge Moyenne des coefficients de sécurité
φ' Alluvions=25° 2,650
Fin de construction sans séisme
φ' Alluvions=30° 2,780
φ' Alluvions=25° 1,900
Fin de construction avec séisme
φ' Alluvions=30° 1,912
Régime permanent (RN) sans φ' Alluvions=25° 2,000
séisme φ' Alluvions=30° 2,150
Régime permanent (RN) avec φ' Alluvions=25° 1,537
séisme de 0,14g φ' Alluvions=30° 1,904
φ' Alluvions=25° 1,354
Vidange rapide sans séisme
φ' Alluvions=30° 1,867
Vidange rapide avec séisme de φ' Alluvions=25° 1,900
0,14g φ' Alluvions=30° 1,912
Tableau 6 : Les résultats de calcul de stabilité du parement en aval sans et avec séisme de 0,14g
Figure 41 : Résultats du calcul de stabilité du parement amont en régime permanent avec séisme (Digue
principale) :
Résultats du calcul de stabilité du parement aval en régime permanent avec séisme (Digue
principale) :
Figure 42 : Résultats du calcul de stabilité du parement aval en régime permanent avec séisme
(Digue principale) :
Page | 67
d. Interprétation
Parement amont : Les valeurs de sécurité correspondant à un fruit initial de 4.8H/1V,
comportent une certaine marge de sécurité, c’est pour cela il est retenu pour la version finale
des plans de construction, vu la taille du barrage étudiée. L’intervalle de sécurité laissé est une
méthode pour se préserver d’éléments imprévisibles qui pourraient apparaître à l’ouverture des
fouilles.
La stabilité du parement en amont est souvent cruciale car il est exposé aux forces
hydrauliques.
Parement aval : En comparant les résultats avec les facteurs de sécurité requis, il s’avère qu’ils
sont tous acceptables voire laissent une marge grande de sécurité,
La stabilité du parement en aval est également importante pour éviter des défaillances
catastrophiques.
Un sol avec un angle de frottement plus élevé a une résistance interne accrue, ce qui peut
contribuer à une meilleure stabilité.
Page | 68
Figure 43 : Génération de maillage dans le barrage en remblai
Page | 69
Figure 45 : L’effet de la pression interstitielle
Poids volumique :
Page | 70
Le poids volumique du sol (ɣ) est un paramètre géotechnique qui est en fonction de la teneur en
eau qui ne varie pas pendant une secousse séismique, alors comme le montre le figure (4-6), le
Cohésion :
Figure 46 : variation du coefficient de sécurité en fonction du poids volumique
La résistance au cisaillement des sols sous l’effet de chargement est dépendue de la cohésion
comme suit : τf = c’+σ’tanυ’
Donc la variation de la cohésion du sol c à une influence sur la résistance au cisaillement donc
une influence sur le coefficient de sécurité, plus le sol est cohérant plus il résiste à l’effort
séismique.
Une variation de la cohésion de 8kpa fait varier le coefficient de sécurité d’une valeur comprise
entre 15% et 17%.
Une contrainte maximale appliquée sur un échantillon de sol qui peut aboutir à ça rupture est
très faible par rapport au module de Young qui correspond à la contrainte nécessaire afin
d’obtenir un allongement relatif maximal =100%=1.
Donc la variation du module de Young E n’a aucune influence sur le coefficient de sécurité.
Coefficient de poisson ʋ :
Etant donné que le coefficient de poisson est le rapport de déformation principale rapporté à la
déformation de sollicitation, donc la variation du coefficient de poisson ʋ n’a pas une influence
sur le coefficient de sécurité.
Page | 72
Du fait que l'angle de frottement interne " υ " d'un sol est approximativement égal à l'angle du
Talus naturel et tandis que les versants naturels sont plus stables avec un angle de pente naturel
petite donc un coefficient de sécurité plus élevé.
Donc le coefficient de sécurité est influencé par la variation de l’angle de frottement interne υ.
Une variation de l’angle de frottement de 5°, fait varier le coefficient de sécurité d’une valeur
comprise entre 12 et 14%.
Influence de la géométrie :
Page | 73
Figure 52 : Variation du coefficient de sécurité en fonction de la hauteur du talus
Le niveau de la nappe de l’état (hw =0) et totalement saturé (hw= 50m) ont un effet sur la
diminution de la valeur de coefficient de sécurité de 20 à 35%.
c. Conclusion :
On a établi un calcul de référence pour un modèle donné : ce calcul n’est pas la solution au
problème, mais une approximation raisonnable de la réalité. On a fait varier individuellement
Page | 74
plusieurs paramètres géotechniques et de la géométrie du modèle de référence et étudier leurs
influences sur le coefficient de sécurité pour des fourchettes de variation raisonnable.
Quelques remarques :
Cette étude ne prétend pas être exhaustive : on a pu faire varier d’autres paramètres, ou
étudier l’influence des variations de paramètres sur d’autre type de résultats.
L’influence combinée de plusieurs paramètres, qui n’a pas été étudiée ici, peut bien sûr être
encore plus importante.
L’influence par une étude en trois dimensions peut être plus importante Cette étude nous permet
de tirer les deux conclusions suivantes :
Le calcul par méthode numérique (MEF) s’avère indispensable pour les calculs de
précision en condition statique que sismique.
Les paramètres qui exercent une influence remarquée sur la stabilité des terrains en
pente dans le cas d’une étude dynamique sous l’effet d’une excitation sismique sont
celle de cisaillement du sol (cohésion et angle de frottement) et le niveau de la nappe.
Page | 75
Chapitre Ⅲ : Etude de stabilité du barrage kharroub
Concepts généraux
Objectif
Sollicitations et approches de calcules
Données Utilisations
Page | 76
Chapitre Ⅳ : Etude de stabilité du barrage Poids
II.Présentation de l’ouvrage
1. Situation géographique
Le barrage de l’Oued Fodda, appelé aussi "barrage Steeg", est le premier des grands
barrages construits en Algérie. L’étude de faisabilité a débuté en 1910 mais sa construction a
commencé en 1926 pour se terminer en 1932. Cet ouvrage est situé à 20 Km environ de la ville
d’Oued Fodda, dans la Wilaya de Chlef au Nord-Ouest de l’Algérie, au confluent d’Oued
Fodda (d’où le nom du barrage). Le cours de cet Oued se fraye un passage dans une gorge
profonde, à travers du massif calcaire connu sous le nom de "Koudiat Larouah".
L'ouvrage construit pour le compte du service central des irrigations en Algérie barre une
vallée profonde composée de calcaires très fissurés qui ont nécessité pour leur
imperméabilisation l’injection à haute pression de produits chimiques et de ciment. Cette
région de l’Oued Fodda a été affectée par de grands séismes et cette nouvelle technique s’avéra
efficace car plusieurs secousses telluriques n’entraînèrent aucune conséquence pour le barrage.
La cuvette du barrage draine les eaux d’un bassin versant de 800 km2. Sa capacité de 225 Hm3
permettait d’irriguer 18 000 Ha sur la rive gauche de la vallée de Chlef en amont du seuil de
Pontéba dans la vallée de l’Oued Fodda.
Page | 77
2. Géométrie du barrage :
Hammam Grouz se trouve près de la marge sud de la Chaîne Tellienne, une zone dans laquelle
il y’a eu des mouvements d’intensités variables presque continus depuis la mise en dépôt des
roches sédimentaires les plus anciennes visibles sur le site (calcaires crétacés). Des phases
variées de deux épisodes majeurs de formation de montagnes (orogènes des Pyrénées et des
Alpes) ont occasionné des déformations successives des strates anciennes, et ont affecté
l’environnement de la sédimentation. En conséquence, il y’a des variations rapides latérales et
verticales de lithologie des roches, et un changement net de sédimentation marine pendant le
Crétacé jusqu’aux dépôts lacustres terrestres pendant le Tertiaire.
Les calcaires crétacés affleurant au site du barrage plongent obliquement vers l’amont avec une
pente modérée (25-30°, ESE/SE). Ils se trouvent sur le flanc sud-est de l’anticlinal du Hammam
Grouz, un pli subsidiaire sur la structure anticlinale majeure du Djebel Grouz. Ce dernier est
une prolongation vers le sud- ouest, un peu décalée par des failles, du massif calcaire du
Frictia- Garnechou, un des trois anticlinaux majeurs tendant vers le NE qui dominent la
structure géologique locale. Les calcaires crétacés sont recouverts en discordance par des
roches sédimentaires plus jeunes, miocènes et pliocènes, qui occupent la cuvette.
C’était pour pallier un besoin crucial en eau potable de la ville de Constantine qu’a été décidé
la mobilisation des ressources en eau par la réalisation du barrage Hammam-Grouz, dans la
partie amont du bassin du Rhummel. Le site de ce barrage, achevé et réceptionné à la fin des
années 1980, se trouve à 1,5 km au sud d’Oued Athmania à 6 km au nord-est de Chelghoum-
Laid et à 45 km au sud-ouest de Constantine. Il est situé à l’amont d’une gorge incisée par
l’Oued-Rhummel dans les calcaires crétacés qui forment Djebel Grouz. La RN5 longe cette
gorge et traverse la cuvette en rive gauche, le barrage boucle l’Oued Hummel à l’entrée de cette
concave et crée un réservoir d’environ 5 km de long sur 2,5 km de large et 37 m de hauteur
avec une longueur de crête de 300 m. Sa capacité initiale est de 45 millions de m3, quant à la
quantité (volume) stockée actuellement elle ne dépasse guère les 15 millions.
Le barrage principal est construit en béton-poids à travers l'Oued Rhummel, dans une gorge
juste à l'amont de la ville d’Oued Athmania. L'eau devant être stockée dans la retenue est
utilisée principalement pour l'alimentation urbaine et éventuellement pour l'irrigation dans le
constantinois.
Page | 78
3. Profil type du barrage
Le barrage de l’Oued Fodda est un barrage de type poids en béton, il possède les
caractéristiques suivantes :
Barrage
24600 0.2 2640
(béton)
Réservoire
2070 0.49 1000
d’eau
Tableau 7 : Propriétés des matériaux du barrage-poids de l’Oued Fodda
Page | 79
5. Failles dans le site du barrage :
Des failles tendant approximativement nord-sud, est-ouest et dans d’autres orientations ont
été cartographiées sur place, découpées dans des sondages et identifiées à partir des
photographies aériennes. La corrélation entre des données de surface et des sondages indique
que quelques failles peuvent avoir des mouvements relatifs de plus de 10m. Un tel mouvement
peut être une sommation de petits rejets sur des surfaces mineures dans une zone faillée plutôt
qu’un grand rejet sur une surface simple de faille.
Souvent des zones faillées sont partiellement ou complètement scellés par des veines de
calcite d’une épaisseur jusqu’à 3m, ou contiennent un remplissage d’argile rouge. Rien
n’apparaît comme étant repérée par une zone majeure de calcaire broyé faible.
6. Séismicité du site :
Les renseignements tectoniques montrent la possibilité de retour assez fréquent des
tremblements de terre au voisinage du site et les rives des tremblements de terre antérieure
montrent la fréquence de tels événements.
Constantine a souffert des tremblements de terre nuisibles dans le passé et n’importe quel point
dans un rayon de 150Km de celle-ci doit être considéré comme étant une origine possible de
tremblements de terre. Une aire d’origine potentielle totale de presque 200 000 Km²,
comprenant deux autres régions séismiques bien définies, dans laquelle la profondeur du foyer
apparaît comme étant superficielle (10 Km) doit être considérée.
Une accélération maximale de 0,3g ayant une période de retour de 100 ans a été déterminée
pour le site de Hammam Grouz, qui correspond à un épicentre de magnitude 5,5 à une distance
de 15 Km du site.
Ce modèle d’éléments finis est créé en utilisant le logiciel ABAQUS qui est l'un des principaux
programmes commerciaux d'éléments finis dans le monde et peut être appliquée à un grand
nombre d'applications dans l'ingénierie. Un modèle d'éléments finis avec plusieurs éléments
solides est utilisé pour modéliser le barrage. Un modèle d'élément finis avec plusieurs éléments
fluides est utilisé pour modéliser le réservoir d’eau.
Page | 80
Les résultats obtenus sont calculés par Abaqus, La modélisation du sol est faite en remplaçant
les encastrements à la base par des ressorts suivant la direction Z, avec caractéristique d’un Sol
rocheux K=60 000 KN /m3
Page | 81
Figure 57 : Modélisation du barrage poids sur ABAQUS
8. Accélérogramme utilisés :
Le calcul dynamique des barrages est synonyme de l’analyse du comportement sismique, car
c’est le seul effet dynamique qu’il peut être amené à supporter. On compte parmi les cas de
charge les forces d’origine sismique mais on les introduit sous forme de charge statique
équivalente obtenue par multiplication de la masse de l’ouvrage par le coefficient sismique
(accélération). La méthode de superposition modale a été utilisée pour la détermination des
déplacements nodaux à la crête du barrage.
Dans le but de mieux évaluer la réponse dynamique du barrage, le système couplé barrage-
réservoir a été soumis à différentes accélérations générées par divers événements sismiques.
Ces essais ont inclus un séisme au Ferndale City Hall le 10/10/38 (Chlef actuellement) avec
une magnitude de 6.60 (MSK), parmi d'autres. La figure ci-dessous illustre les 20 premières
secondes de l'accélérogramme de l'événement en question. Il est important de noter que ces
enregistrements variés ont été réalisés dans le but d'analyser l'influence spécifique de chaque
accélérogramme sur les déplacements observés au niveau de la crête du barrage
Page | 82
Accélérogramme Northwest Calif-02
4.5000E-01
4.0000E-01
3.5000E-01
3.0000E-01
Accélaration m/s2
2.5000E-01
2.0000E-01
1.5000E-01
1.0000E-01
5.0000E-02
0.0000E+00
0.010 0.035 0.048 0.075 0.120 0.180 0.280 0.380 0.500 0.800 1.300 2.000 3.200 4.400 7.00011.000
temps ( sec )
2.0000E-01
1.5000E-01
Accelaration
1.0000E-01
5.0000E-02
0.0000E+00
0.010 0.035 0.048 0.075 0.120 0.180 0.280 0.380 0.500 0.800 1.300 2.000 3.200 4.400 7.00011.000
temps
Page | 83
Accélérogramme Helena, Montana-02
8.0000E-01
7.0000E-01
6.0000E-01
5.0000E-01
Accelaration
4.0000E-01
3.0000E-01
2.0000E-01
1.0000E-01
0.0000E+00
0.010 0.035 0.048 0.075 0.120 0.180 0.280 0.380 0.500 0.800 1.300 2.000 3.200 4.400 7.000 11.000
temps
1.0000E+00
8.0000E-01
Accelaration
6.0000E-01
4.0000E-01
2.0000E-01
0.0000E+00
0.010 0.035 0.048 0.075 0.120 0.180 0.280 0.380 0.500 0.800 1.300 2.000 3.200 4.400 7.000 11.000
temps
9. Résultats et interprétation :
a. Résultats et interprétation sur les contraintes :
A partir des études réalisées dans ce chapitre on peut tirer comme conclusion les points
suivants :
Page | 84
Les résultats indiquent le début de fissures et leur propagation du talon vers l'orteil,
entraînant une réduction de la rigidité de la structure. Ces fissures de traction peuvent
provoquer le glissement et le basculement du barrage avec un déplacement permanent.
Au point de changement de pente ; le point probable de concentration des contraintes, la
contrainte de traction développée sont très importantes et risque de dépasser les
contraintes admissibles
Page | 85
Figure 63 : Concentration des contraintes dans la face amont du barrage
Page | 86
Figure 64 : Concentration des contraintes dans la crête et le pieds du barrage
Le barrage est très rigide cela est dû à la variation croissante (base vers crêtes) de la
fréquence.
Page | 87
Figure 65 : déplacement obtenu par le séisme Northwest Calif-02
Page | 88
Figure 67 : déplacement obtenu par le séisme Helena Montana 01
Interprétation et résultats :
Séisme 1 :
Page | 89
Séisme 2 :
Séisme 3 :
Page | 90
Conclusion :
Page | 91