Analyse Hyd Globale V0
Analyse Hyd Globale V0
Phase :
RAPPORT : HYDROLOGIE
Emetteur :
B
A
0 01/10/2016 1ière diffusion N SOLIMAN S ABOU ZAID F BOUKLOU
Approuvé
Indice Date Modifications - mises à jour Emis par Vérifié par Statut
par
Identification du document :
Code affaire Division OA NOT 001 0
201600290 Emetteur Type Numéro Indice
FEUILLE DE PAGINATION
Indice Indice
Page Observation Page Observation
0 A B C 0 A B C
1 x Page de garde
2 x Pagination
3 x Sommaire
4 x
5 x
6 x
7 x
8 x
CID/PTR/DOA-Rapport M1 Page 1
Etabli le 01.10.2016 D:\Affaires\SAGE_EXPERTISE\Voie Express
CID\MISSION-1\RAPPORT\M1-Analyse hyd [Link]
Indice 0 ED des franchissements interceptés
par le projet de la Voie Express
Guelmim-Laâyoune
SOMMAIRE
8. CONCLUSION.....................................................................................................................55
Phase préparatoire : Cette phase vise à élaborer l’ensemble des données utiles pour
les phases suivantes ;
2. DOCUMENTS DE REFERENCE
Les documents qui ont servi pour la réalisation de cette note sont :
[1] Etude du plan directeur intégré d’aménagement des eaux des bassins Sud
Atlasiques, Sous mission 1.1 – Climatologie - Hydrologie– Volume 1 : Etude de base
Rapport Définitif, Juin 1996, SOGREAH–SCET MAROC, 1996 ;
[2] Etude du plan directeur intégré d’aménagement des eaux des bassins Sud
Atlasiques Sous mission 1.2 –Climatologie - Hydrologie– Volume 2 : Modélisation
des bassins, Rapport Définitif, Juin 1996, SOGREAH–SCET MAROC, 1996 ;
[3] Actualisation des études hydrologiques au niveau des sites de barrages en projet
INGEMA, 1999-2000 ;
[5] Etude d’évaluation des Ressources en Eau de Surface dans Les Bassins De Drâa Et
Guelmim, Mission 1, MAROC INGENOV, ABHSMD, 2011 ;
[8] Dossier d’exécution d’un barrage sur oued Ben khellil Al khibra, 2006
[11] Etude hydrologique et note de calcul hydraulique du Lac Labouir, GCIM 2011
[12] Lac Collinaire Saheb El Harcha, Note Déscriptif et justificative, CID 1989
[14]Etude de conception d’un barrage sur oued Tighmert dans la province de Guelmim
Etude hydrologique et hydraulique, CID 2016
[15] Actualisation de l’étude de l’ouvrage d’art sur oued Sayad situe au pk 1084+000 de
la rn1, Lixus 2015
[16] Etude de reconstruction de l’ouvrage d’art sur l’oued Ben Khlil au Pk 1206 de la RN1,
Projet 7, 2010
[17] Elaboration des standards D’ADM en Assainissement, Mission 2 : Analyse critique &
Méthodologies, LPEE & SAGE 2010
La pluie maximale de durée 24h est de l’ordre de 50mm. Elle ne dépasse pas le
maximum observé. De ce fait, l’évenement de l’année 2014 n’est pas
exceptionnelle. Les dégâts importants constatés serait principalement dûs à une
occupation du sol dans des zones à haut risque d’inondations.
Les figures suivantes présentent les postes pluviométriques d’intérêt pour l’étude. Une
demande a été formulée à l’ABHSMD pour recevoir les données de pluies actualisées.
Néanmoins l’IC a pu construite une base de données de pluies en fonction :
Des données climatologiques publiées sur internet par des sites professionnelles.
Zone Tan-Tan-Laayoune
Les mesures disponibles d’un même lieu ont été regroupé en un seul poste dit Unifié en
adoptant la valeur maximale observée.
Les Pjmax :
Assaka # # # # # # # # # # # # # # # # # # # # # # # # # # # # #
Aoreora # # # # # #
Oued Noun # # # # # # # # # # # # # # # #
AinErrahma # # # # #
# # # # # # # # # # # # #
LAAYOUNE DMN # # # # # # # # # # # # # # # # # # # # # #
# # # # # # # # # # # # # # # # # # #
TANTAN DMN # # # # # # # # # # # # # # # # # # # # # # #
# # # # # # # # # # # # # # # # # # #
Pont Draa # # # # #
# # # # # #
Sidi Ifni # # # # # # # # # # # # # # # # # # # # # # # # # # # # #
Ouijjane # # # # # # # # # # # # # # # # # # # # # # # # # # # # # #
Guelmim
# # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # #
Unifié
Assa # # # # # # # # # # # # # # # # # # # #
Foum El Hassane # # # # # # # # # # # # # # # # # # # # #
Taghjijt Unifié # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # #
Bouizakarne Unifié # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # #
Talmazirte # # # # # # # # # #
Ifrane Anti Atlas # # # # # # # # # # # # # # # # # # # # # # # # # # # # #
Tnine des Ait Erka # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # # #
Jemaa N Tirhirte # # # # # # # # # # # # # # # # # # # # # # # # # # # # # #
Les a et b de Montana
Grâce aux quantiles des Pjmax estimés pour les 19 postes pluviométriques et les
corrélations ci-dessus (r² > 0.9), les paramètres de Montana a, b pour chaque postes et
des iso-valeurs peuvent être approchés :
Les Pan :
Assaka ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ##
AinErrahma ## ## ## ## ## ## ## ## ## ## ##
Sidi Ifni ## ## ## ##
Ouijjane ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ##
LAAYOUNE
## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ##
DMN
GUELMIM
## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ##
DMN
TANTAN
## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ##
DMN
Pont Draa ## ## ## ## ## ## ## ## ## ## ##
Taghjijt ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ##
Bouizakarne EF ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ##
Ifrane de l’Anti Atlas ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ## ##
La zone d’étude présente une grande variablité de la pluie annuelle. De ce fait, il est très
imprudent de se baser uniquement sur une valeur moyenne interannuelle dans les
calculs des débits.
Les figures suivantes présentent les stations hydrométriques d’intérêt pour l’étude. Une
demande a été formulée à l’ABHSMD pour recevoir les données des débits actualisées.
de choisir la formule adéquate pour le calcul des temps de concentration (tc).
En principe, les éléments précédents se basent sur un travail minutieux et fastidieux qui
consiste à dépouiller et découper les hydrogrammes des crues majeures de chaque
station hydrométrique pour aboutir aux valeurs caractéristiques suivantes :
L’étude [7] du barrage Fask a défini les temps de monté des crues majeures au droit de
la station de Taghjijit et Assaka. Le temps de concentration adopté est voisin de 0.75*tc
Kirpich, soit un temps proche de celle du tc Espagnole (dérivée de Kirpich) :
La plus part des barrages collinaires existants sont des petits barrages (Seuils)
affectés par un taux d’envasement très important et par conséquent, ils ont
un impact insignifiant sur le régime des crues des cours d’eau ;
Les barrages, dont l’écrêtement des crues peut affecter les ouvrages de
franchissement de la voie express, sont résumés dans le tableau suivant :
BARRAGE TIGHMART OUED OUARGENNON En phase d'étude ? 61398 224410 1815 4.1 1565 2118 2116 0.1%
BARRAGE TAYERT OUED SIYYAD En phase d'étude Site abondonné 50700 229600 90 3.8 177
BARRAGE GELTA ZERGA OUED DRAA En phase d'étude Pas de visibilité -36450 176730 47000 2.7 3606
BARRAGE BEN KHLIL OUED BEN JLIL En phase d'étude Pas de visibilité -59920 170574 453 4.3 906
BARRAGE SAHB LHARCHA OUED SEHBE EL HARCHA Existant Existant -84300 175200 42
BARRAGE GRARAT RGIBYINE OUED CHEBEIKA En phase d'étude Pas de visibilité -88900 131750 72 119
BARRAGE LABOUIR OUED EL MARMOUTHA En phase d'étude 2017 -240928 75812 956 4.1 1147 1113 3.0%
BARRAGE SAQUIA ALHAMRA OUED SAQUIA ALHAMRA Existant Existant -270775 44383 51 335 400
Nb :
1- La surface complémentaire Barrage Guelta Zergua – Mansour Eddahbi est de l’ordre de 77 686km². Ce dernier
contrôle toute la partie amont Haut Draa et contribue à minimiser le risque de la génération de crues
extraordinaires telle que celle historique de1948 estimée à 6500m3/s !
2- La reprise du calcul de laminage de la crue centennale du Barrage Tighmart (1565m3/s) abouti à un débit de
l’ordre de 1563m3/s. Ce résultat est attendu, vu que le volume de la retenue du barrage (546 198 m3) est très faible
par rapport au volume de la crue estimé à 105 Mm3 !
Les effets de laminage sont très réduits. Seul le barrage Fask peut être pris en
compte ;
Les barrages dont la date du début de réalisation n’est pas fixée ne seront pas pris
en compte.
Jemaa N Tirhirte
Ifrane Anti Atlas
TANTAN DMN
Bouizakarne
AinErrahma
Oued Noun
Talmazirte
Pont Draa
Guelmim
Tifermite
Ouijjane
Aoreora
Sidi Ifni
Taghjijt
Assaka
Pjmax
Pjmax
Unifié
Unifié
Unifié
Assa
Date
Jemaa N Tirhirte
Ifrane Anti Atlas
TANTAN DMN
Bouizakarne
AinErrahma
Oued Noun
Talmazirte
Pont Draa
Guelmim
Tifermite
Ouijjane
Aoreora
Sidi Ifni
Taghjijt
Assaka
Pjmax
Pjmax
Unifié
Unifié
Unifié
Assa
Date
Les séries ci-dessus ont fait l’objet d’un ajustement statistique aux différentes lois usuelles.
Les résultats sont consignés dans les tableaux et les graphiques qui suivent :
Mo Gradex P0
IRE POSTE X Y Z N S Min Max P2 P5 P10 P20 P25 P50 P100
y mm/j mm
1303 Assaka 27700 240600 143 28 22.8 12.7 3.7 60.9 9.9 17.0 21 32 39 47 49 56 63
1096 Aoreora -31100 217600 40 6 20.5 11.1 5.0 38.0 8.7 15.5 19 29 35 42 44 50 56
5856 Oued Noun 20100 243200 115 16 20.5 11.5 8.0 50.0 9.0 15.3 18 29 35 42 44 51 57
247 AinErrahma -400 220850 81 18 16.3 9.8 3.0 40.0 7.6 11.9 15 23 29 35 36 42 47
60033001 Laayoune DMN -273428 45176 63 41 12.9 7.3 0.7 26.4 5.7 9.6 11 18 23 27 28 33 37
60285001 Tan Tan DMN -67515 172665 199 42 16.1 7.6 0.9 29.0 5.9 12.7 15 21 26 30 32 36 41
6260 Pont Draa -43000 183550 25 11 14.3 8.3 6.5 28.2 6.4 10.6 13 20 25 30 31 35 40
4272 Sidi Ifni 36000 273000 59 29 24.9 12.1 8.9 54.8 9.4 19.5 23 34 41 47 49 56 62
5902 Ouijjane 102100 279850 180 30 25.6 10.9 5.2 56.4 8.5 20.7 24 34 41 47 49 54 60
60280001 Guelmim Unifié 45855 229667 300 73 22.9 16.0 3.5 95.0 12.5 15.6 20 34 44 53 56 64 73
1296 Assa 106000 186500 320 20 23.1 13.2 5.2 55.6 10.3 17.1 21 33 41 49 51 58 65
3880 Foum El Hassane 153415 225728 480 21 32.3 20.7 4.0 84.0 16.2 22.9 29 47 59 71 75 86 97
7968 Taghjijt Unifié 109000 236500 572 62 27.3 16.5 0.0 95.0 12.8 19.8 24 38 48 58 62 72 83
2264 Bouizakarne Unifié 81000 250700 646 44 31.6 17.1 7.0 80.0 13.4 23.9 29 44 55 64 67 77 86
7658 Talmazirte 154300 247400 647 10 23.2 18.3 8.0 72.0 14.3 14.9 18 33 45 58 62 76 90
4288 Ifrane Anti Atlas 102110 253701 812 29 34.4 18.4 12.0 76.0 14.4 26.1 31 48 58 69 72 82 92
8704 Tnine des Ait Erka 90000 270000 1050 35 35.4 19.1 7.0 95.0 14.9 26.9 32 50 61 72 75 85 95
4656 Jemaa N Tirhirte 106700 271500 1210 30 36.1 12.0 19.0 60.0 9.4 30.7 34 45 52 59 61 67 74
8208 Tifermite 123600 -288500 1345 22 58.3 23.2 24.0 113.0 18.1 47.8 54 75 89 102 106 118 131
Quantiles des pluies journalières maximales annuelles (mm) pour T=100 ans
Les tableaux suivants récapitulent les données pluviométriques collectées (Pan en mm) :
Nb : la valeur située à l’est du bassin Essaquia Al Hamra a été prise égale à celle de Essmara
4.5.3 Estimation des Pjmax(T), P0, G et Pan pour les bassins versants étudiés.
Etant donné la variabilité dans l’espace des pluies mesurées, le calcul des pluies
de projet devra être fait par le biais d’une intégration :
1
P ( BV ) P( x, y )ds
S
La pluie discrétisée P(x,y) est obtenue par interpolation des valeurs de pluies
affectées au postes pluviométriques.
1 BV-FASK_BGE 112 24.4 13.6 28.9 44.5 55.1 65.6 79.5 90.4
2 BV_AIN-ERRAHMA_STATION 81 14.6 9.3 17.6 27.2 33.9 40.5 49.2 55.7
4 BV_IFNI _STATION 129 21.1 11.1 25.1 37.8 46.6 54.6 64.7 72.2
5 BV-AHRAM 53 12.3 6.5 14.5 21.2 25.3 29.7 35.5 39.8
6 BV-AMMA FATMA 35 9.6 5.6 11.2 16.6 20.3 23.9 28.6 32.2
7 BV-AWDRI 30 9.0 5.6 10.5 15.7 19.2 22.6 27.0 30.3
8 BV-BOUTKEMMART 72 15.5 10.3 18.9 29.9 37.6 45.1 54.9 62.2
9 BV-BEN KHLIL 46 11.1 6.2 13.2 19.5 23.6 27.8 33.3 37.4
10 BV-BOUMALIH 71 13.3 8.1 16.0 24.4 29.9 35.5 42.9 48.3
11 BV_OUED SEYAD 113 21.2 12.8 25.5 43.0 51.0 60.5 72.6 81.6
12 BV-EL OUAAR 33 9.4 5.6 10.9 16.3 20.0 23.5 28.1 31.5
13 BV-NAPPE DE TARFAYA 36 10.2 6.1 12.1 18.8 23.3 27.4 32.6 36.5
14 BV-BOUKILA 80 18.0 11.8 21.8 34.7 43.5 52.2 63.6 72.2
15 BV-CHEBEIKA 33 9.2 5.7 10.6 15.5 18.7 22.0 26.5 29.8
16 BV-DRAA 52 14.4 9.4 17.3 27.0 33.6 40.0 48.4 54.6
17 BV-SEBKHAT AT 'TRAYES 30 8.8 5.7 10.4 16.2 20.0 23.8 28.7 32.3
18 BV-SIDI MED SA'ID 31 9.3 5.8 11.0 17.5 21.8 25.8 31.0 34.9
19 BV-KHAT AL FATER 32 9.4 5.7 11.2 17.8 22.2 26.3 31.6 35.6
20 BV-AS SAQUIYA AL HAMRA 28 7.8 6.0 9.1 13.5 16.5 19.5 23.4 26.3
21 BV-SEBKHAT TAMBAGA 33 9.6 5.7 11.1 16.7 20.6 24.2 28.8 32.3
22 BV-SEHBE EL HARCHA 47 11.4 6.0 13.4 20.0 24.5 28.7 34.2 38.4
23 BV-SERKES 67 14.6 9.4 17.6 27.8 34.1 40.8 49.5 56.0
24 BV BOUSSAFEN+BOUTZEZZARET 69 13.8 8.7 16.7 26.4 31.6 37.7 45.7 51.7
25 BV-SEBAKHT FEDDA 29 8.9 5.7 10.4 15.7 19.3 22.7 27.2 30.5
26 BV-SEBKHAT DAOURA 27 8.4 5.7 9.8 14.7 18.0 21.3 25.6 28.7
27 BV-LAAGUIG 41 10.7 5.8 12.5 18.5 22.7 26.6 31.7 35.6
30 BV-ASSAKA-STATION 102 17.4 11.1 21.2 33.4 42.4 50.5 61.0 68.8
31 BV-OUM LAACHAR 119 21.8 12.75 26.1 41.3 51.9 61.4 73.5 82.4
32 BV-TAGHJIJIT-STATION 110 25 14 29.7 45.5 56.5 67.4 82.2 93.8
33 BV-HOUDA 99.8 18.7 12.05 22.9 36.8 46.3 55.2 66.6 75.0
Les stations ci-dessus sont les seuls points de mesures et de référence qui existent à notre
connaissance. On note un manque significatif d’équipement en stations hydrométriques
de jaugeages des bassins versants, notamment pour ceux qui drainent plus de 1000 km²
et en particulier celui de Saquiat Al Hamra (51 335 Km²) dont la surface total est voisine
de celle de la Moulouya (55 500km2) !. Ce manque est dû principalement aux difficultés
suivantes :
Crues violentes et rapides qui rendent très difficile les opérations de jaugeages ;
Déterrages en permanence des sections des oueds ;
Accès très difficile aux points de mesures ;
Difficulté de faire installer des équipes de jaugeages à proximité des points de
mesures.
Malgré le régime aride qui caractérise la zone du Sud du Maroc, celle-ci posséde un
réseau hydrographique développé qui peut véhiculer des crues violentes générées par
des épisodes pluvieux intenses survenus sur les sous-bassins amont.
Une attention particulière devra être alors, donnée à l’estimation des débits de projet.
L’utilisation directe des formules empiriques qui peuvent être calées dans des contextes
différents, dans le temps et dans l’espace, incite à une grande prudence.
La présente étude tentera, en se basant sur les données hydrologiques disponibles et les
études antérieures effectuées minutieusement dans la zone, de monter une approche
d’estimation des débits au droit des franchissements étudiés, la plus rationnelle possible.
Le bassin de l’oued Assaka draine le bassin de Guelmim, il s’étend sur une superficie de
6896 km² à son embouchure. Le réseau hydrographique y est constitué des 4 principaux
oueds suivants :
Ces cours d’eau sont non pérennes et la plupart du temps à sec. Mais ils peuvent
générer des crues violentes.
IL en découle que les bassins de l’unité de Guelmim présente une bonne aptitude au
ruissellement directe et que seul l’effet d’épandage des crues peut réduire sensiblement
les pointes des crues.
C’est une station qui contrôle l’amont du bassin versant Guelmim, soit environ 1404 km².
On dispose d’un échantillon de 31 valeurs de débits de pointe sur la période 1984-2014.
Le maximum observé est de l’ordre de 1100m3/s, estimé par enquête de crue le 7 janvier
1985. L’évènement de novembre 2014 était de moindre importance avec une pointe
estimée à 602 m3/s.
Les quantiles des débits de pointes par ajustement statistiques peuvent être estimés par
la deuxième droite de tendance soit :
T 10 20 50 100
u Gumbel 2.25 2.97 3.90 4.60
Qp 656 862 1129 1329
Kp 3.44 3.69 3.93 4.07
La méthode consiste à supposer qu’à partir d’un seuil, l’accroissement des volumes
ruisselés (Qj ou R) est proportionnel à celui des pluies abattues. Cette méthode permet
de pallier au manque de données hydrométriques et permet de mette en relief l’effet de
saturation des sols.
L’analyse pluviométrique précédente a permis d’obtenir les données pluviométriques
suivantes :
G P0
ARF T*
S Cp
Qj0
S T0
m m /j mm Km² m 3/s
14 25 0.69 20 1404 3.60 -8.948 76.931 1.82
Le tableau et la figure suivante résume les résultats obtenus en les comparant au valeurs
des débits issus des formules empiriques couramment utilisées :
K a a
2 2 9.88
a N b
20 85 0.742
T Ajustement Gradex Mallet Gauthier FullerII Régional
2 69 780 277
5 383 251 1167 552
10 656 591 424 1460 759
20 862 790 544 1753 967
50 1129 1314 671 2141 1241
C’est une station qui contrôle l’aval du bassin versant Guelmim, soit environ 6557 km².
On dispose d’un échantillon de 32 valeurs de débits de pointe sur la période 1982-2014.
Le maximum observé est de l’ordre de 2200m3/s, estimé par enquête de crue le 7 janvier
1985. L’évènement de novembre 2014 était de moindre importance avec une pointe
estimée à 1796 m3/s.
Les quantiles des débits de pointes par ajustement statistiques peuvent être estimés par
la deuxième droite de tendance soit :
T 10 20 50 100
u Gumbel 2.25 2.97 3.90 4.60
Qp 1099 1586 2216 2688
Kp 2.93 3.31 3.65 3.86
La méthode consiste à supposer qu’à partir d’un seuil, l’accroissement des volumes
ruisselés (Qj ou R) est proportionnel à celui des pluies abattues. Cette méthode permet
de pallier au manque de données hydrométriques et permet de mette en relief l’effet de
saturation des sols.
L’analyse pluviométrique précédente a permis d’obtenir les données pluviométriques
suivantes :
30 BV-ASSAKA-STATION 102 17.4 11.1 21.2 33.4 42.4 50.5 61.0 68.8
G P0
ARF T*
S Cp
Qj0
S T0
PAN
mm/j mm Km² m3/s mm
11.1 17.4 0.56 20 6557 1.90 -153.080 309.79 4.67 102
Qjmax(T)= s x u(T) + Qj0 T>T0
T u Pjm ax(T) Pjmax(A,T) ROJ(T) Qjm ax(T) Qp(T) Kp(T) CN(T) Qgradex
2 0.4 21 12 0.8 62 118 0.6 88.6
5 1.5 34 19 4.1 312 592 2.3 89.5
10 2.3 42 24 7.2 544 1034 2.9 90.0 1034
20 3.0 50 28 10.1 767 1457 3.2 89.9 1457
50 3.9 61 34 15.9 1210 2299 3.7 91.2 2299
Le tableau et la figure suivante résume les résultats obtenus en les comparant aux valeurs
des débits issus des formules empiriques couramment utilisées :
K a a
2 2 9.88
a N b
20 85 0.742
T Ajustement Gradex Mallet Gauthier FullerII Régional
2 118 2445 870
5 592 3660 1731
10 1099 1034 1007 4579 2383
20 1586 1457 1431 5498 3034
50 2216 2299 1847 6712 3895
Etant donnée, la position de la voie expresse, pratiquement à l’aval immédiat des reliefs,
la station Taghjijit sera prise comme station de référence pour l’estimation des débits de
pointes des bassins versants importants, interceptés par la voie de Gulemim au pont
Draa.
C’est une station qui contrôle l’oued Ifni. De part, sa position, elle n’a pas beaucoup
d’intérêt pour la présente étude. Néanmoins nous tenons à examiner ses quantiles, car
c’est la seule station qui jauge un bassin côtier atlantique et dont on dispose des mesures
hydrologiques. Les résultats peuvent être intégrer dans l’analyse régionale qu’on tentera
d’établir par la suite. La station draine un bassin versant de l’ordre de 685 km².
On dispose d’un échantillon de 36 valeurs de débits de pointe sur la période 1979-2014.
Le maximum observé est de l’ordre de 890m3/s (Kp ~ 4.1).
Les quantiles des débits de pointes par ajustement statistiques peuvent être estimés par
la deuxième droite de tendance soit :
T 10 20 50 100
u Gumbel 2.25 2.97 3.90 4.60
Qp 131 396 740 997
Kp 2.48 3.41 3.94 4.19
La méthode consiste à supposer qu’à partir d’un seuil, l’accroissement des volumes
ruisselés (Qj ou R) est proportionnel à celui des pluies abattues. Cette méthode permet
de pallier au manque de données hydrométriques et permet de mette en relief l’effet de
saturation des sols.
L’analyse pluviométrique précédente a permis d’obtenir les données pluviométriques
suivantes :
PAN P0 Gradex P2 P5 P10 P20 P50 P100
N0 Bassins versant
mm mm mm/j mm mm mm mm mm mm
4 BV-Sidi Ifni 129 21.1 11.1 25.1 37.8 46.6 54.6 64.7 72.2
G P0 S Qj0
ARF T* Cp S T0
mm/j mm Km² m3/s
La figure ci-dessus montre que l’évènement de 2014 été exceptionnelle. Il est fort
probable que son occurrence est supérieure à celle que lui confère l’échantillon transmis
par la ABHSMD.
2 23 464 163 4
Le bassin versant de Draa s'étend des montagnes du Haut Atlas Central vers le sud aux
basses collines présahariennes. Il draine une surface totale de l’ordre de 93 000 km².
Compte tenu de ses caractéristiques topographiques, géomorphologiques et
hydrographiques, le bassin de Draa est subdivisé en trois zones :
15 120 km2 pour le Haut Draa ;
14 380 km2 pour la vallée Moyenne du Draa
63 500 km2 pour le sous-bassin du bas Draa.
Ils se répartissent sur les provinces d’Ouarzazate, Zagora, Tata, Assa-Zag et Tan-Tan.
L’oued Draa prend sa source depuis le Haut Atlas aux environs de la ville d’Ouarzazate
jusqu’à son embouchure sur l’océan atlantique à 27 km au Nord-ouest de la ville de Tan-
Tan. Le bassin versant est aménagé par un grand barrage « Mansour Eddahbi » qui
contrôle le haut Draa soit environ 15 000 km². Le bassin versant intermédiaire reste très
important. Il cumul une surface de l’ordre de 77 686 Km² au droit du franchissement
étudié soit la station hydrométrique Pont Draa. Les mesures hydrométriques de cette
station reflètent le comportement du bassin versant de l’Oued Draa à l’état actuel. Elles
seront exploitées pour aboutir, par ajustement statistiques, aux débits de pointes.
En amont immédiat de la voie Express, un site de barrage nommé Gueltat Zergua a été
étudié [10] de façon très détaillée. L’étude a procédé par une modélisation Pluie-débit
sur la base de coefficients d’indexation du sol CN. Cette approche a l’avantage de
mettre en relief l’effet seuil du BV. Ce dernier peut être masqué lors de l’ajustement
statistique à cause d’une taille insuffisante des débits.
L’étude du barrage Guelta Zergua fait ressortir les débits de pointes suivants :
C’est une station qui contrôle pratiquement tout le bassin versant de l’oued Draa et
depuis 1972, elle reflète le régime d’écoulement du bassin intermédiaire (77686 km²).
On dispose d’un échantillon de 31 valeurs de débits de pointe sur la période 1984-2014.
Le maximum observé est de l’ordre de 1100m3/s, estimé par enquête de crue le 7 janvier
1985. L’évènement de novembre 2014 était de moindre importance avec une pointe
estimée à 602 m3/s.
T 10 20 50 100
u Gumbel 2.25 2.97 3.90 4.60
Qp 657 1133 1750 2212
Kp -0.23 0.53 1.13 1.46
C’est une station qui contrôle un « petit afluent » de l’oued Draa qui débouche
directement sur la ville de Tan Tan. Le bassin versant drainé est de l’ordre de 890 km².
On dispose d’un échantillon de 24 valeurs de débits de pointe sur la période 1984-2007.
Le maximum observé est de l’ordre de 1200m3/s survenu le 22/10/2003.
Les quantiles des débits de pointes par ajustement statistiques peuvent être estimés par
la deuxième droite de tendance soit :
T 10 20 50 100
u Gumbel 2.25 2.97 3.90 4.60
Qp 612 820 1089 1291
Kp 3.64 3.89 4.13 4.28
La méthode consiste à supposer qu’à partir d’un seuil, l’accroissement des volumes
ruisselés (Qj ou R) est proportionnel à celui des pluies abattues. Cette méthode permet
de pallier au manque de données hydrométriques et permet de mette en relief l’effet de
saturation des sols.
L’analyse pluviométrique précédente a permis d’obtenir les données pluviométriques
suivantes :
PAN P0 Gradex P2 P5 P10 P20 P50 P100
N0 Bassins versant
mm mm mm/j mm mm mm mm mm mm
9 BV-BEN KHLIL 46 11.1 6.2 13.2 19.5 23.6 27.8 33.3 37.4
G P0
ARF T*
S Cp
Qj0
S T0
m m /j mm Km² m 3/s
6.2 11.1 0.73 10 890 4.30 -340.610 211.63 7.6
Le tableau et la figure suivante résume les résultats obtenus en les comparant aux valeurs
des débits issus des formules empiriques couramment utilisées :
K a a
2 2 9.88
a N b
20 85 0.742
T Ajustement Gradex Mallet Gauthier FullerII Régional
2 119 560 198
5 334 131 838 393
10 612 583 204 1049 541
20 820 728 257 1259 689
50 1089 916 314 1537 885
100 1291 1057 351 1748 1033
Les éléments précédents montrent que le maximum observé a fort probablement une
occurrence plus grande que celle que lui confère l’échantillon étudié. Le débit
centennal est estimé à 1060 m3/s. Néanmoins le débit de projet sera pris égal au
maximum observé soit 1290m3/s.
Une relation régionale préliminaire qui permet d’estimer le débit centennal des
bassins versants importants de la zone d’étude :
Qp100 = 983.98 Ln(S) – 5645.7, avec 685km² <S en km² < 6557km²
Pour le cas étudié, ce modèle est formé principalement par les trois éléments essentiels :
Modélisation pluviométrique :
Cette partie consiste à l’élaboration d’un modèle météorologique pour l’estimation de
la pluie moyenne du bassin, temporaire et spatiale. Cette partie a été traitée dans 3.5.3.
Les Pjmax de chaque bassin versant ont été estimés. Elles seront distribuées selon l’averse
défavorable de type II et multipliées par un coefficient d’abattement spatial (ARF). Ce
coefficient varie en fonction de la surface. Il reflète la non concomitance des Pjmax
d’occurrence T sur tout le bassin versant. Plus bassin est petit ARF tend vers 1.
Les courbes « A, B, C » donnant l’abattement spatial ARF dans Design of Small Dams
concernent une pluie frontale dont le corps à une durée de 6 heures mais dont la durée
totale peut atteindre 12 heures, voire 24 heures et plus. La courbe A correspond à la
cordillère côtière et aux Rocheuses à proximité immédiate du Pacifique. Au Maroc une
telle zone correspondrait éventuellement au Haut Atlas et au Rif dans les versants les plus
concernés par l’effet orographique.
La courbe B correspond à des reliefs moins hauts et à une topographie moins chahutée.
Enfin, la zone C correspond à la transition entre Rocheuses, plateaux intérieurs et
Grandes Plaines de l’Ouest des USA. On voit que le Moyen Atlas et les contreforts des
reliefs participent de la famille B. Enfin, les zones les moins orographiques du Maroc
seraient de la famille C [3].
(P I a )2
Pe
P Ia S
Avec :
o P : Précipitation cumulés au temps t ;
o Ia : Perte initiale ;
o S : Potentiel maximum de rétention (Ia = 0.2S)
( P 0.2 S ) 2
D’où : Pe
P 0. 8 S
Le potentiel "S" et les caractéristiques du bassin versant sont liés par le « Curve Number »
ou le coefficient d’indexation du sol "CN" :
25400 254CN
S (mm)
CN
Le coefficient d’indexation du sol "CN" renseigne sur la capacité du sol à infiltrer l’eau de
pluie. Il existe des tables dans la littérature qui fournissent le "CN" en fonction de
l’occupation du sol et ses états de saturation (sec : (CN(I), humide : CN(II) ou saturé :
CN(III)).
1
Ouvrage Appleid Hydrologie, Ven Te Chow, David R Maidment, Lary W Mays, McGraw-Hill, Inc. 1988
2
Largement utilisé à l’USA, en France et au Maroc il a été adopté par l’étude du plan national de protection contre les
inondations PNI
Pour la présente étude, on retiendra les conclusions de l’étude de l’oued Draa [10] :
Le tableau suivant résume les données de base utilisée par la suite dans la modélisation
pluie débit. Deux états seront examinés :
Etat d’antécédent d’humidité moyen CN(II). Il correspond au CN issus de la
classification précédente et des interceptions initiales qui en découlent (Ia ~ 0.2
S) ;
Etat d’antécédent humide saturé CN(III). Il correspond au CN estimé en fonction du
CN(II) et des Ia qui en découlent.
L
BV S Km² Ia 1 CNI Ia 2 CNII Ia 3 CNIII P100 ARF ARF P100 P% tc Lag
Km
BV1_Oum_Lachar 836.7 21.51 70.25 9.10 84.90 3.93 92.8 82.4 0.75 61.5 60.1 0.9 419.2 251.5
BV0_Houda 43.6 27.27 65.07 11.50 81.60 4.98 91.1 75 0.98 73.5 24.2 1.1 200.8 120.5
BGE FASK 2400.0 15.10 77.08 6.40 88.90 2.76 94.9 90.4 0.65 59.1 112.4 0.7 700.1 420.1
Taghjijit 1404.0 15.56 76.55 6.50 88.60 2.84 94.7 93.8 0.69 64.3 69.9 1.1 449.5 269.7
Sidi Ifni 685.0 11.82 81.13 4.90 91.10 2.16 95.9 72.2 0.78 56.4 23.6 1.1 197.0 118.2
BV2-SEYAD 3708.6 17.12 74.79 7.20 87.60 3.13 94.2 81.6 0.65 53.0 145.2 0.6 883.1 529.8
BV3_Boukila 2158.8 20.18 71.57 8.50 85.70 3.69 93.2 72.2 0.66 47.8 122.9 0.5 813.3 488.0
BV4_Boutkemmart 306.0 16.96 74.97 7.10 87.70 3.10 94.3 62.2 0.89 55.1 34.6 1.1 265.0 159.0
BV5_SER_BOU_BO 580.0 16.34 75.67 6.90 88.10 2.98 94.5 53.85 0.81 43.4 49.6 0.8 366.7 220.0
BV_Ain_Errahma 1635.0 19.36 72.40 8.10 86.20 3.54 93.5 55.7 0.68 37.7 69.5 0.8 480.5 288.3
BV6_BOUMALIH 359.6 16.81 75.14 7.00 87.80 3.07 94.3 48.3 0.90 43.3 29.2 1.6 215.5 129.3
BV7_DRA 77686 26.01 66.13 11.00 82.30 4.75 91.4 54.6 0.65 35.5 600.0 0.07 3885 2331
BV8_AHRAM 80.5 21.34 70.41 8.90 85.00 3.90 92.9 39.8 0.97 38.7 15.7 0.7 155.7 93.4
BV9_BEN KHLIL 890.0 21.51 70.25 9.00 84.90 3.93 92.8 37.4 0.83 31.1 48.6 0.8 362.1 217.2
BV10_SEHB_HARCHA 515.6 35.32 58.99 14.80 77.40 6.45 88.7 38.4 0.87 33.4 52.8 0.4 442.3 265.4
BV11_CHEBIKA 7283.4 26.73 65.52 11.20 81.90 4.88 91.2 29.8 0.77 23.0 183.5 0.2 1265.6 759.4
BV12_LAAGUIG 80.7 26.01 66.13 10.90 82.30 4.75 91.4 35.6 0.97 34.6 17.5 0.8 166.7 100.0
BV13_Amma_Fatma 796.1 27.27 65.07 11.50 81.60 4.98 91.1 32.2 0.84 27.0 84.6 0.4 636.4 381.9
BV14_EL OUAAR 1215.5 36.13 58.44 15.20 77.00 6.60 88.5 31.5 0.81 25.7 88.6 0.3 681.7 409.0
BV16_Awdri 1414.0 34.91 59.27 14.70 77.60 6.38 88.8 30.3 0.81 24.5 84.2 0.3 650.9 390.6
BV22_SEBKHAT AMSEKIR 1051.0 39.25 56.41 16.50 75.50 7.17 87.6 32.3 0.82 26.6 81.2 0.2 683.7 410.2
BV18_Sebkhat Daoura 5147.6 39.04 56.55 16.40 75.60 7.13 87.7 28.7 0.77 22.2 157.4 0.2 1186.9 712.2
BV20_sebkhat at Trayes 1239.0 37.78 57.35 15.90 76.20 6.90 88.0 32.3 0.81 26.3 61.2 0.2 558.6 335.1
BV19_Si SAID 280.1 37.78 57.35 15.90 76.20 6.90 88.0 34.9 0.92 31.9 37.8 0.1 463.8 278.3
BV21_Khat Al Fater 113.7 39.04 56.55 16.40 75.60 7.13 87.7 35.6 0.96 34.2 10.3 0.1 158.6 95.2
BV_23_SAQUIYA_HAMRA 51335 33.92 59.96 14.30 78.10 6.19 89.1 26.3 0.77 20.33 540.0 324.0
BV_ASSAKA 6557.0 24.60 67.38 10.30 83.10 4.49 91.9 68.8 0.65 44.7 184.1 0.5 1072.8 643.7
ST_BEN KHLIL 890.0 21.51 70.25 9.10 84.90 3.93 92.8 37.4 0.83 31.1 48.6 0.8 362.1 217.2
Modèle météorologique de type SCS Storm : Pjmax 100 selon une averse de type II ;
Modèle de perte de type SCS : CN, et Ia aux états d’antécédent d’humidité moyen
et saturé.
Modèle de transfert de type SCS : Lag ~ 0.6*tc.
Les calculs hydrologiques menés précédemment sont récapitulés dans le tableau et les
graphiques suivants :
Q100 Q100
Q10
Q100 Q100 SCS SCS Q100
NBV Bassins versant S Km² ARF P100 P% 0 Max Obs Historique
Régional Gradex Moye Humid Rationnelle
Stat
n e
BV_AIN-ERRAHMA_STATION 1635 38 0.8 1635 1648 601 1159
BV_IFNI _STATION 685 56 1.1 779 1527 1290 1475 1897 890
BV9 BV-BEN KHLIL_STATION 890 31 0.8 1037 1346 1290 220 541 1200
La figure ci-dessus récapitule les résultats obtenus. Elle affiche les éléments suivants :
Les débits de pointe, de la région étudiée, présentent une variation importante qui
dépendrait de :
Pour chaque zone deux états peuvent être considérés : Moyen (CNII) : M et
Saturé (CNIII) : H
Les résultats d’une synthèse régionale basée uniquement sur les stations
hydrométriques disponibles ne peuvent pas être transposés à toute la zone
d’étude. En effet, la variabilité des caractéristiques physiques et pluviométriques
des bassins versants est très importante. Néanmoins et à l’issu de la présente
étude, nous pouvons proposer les courbes enveloppes suivantes :
300 km² < S < 10000 Km² Q = 378.5 Ln(S) - 1648 Q = 860 Ln(S) - 3900
BV2 BV_OUED SEYAD & OUM AACHAR 3709 53 0.6 2441 2583 1836 2763 3168 4.4
BV3 BV-BOUKILA 2159 48 0.5 1909 1991 827 1402 2702 4.5 Q 2014 ~ 1700 m3/s [7]
BV12 BV-LAAGUIG 81 35 0.8 32.5 91.8 126 81 3.3 2 EXUTOIRES, Fonctionnement Simultanés ?
BV13 BV-AMMA FATMA 796 27 0.4 927 1296 56.4 201 224 2.8
BV14 BV-EL OUAAR 1215 26 0.3 1343 1370 30.2 190 308 2.9
BV20 BV-SEBKHAT AT 'TRAYES 1239 26 0.2 1362 1382 31.3 221 313 2.9 Dépressions en amont DU BV19
BV22 BV-SEBKHAT AMSEKIR 1051 27 0.2 1200 1344 22.5 158 275 2.8 Dépressions en amont DU BV19
Le débit sera vérifié en fonction des BV20 et
BV19 BV-SIDI MED SA'ID 280 32 0.1 18.3 89 121 2.9
BV22
BV21 BV-KHAT AL FATER 114 34 0.1 15.5 89 164 88 3.2
Le pont devra être dimensionné pour la
BV23 BV-AS SAQUIYA AL HAMRA 51335 20 5027 520 5797 6000 6000 3.2 capacité du lit mineur qui devra être < au débit
laminé par le Bge projeté