REPUBLIQUE DE CÔTE-D’IVOIRE

Union – Discipline – Travail

MINISTERE D’ETAT, MINISTERE DE L’AGRICULTURE ET DU DEVELOPPEMENT RURAL

Agence pour le Développement de la filière riz (ADERIZ)

PROJET D'AMENAGEMENT HYDRO-AGRICOLE DANS LES REGIONS DE

FOLON ET KABADOUGOU (PAHA– FK)

RAPPORT D’ETUDE D’EXECUTION

BARRAGE DE BIREKO

Société des Grands Travaux Internationaux

SGTI

Version 04

© Janvier 2023
1. Table des matières
CHAPITRE I : ETUDE DE LA RETENUE D’EAU ---------------------------------------------------------------------------- 6
1.1 Topographie du terrain----------------------------------------------------------------------------------------- 6
1.2 Détermination des courbes hauteur-Surface-volume de la cuvette -------------------------------- 6
1.3 Apports d’eau par ruissellement ----------------------------------------------------------------------------- 7
1.4 Détermination de la superficie de bassin versant (S) --------------------------------------------------- 8
1.5 Détermination des paramètres climatiques -------------------------------------------------------------- 8
1.6 Détermination des pertes par dépôt solides-------------------------------------------------------------- 9
1.7 Besoins d’eau de projet ---------------------------------------------------------------------------------------- 9
1.7.1 Besoins agricoles ---------------------------------------------------------------------------------------- 10
1.7.2 Besoins pastoraux --------------------------------------------------------------------------------------- 11
1.8 Pertes par évaporation et infiltration dans la retenue d’eau --------------------------------------- 11
1.8.1 Pertes par infiltration ---------------------------------------------------------------------------------- 11
1.8.2 Pertes par évaporation -------------------------------------------------------------------------------- 11
1.9 Simulation de la retenue d’eau ---------------------------------------------------------------------------- 11
CHAPITRE II : DIMENSIONNEMENT DE LA DIGUE ------------------------------------------------------------------- 13
2.1 Détermination de la cote des Plus Hautes Eaux (P.H.E) ---------------------------------------------- 13
2.2 Détermination du débit crue de projet (Qp) ------------------------------------------------------------ 13
2.2.1 Détermination de la crue décennale (Q10) ------------------------------------------------------- 13
2.2.2 Détermination de coefficient de majoration (C) ------------------------------------------------ 17
2.2.3 Méthode de X0 ou CIEH/EIER ------------------------------------------------------------------------ 18
2.2.4 Laminage de la crue par la méthode de continuité --------------------------------------------- 19
2.3 Détermination de la la cote PHE : ------------------------------------------------------------------------- 22
2.4 Détermination de la cote crête----------------------------------------------------------------------------- 22
2.4.1 Calcul de revanche -------------------------------------------------------------------------------------- 22
2.4.1.1 Calcul de la hauteur des vagues (Hv) -------------------------------------------------------- 23
2.5 Détermination de la cote crête----------------------------------------------------------------------------- 24
2.6 Hauteur de digue (H) ----------------------------------------------------------------------------------------- 24
2.7 II.4) hauteur de la retenue normale (h) ------------------------------------------------------------------ 25
2.8 Largeur de crête (C) ------------------------------------------------------------------------------------------- 25
2.9 Pentes de talus amont et aval ------------------------------------------------------------------------------ 25
2.10 Emprise de la digue (L) --------------------------------------------------------------------------------------- 25
CHAPITRE III : DIMENSIONNEMENT DES COMPOSANTES DE LA DIGUE --------------------------------------- 26

P a g e 2 | 44
 3.1 Tranchée d’ancrage ------------------------------------------------------------------------------------------- 27
3.2 Etude des infiltrations et du dispositif de drainage --------------------------------------------------- 28
3.2.1 Tracé de ligne de saturation -------------------------------------------------------------------------- 29
3.3 Evaluation de débit de fuite (q) ---------------------------------------------------------------------------- 31
3.3.1 Evaluation de débit de fuite à travers le barrage (q1) ------------------------------------------ 31
3.3.2 Evaluation de débit de fuite à travers le sol de fondation (q2) ------------------------------- 31
3.4 Calcul de l’épaisseur e du drain-filtre (e) ---------------------------------------------------------------- 32
3.5 Protection talus aval ------------------------------------------------------------------------------------------ 33
3.6 Protection talus amont --------------------------------------------------------------------------------------- 33
3.7 Etude de stabilité de talus amont ------------------------------------------------------------------------- 34
3.8 Butée de pied --------------------------------------------------------------------------------------------------- 35
3.9 Fossé de pied aval --------------------------------------------------------------------------------------------- 36
CHAPITRE IV : ETUDE DE LA CONDUITE DE VIDANGE -------------------------------------------------------------- 36
4.1 Détermination de temps de vidange --------------------------------------------------------------------- 36
4.2 Dimensionnement de bassin de dissipation type impact -------------------------------------------- 39
CHAPITRE V : DIMENSIONNEMENT DE L’EVACUATEUR DE CRUES --------------------------------------------- 41
5.1 Dimensionnement du chenal ------------------------------------------------------------------------------- 41
5.2 Dimensionnement du coursier ----------------------------------------------------------------------------- 41
5.3 Dimensionnement du convergent ------------------------------------------------------------------------- 42
5.4 Dimensionnement du bassin à ressaut ------------------------------------------------------------------- 42
5.5 Dimensionnement du chenal débouchant dans le drain principal -------------------------------- 43
ANNEXES ---------------------------------------------------------------------------------------------------------------------- 44

P a g e 3 | 44
 FICHE SYNOPTIQUE DU BARRAGE BIREKO

A- LA RETENUE RESULTATS DE L’APD RESULTATS DE L’ETUDE D’EXECUTION

1. Cote du fond du lit mineur : 377,67 m 405,45 m

2. Niveau de retenue normale : 384,50 m 413,55 m
3. Niveau des plus hautes eaux : 384,95 m 414,15 m
4. Volume de retenue normale : 5,19 Mm3 5,18 Mm3
5. Aire en retenue normale : 162,53 Ha 166 Ha

B- LE BARRAGE

2. Couronnement : 385,85 m 415,40 m

3. Longueur en crête : 709 m 668 m
4. Largeur en crête : 5 m 5 m
5. Talus : Aval: 2 H/1V; Amont : 2 H/1V Aval: 2,5 H/1V; Amont: 2 H/1V
6. Hauteur max. sur fondation : 11,89 m 14,09 m
7. Hauteur max. sur TN : 8,2 m 10 m
8. Protection amont de la digue Rip-rap Rip-rap
9. Protection aval de la digue Enherbement Enherbement
10. Tapis-filtrant 20 cm sable + 40 cm gravier + 20 cm sable + 40 cm gravier +
20 cm sable 20 cm sable

C- L’EVACUATEUR DE CRUE+BASSIN DE DISSIPATION

1. Crue de projet (Q100) 13,95 m3/s, débit laminé 3,19 m3/s 54 m3/s, débit laminé 31 m3/s
2. Type : Déversoir latéral à seuil libre Déversoir latéral à seuil libre
3. Longueur de déversement : 4 m 20 m
4. Charge maximale : : 0,4 m 0,6 m
5. Cote seuil : 384,5 m 413,50 m
6. Implantation : Rive droite Rive droite
P a g e 4 | 44
7. Coursier : Béton armé, largeur 2 m, longueur Béton armé, largeur 4 m, longueur 275 m
261,5 m.
8. Restitution : Bassin de dissipation en béton armé Bassin de dissipation en béton armé de type II
de type 2

D- VIDANGE DE FOND

1. Implantation : Rive droite sous la digue Rive droite sous la digue

