Étude de drainage Cité FAD Djibouti
Étude de drainage Cité FAD Djibouti
Août 2024
L’étude a pour but de définir les caractéristiques de l’ouvrage de drainage dans la cité FAD à Djibouti.
Pour mener à bien cette mission, Hydroterra s’est basé sur éléments suivants :
■ L’hydrologie et la climatologie du bassin versant à partir de données météorologiques et de
topographies satellitaires
■ La topographie à l’échelle des ouvrages existants et projetés
■ Une analyse hydraulique de la formation et de l’intensité des écoulements sur la base de l’analyse
hydrologique
La localisation de la cité FAD n’indique pas de risque d’exposition majeur de la cité à des inondations car aucun
bassin versant ne présent en amont de la cité. Toute fois pour éviter la dégradation des voiries et espace bâtis
lors des averses importantes, il est primordial de concevoir un ouvrage de drainage qui permette l’évacuation
des eaux ruisselant de la cité d’une manière contrôlée.
Le présent rapport présente l’ensemble des analyses et méthodes qui ont été menées, à savoir :
Les données climatologiques à l’échelle du pays
La méthode utilisée pour les levées topographiques
Les caractéristiques géomorphologiques des bassins versants de la cité FAD
La méthode pour la détermination du débit maximal à évacuer
Les caractéristiques de l’ouvrage proposé
2.3 INSOLATION
La durée d’insolation est très importante tout au long de l’année. Elle atteint son maximum au mois
d’Octobre et son minimum au mois de Janvier à Djibouti.
Moyenne
J F M A M J J A S O N D
annuelle
Djibouti 10,5 10,0 9,8 9,2 11,2 9,4 8,8 10,4 10,7 10,7 8,7 9,8 9,9
Table 2 – Durée de l’insolation à Djibouti (h)
2.4 VENTS
Comme précisé précédemment les vents sont des alizés du nord-est en hiver et un vent équatorial d’ouest
en été. Le vent est assez fort sur la région de Djibouti avec une moyenne annuelle de 3,7 m/s.
Moyenne
J F M A M J J A S O N D
annuelle
Djibouti 3,7 3,8 3,6 3,5 2,8 3,3 4,9 4,7 3,4 3,2 3,8 3,8 3,7
Table 3 – Vitesse moyenne des vents à Djibouti (m / s)
2.6 EVAPORATION
La moyenne annuelle d’évaporation est estimée à 2192 mm/an. Plus la température est forte et plus
l’évaporation est intense, le maximum est donc mesurée au mois de juillet.
Moyenne
J F M A M J J A S O N D
annuelle
Djibouti 109 95 112 105 140 261 428 378 204 136 114 112 2192
Table 5 - Evapotranspiration à Djibouti (mm)
2.8 PLUVIOMETRIE
Les analyses pluviométriques pour la cité de FAD ont été réalisées à partir des données disponibles des
stations pluviométriques de Djibouti Aérodrome et Djibouti Serpent, les stations les plus proches de la
cité. Le tableau suivant montre la série des pluies maximales de 24 h à Djibouti Serpent (1924-1988) et la
série de pluies journalières maximales à Djibouti Aérodrome (1980 – 2003). Sur la période commune (1980
– 1988) des différences entre les valeurs des deux stations sont présentes : à titre de précaution, sur cette
période le maximum des deux séries de valeurs ont été considérés.
Pour permettre une analyse hydraulique précise et élaborer les dessins des nouveaux ouvrages sur
Autocad, nous avons réalisé une étude en photogrammétrie pour acquérir une couverture topographique
vaste (100 ha) et précise (précision centimétrique).
3.1 GEODESIE
Le système de référence utilisé est le suivant :
Ellipsoïde de référence : IAG GRS 1980
Système géodésique : WGS84
Projection : UTM 38 N
Sytème altimétrique : Local
Figure 2 – Principe schématique de fonctionnement d'acquisition par drone et GNSS RTK Différentiel
L'acquisition par drone permet d'obtenir des images géo Tiff à 2.5 cm de résolution, les différents
traitements permettent d'obtenir un orthophotoplan et un modèle numérique de terrain avec une
extraction très précise des courbes de niveau.
Cette acquisition de la topographie locale sur 100 ha, va permettre de faire des simulations hydrauliques
et de projeter les ouvrages hydrauliques sur Autocad.
Exutoire EST
Exutoire SUD
La situation topographique du terrain naturel de la cité indique deux sens de ruissèlement d’eau pluviale
disposée en plusieurs petits bassins versant, le premier écoulement du Nord au Sud et le second de
l’Ouest à l’Est. Les exutoires de ces écoulements sont tous orientés vers l’oued d’Ambouli situé à moins
d’1 km de la cité.
0.018301
Pente moyenne I (m/m)
BV02
S ou A
Aire du bassin versant 0.07231
(km2)
769.68
Longueur du plus grand chemin L (km)
0.02240758
Pente moyenne I (m/m)
BV03
S ou A
Aire du bassin versant 0.04675
(km2)
647.28
Longueur du plus grand chemin L (km)
0.0183609
Pente moyenne I (m/m)
BV04
S ou A
Aire du bassin versant 0.3019
(km2)
531.96
Longueur du plus grand chemin L (km)
0.02026162
Pente moyenne I (m/m)
BV05
S ou A
Aire du bassin versant 0.10978
(km2)
Longueur du plus grand chemin L (km) 1325.65
0.01769467
Pente moyenne I (m/m)
BV02
S ou A 0.29609
Aire du bassin versant
(km2)
2229.27405
Longueur du plus grand chemin L (km)
0.01586068
Pente moyenne I (m/m)
Afin de réaliser des ouvrages de retenue d’eau d’écoulement superficiel, la présente étude hydrologique
permet d’estimer les apports pluviométriques sur les sites des ouvrages le temps de retour adopté.
[Link] Préambule
L’intensité de précipitation est calculée par l’équation :
𝑖 = 𝑎 ∙ 𝑡𝑛−1
Où les paramètres t et t sont dépendants des paramètres statistiques de l’échantillon moyen µ(x) et
d’un écart quadratique moyen (x) par le biais de l’expression :
1,283
𝛼𝑡 =
𝜎(ℎ𝑡 )
𝜀𝑡 = (ℎ𝑡) − 0,450 ∙ (ℎ𝑡)
En introduisant dans la loi de Gumbel l’expression du temps de retour Tr :
1
𝑇𝑟 =
1 − 𝐹(ℎ 𝑡)
On a l’expression de la valeur de hauteur de précipitation en relation de la durée "t" et du temps de
retour "Tr" :
1 𝑇𝑟
ℎ ( ) = 𝜀 − ln (𝑙𝑛 )
𝑡 𝑟 𝑡
𝛼𝑡 𝑇𝑟 − 1
Table 10 – Djibouti ville – Pluie maximale journalière annuelle (mm) pour différents temps de retour
Ainsi, la pluie maximale journalière avec un temps de retour de 10 ans est estimée à 103.6 mm/24h (pluie
de projet).
[Link] La relation entre la durée de la pluie et la variation de l’intensité de pluie (Courbe IDF)
Comme nous avons vu la relation entre la hauteur de la pluie et la durée de l’événement est représenta
épar la formule :
ℎ = 𝑎 ∙ 𝑡𝑛
La représentation de cette relation dans le plan h, t pour une valeur constante de la fréquence F est
indiquée comme courbe IDF (Intensité-Durée-Fréquence). Rappelons ici que la fréquence d’évènement F
est égale à l’inverse du temps de retour :
𝐹 = 1⁄ .
𝑇𝑟
Les coefficients « a » et « n » sont en fonction du temps de retour Tr, c’est-à-dire la période caractéristique
de chaque valeur de hauteur de pluie durant laquelle on prévoit que cette valeur peut être dépassée au
moins une fois.
Comme indiqué, dans le cadre de la présente étude, le temps de retour a été considéré de 10 ans.
Etant donné que le phénomène présente une croissance plus rapide pour une courte durée de
précipitation, la courbe IDF avec temps de retour 10 ans a été élaborées séparément pendant des durées
inférieures ou supérieures à t=1 h.
Les valeurs des coefficients a et n ont été estimées par deux méthodes différentes :
A. L’intensité des pluies de courte durée a été estimée par l’application de l’abaque basé sur des
14 Cité FAD - Djibouti
données expérimentales par Bell F. (1967) et Goswani A. (1973). Les valeurs reportées dans le
tableau suivant montrent le rapport entre la hauteur de pluie de durée (t) et le total journalier de
précipitation. On voit que presque tout le total de précipitation se mesure pendant les premières
six heures de précipitation.
T 0.25 1h 3h 6h 12 h 24 h
Selon ces hypothèses pour une pluie de 103.6 mm/24h les valeurs des coefficients sont les suivantes :
▶ Pour de temps de pluies inférieurs à une (1) heure :
a = 14,116
n = 0,3266
▶ Pour de temps de pluies majeurs d’une (1) heure :
a = 26,957
n = 0,196
B. L’intensité des pluies de courte durée a été aussi estimée par l’analyse des évènements extrêmes
de type cyclonique encourus à Djibouti le 20 mai 2018 et entre le 21 et le 23 novembre 2019. Les
résultats sont montrés dans le tableau suivant :
h max (mm) 8.6 15.2 27.8 31.4 31.4 31.4 31.4 31.6 31.6 31.6 31.6
21/11/19
0.272 0.481 0.880 0.994 0.994 0.994 0.994 1.000 1.000 1.000 1.000
h max (mm) 10.8 19.8 22.4 35 49.8 50.2 50.2 51.2 51.2 51.2 51.2
22/11/19
0.211 0.387 0.438 0.684 0.973 0.980 0.980 1.000 1.000 1.000 1.000
h max (mm) 10.8 14.4 19.4 44 44.2 44.6 44.6 44.6 50.4 50.4 50.4
23/11/19
0.214 0.286 0.385 0.873 0.877 0.885 0.885 0.885 1.000 1.000 1.000
Max 0.272 0.510 0.880 0.994 0.994 1.000 1.000 1.000 1.000 1.000 1.000
Pour une pluie de 103.6 mm/24h les valeurs des coefficients sont les suivantes :
▶ Pour de temps de pluies inférieurs à une (1) heure :
a = 2,8902
n = 0,8463
▶ Pour de temps de pluies majeurs d’une (1) heure :
a = 83,08
n = 0.0328
Les valeurs des coefficients estimées avec la méthode B (Pluies cycloniques 2018-2019 à Djibouti),
fournissent des valeurs d’intensité pluvieuse légèrement plus élevées que celles estimées avec la
méthode Bell/Goswani.
Exutoire SUD BV 1 BV 2 BV 3 BV 4 BV 5
0.000325 0.02240758 0.02240758 0.02240758 0.02240758
L m 461.80 769.68 647.28 531.96 1325.65
J m/m 0.018301 0.02240758 0.0183609 0.02026162 0.01769467
tc h 0.17 0.23 0.22 0.18 0.39
min 10.25 14.05 13.27 10.99 23.38
Exutoire EST BV 1 BV 2
0.000325 0.000325
L m 1960.4507 2229.27405
J m/m 0.0140205 0.01586068
tc h 0.58 0.61
min 34.57 36.39
En pratique, l’aire drainée est composite et le coefficient résultant (coefficient élémentaire) est la somme
pondérée des coefficients unitaires.
Pour le territoire de Djibouti, et en particulier pour le bassin versant de l’Ambouli, une étude,
commissionnée par le CERD de Djibouti et menée par le groupement SCE-EPTISA en décembre 2013,
indique que pour l’Oued Ambouli le coefficient de ruissellement peut être considéré de l’ordre de 30% à
40 pour les crues courantes assez fortes, et de l’ordre de 45% à 55 % pour des crues fortes exceptionnelles.
Conformément à ces indications et à titre de précaution, pour les bassins versant en amont de la cité FAD,
le coefficient de ruissellement adopté est de 50%.
[Link] Calcul du débit maximal
A partir des résultats obtenus il est possible de calculer le débit maximal par la formule introduite au début
du chapitre :
𝑄 = 0,00278 ∙ 𝑖 ∙ 𝐾 ∙ 𝐶 ∙ 𝑆
Ecoulement SUD BV 1 BV 2 BV 3 BV 4 BV 5
0.00278 0.00278 0.00278 0.00278 0.00278
tc min 10.25 14.05 13.27 10.99 23.38
a 2.8902 2.8902 2.8902 5.0271 5.0271
n 0.8463 0.8463 0.8463 0.8463 0.8463
h mm 20.71 27.05 25.78 105.29 105.29
i mm/h 121.27 115.53 116.54 173.60 173.60
K 1 1 1 0.88 0.88
C 0.5 0.5 0.5 0.5 0.5
S ha 6.185 7.231 4.675 3.019 10.978
Q m3/s 2.04 1.16 0.76 0.64 2.33
Ecoulement EST BV 1 BV 2
0.00278 0.00278
tc min 34.57 36.39
a 2.8902 2.8902
n 0.8463 0.8463
h mm 20.71 27.05
i mm/h 121.27 115.53
K 1 1
C 0.5 0.5
S ha 5.449 5.371
Q m3/s 1.49 1.46
Les caniveaux seront conçus de sorte à collecter en partie ce débit en provenance des eaux des
ruissellements des pluies directement tombée sur la cité, et devront faciliter l’évacuation vers des
exutoires naturels suivant la topographie de la cité.
6.1 DESCRIPTION
Pour faciliter le drainage de l’eau de pluie directement tombé sur la cité, un réseau de caniveaux en béton
armé sera à réaliser. Ces réseaux devront canaliser le débit d’eau en provenance de diffèrent écoulement des sorte à
l’évacuer vers les exutoires naturels. Le débit évacué devra couler vers l’oued d’Ambolie situé à moins d’un kilomètre au Sud-
Est.
Néanmoins il faudra considéré qu’une grande partie des sens des écoulements obtenus avec l’état actuel de terrain sera
modifié à la fin de travaux des terrassements de la cité, le drainage sera donc favorisé par la pente de la voirie qui doit permettre
un écoulements suivant l’exutoire naturels vers l’oued d’Ambouli.
Avantage de la topographie du site : le site du projet est situé en dehors de zone de talweg sur un terrain
topographiquement avantageux, favorisant un ruissellement naturel des cours d’eau directement tombé sur
la cité vers l’oued principal localisé au sud-ouest à environ à 1 kilomètre.
Les débits de ruissèlements et la repartions des bassins versant dans la cité: Sur la surface de 110 ha du cité,
que 8 directions de ruissellent avec des débits de faible intensité sont répertoriés.
Cette situation permet de concevoir un système de drainage peu couté, avec une conception des petits
caniveaux trapézoïdaux sur les rives de voies bitumeux pour permettre un drainage en surface des eaux.
Ou également d’envisagé un écoulement sur les surfaces bitumeux pour les ruissèlements de faible distance
suivant une pente vers des exutoires naturel.
Dans l’idée d’intercepter et casser le débit des eauxde ruissellement en provenance de la partie nord de
la cité un premier caniveau sera conçu sur le grand axe séparant la cité en deux parties nord et sud :
Tandis que les écoulements formés sur la partie sud de la cité seront évacués naturellement sur
la chaussée suivant la pente vers l’oued d’Ambouli.
Le tronçon 1 : A à B e qui intercepte principalement les eaux des ruissèlements en provenance
du nord de la cité.
A A
Les pentes des caniveaux suivront les pentes du terrain de la cité. Le caniveau a été dimensionné en
considérant un remplissage maximal de 95%.
Le débit a été calculé à l’aide la formule suivante :
Devis quantitatif
30 000
88 175
39 686.4
720 972
43 800
SOUS TOTAL (Dollars US) 922 633,4
Détail de caniveaux
Traces de caniveaux