2. Type : Conduite en acier enrobé de béton. Conduite en acier enrobé de béton.
3. Longueur conduites : 44 m 60 m
4. Côte de calage à l’entrée : 377,45 m 408,15 m
5. Restitution : Bassin de type impact puis chenal non Bassin de type impact puis chenal non revêtu
revêtu
6. Pente : 1% 1%
7. Dimension conduite Ø800 mm Ø800 mm
8. Equipement hydromécanique : 1 vanne de garde, Ø800, 1 vanne de garde, Ø800,
1 vanne de vidange type papillon, 1 vanne de vidange type papillon, Ø800
Ø600
9. Prises d’eau : 1 piquage Ø600 mm contrôlé par 1 1 piquage Ø600mm contrôlé 1 vanne papillon
vanne papillon Ø600 Ø600

P a g e 5 | 44
CHAPITRE I : ETUDE DE LA RETENUE D’EAU

Ce chapitre a pour but de vérifier la faisabilité de réalisation du barrage en fonction des

apports et de relief.

1.1 Topographie du terrain

Les travaux topographiques sur le site de BIREKO ont permis de donner l’état de lieux du TN,
de l’emprise de la digue et de la retenue de barrage.

Cette étude permettra de déterminer le volume et la superficie de la cuvette, elle permettra

aussi de dimensionner la digue et des ouvrages connexes.

1.2 Détermination des courbes hauteur-Surface-volume de la cuvette

. A partir du plan topographique de la cuvette, nous déterminons les volumes stockés et les
surfaces correspondant à chaque hauteur.

La méthode de calcul est comme suit :

- Réalisation des courbes des niveaux sur le plan d’état des lieux de la cuvette pour
avoir la surface à chaque hauteur ;

S i  S i 1
- Calcule de la surface moyenne Sm  ;
2
- Calcule le volume partiel Vp  Sm  h ;

- Calcul de volume cumulé VC   V pi ;

i 1

Le tableau ci-dessous présente les résultats :

Cote (m) Volume (m3) Volumes partiels (m3) Surface (m²)

406,5 0 0 0
407 3 604 3 604 14 417
408 108 273 104 669 194 920
409 405 421 297 149 399 377
410 965 494 560 073 720 768
411 1 804 356 838 862 956 956
412 2 901 978 1 097 622 1 238 288
413 4 276 210 1 374 232 1 510 176
413,50 5 181 785 1 523 979 1 660 132
414 5 922 710 1 646 500 1 782 824
415 7 899 394 1 976 684 2 170 543
La représentation de la courbe Hauteur-Volume se présente ci-dessous :
P a g e 6 | 44
 Courbes hauteur-volume/hauteur surface de BIREKO
416
415
414
413
Cotes (m)

412
411
410
409
408
407
406
3 2 1 1 2 3 4 5 6 7 8
Surface(x10Ha) Volume(x10^6m3)
Surface Volume

1.3 Apports d’eau par ruissellement

En l’absence de mesures directes de débits des bassins versants, les écoulements annuels
sont estimés par la méthode de COUTAGNE. Celle-ci est basée sur la superficie du bassin
versant et sur l’analyse des moyennes de paramètres climatiques dont notamment la pluie
et la température moyennes annuelles

V  10 3  Pan  K e  S

Avec :
V : volume des apports en m3 et Ke : coefficient d’écoulement ;
Pan : pluie annuelle en mm et Sb : superficie du bassin versant en km²
Le coefficient d’écoulement s’obtient à partir du déficit d’écoulement donné par la formule

tropicalisée de TURC est: avec

D : déficit d’écoulement en mm ; P : la pluie moyenne annuelle (mm) T : température moyenne
annuelle en °C

D’où
Les apports moyens annuels calculés par cette formule, à l’exutoire constitué par l’axe du
barrage, sont évalués en années moyenne

P a g e 7 | 44
1.4 Détermination de la superficie de bassin versant (S)
La superficie du bassin versant a été déterminée à partir du logiciel Global Mapper
présentées ci-dessous.

Le bassin versant du site de BIREKO couvre une superficie de 63,74 Km² soit 64 Km2

1.5 Détermination des paramètres climatiques

La pluie moyenne annuelle a été déterminée à partir des données de pluies de 1966 à 2021
à la station météorologique de Minignan et la température moyenne annuelle a été évaluée
à partir des données de 1995 à 2002 à la station météorologique d’Odiennée.

Janv Févr Mars Avr Mai Juin Juil Août Sept Oct Nov Déc Total
Pluie moy 2,92 8,2 32 72 111 165 296 353 245 110 27 1,5 1423
(Source : SODEXAM)

Janv Févr Mars Avr Mai Juin Juil Août Sept Oct Nov Déc Total
Température moy 24,9 27 29 29 28 26 25 25 25 25 25 24 26,1
(Source : SODEXAM)

Le calcul des apports permet d’obtenir les résultats suivants :

Valeurs de l'APD Valeurs de l’EE
Apport moyen annuel (m3) 12 781 288 11 648 266

P a g e 8 | 44
1.6 Détermination des pertes par dépôt solides
L’eau de ruissellement entraine avec elle des matériaux solides qu’elle arrache le long de son
parcours, ces matériaux se déposent dès que la vitesse de l’eau devient inferieur à un certain
seuil, ainsi une retenue qui amène l’eau à réduire sa vitesse à une valeur pratiquement nulle
provoque le dépôt de la presque totalité des matériaux solides transportés par les eaux du
ruissellement. Ainsi en Afrique, où les reliefs ne sont pas très accusés mais les précipitations
concentrées, les dépôts solides peuvent parfois avoir des valeurs relativement importantes.

Le volume moyen annuel des apports solides (V) est : V  D S

V : Le volume annuel des apports solides exprimée en (m³/an).

D : dégradation spécifique annuelle exprimée en (m³/km²/an).

S : Superficie du bassin versant exprimée par (km²) ; S = 64km²

 Détermination de la dégradation spécifique annuelle (D)

Nous calculons la dégradation spécifique annuelle (D) par la formule de GRESILLON

2 , 2
P 
D  700   an   S 0,1
 500 

Avec ; D : la dégradation spécifique annuelle, en (m³/km²/an). D= 46,26 m³/km²/an

Pan : Pluie moyenne annuelle exprimée par (mm) ; Pan = 1423 mm

S : Superficie du bassin versant exprimée par (km²) ; S = 64 km²

GRESILLON(APD) GRESILLON (EE)

156 204 152 845

NB : Les pertes par dépôts solides peuvent être réduite par l’auto curage en ouvrant la vanne
de vidange périodiquement.

1.7 Besoins d’eau de projet

Les besoins en eau évalués par l’APD sont de deux types à savoir :
 les besoins agricoles
 les besoins pastoraux

P a g e 9 | 44
1.7.1 Besoins agricoles
Sur la base des données pluviométriques de 1966-2021 de la station de Minignan et des
données de l’ETP de 1994-2002 de la station d’Odienné, le déficit pluviométrique a pu être
évalué et présenté par la courbe suivante.

Déficit pluviométrique à Minignan

Pluie ETP 1/2 ETP
400.0

350.0

300.0
Hauteur de pluie

250.0

200.0

150.0

100.0

50.0

0.0
Janv Févr Mars Avr Mai Juin Juil Août Sept Oct Nov Déc

Cette courbe montre que la pluviométrie à Bireko est monomodale où le pic se situe dans le
mois d’Août avec 353mm de pluie en moyenne.
Le calage des cycles culturaux a été fait suivant celui de l’APD.
Calage 1er cycle Calage 2ième cycle
Spéculations APD EE APD EE
Riz Avril-Juillet Juin-Octobre Août-Nov Déc-Avril
Maraîcher Mai-Avril Oct-Janv Oct-Janv -

Analyse : Le calage des cycles rizicoles a été fait de sorte à profiter au mieux des pluies et
d’éviter les récoltes pendant la saison pluvieuse.
Selon les informations reçues des paysans sur le terrain, le cycle maraicher débute au mois
d’Octobre. Le calage du maraicher a donc été fait de sorte à éviter les périodes de grandes
pluies. De ce fait, un seul cycle a été retenu.
Les besoins en eau ainsi obtenus sont les suivants :
(Calcul des besoins en eau : Annexe 1)

Calage 1er cycle Calage 2ième cycle

Spéculations APD EE APD EE
Riz 9 938 4 035 8 321 19 137
Maraîcher 995 5 876 5 081 -
Total 10 933 9 911 13 402 19 137

P a g e 10 | 44
1.7.2 Besoins pastoraux
Les besoins pastoraux de l’APD utilisés dans le tableau de simulation ont tout simplement
été reconduits dans l’EE.
Sites
Bireko
Besoins
Besoins pastoraux annuels (m3) 252 665

1.8 Pertes par évaporation et infiltration dans la retenue d’eau

En plus des dépôts solides, il y a des pertes par évaporation et infiltration dans la retenue

1.8.1 Pertes par infiltration

L’infiltration moyenne de 1,5 mm/jour soit environ 15 m3/ha/j durant les premières années.
Par la suite avec les dépôts solides, le colmatage du fond de la cuvette réduira de façon
significative ces pertes. Dans le cadre de notre étude nous considèrerons que les infiltrations
sont de 1mm/jour

1.8.2 Pertes par évaporation

- L’évaporation de la période considérée est la différence entre les valeurs de
l'évapotranspiration potentielle (l’ETP) et la pluviométrie mensuelle.

- Les données d’évapotranspiration potentielle (ETP) fournis dans le tableau ci-dessous sont
celles de la station agro météorologique d’Odienné (1994-2002)

Mois Jan Fév Mars Avr Mai Juin Juil Août Sept Oct Nov Déc Total

Pan (mm) 3 8 32 72 111 165 296 353 245 110 27 1 1423

ETP (mm) 119 137 146 149 144 134 122 121 124 128 119 108 1551
Evaporation
116 129 115 77 33 0 0 0 0 18 91 106 686
(mm)
Volume
évapo/infiltr 138371 148312 136850 100573 60235 28200 29140 29140 28200 46409 114073 128991 988493
(m3)
NB : La superficie du plan d’eau pour le calcul des volumes d’évapo/Infiltr est la moyenne du plan
d’eau dans la cuvette. Splan d’eau=94 ha

1.9 Simulation de la retenue d’eau

La méthodologie de conception adoptée consiste à faire varier la hauteur de la digue de sorte
à pouvoir alimenter toute la superficie exploitable du point de vue topographique et
pédologique tout en prenant en compte les prélèvements (besoins agro-pastoraux) et les
pertes (évaporation, infiltration, dépôts solides). La spéculation la plus exigeante du point de
vue hydrique sera privilégiée.

P a g e 11 | 44
Principe de la simulation

L’hypothèse de base est de considérer la cuvette pleine après avoir fixé le niveau de retenue
normale. On établit mensuellement le bilan des apports et des pertes en eau du barrage. Ce
bilan permet en utilisant la courbe hauteur-Volume-Surface, de connaître le volume
disponible dans la cuvette, qui doit varier entre deux valeurs extrêmes :

● Une valeur minimale égale à la tranche d’eau morte au-dessous de laquelle aucun
prélèvement n’est possible et la réalisation du cycle impossible.
● Une valeur maximale au-dessus de laquelle, il y’à écrêtage des apports parce que le
niveau de l’eau à atteint la cote de retenue normale.
Sur une longue série de données agro climatiques consécutives de 1972-2021 (50ans), l’on
simule la réalisation des cycles culturaux. Aussi pour un mois déficitaire, nous considérons le
cycle perdu. La valeur cible adoptée comme taux de réussite satisfaisant est 90%.
Tableaux des simulations : Annexe 2

Désignation Simul.1 Simul.2

Taux de réussite (%) 90 90

Superficie riz (Ha) 241 119,50

Superficie maraichage (Ha) 25 25

Superficie totale exécution (Ha) 266 144,50

Superficie selon APD (Ha) 475 -

Côte PEN (m) 413,55 412,00

Volume de la retenue (m3) 5 181 785 2 901 978

Volume de la retenue selon APD (m3) 5 191 964 -

Côte prise d’eau (m) 408,15 408,15

Hauteur du barrage/Cote PEN (m) 7 5,5

Analyse : La simulation montre que la retenue d’eau domine largement la superficie prévue par
le marché qui est de 115 ha. La capacité de la retenue d’eau optimale pour les 115 ha est d’environ
3 millions de m3 pour une différence de hauteur de digue (à la cote PEN) par rapport à l’APD qui
est de 1,5m. Cependant, durant les saisons sèches, l’eau se fait rare dans la zone, la population
pourrait s’en servir. Aussi, la retenue peut être beaucoup sollicitée pendant cette période à cause
de la non maitrise du nombre exact des bovins transhumants dans la zone.
Il serait donc préférable de maintenir la capacité de la retenue prévue par l’APD pour pallier
aux pertes imprévues.

P a g e 12 | 44
CHAPITRE II : DIMENSIONNEMENT DE LA DIGUE
L’objectif est d’assurer la fonctionnalité et la stabilité de la digue et ses ouvrages connexes.
Dans ce chapitre nous déterminons :

- La cote de Plus Haute Eau (P.H.E).

- La cote crête de la digue.
- La largeur de la crête.
- La hauteur de la digue.
- La largeur de la digue.
- La pente de talus amont et aval.
2.1 Détermination de la cote des Plus Hautes Eaux (P.H.E)
Cette cote est déterminée par la formule suivante :

Cote PHE= Cote RN+ hlam

Avec ;
hlam : tirant d’eau laminée ou la charge maximum sur le seuil au déversoir exprimer en (m)

Plusieurs méthodes permettent de déterminer ce tirant d’eau telles que la méthode X0

méthode X0 de CIEH/ EIER (Comité Inter-Etats d’Etude Hydraulique / Ecole Inter-Etats de
l’Equipement Rural), la méthode de continuité, …

2.2 Détermination du débit crue de projet (Qp)

Le débit de crue centennale ou de projet est obtenu en multipliant le débit de crue
décennale par un coefficient majorateur, qui est fonction des précipitations de même temps
de retour pour le temps de base caractéristique du bassin versant :

Qcp  C  Q10

Avec ; Qcp : la crue de projet est égale à la crue centennal Q100 exprimer par (m³/s).
Q10 : la crue décennale exprimer par (m³/s).
C : coefficient majorateur supérieur à 1 lié à la sécurité de l’ouvrage

 Pour déterminer la crue de projet (Qp) il faut trouver (Q10) et (C).

2.2.1 Détermination de la crue décennale (Q10)

En l’absence de données d’observations des crues sur le site, on exploitera la méthode

d’ORSTOM (Office de Recherche en Science et Technique Outre-Mer) ou de RODIER ET
AUVERAY qui est applicable dans les conditions suivantes :

P a g e 13 | 44
  Pour de bassins versants de superficie inférieur à 120 km².
 S’applique en théorie à toute l’Afrique de l’Ouest.
 La hauteur de précipitation annuelle entre 150 et 1600 mm.

Donc la crue décennale sera déterminée par la formule suivant :

  K  Kr10  P10  Sbv

Q10 
Tb
Q10 : Débit de crue décennale en (m³/s)
 : Coefficient de pointe = 2,5
K : Coefficient d’abattement
Kr10 : Coefficient de ruissellement décennal en (%)
P10 : Pluie maximale journalière décennale en (mm)
S : Superficie du bassin versant = 69 km²
Tb : Temps de base en (s)
 Pour déterminer la crue décennale (Q10) on besoin de trouver les paramètres précédents.
 Le coefficient d’abattement (K)

C’est un paramètre défini par l’équation simplifiée de Vuillaume (1974)

161  0,042  Pan

K 1  log S ; Avec ; Pan : pluie annuelle en (mm) ;
1000
S : superficie du bassin versant en (km²)
 Le coefficient de ruissellement (Kr10)

On détermine le coefficient de ruissellement (Kr10) à partir abaque (voir Annexe 3) qui est
en fonction de la superficie de bassin (S) versant et la classification de la zone (Ri Pi).

 La classe de perméabilité (Pi) et la classe de pente (Ri)

Ces sont des paramètres déterminés suivant la classification de la zone en fonction de la

perméabilité et du relief. À partir de livre des J. RODIER et C. AUVRAY « estimation des debits
de crue decennales pour les bassins versants de superficie inferieur a 200 km² en afrique
occidentale »
Ainsi on définit la classe de perméabilité (P) suivante :

P3 : Bassins assez imperméables comportant des zones perméables d’entendue notable ou

bassins homogènes assez peu perméables.

P a g e 14 | 44
Pour connaître la classe de pente (R) il faut d’abord déterminer les valeurs la pente
longitudinal (il) et la pente transversal (it)

- La pente longitudinale (il) est obtenue en éliminant 20% de part et d’autre de la

longueur du cours d’eau principal, les 60% restant sont rapportés à la différence de
niveau du cours d’eau. (m/km)

H max cd  H min cd
il 
L60%cd
Les paramètres Hmax cd , Hmin cd , L60% cd sont déterminés à partir de logiciel Global
Mapper.
Hmax cd = 411 m ; altitude maximale sur le cours d’eau principal à L20% cd

Hmin cd = 384 m ; altitude minimale sur le cours d’eau principal à L20% cd

Lcd : la longueur totale du cours d’eau principal ; Lcd 11 129 Km

L20% cd : la longueur du cours d’eau principal à 20% ; L20% cd = 2 226Km

L60% cd : la longueur du cours d’eau principal à 60% ; L60% cd = 6 677Km

 il = 0,4%

R2 à des pentes faibles inférieures à 0, 5 : ce sont des terrains de plaines.

 Pluie maximales journalières fréquentielles

L’étude statistique des pluies journalières extrêmes a pour objet de prévenir les crues. Elle est
donc d'un grand intérêt pour anticiper les catastrophes, telles que des inondations, et ainsi
permettre un dimensionnement réfléchi des ouvrages. Dans le cas le plus courant, il n’y a
qu’une averse par 24 heures, de sorte que l’étude statistique des tornades se ramène à l’étude
des averses journalières. Les pluies maximales journalières retenues pour cette étude
couvrent la période de 1972 à 2021 soit une série de cinquante (50) valeurs et sont ajustées par
Gumbel en période humide et en application de la méthode des moments. Les résultats des
données de l’ajustement des pluies maximales journalières est présentée dans le tableau ci-
dessous.

P a g e 15 | 44
Valeurs des quantiles de GUMBEL en période humide à Minignan

Pluie Période de retour Probabilité au non Valeur des quantiles

Ecart type
annuelle (ans) Dépassement f(x) de gauss (mm)
5 0,80 5,48 108
PERIODE 10 0,90 7,40 123
HUMIDE 50 0,98 11,9 155
100 0,99 13,9 169

 Le temps de base (Tb)

On détermine le temps de base (Tb) à partir abaque (voir Annexe 4) qui est en fonction de la
superficie de bassin (S) versant, la classe de pente (Ri), la pluie annuelle Pan et le type de
climat.

 Détermination de temps de montées (Tm)

On détermine le temps de montées (Tm) à partir abaque (voir Annexe 5) qui est en fonction
de la superficie de bassin (S) versant et la classe de pente (Ri), la pluie annuelle Pan et le
type de climat.

P a g e 16 | 44
À partir des paramètres, on peut calculer maintenant la crue décennale (Q10).

VALEURS
CARACTERISTIQUES SYMBOLE UNITE
APD EE
Superficie bv (carte 1/50000 utilisé) Sbv km2 69 64
Coefficient de pointe α
Coefficient d’abattement (vuillaume K 0,81 0,81
Classe de pente Ri R2 R2
Classe de perméabilité Pi P4 P3
Coefficient de ruissellement Kr % 10 31,17
Pluie maximale journalière décennale Pjmax10 mm 130 123
Pluie maximale journalière centennale Pjmax100 mm - 169
Temps de montée Tm min 540 582
Temps de base Tb min 3100 3087
Débit de crue décennale Q10 m3/s 9,82 27
Débit de crue centennale ou débit de projet Qp m3/s 13,95 54

2.2.2 Détermination de coefficient de majoration (C)

Le coefficient de majoration utilisée par l’APD est de 1,42

L'expression donnant le coefficient (C) est :

0 ,12
 Tb 
 
P100  P10  24 
C 1 
P10 K r10
Où
P10 : Précipitation journalière correspondant à une période de retour de 10 ans ;

P100 : Précipitation journalière correspondant à une période de retour de 100 ans

Tb : Temps de base en secondes ; Tb = 51,44 heures= 185 201s

Kr10 : Coefficient de ruissellement de la crue décennale exprimé en fraction et non pas en

pourcentage ; Kr10= 31,17% = 0,32

CALCUL DU COEFFICIENT DE MAJORATION

Q10(mm) Q100(mm) Kr10 Tb10(h) (P100-P10)/P10 ((Tb/24)^0,12)/Kr10 Cp
123 169 0,3117 51,44 0,37398374 3,515555544 2,31476061
Le coefficient retenu est celui proposé par l'ASEER (Association Européenne d'Etudes
d'Equipement Ruraux) qui est de 2. Ce coefficient surestime les crues centenaires
surtout en régime tropical
Le débit de crue centennal ou débit de projet= 2xQ10 = 54m3/s

P a g e 17 | 44
2.2.3 Méthode de X0 ou CIEH/EIER
La méthode de X0 nous permet de déterminer le tirant d’eau laminée (hlam), le débit évacuée
(Qev) la longueur de déversoir (L).

Validité de la méthode
La validité des calculs est limitée dans le cas de très forts laminages. Lorsque β=
Qemax/Qcmax trouvé est inférieur à la valeur Q/Qcmax donnée par l'abaque ci-dessous, il y
a lieu d'être prudent quant à la valeur de Qemax. C'est un cas de figure peu fréquent.

On applique cette méthode par itération sur un tableur Excel :

A partir ce tableau ci-dessus on peut déterminer les résultats suivants :

P a g e 18 | 44
Donc le laminage est valable pour hlam=0,45m, Longueur=72m et Qev=38m3/s
La longueur déversement n’est pas optimale. Dans ces conditions le laminage de la crue sera
évalué par la méthode de continuité comme l’APD.

2.2.4 Laminage de la crue par la méthode de continuité

Pente du radier

Pour un débit évacué donné, la charge dans le barrage sera minimale si l’écoulement sur le
radier est torrentiel. Il faut donner au radier une pente suffisante pour qu’il soit torrentiel.
Ce débit, pour lequel la pente du radier devient critique, est calculé à partir des expressions
suivantes :
Qc2 S 3
  (1) équation du régime critique
g l

 Qc  Ks S I R2 / 3 (2) formule de Manning-Strickler, avec :

 Qc : débit critique
 g : 9,81 m/s-2

P a g e 19 | 44
  S : section mouillée
 l : largeur au miroir
 Ks : coefficient de Manning-Strickler = 70 pour canal en béton
 I : pente longitudinale
 R : rayon hydraulique
En combinant ces 2 équations, on obtient l’expression du tirant d’eau sur le radier :
g 3 (b  2 y 1  m2 ) 4
y( )  (3) avec :
K s2 I (b  my)(b  2my)3
 y : tirant d’eau sur le radier
 m : fruit des rampes d’accès du radier
 b : largeur en base du radier
Cette équation permet de calculer par approximations successives le tirant d’eau sur le radier,
et par suite le débit critique à l’aide de l’expression (1).
Loi de débit de l’évacuateur

En appliquant la conservation de la charge à l’entrée du radier, on a :

Q2
H  y  (1  k ) (4) avec :
2 gS 2
 H : charge dans le barrage
 y : tirant d’eau
 S : section mouillée du radier au passage de la crue
 k : coefficient tenant compte des pertes de charge à l’entrée du radier (k=0,26
pour un rétrécissement brusque).
1
D’où Q  S 2 g ( H  y )  S 2 g ( H  y ) (5) avec   0.89 .
1 k

L’expression (4) permet de calculer la charge correspondant au débit critique du radier. Soit
Hc cette charge. Lorsque la charge dans le barrage varie, on a les 2 cas suivants :
1er cas H  Hc
Dans ce cas, le débit évacué est supérieur ou égal au débit critique et l’écoulement sur le radier
est critique ou torrentiel. La profondeur critique est donc atteinte dès l’entrée du radier et on
Q2 S 3
a:  (6)
g l
En combinant les équations (5) et (6), on obtient l’expression suivante qui donne le tirant
 2 (b  my)
d’eau en fonction de la charge : y  2 H (7)
(2 2  1)b  (4 2  1)my
2è cas H < Hc
P a g e 20 | 44
L’écoulement sur le radier est fluvial et le débit évacué vérifie la formule de Manning-strickler
Q  Ks S I R2 / 3 (8)
En combinant les équations (5) et (8), on obtient l’expression suivante qui donne le tirant
4/3
K 2 I  y(b  my) 
d’eau en fonction de la charge : y  H  s 2   (9)
2 g  b  2 y 1  m2 

Les équations (7) ou (9) permettent de calculer par approximations successives le tirant d’eau
dans le radier correspondant à une charge donnée dans le barrage.
En utilisant cette valeur dans l’expression (5), on obtient le débit évacué, ce qui permet de
définir la loi de débit de l’évacuateur.
Laminage de la crue du projet

On détermine le débit évacué à chaque instant en faisant le bilan en volume des différents
apports pendant l’intervalle de temps donné :
Vs  (Qcm  Qem )t (10) et Vs  Am H (11) avec :

 Vs : variation de volume pendant l’intervalle de temps t

 Qcm : débit moyen de crue pendant l’intervalle de temps
 Qem : débit moyen évacué pendant l’intervalle de temps
 Am : superficie moyenne du plan d’eau pendant l’intervalle de temps
 H : variation de la charge dans le barrage pendant l’intervalle de temps.
Le débit moyen évacué est calculé par approximations successives comme suit :
 on fixe à priori la charge dans le barrage à l’instant considéré
 on détermine à l’aide de la courbe hauteur - surface du barrage la superficie du plan
d’eau à l’instant considéré, puis la superficie moyenne du plan d’eau pendant
l’intervalle de temps
 on calcule à partir des équations (7) ou (9) le tirant d’eau du radier, puis le débit évacué
à l’instant considéré à l’aide de l’équation (5) et le débit moyen évacué pendant
l’intervalle de temps
 on calcule la variation de volume à l’aide des équations (10) et (11)
 on corrige éventuellement la charge estimée et on reprend le processus jusqu’à la
convergence des variations de volume calculées précédemment.
Ce processus se programme sur tableur. Au passage de la crue du projet, nous obtenons le
débit maximum évacué et la charge maximale sur le seuil.

P a g e 21 | 44
Données d’entrée
Données sur le barrage
Crue (m3/s) 54,0
Temps de base (min) 3 087
Temps de montée (min) 582
Cote PEN (m) 413,55
Caractéristiques du radier
base(m) 20,00
fruit 10,00
Ks 70
Coefficient débit 0,82

Tableau de calcul du laminage : Annexe 6

Conclusion :
La hauteur d’eau déversant est de hlam=0,60m, Longueur=20m et Qev=31m3/s

2.3 Détermination de la la cote PHE :

Cote PHE= Cote RN+ hlam= 414,15m

2.4 Détermination de la cote crête

La cote crête est définie par :
Cote crête = Cote PHE + Revanche (R)

2.4.1 Calcul de revanche

La revanche R est la différence de cote entre le PHE et la crête.
La connaissance du fètch et de la vitesse des vents dominants dans la zone permet
l’utilisation du calcul de la revanche avec la formule de Gaillard

V2
R  0,75H V  et V  1,5  2  H V
2g

P a g e 22 | 44
Avec ; R : Revanche selon la formule de Davis (m)
V : Vitesse de propagation des vagues (m/s) ;
Hv : Hauteur des vagues (m)
g : Accélération de la pesanteur (m/s²) ; g = 9,81 m/s²

 Pour déterminer la revanche (R) on va suivre les étapes suivantes.

1) On détermine la hauteur des vagues (Hv).
2) On détermine la vitesse de propagation des vagues (V).
2.4.1.1 Calcul de la hauteur des vagues (Hv)
La hauteur des vagues dépend de plusieurs facteurs et principalement de la force du vent, de
sa durée et de l’étendue sur laquelle il souffle avec une direction à peu près constante (fetch).
L’action des vagues dépend essentiellement de la dimension de la retenue et de sa localisation
géographique
Pour le calcul de la hauteur des vagues, on utilise les formules suivantes :
 Formule de Stevenson (Thomas, 1976) :
Cette méthode dépend de la valeur du fetch.

 Pour f  30 km : HV  0,76  (0,34  f )  (0,26  f 1/ 4 )

 Pour f  30 km : HV  0,34  f
Avec,
Hv : Hauteur des vagues en (m) ;
f : Longueur du fetch en (km) ; f = 2,58 km (lecture directe sur la cuvette à l’aide de logiciel
covadis)
 Formule de Molitor (1935) :
Dans cette formule empirique, en fait intervenir la vitesse du vent " w " exprimée en km/h et
la longueur du fetch « f » en Km :

 Pour f  30 km : HV  0,76  (0,032  w  f )  (0,26  f 1/ 4 )

 Pour f  30 km : H V  0,34  ( w  f )1 / 4

P a g e 23 | 44
Avec,
Hv : Hauteur des vagues en (m) ;
f : Longueur du fetch en (km) ; f= 2,58 km

w : Vitesse du vent en (km/h) ; déterminé en ligne sur «WeatherSpark »

(https://fr.weatherspark.com) qui donne l’historique de vent mensuelle de 2015 jusqu’à 2021
pour Odienné ce représente dans le tableau suivant :

Vitesse de vent maximale journalière annuelle d’Odienné de 2015 à 2021

Années 2015 2016 2017 2018 2019 2020 2021
Vitesse de vent (km/h) 51 96 75 145 47 165 37

La vitesse moyenne est de 88 km/h , mais la valeur retenue est de 100 km/h.
Les résultats obtenus sont consignés dans le tableau ci-dessous.
Hauteur de vague (h) Revanche_GAILLARD
Fètch(F) Vitesse vent (w)
MOLITOR STEVENSON MOLITOR STEVENSON
Bireko 2,58 100 0,94 0,84 1,29 1,15

Conclusion :
La valeur de la revanche retenue sera de R=1,3m en tenant compte des tassements
éventuels de la digue. Cette valeur est supérieure à celle de l’APD (0,95m) de 0,35m.

2.5 Détermination de la cote crête

La valeur de la cote crête est donnée par l’expression suivante :
Cote crête = Cote PHE + Revanche (R)

COTE CRETE = 414,15+ 1,30  COTE CRÊTE = 415,45 m

Conclusion :
La digue de BIREKO est caractérisée les cotes suivantes :

- Cote R.N : 413,55 m

- Cote PHE : 414,15 m.
- Cote crête : 415,45 m
2.6 Hauteur de digue (H)
.
H= cote crête – cote lit mineur

H= 415,45 – 405,45  H= 10 m

P a g e 24 | 44
2.7 Hauteur de la retenue normale (h)
h= cote RN – cote lit mineur

h= 413,55 - 405,45  h= 8,10 m

2.8 Largeur de crête (C)

La largeur de crête devra être suffisante pour permettre la circulation, l’accès et le contrôle
des ouvrages. Sa détermination est faite à partir de deux formules suivantes :
 Selon Knappen : 𝐶 = 1,65√𝐻 ==> 𝐶 = 5,21𝑚

 Selon PREECE : 𝐶 = 1.1√𝐻 + 1 ==> 𝐶 = 4,48 𝑚

On prend une largeur de crête 5 m.

2.9 Pentes de talus amont et aval

 Talus amont
La pente de talus de l’APD est adoptée (1V/2H)

 Talus aval
Selon nos expériences en construction de barrages en terre, les ravinements apparaissent
généralement sur les talus aval pris à (1V/2H) pour les digues d’hauteur avoisinant les 10m.
De ce fait, nous proposons un talus de 1V/2,5H afin d’assurer une pente plus douce et
minimiser les ravinements sur le talus due aux précipitations.

2.10 Emprise de la digue (L)

On détermine l’emprise de la digue (L) géométriquement :

on a : L  X1  C  X 2
H H
et : tan(  amont )  ; tan(  aval ) 
X1 X2
P a g e 25 | 44
 H H
 X1  et X2 
tan(  amont ) tan(  aval )
H H
 L C 
tan(  amont ) tan(  aval )

 1 1 
L  C  H    
 tan(  amont ) tan(  aval 
)

Avec ;

L : largeur d’emprise en (m).

C : largeur de crête en (m) ; C = 5 m.

H : Hauteur de la digue en (m) ; H=10 m

 amont : l’angle de talus amont (1V/2H) ; tan(  amont ) 

1
2

 aval : l’angle de talus aval (1V/2,5H) ; tan(  aval ) 

1
2,5

L=50 m

CHAPITRE III : DIMENSIONNEMENT DES COMPOSANTES DE LA DIGUE

L’objectif de cette étude est de vérifier les dimensions de chaque élément pour assurer le
fonctionnement et la stabilité de la digue

Dans ce chapitre on vérifie :

- Les dimensions de la tranchée d’ancrage.

- L’étudie du tracer de la courbe des infiltrations dans la digue.
- L’évaluation de débit des fuites.
- Le Dimensionnement de l’épaisseur de filtre.
- La détermination de la protection de talus amont de la digue.
- La détermination de la protection de talus aval de la digue.
- L’évaluation de la stabilité de pente amont de la digue

P a g e 26 | 44
3.1 Tranchée d’ancrage
La tranchée d’ancrage a pour objectif d’assurer la fondation et l’étanchéité du barrage. Elle
allonge les lignes de fuite et évite le phénomène d’érosion interne régressive qui s’auto
accélère au fur et à mesure de l’augmentation du gradient hydraulique h/L (h : hauteur d’eau
à la retenue normale et L : Emprise de la digue).

La hauteur d’eau (h) à l’amont du barrage étant imposée, on déterminera la longueur (Lv) des
lignes de fuite dans la fondation qui est présentée dans le schéma suivant :

On utilise la règle de LANE qui proposé par la condition suivante :

1 / 3  L  LV
C
h

Où ; L : Emprise de la digue en (m) ; L = 50 m.

Lv : Longueur des cheminements verticaux en (m).

h : : hauteur d’eau (m) ; h= 8,70 m

C : Coefficient de LANE qui dépend de la nature du terrain. Déterminé par le tableau suivant

P a g e 27 | 44
 COEFFICIENT DE LANE

N° NATURE DU TERRAIN C
1 Sables fins et limons 8,5
2 Sables fins 7
3 Sables moyens 6
4 Gros sables 5
5 Petits graviers 4
6 Gros graviers 5
7 Mélange de graviers et de gros galets 2,5
8 Argile plastique 3
9 Argile consistante 2
10 Argile dure 1,8

 On prend la coefficient C = 2 => Argile consistante déterminée à partir les résultats

géotechniques du sondage sur l’axe de la digue.

1 / 3  L  LV 16,67  LV
C   2  LV  0,73 m
h 8,70

Donc la profondeur de tranchée d’ancrage doit respecter la profondeur proposée par la

compagne de sondage.

 •La tranchée sera exécutée le long de l’axe de la digue, la profondeur doit atteindre le
bon sol approuvé par l’étude de résultats des sondages contradictoires faite
précédemment.
 •La largeur de base est de 4 m. Elle assure la facilité des travaux de terrassement et
compactage.
 La pente de talus sera de 1V/1H assurant la stabilité des parois

3.2 Etude des infiltrations et du dispositif de drainage

La perméabilité dans un barrage en terre existe toujours, même si les pertes d’eau par
infiltration sont faibles et pour déterminer ces infiltrations, il est impératif d’évaluer les
éléments suivants :

- La ligne de saturation ou phréatique,

- Le débit de fuite de l’eau par infiltration,

P a g e 28 | 44
 - L’épaisseur du drain.

3.2.1 Tracé de ligne de saturation

Dans le cas d’une digue homogène drainée, reposant sur une assise imperméable, KOSENY a
montré que la ligne de saturation (ligne phréatique) est une parabole d’axe horizontal, dont
le foyer est l’extrémité amont du drain, auquel se raccorde la ligne phréatique

Elle a pour équation : y2 -y02- 2xy0 = 0

Avec : Y0 = h2  d 2 - d

Les paramètres d, h, Y0 sont représentés sur la figure suivante

A : ligne phréatique (- - - - - -)
B : parabole de KOSENY ( )
La parabole de base coupe le plan d’eau au point A à une distance (AB = 0.3b). La ligne de
saturation dans le massif se confond sensiblement avec la parabole de KOSENY.

Pour tracer cette ligne il faut raccorder la parabole de KOSENY au point B du plan d’eau par
une courbe normale au parement amont en B [figure (7-1)].

HYPOTHESES

L : Emprise du barrage (L = 50 m)

C : largeur de la crête du barrage (C = 5 m) ;

P a g e 29 | 44
 H : hauteur totale de barrage (H= 10 m) ;

hlam : charge maximale sur déversoir (hlam = 0,6 m)

h : hauteur de la retenue normale (h = 8,10 m)

Pentes des talus tg αamont = 1/ 2 et tg αaval = 1/ 2,5

A partir du graphique ci-dessus on détermine le traçage des courbes de saturation

 Calcul de la longueur du drain (Ld)

La longueur de drain et égale à (1/3) de largeur de la digue elle est estimée par la formule
1
suivante : Ld =  L ; avec : L : Emprise de la digue en (m) ; L = 50 m
3

Soit Ld =16,67 m

 Calcul de (b) :

C’est la projection horizontale de la ligne mouillée du talus amont. Elle est définie par la
formule suivante :

b = H  R 
1
tg amont

Avec ; H – R = h + hlam (hauteur d’eau à l’amont du barrage)

tg αamont = 1/ 2 : pente du talus amont.

 b = 17,4 m

 Calcul de (d ) :

d = L – l – 0.7 b  d = 21,15 m

 Calcul de (Y0 ) :

Y0 = (h  h max) 2  d 2 - d  Y0 = 1,72 m

D’où l’équation de la courbe de débit de fuite :

y 2  y 02 𝑦 2 −(1,71)2
x = 2×1,71
 x = 0,291 y² - 0,859
2 yo

P a g e 30 | 44
3.3 Evaluation de débit de fuite (q)
Le débit de fuite à travers une digue se décompose en deux types :

- Débit de fuite à travers le barrage (q1)

- Débit de fuite à travers le sol de fondation (q2)

3.3.1 Evaluation de débit de fuite à travers le barrage (q1)

En appliquant la loi de DARCY, on évalue approximativement le débit de fuite à travers le
corps du barrage. Pour un mètre linéaire du barrage, la loi de DARCY s’écrit comme suit :

h
q1 = K1 * S * i = K1 * S *
L1

Avec; q1 : débit par mètre linéaire du barrage ;

S : longueur mouillée du parement amont du barrage,

K1 : perméabilité horizontale du corps de digue en graveleux latéritiques argileux compactés :

K1 = 1,8.10-8 m/s (d’après l’APD).

h : la hauteur d’eau à l’amont du barrage, h = 8,65 m.

L1 : longueur moyenne des lignes de courant. Elle est prise comme étant la distance séparant
l’extrémité amont du drain horizontal et le centre de la digue mouillée du parement amont.

 Calculons L1 et S

2
h
L1  d  0.2b    
2
et S b2  h2
2

 L1= 25,13 m et S = 19,99 m

Le débit par mètre linéaire du barrage q1 est : q1 = 1,245 x 10-7 m3/s/m

3.3.2 Evaluation de débit de fuite à travers le sol de fondation (q2)

Le débit q2 à travers la fondation à la charge maximale est donnée par la formule suivante :

q2 = K2 (T – F) * h / (B+ 2F)

P a g e 31 | 44
Où K2 : Coefficient de perméabilité des fondations K2 = 2,4.10-6m/s (d’après l’APD).

T: Profondeur de la couche perméable, T=5 m

F: Profondeur de la para-fouille, F = 4,73 m

h: Hauteur d’eau dans le barrage à la retenue normale, h= 8.6 m

B: Largeur de base de la para-fouille, B = 4 m

Le calcul donne q2 = 4,19 10-7m3/s/m

 Evaluation du débit total q :

q = q1 + q2 = 5,43 10-7 m3/s/m.

3.4 Calcul de l’épaisseur e du drain-filtre (e)

On prend un débit q égal au double du débit de fuite escompté. Ainsi, l'épaisseur (e) du filtre
à prévoir sera :

q  Ld
e  2
K
Avec ; ld : Longueur du drain , ld = 16,67 m

K : Coefficient de perméabilité du sable ; K = 10-4 m/s

On obtient environ e = 0,60 m. La valeur de l’APD est reconduite e= 0,8m.

Le tapis sera en couche interne de gravier enveloppé par du sable composé de :

P a g e 32 | 44
 - 20 cm sable 0/4
- 40 cm gravier 4/15
- 20 cm sable 0/4

3.5 Protection talus aval

La protection proposée est la réalisation d’un gazon semer sur une couche de terre végétale
de 20 cm d’épaisseur. Cette protection permet :

- Protéger le corps de la digue des variations climatiques

- Éviter les ravinements

3.6 Protection talus amont

La protection du talus amont protège la digue contre le batillage ou action érosive des vagues.
Le choix proposé par l’APD est maintenu. L’opération consiste à mettre en place un
enrochement en rip-rap avec une couche de transition en graveleux latéritique suffisamment
grossier d’épaisseur 20 cm.

 Détermination d’épaisseur de couche d’enrochement avec granulométrie

Les dimensions minimales des blocs sont en fonction de la hauteur des vagues déterminée à
partir de l'abaque de l'U.S Army corps of Engineers suivant :

EPAISSEUR MINIMAL DE LA COUCHE

HAUTEUR DES VAGUES (m) D50 MINIMUM (m)
D'ENROCHEMENT (m)

0 à 0,3 0,3 0,2

0,3 à 0,6 0,4 0,25
0,6 à 1,2 0,45 0,3
1,2 à 1,8 0,55 0,4
1,8 à 2,4 0,77 0,45
2,4 à 3 0,8 0,55
Abaque de l'U.S Army corps of Engineers

La hauteur des vagues (Hv) telle que déterminée dans la partie calcul de la revanche se situe
dans l’intervalle suivant : 0,6 < HV < 1,2  la couche minimale d’enrochement est 0,45 m et
les blocs auront au moins 30 cm de taille.

P a g e 33 | 44
3.7 Etude de stabilité de talus amont
Notre étude se fera sur la base de la méthode de Bishop. A l’aide du logiciel Talren, la stabilité
des talus a pu être vérifiée avec un coefficient de sécurité F ≥ 1.5.
Données d’entrée du projet

Stabilité du talus aval

P a g e 34 | 44
Stabilité du talus amont

Conclusion : Les pentes de talus sont stables

3.8 Butée de pied

 Butée de pied amont

La butée de pied a pour fonction d’ancrer la protection amont de la digue et de prendre une
partie la charge des poussées engendrées par le poids propre des pierres. La hauteur
d’ancrage sera de 50 cm au minimum.

 Butée de pied aval

Elle permet la continuité du drainage et le bloque de la désagrégation du tapis filtre. La

granulométrie est de 5mm – 400 mm. Elle sera une continuité granulométrique du gravier du
tapis filtre. La posé sera manuelle. La hauteur est de 1m supérieur à la hauteur du tapis filtre.
Il est préférable de poser un géotextile qui enveloppe la butée pour éviter la migration des
particules solides.

P a g e 35 | 44
3.9 Fossé de pied aval
Le fossé de pied aval permet de canaliser et évacuer les eaux de drainage et de ruissellement.
Ce fossé de section trapézoïdale de talus 1/1 sera en béton légèrement armé avec une base
de 50 cm, et une profondeur variable.

CHAPITRE IV : ETUDE DE LA CONDUITE DE VIDANGE

Dans ce chapitre nous étudions :

- Le temps de vidange.

- Le dimensionnement du bassin de dissipation type impact.

4.1 Détermination de temps de vidange

On retient le choix de l’APD pour une conduite de diamètre 800 mm.

 Evaluation du temps de vidange

𝑉
𝑇=
𝑄

Avec : T: temps de vidange partiel en (secondes)

V : volume d'eau écoulée entre deux niveaux en (m3)
Q : débits moyens entre deux niveaux en (m3/s)

- Le volume d'eau écoulée entre deux niveaux est déterminé à partir la de courbe
hauteur-surface.
- Le débit moyen entre deux niveaux est exprimé par :

𝑄 = 𝑉𝑐 × 𝑆

Où ; 𝑉𝑐 : vitesse de la conduite en (m/s)

𝜋D²
S : Section de la conduite ; 𝑆 = ; avec D = 800 mm
4

 On va déterminer par la suite la formule de vitesse de la conduite.

Les conditions de vidange doivent tenir compte des pertes de charge linéaires et singulières
VC2
suivantes : J L  LC  2
K c  RH4 / 3
Pertes de charge linéaire :

P a g e 36 | 44
Avec ;
J L : pertes de charge linéaire (en m) ;
Lc : longueur de la conduite (m) ;
Vc : Vitesse de la conduite (m/s) ;
Kc : Coefficient de rugosité de Strickler (ml/3/s)
RH : rayon hydraulique (m) ; pour les conduites ayant une section circulaire. R  D
H
4

VC2
Pertes de charge singulière : J S  (0,55  C ) 
2g

Avec ;
JC = pertes de charge singulière (m) ;

Vc : Vitesse de la conduite (m/s) ;

C : Coefficient de débit, dépendant du degré d’ouverture de la vanne, fourni dans le tableau

suivant :

% Ouverture de la vanne 30 50 70 80 100

C 120 19 3,5 1,5 0,2
Source : (Touihri, 2004)

La vidange tient compte aussi des contraintes suivantes :

La vitesse dans la conduite de vidange ne doit pas être ni très élevée ni trop faible. On admet
des vitesses de 3 à 7 m/s, et on veillera à ce qu’en fin de vidange, on puisse garder une vitesse
suffisante pour éviter le colmatage par les dépôts solides. La vitesse de 0,4 m/s semble être
un minimum. Mais dans notre étude, nous admettons que la vitesse maximale est de 7 m/s et
une vitesse minimale est z de 2 m/s permettant l’auto curage de la conduite en faible charge.

En appliquant le théorème de BERNOUILLI, assimilant la retenue à un déversoir alimentant

une conduite qui débouche à l’air libre, l’on a :
2
VRN V2
Z RN  PRN   Z C  PC  C  J L  J S (1)
2g 2g

PRN : la pression à la surface de la retenue ; PRN = Patm

Pc : la pression en aval de la retenue normale ; Pc = Patm

P a g e 37 | 44
Zc : Cote aval de la conduite de vidange ;
ZRN : cote de la retenue normale ;
JL : perte de charge linéaire (m) ;
Ks : Coefficient de MANNING-STRICKLER ;
LC : longueur de la conduite (m) ;
D : diamètre de la conduite (m) ;
Vc : vitesse dans la conduite (m/s) ;
Js : Pertes de charge singulières (m) ;

VC2
Z RN  Z C   JL  JS (2)
La relation (1) devient : 2g

Alors la vitesse Vc est donnée par la formule suivante :

2g  Z
La relation (2) devient : Vc  D
2 g  LC et ; Z  Z RN  Z C ; RH 
1  (0,55  C )  4
KC2  RH4 3
Avec :

𝑉𝑐 : vitesse de la conduite (m/s) ;

Z : différence des niveaux d’eaux en amont et en aval ;(m)

Zc : Cote aval de la conduite de vidange ;

ZRN : cote de la retenue normale ;
g : accélération de la gravité en (m/s²) ; g = 9,81 m/s2
C : Coefficient de débit dépendant du degré d’ouverture de la vanne; C = 0,2
Kc : Coefficient de rugosité de Strickler en (ml/3/s) ; KC = 70

Lc : longueur de la conduite en (m) ; LC = 60 m

RH : rayon hydraulique en (m) ; RH = 0,2 m

D : diamètre de la conduite en (m) ; D = 0,8 m

Le calcule a été effectués pour une conduite de diamètre 800 mm et une longueur de 60 m.
Le tableau suivant donne les volumes partiels, les vitesses, les débits et les temps de vidange.

Evaluation du temps de vidange

P a g e 38 | 44
 Cote moy Z (m) Q (m³/s) V (m³) T (s) T (jours)
413,55 6,75 2,8344 1 510 366 532 862 6
413 6 2,6723 1 374 232 514 244 6
412 5 2,4395 1 097 622 449 938 5
411 4 2,1820 838 862 384 455 4
410 3 1,8896 560 073 296 393 3
409 2 1,5429 297 149 192 594 2
408 1 1,0910 104 669 95 940 1
407 0,25 0,5455 3 604 6 607 0
406,5 0 0 0 0 0
Total de temps en jours 29

4.2 Dimensionnement de bassin de dissipation type impact

Dans ce type de bassin la dissipation est réalisée par le choc du jet incident sur un écran
vertical qui est en général une poutre traversant le bassin de part en part, Le jet réfléchi vers
l'arrière, heurte le jet incident et crée des tourbillons qui absorbent l'énergie,
Ce type d'ouvrage est extrêmement efficace,
Les dimensions principales sont en relation avec le débit. Elle se déduisent toutes de la largeur
(l) du bassin que l'on détermine à l'aide de l'abaque suivant :

P a g e 39 | 44
A partir l’abaque on détermine (l) la largeur du bassin pour un débit débit maximal donné par
le tableau ci-dessus Q = 3 m3/s  l = 3m

Les autres dimensions du bassin sont déduites à partir de la valeur de l=3m et l’abaque.

P a g e 40 | 44
CHAPITRE V : DIMENSIONNEMENT DE L’EVACUATEUR DE CRUES
L’objectif dans ce chapitre est axé sur le dimensionnement des éléments en gras.

L’évacuateur de crues est constitué de

- chenal d’amenée.
- radier submersible avec des rampes d’accès à 10% dont la largeur en base représente
la longueur de déversement.
- chenal
- coursier.
- convergence.
- bassin de dissipation.
- chenal évacuateur.
5.1 Dimensionnement du chenal

La pente du chenal doit être telle que le régime d’écoulement y soit fluvial. Il est de forme
rectangulaire.
Q²
L’écoulement fluvial est obtenu pour Yn > Yc avec : Yc  3 et Ic=0,25%.
gL ²

La profondeur normale se déduit de la formule de Manning Strickler : Q  K .S . I .R . Avec

2/3

K=70
Le tracé du chenal devra se faire avec une pente inférieure à 0,25% pour que l’écoulement y
soit fluvial. Ainsi, une pente I = 0,20% sera adoptée dans le chenal, on obtient :
Désignation Q (m3/s) L(m) Ic(%) Yc (m) Yn (m)
Valeurs 31 20 0,25 0,63 0,67
Yn est bien supérieur à la profondeur critique. Le régime y est donc fluvial pour I=0,20%.

5.2 Dimensionnement du coursier

La pente du coursier doit être telle que le régime d’écoulement y soit torrentiel. Il est de forme
rectangulaire. L’écoulement torrentiel est obtenu pour I > Ic avec : Ic=((Q/K.R2/3.S)² et Yc>Y .
La largeur du coursier est déterminée par l= 2yc. On procédera par approximations successives
à partir d’une première valeur de l= 2yc, puis en affinant progressivement de sorte à obtenir
une hauteur de coursier adapter à notre terrain et qui tient compte de la revanche à respecter.
Les différentes valeurs obtenues se résument dans le tableau ci-après.
Conditions pour un
Données Conditions réelles
régime torrentiel
Longueur
Désignation Q (m3/s) l (m) IC(%) Yc(m) I(%) Y(m)
(m)
Valeurs 31 4 0,40 1,83 2,5 0,94 275m

P a g e 41 | 44
En ce qui concerne la profondeur du coursier, il convient d’ajouter une revanche pour éviter
les risques de débordement de l’eau. ‘‘ Design of Small Dams’’ propose la formule suivante :
Revanche=0,6+0,05.V 3√𝑦
Avec : V = vitesse en m/s et y =tirant d’eau en m
Profondeur Profondeur Epaisseur
y (m) v (m/s) R (m) retenue (m) (m)
Désignations calculée (m)
Valeurs 0,94 8,24 1,00 1,94 2 0,20*
*L’épaisseur des voiles du coursier est prise à H/12 soit 2/12=0,17m. On retient 20cm.

5.3 Dimensionnement du convergent

Un convergent faisant passer la section de la longueur de déversement à la largeur du coursier,
il est donné par la formule suivante :Lconvergent= 2,5(Ldéversant - lcoursier);
Longueur convergent
Sites
Longueur déversant (m) Largeur coursier (m) (m)
Bireko 20 4 40

5.4 Dimensionnement du bassin à ressaut

La résolution du système d’équation ci-dessous permet de calculer le nombre de Froude dont
dépendent directement les caractéristiques du ressaut.
V  2 g (0,8( H  h)  y )
V
Q F
y gy
L 'V
H et h respectivement la hauteur de chute (m) et la charge sur le déversoir(m)
V et y1 sont la vitesse (m/s) et le tirant d’eau(m) au pied du coursier
L’: largeur du coursier à sa jonction avec le bassin
Q : le débit laminé (m3/s)
F : le nombre de Froude
Désignation H(m) h(m) L'(m) Q (m3/s) V (m/s) Y1(m) F
Valeurs 9,5 0,60 6,0 31 12,38 0,37 6,10
Ces valeurs de F, nous indique que le bassin à adopter pour la dissipation de l’énergie de la
chute est le bassin USBR de type II. Les différents paramètres de dimensionnement sont
déterminés graphiquement par abaque.
Les caractéristiques du bassin de dissipation suivant les sites sont :
Désignation hauteur blocs chicanes Hauteur seuil terminal
Y2(m) Longueur(m) Largeur(m)
(m) (m)
Valeurs 3,34 8 6 0,70 0,60

P a g e 42 | 44
5.5 Dimensionnement du chenal débouchant dans le drain principal
Le canal de restitution est la partie de l’évacuateur de crues qui raccorde le bassin
d'amortissement avec le drain en aval.

Cette partie est de forme trapézoïdale avec un fruit de 1/1, a une pente du canal de I =
0,35% et il est en enrochement.

Yn est déterminé à partir de la formule de Manning Strickler

2/3
𝑄𝑙𝑎𝑚 = 𝐾. 𝑆𝑚 . 𝑅𝐻 . 𝐼1/2

𝑆𝑚
Où 𝑆𝑚 = 𝑌𝑎 × (𝐿 + 𝑌𝑎 ) ; 𝑅𝐻 = ; 𝑃𝑚 = 𝐿 + 2√2 × 𝑌𝑎
𝑃𝑚

Qlam : Débit laminé en (m3/s) ; Qlam= 31 m3/s

K : coefficient de Manning Strickler ; K = 40 (Pour rivière)

Sm : Section mouillée en (m²) ;

Pm : Périmètre mouillé en (m) ;

Rh : Rayon hydraulique en (m) ;

L : Largeur de chenal d’évacuation en (m) ; L= 8 m ;

I : pente en (m/m) ; I = 0,35 % ;

m : fruit de berge ; m=1

Par un calcul itératif, la profondeur yn obtenue est de 1,3m. La profondeur du chenal

d’évacuation peut être prise à 1,5m en tenant compte d’une revanche de 20cm.

La longueur de protection du lit en aval du bassin de dissipation est donnée par la formule
suivante :

Lp ≥ 2L

Avec L : Longueur de bassin de dissipation ; L = 8 m

On obtient : Lp ≥ 2 x 8 = 16 m

P a g e 43 | 44
ANNEXES

P a g e 44 | 44