M HDR - Urb 2018 09

République Algérienne Démocratique et Populaire
Ministère de l’Enseignement Supérieur et de la Recherche Scientifique

Université Mohamed Seddik Benyahia Jijel
Faculté des sciences et de la Technologie
Département de Génie Civil et Hydraulique
MÉMOIRE DE FIN D’ÉTUDES

En vue de l’obtention du diplôme de Master en Hydraulique
Option : Hydraulique Urbaine
Présenté par :
DOUAS Fatima-Zohra
KHEDIMALLAH Mariama
Thème
ETUDE DE LA PROBLEMATIQUE DU RISQUE

D'INONDATION EN MILIEU RURAL
Dirigé par :
- Mr. ABDI Ishak
Soutenu publiquement le 28/ 06 / 2018 devant le jury composé de :

- Président : RATIAT. A
- Examinateur :KRIMIL .F
Promotion 2018
Les problèmes d'estimation des crues constituent un souci majeur aux techniciens chargés de
dimensionner les ouvrages hydrotechniques. Dans ce travail, nous avons quantifié le débit de crue
centennale par différentes approches de modélisation hydrologique dans le bassin versant de l'Oued
Djendjen situé à Texenna au Nord-Est de la wilaya de Jijel. Dans ce cadre nous avons appliqué des
méthodes statistiques, des méthodes empiriques et des méthodes basés sur la relation pluie-débit.
Dans le cas des méthodes statistiques, l’ajustement par la méthode de Gumbel à base de données
hydrométrique et la méthode du Gradex à base de données hydro-pluviométrique, ont été
appliquées. Pour la modélisation pluie-débit nous avons utilisé deux modèles. Le premier est le
modèle SCS-CN retenu pour la fonction de production et le modèle NASH pour la fonction de
transfert.
L'étude hydrologique, par l'application de ces méthodes d'estimation des crues, a révélé des
différences dans les résultats obtenus. Ceci est expliqué par le manque de fiabilité de certaines de ces
méthodes.
Enfin, nous avons effectué une étude hydraulique afin de localiser les zones susceptibles d'être
inondés. Nous avons choisi de mener cette étude d'inondabilité avec les logiciels: ArcGis et HEC-RAS.
Mot clés
Débit, Gumbel, Gradex, SCS-CN, hydrogramme unitaire, HEC-RAS,ARC-GIS, crue, inondation
i
Flood estimation problems are a major concern in the hydro technical structures design field. In this
work, we quantified the 100-year flood discharge by different hydrological modeling approaches
applied to the Oued Djendjen watershed located in Texenna in the north-east of the Jijel
department. In the scope of this work we have applied statistical methods, empirical methods and
methods based on the rainfall-runoff relationship. In the case of statistical methods, the
hydrometric-based Gumbel adjustment and the rainfall-based Gradex method were applied. For the
rainfall-runoff modeling we used two models. The first is the SCS-CN model selected for the
production function and the NASH model for the transfer function.
The hydrological study, by the application of these methods of estimating floods, revealed
significant variations in the obtained results. This may be explained by the unreliability of some of
these methods.
Finally, we carried out a hydraulic study to locate the areas likely to be flooded. We conducted this
flood risk study using the softwares: ArcGis and HEC-RAS.
Keywords
Flow, Gumbel, Gradex, SCS-CN, unit hydrograph, HEC-RAS, ARC-GIS, flood
ii
‫مـلـخـص‬
‫يعتبر تقدير مخاطر الفيضانات مصدر قلق في مجال تصميم المنشئات الهيدروتقنية‪ .‬في هذا العمل ‪،‬‬
‫قمنا بتقدير كمية الفيضانات لمدة ‪ 100‬عام من خالل طرق نمذجة هيدرولوجية مختلفة و التي‬
‫قمنا بتطبيقها على الحوض المتدفق لواد جن جن المتواجد في تاكسنة في الشمال الشرقي من والية‬
‫جيجل‪ .‬في إطار هذا العمل قمنا بتطبيق طرق إحصائية‪ ،‬طرق تجريبية و طرق مبنية على أساس‬
‫‪ Gumbel‬المبني على‬ ‫العالقة "تساقط ‪ -‬تدفق"‪ .‬في حالة األساليب اإلحصائية ‪ ،‬تم تطبيق معادلة‬
‫القياس الهيدرولوجي وطريقة ‪ Gradex‬المبنية على قياس التساقط‪ .‬بالنسبة للنمذجة "تساقط ‪-‬‬
‫تدفق" ا ستخدمنا نموذجين‪ .‬األول هو نموذج ‪ SCS-CN‬المحدد لوظيفة اإلنتاج ونموذج ‪NASH‬‬
‫لوظيفة النقل‪.‬‬
‫كشفت الدراسة الهيدرولوجية المتبعة‪ ،‬عن طريق تطبيق هذه الطرق لتقديرمخاطر الفيضانات‪ ،‬فروقا‬
‫في النتائج التي تم التوصل إليها‪ .‬يمكننا تفسير ذلك بعدم نجاعة بعض هذه الطرق ‪.‬‬
‫أخيرًا ‪ ،‬أجرينا دراسة هيدروليكية لتحديد المناطق التي من المرجح أن تغمرها المياه‪ .‬تم إجراء‬
‫دراسة مخاطر الفيضانات هذه باستخدام برنامجي ‪ ArcGis :‬و ‪.HEC-RAS‬‬
‫كلمات داللية‬
‫تدفق‪ ،SCS-CN ،Gradex ،Gumbel ،‬هيدروغرام وحدوي‪ ،ArcGis HEC-RAS ،‬فيضان‬
‫‪iii‬‬
Tout travail réussi dans la vie nécessite en premier lieu la bénédiction de dieu, et en suite l’aide et le
support de plusieurs personnes.
Nos remerciements s'adressent à notre encadrant, Monsieur ABDI Ishak, pour avoir accepté de
diriger ce travail. Sa rigueur scientifique, son soutien et ses conseils nous ont été d’une aide
inestimable.
Nous tenons aussi à remercier les membres du jury pour l’honneur qu’ils nous font en acceptant
d'évaluer ce travail.
Nous n'oubliant pas de remercier le staff qui nous encadrer lors de notre stage à l'URBACO.
Nous remercions nos collègues de la première promotion de Master en hydraulique à l’université de

Jijel ainsi que tous les enseignants et les membres du staff administrative et en particulier ceux du
département de « Génie Civil et Hydraulique »
Enfin, Nous tenons à exprimer notre reconnaissance à toutes les personnes qui ont contribué, depuis
notre enfonce, dans notre éducation morale et scientifique.
iv
A mes parents qui me sont très chers.
A tous mes frères et sœurs.
A toute ma famille.
A ceux qui m’aime et m’ont aidé sans qu’ils le sachent.
Aux chères amies Nadjat, Habiba, Ilham, Meriame qui m’ont toujours soutenu.
Fatima Zohra
v
À ma mère,
« Tu m’as donné la vie, la tendresse et le courage pour réussir.
Tout ce que je peux t’offrir ne pourra exprimer l’amour et la reconnaissance que je te porte.
En témoignage, je t’offre ce modeste travail pour te remercier pour l’affection dont tu m’as toujours entourée »
À mon père,
« L’épaule solide, l’œil attentif compréhensif et la personne la plus digne de mon estime et de mon respecte
Aucun dédicace ne saurait exprimer sentiments, que dieu te préserve et te procure santé et longue vie »
À ma chère amie Amel
« Notre passé se ressemble,
Notre présent nous assemble
Notre futur c'est ensemble. »
À mes frère Abd-El Basset et Abd-Allah.
À mes sœurs
À mes familles.
À tout mes amis.
Mariama.
vi
RESUME ................................................................................................................................................................. I
ABSTRACT ............................................................................................................................................................. II
‫ مـلـخـص‬................................................................................................................................................................ III
REMERCIEMENTS ................................................................................................................................................. IV
DEDICACES ............................................................................................................................................................ V
DEDICACES ........................................................................................................................................................... VI
TABLE DES MATIERES .......................................................................................................................................... VII
TABLE DES FIGURES ............................................................................................................................................... X
TABLE DES TABLEAUX......................................................................................................................................... XIII
INTRODUCTION GENERALE ................................................................................................................................. XV
I. NUMERO DU CHAPITRE ................................................................................................................................ 1
CHAPITRE I ............................................................................................................................................................ 1
PROBLEMATIQUE &GENERALITES SUR LES INONDATIONS .................................................................................... 1
I.1. CYCLE HYDROLOGIQUE ....................................................................................................................................... 1

I.1.1. Phases principales ................................................................................................................................... 2
I.2. LES PARAMETRES FONDAMENTAUX DU COURS D'EAU ............................................................................................... 5
I.2.1. Lit de la rivière ......................................................................................................................................... 5
I.2.2. Rive / berge ............................................................................................................................................. 6
I.2.3. Ripisylve ................................................................................................................................................... 7
I.2.4. Alluvions et substratum ........................................................................................................................... 7
I.3. LES EVENEMENTS-EXTREMES .............................................................................................................................. 7
I.3.1. Etiage et tarissement .............................................................................................................................. 7
I.3.2. Crue ......................................................................................................................................................... 8
I.3.3. Type de crue .......................................................................................................................................... 14
I.3.4. Inondation ............................................................................................................................................. 15
I.4. DIFFERENCE ENTRE INONDATION ET CRUE............................................................................................................ 19
I.5. RISQUE DE L'INONDATION ................................................................................................................................ 19
I.5.1. Risques majeurs-inondations ................................................................................................................ 20
I.5.2. Dégât éventuel ...................................................................................................................................... 20
I.5.3. Dans le monde ....................................................................................................................................... 21
I.5.4. Contexte l'Algérien ................................................................................................................................ 23
I.6. CONCLUSION ............................................................................................................................................. 25
II. NUMERO DU CHAPITRE ...............................................................................................................................27
CHAPITRE II ..........................................................................................................................................................27
II. ETAT DE L'ART SUR LA MODELISATION DES INONDATION ...........................................................................27
II.1. METHODE STATISTIQUE ................................................................................................................................... 27

II.1.1. L'analyse fréquentielle ...................................................................................................................... 27
II.1.1.3. Les courbes " débit- durée-fréquence ............................................................................................... 31
II.2. LA METHODE DU GRADEX ................................................................................................................................. 31
vii
II.3. LA METHODE RATIONNELLE ............................................................................................................................... 33
II.3.1. Principe de la méthode ..................................................................................................................... 33
II.4. LES METHODES EMPIRIQUES.............................................................................................................................. 37
II.5. TRANSFORMATION PLUIE-DEBIT ........................................................................................................................ 41
II.5.1. Fonction de production ..................................................................................................................... 41
II.5.2. Méthodes à base physique :.............................................................................................................. 48
II.5.3. Fonction de transfert ........................................................................................................................ 52
II.5.4. L'Hydrogramme Unitaire Synthétique Régional de SCS .................................................................... 59
II.6. CONCLUSION ................................................................................................................................................. 61
III. NUMERO DU CHAPITRE ...............................................................................................................................62
CHAPITRE III .........................................................................................................................................................62
III. ZONE D'ETUDE ET CARACTERISTIQUE MORPHOMETRIQUES ........................................................................62
III.1. PRESENTATION DE LA ZONE D'ETUDE .................................................................................................................. 62

III.1.1. Présentation de la Wilaya ................................................................................................................. 62
III.1.2. Présentation du bassin versant ......................................................................................................... 66
III.2. ETUDE DES CARACTERISTIQUES MORPHOLOGIQUE DU BASSIN VERSANT DE L'OUED DJENJEN ........................................ 70
III.2.1. Caractéristiques Géométriques ........................................................................................................ 70
III.2.2. Caractéristiques de forme ................................................................................................................. 71
III.2.3. Caractéristiques de relief .................................................................................................................. 72
III.3. CONCLUSION ................................................................................................................................................. 81
IV. NUMERO DU CHAPITRE ...........................................................................................................................82
CHAPITRE IV .........................................................................................................................................................82
IV. ESTIMATION DU DEBIT DE POINTE...........................................................................................................82
IV.1. ANALYSE ET TRAITEMENT DES DONNEES .............................................................................................................. 82

IV.2. CALCUL DE DEBIT DE POINTE AU DROIT DE LA STATION MISSA .................................................................................. 83
IV.2.1. Ajustement statistique ...................................................................................................................... 84
IV.2.2. Teste d'adéquation par la loi de Gumbel .......................................................................................... 86
IV.2.3. La méthode Gradex ........................................................................................................................... 89
IV.2.4. Calcul des quantiles pour les différentes périodes de retour ............................................................ 91
IV.2.5. Application des méthodes empiriques .............................................................................................. 94
IV.2.6. Modélisation pluie-débit ................................................................................................................... 94
IV.3. PLUIE DE PROJET........................................................................................................................................... 100
IV.3.1. Hyétogramme des pluie synthétiques ............................................................................................. 100
IV.3.2. Méthode de la pluie de projet double triangle symétrique ............................................................. 101
IV.4. APPLICATION DE LA METHODE SCS-CN ET HUI DU MODELE DE NASH SUR LES PLUIES DE PROJET ............................. 103
IV.5. INTERPRETATION DES RESULTATS ..................................................................................................................... 105
IV.6. CONCLUSION ............................................................................................................................................... 108
V. NUMERO DU CHAPITRE .............................................................................................................................109
CHAPITRE V ........................................................................................................................................................109
V. ETUDE HYDRAULIQUE ...............................................................................................................................109
V.1. L’OBJECTIF DE L’ETUDE HYDRAULIQUE .............................................................................................................. 109

V.2. PRESENTATION D'OUTILLES D'INFORMATION UTILISEES ......................................................................................... 110
V.2.1. Présentation du modèle HEC-GEORAS ............................................................................................ 110
V.3. ETAPE DE SIMULATION PAR HEC-RAS ............................................................................................................. 113
V.4. PRESENTATION GENERALE .............................................................................................................................. 113
V.5. ELABORATION DU FICHIER DE LA GEOMETRIE ...................................................................................................... 115
viii
V.6. SIMULATION HYDRAULIQUE ............................................................................................................................ 115
V.6.1. Le choix du coefficient de la rugosité (coeff. de Manning Strickler) ................................................ 115
V.6.2. Donné hydrométrique et Les condition au limite ............................................................................ 116
V.6.3. Les profils en long ........................................................................................................................... 117
V.6.4. Les profils en travers ....................................................................................................................... 119
V.6.5. Vitesse d'écoulement ...................................................................................................................... 120
V.6.6. la courbe de tarage ......................................................................................................................... 120
V.6.7. principaux paramètres hydrauliques .............................................................................................. 121
V.7. CONCLUSION ............................................................................................................................................... 122
CONCLUSION GENERALE ....................................................................................................................................123
BIBLIOGRAPHIE ..................................................................................................................................................125
ANNEXE A ..........................................................................................................................................................129
ix
Figure ‎I-1: Conceptualisation du cycle de l’eau(Raghunath,2006) ......................................................................... 2
Figure ‎I-2: Cycle de l'eau(Laborde,2000) ................................................................................................................ 2
Figure ‎I-3: Lit mineur d'un cours d'eau (Bachi, 2011)............................................................................................. 5
Figure ‎I-4: Lit majeur d’un cours d’eau (Bachi, 2011)............................................................................................. 6
Figure ‎I-5: Lit moyen du cours d'eau(Bachi, 2011) ................................................................................................. 6
Figure ‎I-6: Représentation de termes géographique usuels (Degoutte, 2012) ...................................................... 7
Figure ‎I-7: Exemple de courbe de tarissement (Musy, 1998)................................................................................. 8
Figure ‎I-8:Hydrogramme de crue. (Bennis, 2009) .................................................................................................. 9
Figure ‎I-9: Influence de la forme de bassin versant sur l'hydrogramme de crue (musy et higy, 2003) ............... 12
Figure ‎I-10: Inondation de Plaine. (Merabet ,2006) ............................................................................................. 15
Figure ‎I-11: Inondation par remontées des nappes phréatiques. (Merabet, 2006) ............................................ 16
Figure ‎I-12: Inondation par crues torrentielles. (Merabat, 2006) ........................................................................ 17
Figure ‎I-13: Inondation par ruissellement en secteur urbain. (Merabat, 2006) .................................................. 18
Figure ‎I-14: Types des inondations ....................................................................................................................... 19
Figure ‎I-15 : En chine, Huang He, le fleuve Jaune (Photo: Wikimédia Commons/Leruswing). ............................ 22
Figure ‎I-16: Les inondations aux abords de Houston. (David J. Phillip - figaro.fr) ................................................ 22
Figure ‎I-17: cas significatifs des inondations survenues en Algérie. (Merabet, 2006) ......................................... 25
Figure ‎I-18: Inondations de bab el Oued - 2011 ................................................................................................... 25
Figure ‎I-19: Inondations de Ghardaïa - 2008 ........................................................................................................ 25
Figure ‎I-20: Inondations de Skikda (Les allées) - 1984.......................................................................................... 25
Figure ‎I-21: Inondations de Jijel (Chekfa) ............................................................................................................. 25
Figure ‎II-1:Courbe QDF(Hingray et al ,2009) ........................................................................................................ 31
Figure ‎II-2: Principe de la méthode du basée sur une distribution de Gumel (Musy, 1998). .............................. 33
x
Figure ‎II-3: Définition les lignes isochrones.(Satin et al,1999) ............................................................................. 33
Figure ‎II-4:Influence d'un élément de pluie uniforme d'intensité i à l'exutoire.(Satin et al,1999) ...................... 34
Figure ‎II-5: Distribution de la Pluit nette par la méthode de horton.(Musy,1998) .............................................. 43
Figure ‎II-6:Méthode de l'indice φ.(Musy,1998).................................................................................................... 44
Figure ‎II-7:Méthode du SCS.(Musy,1998) ............................................................................................................. 46
Figure ‎II-8: Schématisation du processus infiltration selon Grenne et Ampt. (Musy, 1998) ............................... 49
Figure ‎II-9:Hypothèse de linéarité de l'Hydrogramme Unitaire (Musy, 1998) ..................................................... 54
Figure ‎II-10:Définitions des caractéristiques de l'Hydrogramme Unitaire. (Hingray et al, 2009) ........................ 55
Figure ‎II-11: La réponse impulsionnelle.(Hingray, 2009) ...................................................................................... 57
Figure ‎II-12:Construction de la courbe en S (Hingray et al, 2009)........................................................................ 58
Figure ‎II-13: Hydrogramme unitaire synthétique sans dimension - SCS, et son approximation triangulaire. ..... 59
Figure ‎II-14 : Schématisation d'un bassin versant du 3𝑒 ordre et probabilités de transition d'un cours d'eau
d'ordre 𝑖à un cours d'eau d'ordre 𝑗 > 𝑖 (Hingray, 2009) ........................................................................... 61
Figure ‎III-1: Situation géographique de la région d’étude (source : google map) ................................................ 63
Figure ‎III-2: Variations mensuelles des précipitations. Boudjoudjou, 2010)........................................................ 64
Figure ‎III-3: Régime saisonnier des précipitations. (Boudjoudjou, 2010) ............................................................ 65
Figure ‎III-4: : Carte hydrographique de la région de Jijel au 1/200 000 (Boudjoudjou, 2010) ............................. 66
Figure ‎III-5: Bassin versant de l’Oued Djendjen (Fatih, 2009) .............................................................................. 67
Figure ‎III-6 : Carte d'occupation de sol "bassin versant de Oued Djendjen" ....................................................... 69
Figure ‎III-7: Carte des limites de la zone d'étude. .................................................................................................... 70
Figure ‎III-8: Rectangle équivalant du bassin versant d'oued Djenjen. ................................................................. 72
Figure ‎III-9: Courbe hypsométrique du bassin versant ........................................................................................ 73
Figure ‎III-10: Carte caractéristique des altitudes du bassin versant d'Oued Djendjen ............................................ 74
Figure ‎III-11: Les ordres de cours d'Oued Djendjen selon la classification de Horton. ........................................ 80
Figure ‎IV-1: Diagramme explicative de la méthodologie d'extraction des données débit. ................................ 83
xi
Figure ‎IV-2:Ajustement statistique de débit maximum journalier de la station Missa. ....................................... 84
Figure ‎IV-3: Ajustement statistique de débit de pointe maximale journalier de la station Missa. ..................... 85
Figure ‎IV-4:Ajustement statistique de débit de point horaire de la station Missa. ............................................ 85
Figure ‎IV-5: Ajustement par la loi de Gumbel ...................................................................................................... 88
Figure ‎IV-6:Organigramme explicatif les différentes étapes appliqués de la méthode de Gradex ..................... 89
Figure ‎IV-7: Ajustement statistique des pluies maximales journalières de la station Texenna ........................... 91
Figure ‎IV-8Principe de méthode Gradex (les débits maximaux journaliers) ........................................................ 93
Figure ‎IV-9:Principe de méthode Gradex (les débits de pointe journaliers) ........................................................ 93
Figure ‎IV-10: Les composantes de l’écoulement( Cour Musy,2005) .................................................................... 95
Figure ‎IV-11 :Ruissèlement simulé vers sue ruissèlement observé ..................................................................... 98
Figure ‎IV-12: Calibration modèle cascade de NASH ............................................................................................. 99
Figure ‎IV-13:Validation du modèle de NASH ...................................................................................................... 100
Figure ‎IV-14:la pluie de projet double triangle (Kwan Tun lee,2008) ................................................................ 102
Figure ‎IV-15: Organigramme de l'application des modèles P-Q ........................................................................ 104
Figure ‎IV-16: Hydrogrammes de crue obtenu par la modélisation P-D pour T=100ans .................................... 105
Figure ‎V-1: Représentation schématique du dispositif informatique utilisé ..................................................... 110
Figure ‎V-2: Organigramme récapitulative les tâche du logiciel HEC-RAS (US ACE, 2010)................................ 111
Figure ‎V-3: Représentation des termes de l’équation d’énergie(US ACE, 2002) ............................................... 112
Figure ‎V-4: L’interface du logiciel HEC-RAS 5.2.0 ............................................................................................... 114
Figure ‎V-5: Schéma illustrant les fichiers d'un projet(HEC-RAS users mnual).................................................... 114
Figure ‎V-6: L’introduction des conditions aux limites pour l’écoulement ......................................................... 117
Figure ‎V-7: le profil en long du tronçon modélisé d’Oued Djendjen .................................................................. 118
Figure ‎V-8: Profile en traverse de quelque sections........................................................................................... 119
Figure ‎V-9: Profil du vitesse ................................................................................................................................ 120
Figure ‎V-10: courbe de tarage ............................................................................................................................ 120
xii
Tableau ‎I-1: Divers inondations au monde (Benamia, 2012) ............................................................................... 22
Tableau ‎I-2: Diverses inondations en Algérie (Boulghora, 2007; Yahiaoui, 2012). .............................................. 24
Tableau ‎II-1: loi statistique (Dubreuil, 1974) ........................................................................................................ 28
Tableau ‎II-2: Exemple de formule pour l'estimation du temps de concentration 𝑡𝑐(Anctil et al, 2012) ............. 36
Tableau ‎II-3:Formules utilisant uniquement la superficie du bassin (Ouarda et al ,2001) .................................. 38
Tableau ‎II-4:Formules faisant intervenir la précipitation. (Ouarda et al ,2001)................................................... 39
Tableau ‎II-5:Formules faisant intervenir la période de retour(Ouarda et al ,2001)............................................. 39
Tableau ‎II-6: Paramètres de l’équation de Horton pour différents types de sol (Musy, 1998) ........................... 42
Tableau ‎II-7: Capacités initiale 𝒊𝟎, et finale 𝒊𝒇 d’infiltration en fonction des classes de sols définies par le SCS
(Musy,1998) ............................................................................................................................................... 43
Tableau ‎II-8:CN pour des conditions normales d’humidité du sol (Chow, 1988) ................................................ 48
Tableau ‎II-9: Les principales fonctions de transformation pluie-débit(Hingray et al,2009)................................. 52
Tableau ‎III-1: Le taux d’occupation du sol du bassin versant de l’oued Djendjen. .............................................. 69
Tableau ‎III-2: La répartition des altitudes en fonction des surfaces partielles. ................................................... 73
Tableau ‎III-3: Classification du relief par l'ORSTOM selon Ds .............................................................................. 75
Tableau ‎III-4: Différente formule du temps de concentration ............................................................................. 77
Tableau ‎III-5:Résultats de temps de concentration.............................................................................................. 78
Tableau ‎III-6: Formule des rapports de confluence (Chow, 1988) ....................................................................... 79
Tableau ‎III-7:Résultats de rapport de confluence ................................................................................................ 79
Tableau ‎III-8 :Caractéristiques géomorphologiques de bassin versant d'oued Djendjen .................................... 80
Tableau ‎IV-1:Teste khi2 ........................................................................................................................................ 86
Tableau ‎IV-2: Les paramètres obtenus par la loi de Gumbel ............................................................................... 87
Tableau ‎IV-3:Les quantiles obtenus par l'ajustement des débits de pointe maximale journalière ................... 87
xiii
Tableau ‎IV-4:Les quantiles obtenus par l'ajustement des débits maximaux journaliers. .................................... 87
Tableau ‎IV-5:Les quantiles obtenus par l'ajustement des débits de pointe horaire ............................................ 88
Tableau ‎IV-6:les paramètres obtenus par la loi de Gumbel pour les Pmax, j ..................................................... 90
Tableau ‎IV-7: les résultats d'ajustement des pluies maximales journalières par la loi de Gumbel ..................... 90
Tableau ‎IV-8: paramètres de la loi Gradex ........................................................................................................... 91
Tableau ‎IV-9: Les quantiles obtenus par les débits maximaux journaliers. ......................................................... 92
Tableau ‎IV-10: Les quantiles obtenus par les débits de pointe maximaux journaliers. ...................................... 92
Tableau ‎IV-11 : les résultats obtenus de débit de crues basé sur les méthodes empiriques .............................. 94
Tableau ‎IV-12 : Fonction objective ....................................................................................................................... 97
Tableau ‎IV-13 : calibration et validation de la fonction de production ................................................................ 97
Tableau ‎IV-14: Les paramètres du model pluie-débit ........................................................................................ 103
Tableau ‎IV-15: Les débits obtenu par la modélisation P-D pour T=100ans ....................................................... 104
Tableau ‎IV-16: tableau récapitulatif des débits centennaux selon plusieurs méthodes ................................... 106
Tableau ‎V-1: Coefficients de Manning ............................................................................................................... 116
Tableau ‎V-2: les principaux paramètre hydraulique calculé pour chaque ....................................................... 121
xiv
Les crues désignent une période de fort débit d'un fleuve ou d’une rivière ou le déparquement
d'une certaine hauteur d'eau, il peut y en avoir plusieurs au cours de l’année. Cette hauteur d’eau
peut être qualifiée en fonction de sa durée de retour. En fait, les épisodes de pluie ou d’orages
intenses, peuvent engendrer une catastrophe hydrologique qui s'appelle l'inondation. Ce terme est
d’acceptation plutôt topographique ou géomorphologique, il faut que le débit de crue dépasse un
seuil donné, au-delà duquel la rivière sort de son « lit mineur », c’est-à-dire son tracé normal, puis
déborde dans son lit dit « majeur », mais aussi le remplissage d’une cuvette topographique par les
eaux de ruissellement local.
Les inondations sont parmi les catastrophes naturelles qui causent des pertes et des dommages
matériels importantes dans le monde. Elles constituent le risque le mieux distribué sur la planète
(White, 1999).Comme beaucoup d’autres pays du monde, l'Algérie est exposée aux inondations
dévastatrices. Si par le passé, ces inondations étaient surtout connues suite aux débordements des
grands cours d’eau dans les grandes plaines agricoles, depuis une vingtaine d’années, elles ont surtout
touché les grandes villes et centres urbains, où elles ont engendré de nombreuses pertes en vies
humaines et des dégâts matériels considérables.
Afin d'évaluer le danger posé par cette catastrophe naturelle dans le territoire, il faut bien estimer le
débit de crue pour une période de retour désirée. Pour ce fait, nous allons mener une étude qui traite
un sujet liée à ces phénomènes à savoir, les inondations dans le bassin d'oued Djendjen situé à la
région Texenna (Wilaya de Jijel).
Dans ce travail, on s’est donc intéressé à l’étude des crues extrêmes pour plusieurs périodes de
retour, en particulier le débit de crue centennal, pour un bassin versant jaugé.
Au fil des années, de nombreuses méthodes et modèles ont été élaborés pour permettre la
détermination des débits de crue. Dans la pratique, le choix d'une méthode dépend essentiellement
des données à disposition sur le bassin versant étudié. Principalement, les données de débits qui
constituent un facteur important à l'étude hydrologique.
Ce mémoire a pour objectif de présenter et appliquer quelques approches utilisées à savoir celles
basées sur les méthodes statistiques (méthode de Gradex, méthode de Gumbel), celles basées sur les
méthodes empiriques, et celles basées sur la modélisation pluies-débits.
xv
La multiplicité de ces méthodes nous amène à interroger, les quelles de ces méthodes sont fiables et
donnent des résultats cohérentes?
Afin de résoudre notre problématique nous menons une étude comparative est évaluative pour
vérifier la validité des méthodes prévues, par la détermination des débits de cours pour une période
de 100 ans.
Notre travail s’intéresse à l’emploi d’un couplage de modèles hydrologique et hydraulique pour la
modélisation des crues et des inondations au bassin versant d'Oud Djendjen.
D'une part, ce débit est utilisé notamment dans la perspective de la modélisation hydraulique qui
permet d'obtenir par le calcul, dans tous les secteurs exposés du cours d'eau, les caractéristiques
hydrauliques des écoulements telle que la hauteur d'eau, la vitesses du courant, la surface mouillé, le
nombre de Froude etc.., atteintes pour la gamme de débit simulées, à l'aide des outils comme ArcGis
et HEC-RAS.
Les résultats issus de l’étude hydrologique ainsi que la géométrie, la pente et la rugosité de l’Oued
Djendjen sont utilisés comme données d’entrée dans le modèle hydraulique développé à laide de
HEC-RAS. En sortie, on aura les lames d’eaux correspondantes et quelques variables
caractéristiques d’un écoulement fluvial.
Enfin, et d'après la simulation hydraulique nous pouvons déterminer les zones susceptibles d'être
inondable à la région étudiée.
xvi
I. Numéro du chapitre
— Aristotle
L es inondations sont des phénomènes complexes. Elles représentent la

source de catastrophes naturelles la plus fréquente dans plusieurs
régions du monde et sont responsables de la majorité des dégâts
provoqués par l'ensemble de catastrophes naturelles. Ce chapitre présente des
notions et concepts fondamentaux liés à ce phénomène, ainsi que les différentes
causes qui l'engendre, leurs types, et les risques associés. A la fin du chapitre,
nous passons en revue certaines des inondations qui ont eu lieu en Algérie et
par ailleurs dans le monde
I.1. Cycle hydrologique

Comme point de départ pour l'étude de l'hydrologie, il est utile de considérer le cycle
hydrologique. C'est un modèle conceptuel de la façon dont laquelle l'eau transite entre la
terre et l'atmosphère dans différents états (gaz, liquide ou solide). (Davie, 2002)
La notion de cycle de l’eau a donné lieu, par le passé, à de nombreuses spéculations. Comme
le souligne Roche (1969), « l’idée simple que la pluie pouvait être à l’origine de toutes les
eaux coulant dans les rivières a mis longtemps à faire son chemin ».
En peut définir les mouvements de l’eau sur la terre comme des flux dans un système à
circuit fermé, où les pertes dans l’espace et les apports depuis le noyau de la terre sont
négligeable. Une représentation simple de ce circuit implique des mouvements d’eau des
masses océaniques vers l’atmosphère, de l’atmosphère vers les masses continentales, puis des
masses continentales vers les océans (figure I-1).
1
Chapitre I — Problématique &généralités sur les inondations
Figure ‎I-1: Conceptualisation du cycle de l’eau(Raghunath,2006)
Ce schéma peut être résumé par l'organigramme suivant:
Atmosphère
Précipitation
Atteint la Atteint les plans

S'évapotranspire S'évapore
végétation d'eau libre
S'évapore Atteint le sol Ruisselle
S'évapotranspire S'infiltre S'écoule Atteint les océans
Figure ‎I-2: Cycle de l'eau(Laborde,2000)
I.1.1. Phases principales

Le cycle de l'eau est définit par l'ensemble de tous les processus de transformation de
l'eau sur la terre dont les phases principales sont :
2
I.1.1.1. Précipitation
Les précipitations, ou, plus justement, les météores, sont constituées par toutes les eaux
météorologiques qui tombent sur la surface de la terre sous forme liquide (pluie, bruine) ou
sous forme solide (neige, grêle, grésil) ou encore qui se déposent sur le sol (rosée, gelée
blanche, givre). (Mazzouli, 2012)
Elles sont provoquées par un changement de température ou de pression. La vapeur d'eau

de l'atmosphère se transforme en liquide lorsqu'elle atteint le point de rosée par
refroidissement ou augmentation de pression. Pour produire la condensation, il faut
également la présence de certains noyaux microscopiques, autour desquels se forment des
gouttes d'eau condensées. La source de ces noyaux peut être océanique, continentale ou
cosmiques. Le déclenchement des précipitations est favorisé par la coalescence des gouttes
d'eau. L'accroissement de poids leur confère une force de gravité suffisante pour vaincre les
courants ascendants et la turbulence de l'air, et atteindre le sol.
Les précipitations sont exprimées en intensité (mm/h) ou en lame d'eau précipitée (mm)
(rapport de la quantité d'eau précipitée uniformément répartie sur une surface). (Musy,2005)
On distingue trois types de précipitations:
 Les précipitations convectives : sont des précipitations de courte durée, très intenses,
de type orageux avec une faible extension spatiale. Ce type de précipitation dû à une
ascension rapide de masses d'air dans l'atmosphère;
 Les précipitations orographiques, se produisant au niveau des massifs montagneux. Il y

a rencontre entre la masse d'air chaud et humide et un relief topographique. L'intensité
et la fréquence des précipitations assez régulières dépendent de l'altitude, de la pente et
de son orientation;
 Les précipitations frontales (type cyclonique) se produisant quand il y a la rencontre

entre deux masses d'air de température et de vitesse différentes. Le glissement d'une
masse d'air chaud sur une masse d'air plus froid donne naissance à des précipitations
peu intenses mais longues et étendues. (Mazzouli, 2012)
I.1.1.2. Transport par les vents et les courants
Les nuages formés par évaporation peuvent être transportés par les vents et les
courants. Ces mouvements d'aire sont générés par le gradient de pression qui existe entre ces
3
centres de haute et basse pressions. L'existence de ces centres est directement reliée au
gradient de température entre des endroits exposés différemment au soleil.(Bennis, 2009)
I.1.1.3. Evaporation
Le processus d'évaporation dépond de la température, du vent, des rapports

énergétiques entre la surface et l'atmosphère et de facteur physique caractérisant la nature et
les propriétés de la surface. L'évaporation à partir des surfaces d'eau libre dépond de la
surface et de la salinité de l'eau. L'évaporation à partir d'un sol nu dépond de la surface de la
teneur en eau du sol, de sa capillarité, de son albédo. Le flux évaporation à partir d'un
couvert végétal dépond de la résistance aérodynamique au vent et de la résistance de la
canopée ou de la surface au moment de l'eau à travers les stomates du végétal. (Mazoulli,
2012)
I.1.1.4. Transpiration
La transpiration est le processus par lequel la vapeur d'eau s'échappe des feuilles de la
plante, la quantité d'eau transpirée par les plantes dépend des mêmes facteurs
météorologiques influençant l'évaporation physique mais aussi l'humidité du sol dans la zone
racinaire, de l'espace de la plante et de son âge, du développement de son feuillage et de la
profondeur des racines. (Mazoulli, 2012 ; Reghounath, 2006)
I.1.1.5. Evapotranspiration
L'évapotranspiration regroupe le flux de transpiration et le flux d'évaporation.

Ressource indispensable à la croissance des plantes des plantes. L'eau est un élément
constitutif majeur de la matière végétale, mais aussi une source d'hydrogène et d'oxygène
pour la photosynthèse (hydrolyse). L'eau consommée par les plantes n'est cependant pas
essentiellement consacrée à la formation de matière végétale mais à la transpiration, moteur
de la circulation de la sève depuis les racines vers les feuilles. (Roche et al, 2012)
I.1.1.6. Infiltration/ Percolation
 L'infiltration désigne le mouvement de l'eau pénétrant dans les couches superficielles

du sol et l'écoulement de cette eau dans le sol et le sous-sol, sous l'action de la gravité
et des effets de pression.
 La percolation représente plutôt l'infiltration profonde dans le sol, en direction de la

nappe phréatique. Le taux d'infiltration est donné par la tranche ou le volume d'eau qui
4
s'infiltre par unité de temps (mm/h ou m3/s). La capacité d'infiltration ou l'infiltrabilité

est la tranche d'eau maximale qui peut s'infiltrer par unité de temps dans le sol et dans
des conditions données. L'infiltration est nécessaire pour renouveler le stock d'eau du
sol, alimenter les eaux souterraines et reconstituer les réserves aquifères. De plus, en
absorbant une partie des eaux de précipitation, l'infiltration peut réduire les débits de
ruissellement. (Musy, 1999)
I.1.1.7. Ruissellement
L'excédent des précipitations qui ne s'est pas infiltré ou évaporé ou n'a pas été
intercepté par la végétation, va s'écouler selon la pente du terrain. C'est le ruissellement de
surface qui alimente les rivières et les fleuves se déchargeant dans les mers et océan sous
formes ruissellement.(Bennis, 2009)
I.1.1.8. Interception
Une partie de l'eau des précipitations n'arrive pas au sol. Une fraction est interceptée
verticalement par le couvert végétal et peut s'évaporer dans l'atmosphère ou retourner au sol
par égouttement de la végétation ou par ruissellement le long des troncs dans le cas des
forêts. (Mazzouli, 2012)
I.2. les paramètres fondamentaux du cours d'eau

I.2.1. Lit de la rivière
Le lit d'un cours d'eau est l'espace qu'il occupe en permanence ou temporairement. On
distingue:
 Le lit mineur, est l'espace occupé par l'écoulement pour les crues courantes délimité par
les berges et occupe en période de pluviosité normal; (Mazzouli, 2012)
Figure ‎I-3: Lit mineur d'un cours d'eau (Bachi, 2011).
5
 le lit majeur est plaine inondable. Il est limité par les plus hautes eaux. les bordures du
lit majeur ne sont mises en eau que les crues extrêmes avec une hauteur d'eau assez
faible. Les vitesses d'écoulement y sont faibles et les particules les plus fines (limons,
argiles) s'y déposent par sédimentation, (Mazzouli, 2012)
Figure ‎I-4: Lit majeur d’un cours d’eau (Bachi, 2011)
 Le lit d'étiage, occupé pendant la période de sécheresse et correspondant au niveau des

basses eaux normales
 Lit moyen : Il correspond à l’espace fluvial ordinairement occupé par la ripisylve.
Figure ‎I-5: Lit moyen du cours d'eau(Bachi, 2011)
I.2.2. Rive / berge

Rive et berges sont souvent confondues à tort. La berge est le talus incliné qui sépare le
lit mineur et le reste du lit majeur. Sa localisation est donc assez précise.
La rive est milieu géographique qui sépare les milieux aquatique et terrestre. Elle démarre au
sommet de la berge et constitue une plate plus ou moins étendue qui reste sous l'influence
du milieu aquatique. (Degoutte, 2012)
6
I.2.3. Ripisylve
La ripisylve est la formation végétale naturelle située sur la rive. Elle peut être limitée à
un cordon arboré étroit qui souligne le bord du lit mineur de la rivière ou bien elle est un
véritable foret alluvial s’étendant sur plusieurs dizaines ou centaine de mètre de pat et d’autre
du lit mineur.
Cette foret occupe tout ou partie du lit majeur. C’est un milieu inféodé à la rivière,
particulièrement riche en termes de diversité floristique. Il comporte des strates herbacées
souvent très diversifiées et des strates buissonnantes, arbustives et arborescentes composées
d’un nombre restreint d’espèces.
La végétation des bancs, des berges et des rives doit être considérée comme un continuum.
(Degoutte, 2012)
I.2.4. Alluvions et substratum

Une rivière coule généralement sur ses alluvions. Les alluvions sont les grains fins ou
grossiers alternativement déposés ou repris par le courant. Les alluvions recouvrent le
substratum rocheux formé d'un roche dure ou bien d'un roche plus moins tendre (schistes,
grés, marnes). (Degoutte, 2012)
Figure ‎I-6: Représentation de termes géographique usuels (Degoutte, 2012)
I.3. Les événements-extrêmes

I.3.1. Etiage et tarissement
Les basses eaux correspondent aux plus petites hauteurs d'eau observées sur une rivière.
D'une façon plus générale, le terme "étiage" les plus petits débits ou les plus petites hauteurs
7
d'eau observés sur un cours d'eau. La période de basses eaux désigne quant à elle la (les)
période(s) de l'année où les débits ou les hauteurs d'eau sont les plus faibles. (Musy, 1998)
Le tarissement d'un cours d'eau constitue le processus principal durant les étiages.(Lang,
2007) Il est défini comme la décroissance régulière des débits des basses eaux en fonction du
temps. En absence de pluies, l'évaporation et la transpiration végétale épuisent
progressivement les réserves en eau souterraine du bassin versant. Les débits décroissent
alors régulièrement. Le tarissement est non influencé en absence de pluie et
d'évapotranspiration (FigureI-7.a) ou influencé dans le cas où de nouvelles pluies viendraient
alimenter les réserves (Figure I-7.b). (Musy, 1998)
Temps (J) Temps (J)

a) Non influencé b) Influencé
Figure ‎I-7: Exemple de courbe de tarissement (Musy, 1998)
I.3.2. Crue
La crue correspond à l’augmentation de la quantité d’eau qui s’écoule dans la rivière et
peut concerner l’ensemble du lit majeur de la rivière. De nombreux géographes et
hydrologues ont adopté le critère qu'une rivière est en crue lorsque son débit est trois à cinq
fois supérieur à son débit moyen. De façon plus pratique, on admet qu'une rivière est en
crue lorsqu'elle déborde des limites de son lit mineur. Il s'agit d'un phénomène naturel
périodique qui n'est exceptionnel que lorsque les débits deviennent considérables par
rapport à son module. On parle alors de crue critique, laquelle peut engendrer une
inondation sur les zones riveraines. Une crue se caractérise par son hydro-gramme graphique
qui représente les variations de débit en fonction du temps. Plus précisément, c'est la partie
montante de Cette hydrogramme qui est appelé « crue », la partie descendante étant « la
décrue ». Une crue se définit par différents critères : sa genèse, sa durée, sa fréquence, son
débit de pointe et son volume (Salomon, 1997).
8
I.3.2.1. Hydrogramme de crue
L'hydrogramme est la série chronologie des débits successifs enregistrés à un endroit

donné d'une rivière à la suite de précipitations ou de fonte des neiges.
La figure (I.8) montre un schéma d'hydrogramme type mesuré à l'exutoire d'un bassin
versant.la nomenclature utilisée est suivant:
 Le temps de montées 𝑇𝑀 : est le temps écoulé entre les points B et C qui délimitent la
courbe de concentration;
 Le temps de base𝑇𝐵 :est le temps écoulé entre les points B et D; il représente la durée
de ruissellement de surface,
 Le temps de réponse 𝐿𝑎𝑔 est le temps écoulé entre F et C, F représentant le centre de

l'hyétogramme;
 Le débit de pointe 𝑄𝑃 : est le débit maximal atteint en période de cure. C'est ce débit qui
est normalement utilisé dans la conception des ouvrages hydrauliques(ponceau,
conduite d'égout, évacuateur de crues.....) (Bennis, 2009)
Figure ‎I-8:Hydrogramme de crue. (Bennis, 2009)
9
I.3.2.2. Genèse de crue
Il est relativement aisé d'imaginer la genèse d'une crue à partie de la pluie. Sur une
surface réduite, la fraction de la pluie qui n'a pas été interceptée et ne s'est pas infiltrée,
ruisselle sur le terrain. Cette eau de ruissellement chemine sur le sol par des voies diverses,
une partie s'infiltre éventuellement et le reste rejoint le réseau hydrographique. Une
multitude de filets d'eau courent ainsi sur le sol, chacun n'obéissant apparemment qu'à son
propre caprice; ces filets se rejoignent, se subdivisent, se mélangent et parviennent à un
collecteur plus important.
La course erratique du filet d'eau est en réalité soumise à des lois physiques très strictes mais
la complexité de leur application dans la nature est telle, chaque brin d'herbe, chaque caillou
ayant une action personnelle sur l'écoulement, que seul le résultat final est accessible à une
observation quantitative utile. Les caprices du parcours ont battu les cartes et on se trouve
devant un fait statistique. Au cours de ce brassage, chaque écoulement élémentaire, pris à un
instant donné, dépend de la nature physique du bassin à cet emplacement et des conditions
de l'averse. Il est donc naturel que l'influence des caractéristiques du terrain se manifeste
également de façon statistique sur le ruissellement (principalement: capacité apparente
moyenne d'infiltration, pente, couverture végétale). (Roche, 1963)
I.3.2.3. Processus de la genèse des crues et des inondations
Comprendre ce processus de formation des crues revient à analyser les différents

facteurs concourant à la formation et à l’augmentation temporaire des débits d’un cours
d’eau, on distingue :
I.3.2.3.1. Abondance des précipitations
C'est le principal processus explicatif des crues et des inondations. On peut distinguer
plusieurs types :
 Les averses brutales de type orageux où les orages porteurs de précipitations ont des
lames d'eau de 200 à 300 mm en quelques heures.
 Les averses durables qui sont moins brutales car elles comportent en fait des séries
d’averses successive s'entrecoupées de périodes de tranquillité, mais elles sont tout
aussi dévastatrices. (Merabat,2006)
Les pluies de type orographique et cyclonique, En règle générale, la puissance de la crue et

les dégâts qu'elle occasionne sont d'autant plus forts que les précipitations sur le bassin
10
versant ont été fortes, intenses et durables, que les quotients d'écoulement sont plus élevés et
que la morphologie et la nature des matériaux se prêtent à l'érosion. (Salomon,1997)
I.3.2.3.2. Le ruissellement
Représenter la part de l’eau qui n’a pu s’infiltrer dans le sol. Il dépend de la nature du
sol, de sa pente, de son occupation de surface et de l’intensité de l’épisode pluvieux.
(Merabat, 2006)
I.3.2.3.3. Fusion brutale des neiges
En effet, en fin d'hiver l'arrivée de front chaud (avec vent et pluie) va fusionner
rapidement la neige qui couvre les bassins versants montagneux (Alpes, Pyrénées en France
par exemple) ce qui provoque de grandes crues particulièrement importantes sur un sol
parfois gelé et sans végétation. Mais le plus souvent l'alimentation d'origine neigeuse se
combine à celle des averses. (Salomon, 1997)
I.3.2.3.4. Débâcle des glaces fluviales
Celle-ci intervient à la suite de la rupture de barrages naturels de glace que ce soit en

haute montagne (Alpes, Pyrénées) soit dans les pays des hautes latitudes (cas de grands
fleuves sibériens ou canadiens). (Merabat, 2006)
I.3.2.3.5. Accidents
Rupture accidentel de barrage naturel ou artificiel. Parfois cet accident est lui-même lié à
un autre phénomène exceptionnel (glissement de terrain, lié à des précipitations importantes,
volcanisme, séisme…). (Merabat, 2006)
I.3.2.3.6. La propagation de l’onde de crue
Elle est fonction de la structure du lit et de la vallée alluviale, notamment de la pente et

des caractéristiques du champ d’inondation. (Asnouni, 2014)
I.3.2.4. Facteur qui influe sur les crues
Un grand nombre de facteurs autres que celles dictés précédemment intervient, souvent
de façon aggravante, dans les dégâts occasionnés par les inondations. Parmi ces derniers, on
peut rappeler l’influence des facteurs naturels et anthropiques.(réf)
11
I.3.2.4.1. Influence des facteurs naturels
 Surface et forme du bassin versant: pour une même surface, l'allure de l'hydrogramme
de crue résultant d'une pluie donnée, est très différente suivant la forme du bassin
versant ; un bassin très allongé ne réagit pas comme un bassin de forme ramassée ce
qui a amené l'utilisation d'un indice « coefficient de compacité de Gravelius » pour
caractériser la forme du bassin versant. Plus l’indice est faible, plus la concentration des
eaux apportées par les affluents est rapide et plus les crues risquent d'être brutales et
bien différenciées (Roche, 1963).
Figure ‎I-9: Influence de la forme de bassin versant sur l'hydrogramme de crue (musy et higy,
2003)
 Configuration du lit du cours d'eau : la pente exerce une influence directe sur la rapidité
de l'écoulement et donc sur la puissance de la crue. Si le profil en long du cours d'eau
est assimilable à une suite de segments plus ou moins pentus, il exerce une action
visible sur la crue. Le profil en travers est aussi important ; pour un même débit de
crue, un cours d’eau encaissé verra sa hauteur d'eau monter beaucoup plus vite qu'un
cours d'eau à profil plus évasé. En contrepartie, ce dernier a de plus grands risques de
débordements (Cosandy &Robinson, 2012).
 Densité des cours d’eaux et perméabilité du bassin versant : La densité des cours d'eau
est fonction de la nature du sol. Le réseau est d'autant plus développé et complexe que
le terrain est moins perméable. Plus l’imperméabilité est forte et plus les eaux de pluies
sont disponibles pour le ruissellement de surface. Ainsi, les cours d'eau drainant les
12
régions imperméables ont une plus forte probabilité à développer des crues
dangereuses (Asnouni, 2014).
 Sol: le sol peut stocker de l’eau en remplissant partiellement (présence d’air) ou

totalement (absence d’air) sa porosité. Cette capacité de stockage dépend de la
composition granulométrique des éléments solides du sol (texture), du mode
d’arrangement de ces éléments avec la matière organique (structure) et de leur nature
minéralogique. Pour une masse de sol sec donné, le stock en eau d’un sol argileux et
bien structuré sera plus élevé que celui d’un sol sableux et compact. Le sol constitue
une réserve disponible d’eau pour le développement de tous les organismes vivants,
notamment les plantes. L’eau non stockée s’écoule en profondeur (verticalement ou
latéralement) pour alimenter les nappes souterraines, les lacs, les rivières et la mer. La
porosité d’un sol conditionne l’infiltration verticale de l’eau depuis les parties
supérieures vers les parties inférieures. Tout ralentissement de l’infiltration entraînera
une stagnation et une accumulation de l’eau dans la porosité. Si ce phénomène se
produit dès la surface, l’eau ne pourra pas s’infiltrer et s’écoulera par ruissellement en
suivant la pente topographique. Au-dessus d’une nappe phréatique plus ou moins
mobile latéralement, l’eau qui s’infiltre depuis la surface remplit progressivement la
porosité et sature le sol en remontant vers la surface. Le manque d’absorption de l’eau
par les sols est une des principales causes de la genèse des inondations. Lorsque
l’intensité des pluies est supérieure à la capacité d’infiltration de la surface du sol, on
parlera de ruissellement « hotonien ». En revanche, un sol, partiellement ou totalement
saturée en eau par une nappe souterraine superficielle, génèrera du ruissellement « par
saturation ». Dans des cas bien particuliers, ces deux types de ruissellement peuvent se
produire en même temps et produire subitement des inondations catastrophiques.
(Montoroi, 2012)
 Couvert végétal: l'on sait depuis longtemps qu'une couverture végétale dense réduit et
ralenti considérablement l'écoulement (Cosandy & Robinson, 2012), elle agit sur le
ruissellement, retient une quantité de pluie, évapore l’eau, donc elle est considérée
comme un facteur indispensable dans la protection contre les crues.
I.3.2.4.2. Facteurs anthropiques
Les actions de l’homme se résument principalement à :
13
 Occupation de zones riveraines : la concentration des personnes et l’accumulation des

biens dans les champs d’inondation génère le risque d’inondation par accroissement
des enjeux, augmentation des vulnérabilités et absence de prévention-précaution;
 L’imperméabilisation des sols en conséquences d'aménagements urbains ou agricoles;
 Défaillance des dispositifs de protection comme les digues et déversoirs;
 L’aménagement parfois hasardeux des cours d’eau Beaucoup de rivières ont été
modifiées localement sans se soucier des conséquences en amont ou en aval. Ces
aménagements peuvent avoir pour conséquences préjudiciables l'accélération des crues
en aval et l'altération du milieu naturel. (Salomon, 1997 ; Bachi, 2011)
I.3.3. Type de crue

La typologie la plus simples et la plus classique distingue les crues rapides et les crues
lents.
 Les crues rapides: se produisent à la suite d'épisodes de pluie intense et brutaux qui
intéressent en général des surfaces restreintes. La conséquence la plus notable pour les
riverains et la rapidité de la montée des eaux.
 Les crues lentes: concernent des bassins versants de grande dimension qui sont
rarement entièrement soumis à des épisodes pluvieux simultanés.
Les mouvements des grandes masse d'air donnent lieu le plus souvent à des pluies longues
mais d'intensité modérée et intéressant de grandes étendues. Elles provoquent des crues
lentes.au contraire, des épisodes convectifs intéressant des zones de plus faible superficie
que les précédents peuvent provoque des crues très rapides. (Degoutte, 2012)
La typologie ci-dessus est particulièrement simple car elle porte seulement sur la forme de
l'hydrogramme de crue, pointu ou étalé. Pour aller plus loin, il faut tenir compte de ce qui
caractérise l'événement. L'événement considéré n'est pas une pluie mais une crue, c'est-à-dire
une transformation de la pluie en débit par le bassin versant.
En parlant de crue de plaine ou torrentielle, on voit que l'on est tenté de parler de la zone
affectée par la crue et que l'on se rapproche plus du complexes, et c'est la raison que de celui
de crue. Toutes ces distinctions sont complexes, et c'est la raison pour laquelle il n'existe pas
de classification universelle à la fois simples et complète. (Degoutte, 2012)
14
Nous en proposons une, mais d'autres seraient tout aussi pertinentes sinon plus, selon
l'objectif que l'on poursuit. C'est en fait une classification des crues, c'est-à-dire des causes.
Nous constaterons sur divers exemples qu'un même épisode peut être classé de plusieurs
façons, car les frontière entre les phénomènes sont indicative, car les origines météologiques
peuvent être mixtes ou car un même épisode peut intéresser la montagne et la
plaine.(Degoutte, 2012)
I.3.4. Inondation
L'inondation est une submersion (rapide ou lente) d'une zone habituellement hors de
l'eau. Les inondations comprennent les débordements hors d'une cour d'eau, les remontées
des nappes; les ruissellements résultant de fortes pluies d'orage, l'inondation par le
débordement des cours d'eau est la plus fréquente dans le monde. La crue est le facteur
majeur induisant à l'inondation. C'est une augmentation du débit d'eau qui s'écoule dans le
cours d'eau. (Remini, 2005)
I.3.4.1. Type d'inondation
On distingue plusieurs types d’inondations provoquées par le débordement des cours

d’eau et des rivières :
I.3.4.1.1. Inondations de plaines
Les inondations de plaine se produisent à la suite d'épisodes pluvieux océaniques

prolongés mais d'intensités modérée, s'abattant sur des sols où le ruissellement est long à
déclencher, sur des bassins versants moyens à grands (supérieur à 500 km2). Le cours d'eau
sort lentement de son lit ordinaire pour occuper son lit majeur et inonder la plaine pendant
une période relativement longue. Ces phénomènes concernent particulièrement les terrains
bas ou mal drainés. Sa dynamique lente perdure plusieurs semaines. (Merabet, 2006)
Figure ‎I-10: Inondation de Plaine. (Merabet ,2006)
15
I.3.4.1.2. Inondations par remontées des nappes phréatiques
Elles correspondent à une inondation par débordement indirect. Elles se manifestent

par la remontée de la nappe phréatique qui affleure en surface et/ou par l’intrusion d’eau
dans les différents réseaux d’assainissement. (Merabet, 2006)
Les désordres liés à ce type d’inondation se traduisent par des remontées sous bâtisses,
l'ennoyage des ouvrages souterrains, la déstabilisation des pentes, des tassements de remblai,
des glissements de terrain, des perturbations sur les réseaux publiques, des pressions sous les
constructions, etc. (Boubchir, 2007)
Figure ‎I-11: Inondation par remontées des nappes phréatiques. (Merabet, 2006)
Les désordres liés à ce type d’inondation se traduisent par des remontées sous Bâtisses,
l'ennoyage des ouvrages souterrains, la déstabilisation Tassements de remblai, des
glissements de terrain, des perturbations sur les réseaux Publiques, des pressions sous les
constructions, etc.
I.3.4.1.3. Inondations par crues torrentielles
Les inondations torrentielles sont le résultat de crues débordantes des torrents et des
rivières torrentielles. Elles se caractérisent par une montée rapide des eaux, une vitesse
d’écoulement très élevée et une durée de submersion très courte. Ces crues demeurent,
exceptionnelles et constituent un grave danger pour les populations. (Boubchir, 2007)
16
Figure ‎I-12: Inondation par crues torrentielles. (Merabat, 2006)
Ces inondations touchent principalement les zones de montagne et les cours d'eau du
pourtour méditerranéen, elles ont des vitesses d'écoulement importantes, même dans le lit
majeur. Les bassins versants qui peuvent être affectés par ces phénomènes sont celles qui
ont de petite à moyenne superficie et une forte pente moyenne. Les crues torrentielles se
caractérisent par un très fort transport solide et une profonde modification du lit à l'occasion
de l'événement.
Les dommages imputables à ces phénomènes sont avant tout liés à la vitesse du courant,
renforcés par les matériaux que peuvent charrier les rivières générant de telles inondations.
(Bachi, 2011)
Les inondations torrentielles se caractérisent par un très fort transport solide et une
profonde modification du lit à l'occasion de l'événement. Les dommages imputables à ces
phénomènes sont avant tout liés à la vitesse du courant, renforcés par les matériaux que
peuvent charrier les rivières générant de telles crues. (Bachi, 2006)
I.3.4.1.4. Inondations par ruissellement superficiel
Ce type d'inondation ne se produit pas par débordement d'un cours d'eau, mais par
ruissellement intense généralisé. En milieu urbanisé, la capacité d'absorption par le réseau
d'assainissement superficiel ou souterrain est dépassée. En milieu rural, ce sont des vallons
secs qui se mettent à couler. L'inondation provient à la fois de l'amont, par ruissellement
d'eau sur les fondes supérieurs et de l'aval, par saturation des réseaux d'assainissement. Les
crues peuvent être transformées en véritables torrents et emporter des personnes ou de
voitures. (Degoutte, 2012)
17
Figure ‎I-13: Inondation par ruissellement en secteur urbain. (Merabat, 2006)
I.3.4.1.5. Inondations par rupture d’ouvrage ou d’embâcle
Dans le cas de rivières endiguées, l’inondation survient brutalement soit par

débordement au-dessus de la digue, soit par rupture de la digue. Le phénomène peut être
très brutal et d’autant plus dommageable que le site est proche de la digue. Le fait d'être
derrière un ouvrage de protection dimensionné pour un certain niveau de crue peut donc
rehausser le seuil de risque en cas de rupture ou dépassement de ce dernier. Des secteurs
habituellement hors de l'eau peuvent se trouver brutalement inondés. Un embâcle consiste
en l’obturation d’un cours d’eau par la constitution d’une digue naturelle entraînant une
retenue d’eau importante. La digue peut être constituée par des éléments solides arrachés à
l’amont et charriés par le cours d’eau ou par un glissement de terrain. La rupture d’embâcle
peut se produire plusieurs jours après une période de pluies exceptionnelles ou l’apparition
d’un mouvement de terrain. (SPPPI, 2004 ; Bachi, 2006)
I.3.4.1.6. Inondations marines
Les submersions marines sont des inondations temporaires de la zone côtière par la mer
dans des conditions météorologiques (fortes dépressions et vents de mer) et forts
coefficients de marée. Elles se traduisent par l'invasion par des eaux salées particulièrement
agressives. Elles se manifestent soit lors d’un raz de marée ou de tsunami (occurrence très
faible, mais phénomène dévastateur), soit lors d’une tempête (surcote marine, vents et
précipitations importants) ou en cas de rupture des défenses contre la mer (dans ce cas, les
risques se concentrent le long du boulevard de littoral dans les secteurs dépressionnaires).
(SPPPI, 2004)
La Figure ‎I-14 représenté un organigramme récapitulatif des types d'inondations que nous
avons vus.
18
Inondations par Inondations par

crues torrentielles ruissellement
superficiel
Inondations par
remontées des Inondations
nappes marines
phréatiques
Inondations par
Inondation de rupture
plaines d’ouvrage ou
d’embacle
Figure ‎I-14: Types des inondations
I.4. Différence entre inondation et crue

Sous les termes de crues et d’inondations se cache un phénomène multiforme.
Selon Veyret&Meschinet de Richemond, 2003 : « La crue définit les hautes eaux qui peuvent
demeurer dans le lit mineur du cours d’eau. Dès que celui-ci ne suffit plus à contenir
l’écoulement, l’eau déborde et se répand dans le lit majeur provoquant une inondation. Il
peut donc y avoir une crue sans inondation » La crue se caractérise donc comme une
augmentation des eaux d’un cours d’eau tandis que l’inondation induit la notion de
débordement et d’invasion d’un territoire. Les inondations sont considérées comme des
phénomènes d’origine climatique et hydrologique.
I.5. Risque de l'inondation

Tout risque ne peut être compris ou expliqué sans une étude spatio-temporelle précise.
Comme le souligne Charles-Pierre Péguy (1989) : « On ne peut négliger un facteur spatial : la
repense diffère selon que l’on considère une petite région ou un continent Prenons le cas
des cyclones : l’occurrence de ceux-ci en un point ou un autre ou plutôt selon une trajectoire
ou une autre, croit très vite avec la surface considérée ».
19
Le risque d’inondation reste aujourd’hui le risque le plus répandu, faisant le plus de victimes
et de dégâts. Il contribue pour 58% dans le nombre des victimes des catastrophes naturelles
et pour 31% dans le montant des pertes économiques. La réaction face à ce risque
d’inondation a consisté pendant plusieurs décennies à lutter contre les crues qui signifiait
d'abord que l'on cherchait à agir sur le phénomène naturel la crue et uniquement sur lui, c'est
-à-dire sur une seule composante du risque « l'aléa » et non sur les enjeux ou leur
vulnérabilité. (Bachi, 2011).
I.5.1. Risques majeurs-inondations

Le risque d’inondation est la combinaison de la probabilité d’une inondation et des
conséquences négatives potentielles pour la santé humaine, l’environnement, le patrimoine
culturel et l’activité économique associées à une inondation. (Journal Officiel de l’Union
Européenne, 2007)
Les types de risques auxquels chacun de nous peut être exposé sont regroupés en cinq
familles:
 Risques naturels : avalanche, feu de forêt, inondation, mouvement de terrain, cyclone,

tempête, Séisme et éruption volcanique ;
 Risques technologiques : d'origine anthropique, ils regroupent les risques industriels,

nucléaire, biologique, rupture de barrage... ;
 Risques de transports collectifs (personnes, matières dangereuses) : sont des risques

technologiques, On en fait un cas particulier car les enjeux varient en fonction de
l'endroit où se développe l'accident ;
 Risques de la vie quotidienne (accidents domestiques, accidents de la route...) ;
 Risques liés aux conflits. (Merabet, 2006)
I.5.2. Dégât éventuel

I.5.2.1. Dommages directs Souvent divisés en dommages
 à la bâtisse
 au contenu de la bâtisse
Dommages aux structures, infrastructures et immeubles privés ou publics; certains sont

facilement monnayables comme les infrastructures et les résidences alors que d'autres le sont
20
beaucoup plus difficilement comme ceux des secteurs commerciaux ou industriels. (Bachi,
2011)
I.5.2.2. Dommages indirects
Dommages généralement économiques résultant de l’inondation. Les exemples

classiques sont la perte de roulement des commerces et industries ou encore les coûts des
secours, des détournements routiers et ferroviaires, des délais et autres phénomènes causés
par la crue. (Bachi, 2011)
I.5.2.3. Dommages secondaires
Dommages surtout économiques que subissent ceux qui ne sont pas directement
inondés. Par exemple, les pertes d'approvisionnement ou d'achalandage par les commerces
et industries. (Bachi, 2011)
I.5.2.4. Dommages d'incertitude
Dommages liés à I’ incertitude des victimes quant au moment et à la probabilité de la

prochaine inondation. Se calcule par le montant que les propriétaires seraient prêts à
débourser pour être délivrés du problème. La perception joue un rôle. (Bachi, 2011)
I.5.3. Dans le monde

Dans le monde, les inondations sont au rang de premier risque naturel et représentent
en France plus de la moitié des risques naturels avec des enjeux humains. Une commune
française sur trois est concernée par les inondations (dont 300 grandes agglomérations).
Pour 160 000 km de cours d'eau, une surface de 22 000 km2 est reconnue particulièrement
inondable : deux millions de riverains sont concernés. Les dégâts causés par les inondations
représentent environ 80 % du coût des dommages imputables aux risques naturels, soit en
moyenne 600 M€ - coût annuel calculé sur la période 1980-2000. (Benamia, 2012)
Entre 1994 et 2004, les catastrophes naturelles ont touchées 2,5 milliards de personnes et
ont causé le décès de 478 100 personnes, et des pertes économiques estimées à 690 milliards
de dollars. Parmi ces catastrophes, les inondations demeurent l’une des catastrophes
naturelles les plus dévastatrices, touchant plusieurs pays à travers le monde : (Boulghora,
2006)
21
Tableau ‎I-1: Divers inondations au monde (Benamia, 2012)
Lieu La date Les sinistres

Hambourg - Allemagne 17 février 1962 300 décès.
Marrakech - Maroc 17 août 1995 150 décès
Haïti Septembre 2004 1 300 décès, 1 100 disparus et 3 000 blessés.
Malaisie, Thaïlande, Sri Décembre 2004 290 000 décès et 30 000 disparus.
Lanka et Inde
Maharastra – Inde Juillet 2005 inondations et glissements de terrains, 969 décès
Europe Juillet 2005 70 décès
Figure ‎I-15 : En chine, Huang He, le fleuve Jaune a vu son cours fluctuer au fil du temps, de façon
parfois brutale.(Photo: Wikimédia Commons/Leruswing).
Figure ‎I-16: Les inondations aux abords de Houston. (David J. Phillip - figaro.fr)
22
I.5.4. Contexte l'Algérien

L’Algérie est confrontée aux phénomènes de crues et inondations, qui se manifestent de
façon catastrophique constituant ainsi une contrainte majeure pour le développement
économique et sociale (Lahlah, 2001).
Plusieurs régions du pays sont régulièrement menacées par ces catastrophes naturelles
dont les effets sont souvent intensifiés par d'autres facteurs qui aggravent les effets des
crues, ces événements dramatiques engendrant souvent des bilans lourds de dégâts humains
et matériels ne sont pas toujours liées à des situations météorologiques exceptionnelles et
surviennent dans beaucoup de régions suite à des épisodes pluvieux saisonniers et n'ayant
rien d'exceptionnel. L'inventaire des inondations catastrophiques à travers le pays établi pour
les 30 dernières années montre qu'il n'existe pas de régions prémunies contre ce risque et
que ces événements sont imprévisibles dans le temps et dans l'espace. Cet inventaire fait
ressortir les grandes inondations engendrées par des pluies exceptionnelles généralisées sur
des grands bassins versants et pouvant toucher plusieurs régions atteignant parfois l'ampleur
d'une catastrophe nationale. La figure I-17 représente l'ampleur des inondations en Algérie.
(Lahlah, 2004)
Le tableau suivant cite les inondations les plus importantes qu'a connues l'Algérie ces
dernières décennies:
23
Tableau ‎I-2: Diverses inondations en Algérie (Boulghora, 2007; Yahiaoui, 2012).

Lieu La date Les sinistres
Au moins 11 personnes sont mortes, et près de

EL BAYADH 01 octobre 2011 170 maisons détruites, et deux ponts dans des
inondations très importantes
mercredi et jeudi
Bechar 15 et 16 octobre au moins de 08 personnes sont mortes noyées
2008
plus de 40 morts et des dégâts matériels très
Ghardaïa Octobre 2008
lourds
Sidi Bel Abbes 24 Octobre 2000 02 décès et d'importants dégâts matériels
Ouest algérien 22 octobre 2000 plus de 24 décès
Bab El Oued – 10-11 novembre 733 victimes, 30 000 sans-abris et importants

Alger 2001 dégâts matériels
Plusieurs régions 60 décès et des dizaines de disparus au cours de
Octobre 1994
du pays dix jours d’inondations.
Ghilizane 20 octobre 1993 22 décès et 14 blessés à Oued Rhiou.
20 morts, 500 têtes de bétail emportées et dégâts

Jijel
évalués à 50 millions de dinars
1140 familles sinistrées, 200 hectares de cultures
Constantine
détruites.
Skikda 8000 habitations envahies.
03 février 1984 03 disparus, 1957 familles sinistrées, 02 ponts

Guelma détruits et canalisations d’AEP détruites à 100
%.
777 familles sinistrées et pertes importantes en
Khenchela
Cheptel.
Oum – El –
144 maisons détruites et 429 familles sinistrées
Bouaghi
El – Eulma 01 septembre 1980 44 morts, 50 blessés et 365 familles sinistrées.
Azzazga 12 octobre 1971 40 morts et des centaines d’habitations détruites.
24
Figure ‎I-17: cas significatifs des inondations survenues en Algérie. (Merabet, 2006)
Figure ‎I-18: Inondations de bab el Oued - Figure ‎I-19: Inondations de Ghardaïa - 2008
2011
ffffffffffffffffff
ffffffffffffffffff
Figure ‎I-20: Inondations de Skikda (Les Figure ‎I-21: Inondations de Jijel (Chekfa)
allées) - 1984
I.6. CONCLUSION ffffffffffffffffff
ffffffffffffffffff
La compréhension profonde des phénomènes hydrologiques est indispensable pour
toute étude dans le domaine. Nous avons abordé les principales notions liées aux
inondations. Ces dernières sont engendrés par un autre phénomène hydrologique causé par
25
les fortes précipitations continues pendant une longue duré dites "crue". Il est important de
mener une étude détaillée sur les éléments que nous avons présenté dans ce chapitre pour
mieux comprendre la genèse de ce phénomène. Cela nous sera utile pour pouvoir examiner
et analyser correctement le problème de l'inondation.
26
II. Numéro du chapitre
— Confucius
N ous dirons qu'une rivière est en crue lorsque son débit dépasse une
valeur de faible probabilité; il faut donc distinguer les crues,
phénomène peu fréquent, des hautes eaux. La question que l'on peut
se poser donc à propos des crues est la valeur du débit de pointe de probabilité
donnée (par exemple pour une ville : crue centennale, Un chantier : crues
décennale, Evacuateur de crue d'un barrage : crues millénaire)
Des nombreuses méthodes approuvées à l'estimation de la valeur du débit de
crue (centennale et décennale) d'un bassin versant, sont décrites dans le
présent chapitre.
II.1. Méthode statistique

L'analyse statistique des débits de crue a pour l'objective soit l'estimation de la période
de retour des événements observés soit d'estimer le débit qui correspond à une période de
retour choisie. (Hingrayet al, 2009). Elle s'applique lorsqu'on dispose de longues séries de
données sur lesquelles il est possible d'envisager des traitements statistiques appropriés.
Parmi ces méthodes, on privilégiera l'analyse fréquentielle qui est chère aux hydrologues et
qui permet la détermination de certaines caractéristiques des crues. (Musy, 1998)
II.1.1. L'analyse fréquentielle

L'analyse fréquentielle est une méthode statistique de prédiction consistant à étudier les
évènements passés, caractéristique d'un processus donné (hydrologique ou autre), afin d'en
définir les probabilités d'apparition future. (Meylan et al, 1999)
27
Chapitre II — Etat de l'art sur la modélisation des inondation
Cette méthode se base sur le principe de l'ajustement d'une série d'observations à une
loi de distribution statistique connue. Après quoi, il convient d'en vérifier son adéquation. Il
est alors possible d'en extraire des résultats relatifs à des phénomènes rares et très rarement
observé. (Musy, 1998)
II.1.1.1. Certaines lois adoptées à l'étude statistique
Pour les débits des crues, les fonctions de distribution les plus couramment utilisées
sont les suivant :
Tableau ‎II-1: loi statistique (Dubreuil, 1974)

Fonctions de
Lois les paramètres de la loi
repartitions
µ
1 1 2 La moyenne X
Gauss 𝐹 𝑥 = 𝑒 −2µ 𝑑µ
2𝜋 l'écart typeδx
−∞
µ
1 1 2
𝐹 𝑥 = 𝑒 −2µ 𝑑µ
2𝜋 La moyenne y
Galton −∞
𝑦 −𝑦 l'écart typeδy
µ=
avec 𝛿
𝑦 = 𝑙𝑜𝑔𝑥
6
α= δ
Gumbel 𝐹 𝑥 = 𝑒 −𝑒 𝑎𝑣𝑒𝑐
𝑦 π x
x0 = x − γα
(γ ∶ nombre d′ euler = 0,577)
II.1.1.2. Loi de Gumel-double exponentielle
Cette loi permet l'étude de la distribution des fréquences de valeurs extrêmes

(maximums ou minimums annuels par exemple).
La fréquence cumulative de la distribution de Gumel pour les valeurs extrêmes maximales à

la forme de l'exponentielle double suivante: (Bennis, 2014)
𝑥 −𝑥 0
𝐹 𝑥 =𝑒 −𝑒 𝛼
− ∞ ≤ 𝑥 ≤ +∞ II. 1
𝛼,𝑥0 : Coefficient d'ajustement;
𝛼: Paramètre d'échelle appel Gradex;
28
𝑥−𝑥 0
: La variable réduite de Gumel.
𝛼
Par changement de variable:
𝑥 − 𝑥0
𝑦= II. 2
𝛼
Donc la fonction de répartition de Gumel est :
𝑦
𝐹 𝑥 = 𝑒 −𝑒 − ∞ ≤ 𝑥 ≤ +∞ II. 3
𝑥−𝑥 0
l'équation 𝑦 = , présentée sous forme 𝑥 = 𝛼𝑦 − 𝑥0 est l'équation d'un droite qui
𝛼
représente la loi de Gumbel sur papier à probabilité Gumel. (Touaïbia, 2015 ; Laborde 2000)
A. Les paramètres de la loi
6
𝛼= ∙ 𝛿𝑥 II. 4
𝜋
𝑥0 = 𝑥 − 𝛾 ∙ 𝛿𝑥
Avec: II. 5
𝛾 ∶ 𝑛𝑜𝑚𝑏𝑟𝑒 𝑑′ 𝑒𝑢𝑙𝑒𝑟 0,577

B. Procédé d'ajustement
 Classement des valeurs par ordre croissant en leur affectant un numéro d'ordre
 Calculer la fréquence expérimental;
 Calcul les caractéristique empirique de la série;
 Calculer la variable de Gumel pour chaque valeur observée;
y = − ln − lnF x II. 6
 Reporter les valeurs observées sur papier Gumel ;
 Calculer Le coefficient de corrélation entre les valeurs observées et la variable de

Gumel dont la formule général est:
𝑛
𝑖=1 𝑋𝑖 − 𝑋 𝑌𝑖 − 𝑌
𝑟=
𝑛 𝑛
II. 7
𝑋𝑖 − 𝑋 2 𝑌𝑖 − 𝑌 2
𝑖=1 𝑖=1
𝑋 et 𝑌sont respectivement les moyennes arithmétiques des variables 𝑋et 𝑌.
29
 Si la corrélation est bonne, calculer les paramètres d'ajustement de la droite de Gumel.

la droite de régression ou droite de Gumel est(Touaïbia, 2015) :
𝑥0 = x − 0.45 ∙ 𝛿𝑥 II. 8
où
𝛼 est la pente de la droite
𝑥0 : l′ ordonnée à l′ origine
y ∶ la variable de Gumbel pour une probabilité donnée
C. Quantile
𝑋𝐹 = 𝛼 ∙ 𝑌 + 𝑋0 II. 9
D. Intervalles de confiance d'un quantile
Le CTGREF(1987) propose de construire intervalle de confiance à partir de la relation

suivant pour loi de Gumel (Meylan et al, 2008)
u  u2 
1 1
2
 1  1.396  K  1.1  K 
2 2
 (1.1  K  0.5772)
N N
T1,2 
u2  II. 10
1
1  1.1  2
𝑢1−𝛼 : le quantile de la variable normale centrée réduite correspondant au niveau (1 − 𝛼)

2
choisi pour l'intervalle de confiance
K : le facteur de fréquence pour la loi de Gumel, correspondant au temps de retour choisi:
K  0.7797  ( y p  0.5772) II. 11
L'intervalle de confiance pour les quantiles est alors défini par:
XF − T2 δx < X F < XF + T1 δx II. 12
30
II.1.1.3. Les courbes " débit- durée-fréquence"
L'analyse fréquentielle des débits de crue peut s'effectuer sur les débits moyens maximums
correspondant à différents durées. Ceci conduit alors à l'élaboration de courbes QDF
(courbe débit Q, durée D, fréquence F de non-dépassement).
Les courbes QDF synthétisent les relations statistiques existant entre les caractéristiques
principales de l'hydrogramme de crue, c'est-à-dire le débit moyen maximal Q, la durée d et la
fréquence F, d'où leur nom «QDF». En d'autres termes, elles donnent la valeur du débit
moyen maximal de crue, estimé sur une durée d, correspondant à une période de retour T
donnée (Figure II-2). (Hingray et al, 2009)
Figure ‎II-1:Courbe QDF(Hingray et al ,2009)
II.2. La méthode du Gradex

La méthode du Gradex est une approche hydro-pluviométrique probabiliste du calcul
des débits de crues extrêmes d'une rivière, pour des fréquences de dépassement comprises
entre 102 et 104, soit encore des événements de durée de retour située entre cent et dix mille
ans.
Elle s'applique à des bassins versants de quelques dizaines à plusieurs milliers de kilomètres
carrés recevant des pluies relativement homogènes du point de vue spatial. Pour ces bassins,
le phénomène prépondérant dans la genèse des crues est la transformation de la pluie
météorique en ruissellement direct de surface.
31
La méthode est basée sur le processus physique global (CFGB, 1994):
Pluie brute  Pluie efficace  Fonction de transfert Volume, pointe et forme de la crue
La méthode de s’appuie sur les hypothèses suivantes:
 Les débits maximum recherchés sont provoqués uniquement par des pluies maximales,
uniformément réparties sur le bassin. Il n’y a donc pas diverses origines de formation
des crues (neige, glace), tout comme une grande variabilité spatio-temporelle de la pluie
qui provoque l’écoulement.
 Les pluies maximales et les débits correspondants (débits maximaux) suivent une
même loi de distribution statistique, dite des "extrêmes" en raison de la nature du
phénomène recherché (crues rares). La loi de Gumel est souvent utilisée dans ce but et,
dans ce cas uniquement, le caractère exponentiel de cette distribution est décrit par la
pente de la droite d’ajustement des pluies observées, mesurée sur un diagramme de
probabilité adéquat. La pente de cette droite est le gradient de cette distribution
Exponentielle, d’où le nom de la méthode 𝑡𝑔(𝛼) sur la (Figure II-2).
 A partir d’une certaine valeur de pluie, correspondant à un état de saturation en eau du

bassin, tout excédent de pluie provoque le même excédent du débit (les pertes à
l’écoulement dues à la rétention, l’évaporation, l’infiltration tendent vers une
constante). En terme statistique, ceci revient à dire qu’à partir de cette valeur, le
comportement asymptotique du débit sera identique à celui des pluies. En raisonnant
sur une loi de Gumel et sur un papier de probabilité ad-hoc, la droite de distribution
des débits est alors parallèle à partir de ce seuil à celle des pluies (Figure ‎II-2) est alors
possible de déterminer par la suite la valeur des débits correspondant à des pluies
maximales par simple extrapolation statistique.
32
Figure ‎II-2: Principe de la méthode du basée sur une distribution de Gumel (Musy, 1998).
II.3. la méthode rationnelle
II.3.1. Principe de la méthode
La méthode rationnelle est l'une des plus vieilles méthodes de prévision des crues
(milieu du XIXesiècle).elle est un excellent outil pour la détermination d'un hydrogramme.
En effet, permet tout au long du développement du calcul, de rationaliser les résultats et de
dégager ainsi les meilleures caractéristiques du projet à retenir.
Figure ‎II-3: Définition les lignes isochrones.(Satin et al,1999)
Dans la pratique, la démarche consiste à estimer les débits produits à l’aval de secteurs
élémentaires 𝐴1 , 𝐴2 , … , 𝐴𝑛 délimités par des lignes isochrones de ruissellement (Figure II-3)
Ainsi, la zone 𝐴𝑗 est délimitée par les isochrones 𝑡j−1 et 𝑡j .
Si l’on pose 𝑡 = ∆𝑡, 𝑡2 = 2∆𝑡, . . . , 𝑡𝑛 = 𝑛 ∆𝑡 et qu’il tombe sur cette aire une averse
uniforme d’intensité 𝑖 et de durée𝑡𝑟 , la valeur 𝐶𝑗𝑖 sera la partie de l’averse qui ruisselle sur
cette aire élémentaire 𝐴𝑗 .
33
Dans cette hypothèse, l’évolution du débit mesuré à l’exutoire pour une durée 𝑡𝑟 =
∆𝑡, … , se présenteainsi :
𝑄 (𝑜) = 0
𝑄 (∆𝑡) = 𝐶1 𝑖𝐴1
𝑄 (2∆𝑡) = 𝐶 2 𝑖 𝐴2
........
𝑄 (𝑛∆𝑡) = 𝐶𝑛 𝑖𝐴𝑛
et la somme des temps 𝐴𝑡 représentera la valeur finale du temps que mettra une goutte
d’eau, la plus éloignée hydrauliquement, pour parvenir à l’exutoire considéré. C’est cette
valeur finale qu’on appelle « temps de concentration».
Le pas étant délimité par les lignes isochrones (∆𝑡), si l’averse d’intensité 𝑖 dure 2∆𝑡, 3∆𝑡, . . .,
on obtient l'hydrogramme d’écoulement en ajoutant successivement les hydrogrammes
élémentaires correspondant aux intervalles ∆𝑡, 2∆𝑡, 3∆𝑡, . . ., en allant de gauche à droite sur
le graphique (Figure II-4)
Figure ‎II-4:Influence d'un élément de pluie uniforme d'intensité i à l'exutoire.(Satin et al,1999)

Ainsi, les débits à l’exutoire auront les valeurs suivantes :
𝑄 (0) = 0
𝑄 (∆𝑡) = 𝐶1 𝑖𝐴1
𝑄 (2∆𝑡) = 𝐶1 , 𝑖𝐴1 + 𝐶2 𝑖𝐴2
𝑄 (𝑛∆𝑡) = 𝐶 𝑛 −1 𝑖𝐴 𝑛 −1 ) + 𝐶𝑛 𝑖𝐴𝑛
34
On conviendra alors que, pour une averse uniforme d’intensité 𝑖 et de durée 𝑛∆𝑡 = 𝑡𝑐 tout
le bassin contribue à la formation du débit 𝑄 (𝑛∆𝑡) à l’exutoire. Par ailleurs, si l’averse durait
plus longtemps que𝑛∆𝑡 avec la même intensité𝑖, la sommation d’hydrogrammes identiques,
décalés dans le temps, montrerait que le débit de pointe de l’hydrogramme demeurerait égal
à la valeur de Q (𝑛∆𝑡).
En conclusion, lors d'une averse d'intensité 𝑖 uniforme, la contribution maximale du

ruissellement à l’exutoire sera obtenue lorsque la durée de l’averse sera au moins égale au
temps de concentration 𝑡𝑐 du bassin. (Satin et al,1999)
II.3.1.1. Formule générale
On obtient finalement l’expression de la formule rationnelle :
𝑄𝑝 = K . C . i. A II. 13
Dans laquelle :
𝑄𝑝 : Débit de pointe de l’hydrogramme en 𝑚3 /𝑠 ;
K : Coefficient d’homogénéité se rapporter aux unités1/360;
C:Coefficient de ruissellement dans la limite 𝑜 < 𝑐 <1 ;
i : Intensité en mm/h ;
A :Superficie du bassin en hectares. (Satin et al,1999)
II.3.1.2. Intensité
Sous les hypothèses de la méthode rationnelle, la durée critique d'une pluie d'intensité
moyenne maximale de période de retour 𝑇 est égale au temps de concentration 𝑡𝑐 du bassin.
L'intensité de pluie correspondante est obtenus à partir des courbes IDF étables pour le
bassin étudie ou la région concernée. (Hingrayet al, 2004)
35
‎II-5:Courbe IDF pour différents temps de retour(échelles logarithmique) (Maylan et al,1999)

II.3.1.2.1. Le temps de concentration
Le temps de concentration d’un bassin versant est défini comme le temps mis par l’eau
pour parcourir la distance entre le point le plus éloigné hydrauliquement de l’exutoire et ce
dernier. (Bénis, 2007)
Il existe dans la littérature technique plusieurs formules pour estimer le temps de

concentration 𝑡𝑐 . Nous citons ici les formules les plus utilisées dans le tableau ci-dessous:
Tableau ‎II-2: Exemple de formule pour l'estimation du temps de concentration 𝑡𝑐 (Anctil et al,
2012)
Méthode Expression
𝑡𝑐 = 𝑡𝑕 + 0,00195. 𝐿0,77 𝐼 −0,385
Où
𝑡𝑕 : temps d'humectation [min];
Kripiche 𝟏𝟗𝟒𝟎 L:longueur de la ligne d'eau principale [m];
I:pente moyenne entre l'exutoire et le point du plus
éloigné [m/m].
𝑡𝑐 =0,000325𝐿0,77 𝐼 −0,385 Où
Chow(1962) L : La longueur maximal du cours d'eau[m];
I: est la pente moyenne du bassin (m/m).
36
Formule de la FAA (Fédéral Aviation 𝐿

Agency, 1970) 𝑡𝑐 = 𝑘. 1,1 − 𝐶 .
𝐼 0,333
Où
𝑘 = 3,26 c'est le coefficient le ruissellement;
L: la longueur de parcours 𝑚 ;
I: est la pente représentative, en %.
𝐿0,6 𝑛0,6
𝑡𝑐 = 𝐾
𝑖 0,4 𝐼 0,3
Où
Formule de Ponde cinématique 𝑘 = 6,92
(Viessman et Lewis, 1956) L : est la longueur de parcours en (m);
I : est l'intensité de l'averse en (mm/h);
n : le coefficient de Manning.
I : est la pente moyenne du bassin (m/m).
II.3.1.2.2. Coefficient de ruissellement C
C'est le paramètre de la méthode rationnelle le plus difficile à déterminer. Ce paramètre

englobe le degré d’imperméabilité, les pertes dues à l’infiltration, le stockage dans les
dépressions, la pente du terrain, l’intensité de l’averse et les conditions d’humidité du sol.
Ce paramètre doit être choisi judicieusement car une surestimation de C peut conduire à un
surdimensionnement des ouvrages et donc à des coûts élevés. D’un autre côté, une sous-
estimation de C peut remettre en cause la sécurité d’exploitation des ouvrages.(Bénis,
2014)Le coefficient de ruissellement pondéré se calcule par la relation suivante pour n zones:
𝑛
𝑖=1 𝐶𝑖 𝐴𝑖 II. 14
𝐶=
𝐴
II.4. Les méthodes empiriques

Contrairement aux méthodes statistiques, les modèles empiriques ont généralement une
structure imposée résultant de l'expérience d'une certaine conceptualisation du
comportement hydrologique du bassin versant. Si le degré de conceptualisation est
important, le modèle est dit pseudo-empirique.
37
L'avantage des méthodes empiriques est que le nombre d'observation nécessaires à leur par
métrisation est généralement réduit par rapport à celui nécessaire pour la mise en place d'un
modèle statistique. Elles permettent aussi en général de mieux tenir compte des spécificités
locales du comportement hydrologique considéré. (Hingray et al, 2009)
Il existe un grand nombre de formules empiriques, en cas général, mises au point pour une
région donnée ; leur utilisation doit être faite avec une prudence extrême. Elles peuvent
servir toutefois à caractériser spatialement le phénomène crues en étudiant la répartition
spatiale des index figurant dans ces formules dans les tableaux ci-dessous:(Bois, 2000)
Tableau ‎II-3:Formules utilisant uniquement la superficie du bassin (Ouarda et al ,2001)
Méthode Formule générale

𝑘
𝑙
𝑄 𝐴 10
Francou-Rodier =
𝑄0 𝐴0
𝑄: débit maximum en m3/s;
𝑄0 : débit maximum d'une crue observée dans un
bassin de superficie 𝐴0 (en𝐾𝑚2 );
𝐴 : superficie du bassin versant en 𝐾𝑚2 ;
𝑘 : coefficient de Francou-Rodier.
𝑥
𝑄1 = 𝑄2 𝐴1 𝐴2
𝑄1 : débit de pointe ou débit quotidien moyen
𝑚3 𝑠à déterminer pour non jaugé;
Calcul proportionnel
𝑄2 : débit de pointe ou débit quotidien moyen
𝑚3 /𝑠 au site jaugé, tel que déterminé par une
analyse fréquence de crue;
𝐴1 : superficie en 𝐾𝑚2 du bassin
hydrographiqueen amont de l'endroit où l'on
cherche à établir le débit de crue;
𝐴2 : superficieen Km2 du bassin
hydrographique convergeant vers la station de
jaugeage;
x : exposant habituellement pris égale à 1.
38
Tableau ‎II-4:Formules faisant intervenir la précipitation. (Ouarda et al ,2001)
Q = λ ∙ m ∙ Hm ∙ A
𝑄: débit de crueen 𝑚3 /𝑠 ;
Iskowski 𝐴: superficie du bassinen 𝑘𝑚2 ;

𝐻𝑚 : hauteur moyenne annuelle des
précipitations sur le bassinen 𝑚𝑚 ;
λet 𝑚: des coefficients à déterminer.
Cr ∙ H ∙ A
Q=
3.6 ∙ t c
𝑄 : débit de crue maximum en m3/s;

𝐶𝑟 ∶ coefficients de ruissellement du bassin
pour la crue considérée;
Tura𝐳𝐳a
𝐻 : hauteur totale maximale de précipitation
correspondant au temps de concentration du
bassin, en 𝑚𝑚 ;
𝑡𝑐 : temps de concentration en heures;
𝐴 : aire du bassin versant en 𝐾𝑚2 .
Qm = Hm ∙ A ∙ log A + 2
Qm : débit maximum annuel moyen en [𝑚3 /s];

Fersi A : superficie du bassin versant en 𝑘𝑚3
Hm : ruissellement annuel moyen en 𝑚𝑚
Tableau ‎II-5:Formules faisant intervenir la période de retour(Ouarda et al ,2001)
39
𝑄𝑟
𝐴
= 2𝑘 𝑙𝑜𝑔 +𝑎𝐻𝑚 ( ) 1 + 4𝑙𝑜𝑔 𝑇 − 𝑙𝑜𝑔 𝐴
𝐿
𝑄𝑟 : débit de pointeen [m3/s] de période de

retour T;
𝐴 : superficie du bassin versant en 𝑘𝑚2 ;
𝐿 : longueur de talweg principal en 𝑘𝑚 ;
𝐻𝑚 : hauteur moyenne de la précipitation en
Maillet et Gauthler
𝑚;
𝑘: paramètre (facteur pluviométrique considéré
dans sa variabilité);
𝑎: paramètre (facteur pluviométrique considéré
moyennement).
Qr = q 1 + 0.8 log T
Qr : débit de crueen m3/s de temps de retour T (en
années);
𝑞:moyenne des débits maxima journaliers de chaque
Fuller
année; calculée d'après les données disponibles (il
s'agit en quelque sorte d'une crue "inter-annuelle").
𝑄𝑟 = 𝐾 ∙ 𝑃 ∙ 𝐴0.58 ∙ 𝐼 0.42
Qr : débit de pointeen 𝑚3 /𝑠 de période de

retour T années;
𝑘 : coefficient dépendant de la nature du
Mac Math
bassin versant (couvert végétal, topographie);
𝑃 : pluie maximale sur 24 heures en 𝑚 de
période deretour T années;
𝐴 : superficie du bassin versant en 𝑘𝑚2 ;
I : pente moyenne du bassin versant.
40
II.5. Transformation pluie-débit

Prévoir le futur aide à mieux gérer le présent. Ceci est vrai aussi dans le domaine de
l’eau où le gestionnaire d’un système hydrique est intéressé à connaître les débits futurs qui
alimentent les réservoirs. Comme les données de base sont les variables météorologiques, il
faut connaître les relations entre ces derrières et le débit dans la rivière. Ces transformations
permettent de relier les précipitations enregistrées ou prévues, aux débits d’eau observés ou
prévus à certains endroits stratégiques dans une rivière ou dans un réseau d’égout.
(Bennis,2007)
Les relations pluie-débit évoluent lors des différentes phases de l’averse :
 L’imbibition est caractérisée par une infiltration qui dépend de la nature et de

l’occupation des sols, du relief, du degré de saturation. La vitesse d’infiltration d’une
pluie est obtenue en se référant au coefficient 𝐾 de la formule de Darcy.
 La phase transitoire correspond à la constitution d’un stock d’eau, qui peut être
important sur les terrains plats 𝑝𝑒𝑛𝑡𝑒 < 1 % , avant d’obtenir la mise en pression
nécessaire au ruissellement.
 Le régime permanent est obtenu par une intensité, un ruissellement et un écoulement

constants.
 La vidange intervient à la fin de l’averse par un prolongement dans le temps des

apports d’eau décroissants. Plus le parcours dans le bassin versant ne sera long et son
relief faible de vidange sera long. (Satin et al,1999)
II.5.1. Fonction de production

La fonction de production permet de définir l'intensité des pluies efficaces et donc du
coefficient d'écoulement à l'appliquer aux intensités de pluies brutes à chaque pas de temps
de calcul.(Roche et al, 2012)
Parmi l'ensemble des modèles existants, on peut retenir deux grandes approches, à savoir:
 Modèles basés sur des relations empiriques;
 Modèles à base physique.
41
II.5.1.1. Modèles empiriques
Les relations empiriques expriment une décroissance de l'infiltration en fonction du

temps à partir d'une valeur initial, on citons :
II.5.1.1.1. Le modèles de horton
𝑖 𝑡 = 𝑖𝑓 + 𝑖0 − 𝑖𝑓 . 𝑒 −𝛾.𝑡 II. 15
Où:
𝑖 𝑡 : Capacité d’infiltration au temps t [mm/h];
𝑖0 , 𝑖𝑓 : Capacité d’infiltration respectivement initiale et finale, dépendant surtout du type

de sol [mm/h];
t : Temps écoulé depuis le début de l’averse [h];
𝛾 : Constante empirique, fonction de la nature du sol [𝑚𝑖𝑛−1 ].
L'utilisation de ce type d'équations, quoique répandue, reste limitée, car la détermination des
paramètres 𝑖0 , 𝑖𝑓 , 𝛾 et présente certaines difficultés. Ces paramètres peuvent cependant être
estimés entre autre par régression à partir de données d'infiltration cumulative.
Le paramètre 𝑖0 dépend du type de sol et des conditions initiales d’humidité le paramètre

𝑖f est une caractéristique constante d’un type de sol donné et le paramètre y dépend
également de la direction de l'écoulement.
Les tableaux ‎II-6 et II-7 précisent quelques valeurs des paramètres des formules de Horton
pour différents types de sols ainsi que pour les quatre classes de sols définies par le Soil
Conservation service américain. (Musy, 1998)
Tableau ‎II-6: Paramètres de l’équation de Horton pour différents types de sol (Musy, 1998)
Type de sol iO [mm/h] if [mm/h] 𝛾 [min−1 ]
42
Standard nu 280 6-220 1.6

Agricole avec tourbe 900 20 - 290 0.8
Tourbe 325 2-20 1.8
Sable fin nu 210 2-25 2.0
Argile avec tourbe 670 10-30 1.4
Tableau ‎II-7: Capacités initiale 𝒊𝟎 , et finale 𝒊𝒇 d’infiltration en fonction des classes de sols définies
par le SCS (Musy,1998)
Type de sol Groupe 𝑖0 , [mm/h] if [mm/h]
Sable, silt A 250 12-8
Limon sableux B 200 8-4
Limon argileux C 130 4-1
Argiles, sols salins D 75 1-0
Figure ‎II-5: Distribution de la Pluit nette par la méthode de horton.(Musy,1998)
La méthode de Horton fait l'hypothèse que le ruissellement de surface apparaît dès que
l'intensité de précipitation dépasse la capacité d'infiltration du sol. Cette capacité vari au
cours du temps, comme l'illustre la (Figure ‎II-5).
43
II.5.1.1.2. Méthode de 𝝋 constant
L'indice 𝜑 est défini comme une capacité moyenne d'infiltration constante dans le
temps pour une averse donnée, au-dessus de laquelle tout excédent de précipitation se
retrouve sous forme de ruissellement de surface à l'exutoire.
L'indice est obtenu en traçant une parallèle à l'axe du temps de telle sorte que la surface du
hyétogramme située au-delà de cette droite représente la hauteur de lame ruisselée, c'est-à-
dire la hauteur de pluie nette (Figure ‎IIII-7).
Cet indice d'infiltration intègre de manière très simple les divers processus d'interception, de
rétention de surface, d'évaporation et d'infiltration. Pour de pluies importantes (événements
rares), l'interception et le stockage de surface peuvent être négligés par rapport à la
précipitation totale.
La variation de la capacité d'infiltration durant une averse restant difficile à connaître,

l'hydrologue est fréquemment tenté d’utiliser cette méthode. Cependant, l'hypothèse de
pertes par infiltration constantes dans le temps peut conduire à des sous-estimations des
volumes ruisselés, en particulier pour des averses de longue durée où le sol peut atteindre un
état de saturation important.(Musy, 1998)
Figure ‎II-6:Méthode de l'indice φ.(Musy,1998)

II.5.1.1.3. Méthode du SCS-CN
Apparue vers 1950, la méthode dite du SCS est le fruit de plus de vingt années
d'analyses des relations pluies-débits sur des bassins versants ruraux de petite taille. Cette
méthode fut développée par trois consultants privés, Homer, Horton et Sherman qui
44
analysèrent les résultats obtenus à l'aide d’infiltromètres. Sherman sera alors le premier à
proposer une représentation graphique du ruissellement en fonction de la pluie.
Cette représentation permet d'effectuer le constat que les points représentant les volumes
ruisselés se situent au-dessous d'une droite de pente égale à l'unité. A partir de cette
constatation, Mockus suggéra que le ruissellement peut en être déduit si l'on prend en
considération le type de sol et sa couverture, les précipitations antécédentes, la durée de la
précipitation et la période de l'année durant laquelle l'événement arrive. Fort de ces
remarques, Mockus proposa alors une relation faisant intervenir ces divers paramètres sous
la forme d'indices. (Musy,1998)
Développements de la méthode
On définit tout d'abord les grandeurs suivantes:
P : Précipitation totale (Pluie brute);
Q : Pluie nette (Précipitation participant au ruissellement);
P - Q : Pertes additionnelles. Ce sont les précipitations infiltrées après le début du

ruissellement;
S: Pertes maximales potentielles;
𝐼𝑎 : Pertes initiales.
On a de suite la situation suivante (Figure ‎II-7):
On fait ensuite l'hypothèse du SCS qui consiste à admettre que la somme des courbes
représentant les précipitations cumulées et celle représentant l'infiltration cumulative est
constante. Les deux courbes sont donc symétriques par rapport à un axe parallèle à l'axe des
abscisses passant par le point d'intersection o sur la partie de droite de la figure ci-dessus:
45
Figure ‎II-7:Méthode du SCS.(Musy,1998)
En appliquant les règles de la géométrie élémentaire relatives aux triangles semblables issus
de cette figure, on peut écrire simplement:
𝑄 𝑃−𝑄
= II. 16
𝑃 𝑆
Et le ruissellement 𝑄s'exprime par :
𝑃2
𝑄= II. 17
𝑃+𝑠
En introduisant les pertes initiales 𝐼𝑎 ,le terme P devient 𝑃 − 𝐼𝑎 et l'éqaution précédente

s'écrit :
𝑃 − 𝐼𝑎
𝑄= II. 18
𝑃 − 𝐼𝑎 + 𝑆
Après avoir effectué de nombreuses expériences, le SCS a proposé une relation empirique
supplémentaire liant les pertes initiales aux pertes maximales potentielles:
𝐼𝑎 = 0.2 ∙ 𝑆 II. 19
Et l'expression de la pluie nette s'écrit:
(𝑃 − 0,2 ∙ 𝑆)2
𝑄= II. 20
𝑃 + 0.8 ∙ 𝑆
Afin de standardiser les courbes représentant la pluie nette en fonction de la pluie brute, le
SCS a introduit un paramètre sans dimension se nommant "Curve Number" (CN). Ce
46
paramètre est normé tel que: 0 < CN <100 et rapporté en unités métriques, il s'exprime dès
lors par:
25400
𝐶𝑁 = II. 21
𝑆 + 254
Compte tenu des expériences réalisées, il est possible de définir plusieurs paramètres CN
suivant les conditions antécédentes d'humidité. On a ainsi:
25 400
 Conditions normales: 𝐶𝑁 𝐼𝐼 = 𝑆+254
 Condition sèches:𝐶𝑁(𝐼) = 𝑓𝑐𝑡 ((𝐼𝐼))
 Condition humides: 𝐶𝑁(𝐼𝐼𝐼) = 𝑓𝑐𝑡 ((𝐼𝐼))
En fin, le Soil conservation service a établi des relations entre les valeurs du CurveNumber
en fonction du type et de l'utilisation de différents sols. Tableau ‎II-8 donne quelques
exemples pour des conditions normales d'humidité (CNII) et pour des pertes initiales
𝐼𝑎 = 0.2 ∙ 𝑆 II. 20
47
Tableau ‎II-8:CN pour des conditions normales d’humidité du sol (Chow, 1988)
classes de sols
Occupation du sol
A B C D
Sols cultivés bien entretenus 72 81 88 91
Sols cultivés mal entretenus 62 71 78 81
Pâturage dans de bonnes conditions 68 79 86 89
Pâturage dans de mauvaises conditions 39 61 74 80
Prés 30 58 71 78
Forêts à faible taux de couverture 30 58 71 78
Forêts à taux de couverture élevé 25 55 70 77
Place de jeux, parcs, golf, cimetières:
- plus 75 % de la surface est enherbée 39 61 74 80
- 50 - 75 % de la surface est enherbée 49 69 79 84
Zones commerciales (85 % de surfaces 89 92 94 95
imperméables)
Zones industrielles (72 % de surfaces 81 88 91 92
imperméables)
Zones résidentielles:
77 85 90 91
Taux de surface imperméables: 65%
61 75 83 87
38%
57 72 81 86
30%
54 70 80 85
25%
51 68 79 84
20%
Parkings, toits etc.. 98 98 98 98
Routes et chemins: Revêtement bitumineux 98 98 98 98
Revêtement gravelé 76 85 89 91
Revêtement en mauvais état 71 82 87 89
II.5.2. Méthodes à base physique :

Ces modèles dérivent d'une manière simplifiée le mouvement de l'eau dans le sol, en
particulier au niveau du font d'humidification, en fonction de certains paramètres
physique.(Musy, 1998)
II.5.2.1. Modèle de Philip
Philip a proposé une méthode de résolution de l'équation de l'infiltration verticale pour

certaines conditions initiales et limites. L'infiltrabilité, exprimée sous la forme de séries de
puissances, peut être simplifiée comme suit:
1 II. 22
𝑖 𝑡 = ∙ 𝑆 ∙ 𝑡 −0,5 + 𝐴𝑡
2
48
Où: 𝑖 𝑡 : Taux d'infiltration [cm-s-1];
𝑆 :Sorptivité [cm-s-1/2],
t : Temps [s];
𝐴𝑡 : Composante gravitaire, fonction de la conductivité hydraulique à saturation [cm/s].
La sorptivité représente la capacité d'un sol à absorber l'eau lorsque l'écoulement se produit
uniquement sous l'action du gradient de pression.
Pour l'infiltration horizontale, le gradient gravitationnel n'est pas impliqué. L'infiltrabilité

s'exprime donc ainsi, selon Philip:
1
𝑖 𝑡 = ∙ 𝑆 ∙ 𝑡 −0.5 II. 23
2
La sorptivité est alors définie par la lame infiltrée I en écoulement horizontal:
𝐼
𝑆= II. 24
𝑖
La sorptivité dépend des conditions initiales et des conditions aux limites du système. Elle
est fonction des teneurs en eau initiale du sol 𝜃𝑖 , et imposée ensurface 𝜃0 .(Musy, 1998;
Roche et al, 2012)
II.5.2.2. Modèle Grène et Ampt
Ce modèle repose sur des hypothèses simplificatrices qui impliquent une schématisation
du processus d'infiltration selon la loi de Darcy (figure II-8).Il s'avère satisfaisant dans le cas
de son application à un sol dont la texture est grossière.
Figure ‎II-8: Schématisation du processus infiltration selon Grenne et Ampt. (Musy, 1998)
49
Les hypothèses principales de la méthode sont:
 Un front d'humidification parfaitement défini;
 Une zone de transmission dont la teneur en eau et la conductivité hydraulique sont

constantes dans le temps et l'espace;
 Les forces de succion au front d'humidification restent constantes quel soit le temps
et la position considérée.
Basé sur la loi de Darcy, le modèle inclut les paramètres hydrodynamiques du sol:
𝜕𝐼 𝐻0 − 𝐻𝑓 𝑡 𝑕0 − 𝑕𝑓 − 𝑍𝑓 𝑡
𝐼 𝑡 = = 𝐾𝑠 ∙ = 𝐾𝑠 ∙ II. 25
𝜕𝑡 𝑧𝑓 𝑡 𝑍𝑓 𝑡
Avec:
𝑖 𝑡 : Taux d'infiltration [mm/h];
I: Infiltration cumulative [mm]; temps [h];
𝐾𝑠 : Conductivité hydraulique à saturation [mm/h];
𝐻0 : Charge hydraulique totale en surface [mm];
𝐻𝑓 𝑡 : Charge hydraulique totale au niveau du front d'humidification [mm];
𝑍𝑓 𝑡 : Profondeur atteinte par le front d'humidification (autrement dit hauteur d’eau

infiltrée depuis le début de l’alimentation) [mm];
𝑕0 : Charge de pression en surface [mm];
𝑕𝑓 : Charge de pression au front d'humidification [mm].
Une des hypothèses du modèle de Green et Ampt stipule que la teneur en eau de la zone de
transmission est uniforme. L'infiltration cumulative I(t) résulte alors du produit de la
variation de teneur en eau et de la profondeur du front d'humidification.
Où:
𝐼 𝑡 = 𝜃0 − 𝜃𝑖 ∙ 𝑍𝑓 𝑡 II. 26
50
I(t) : infiltration cumulative;
𝜃0 : Teneur en eau imposée en surface;
𝜃𝑖 : Teneur en eau initiale du sol;
𝑍𝑓 𝑡 : Profondeur atteinte par le front d'humidification.
En combinant ces expressions, on obtient:
𝑕0 − 𝑕𝑓 − 𝑍𝑓 𝑡
𝐼 𝑡 = − 𝜃0 − 𝜃𝑖 ∙ 𝑍𝑓 𝑡 = 𝐾𝑠 ∙ II. 27
𝑍𝑓 𝑡
Pour un écoulement horizontal, l'infiltrabilité varie de façon inverse par rapport à

l'infiltration cumulative:
𝑕0 − 𝑕𝑓 ∙ 𝜃0 − 𝜃𝑖 II. 28
𝑖 𝑡 = 𝐾𝑠 ∙
𝐼 𝑡
Pour un écoulement vertical, l'infiltrabilité devient alors par approximation, pour des valeurs
de t élevées:
𝑖(𝑡) = −𝐾𝑠 II. 29

La méthode de Green et Ampt s'avère particulièrement satisfaisante dans le cas d'infiltration
dans un sol sec à texture grossière.
Cette méthode reste cependant empirique puisqu'elle nécessite la détermination

expérimentale de la valeur de la charge de pression au front d'humidification. En pratique,
cette valeur est obtenue empiriquement d'après la teneur en eau du sol.(Musy, 1998)
51
Le tableau ci-dessous résume les principales fonctions de transformation pluie-débit:
Tableau ‎II-9: Les principales fonctions de transformation pluie-débit(Hingray et al,2009)

Modèle Auteur Fonction Paramètre
𝑕0 et 𝑕𝑓 charges de pression en
Green - 𝑕0 − 𝑕𝑓 surface et au front d'humidité;
𝑖 𝑡 = 𝐾𝑠 1 +
Ampt 𝑍 𝑡 𝑍 𝑡 : profondeur du d'humectation
A base atteinte au temps t .
1
physique 𝑆 : Sorptivité [𝑐𝑚 − 𝑠 − 2 ],
1 A: composante gravitaire, fonction

Philip 𝑖 𝑡 = ∙ 𝑆 ∙ 𝑡 −0,5 + 𝐴𝑡
2 de la conductivité hydraulique à
saturation [𝑐𝑚 ∙ 𝑠 −1 ].
𝑖0 , 𝑖𝑓 : capacité d’infiltration
respectivement initiale et finale,
𝑖 𝑡 dépendant surtout du type de sol
Horton
= 𝑖𝑓 + 𝑖0 − 𝑖𝑓 . 𝑒 −𝛾.𝑡 [𝑚𝑚/𝑕];
𝛾 : constante empirique, fonction de
Empirique la nature du sol [𝑚𝑖𝑛−1 ],

(𝑃 − 0,2 ∙ 𝑆)2
𝑄=
𝑃 + 0.8 ∙ 𝑆
P : Précipitation totale (Pluie brute);
SCS- AVEC
S : Pertes maximales potentielles.
CN 25400
𝐶𝑁 =
𝑆 + 254
II.5.3. Fonction de transfert

Les pluies nettes ou efficaces étant calculées à chaque pas de temps, il reste à tenir
compte de leurs délais de transfert au sein du réseau hydrographique du bassin versant.
(Roche et al,2012)
Dans le suivant on va présenter les concepts est les modèles utilisés pour le transfert des
ruissellements directs issus de la pluie nette (Hingray et al,2009)
52
II.5.3.1. Fonction conceptuelle
II.5.3.1.1. Le modèle à réservoir linéaire
Le modèle à réservoir linéaire est de loin le plus utilisé pour modéliser la fonction de
transfert. Le comportement du réservoir est régi par seulement deux équations: l’équation de
continuité et l’équation de vidange. La linéarité réside dans l’équation de vidange, le débit
sortant étant supposé une fonction linéaire de la hauteur d’eau (ou du volume d’eau) stocké
dans le réservoir:
𝑑𝑆 𝑡 II. 30
Equation de continuité : =𝑖 𝑡 −𝑞 𝑡
𝑑𝑡
Equation de vidange: 𝑞 𝑡 =𝑆 𝑡 𝐾 II. 31
Où
𝑆 𝑡 : Le volume stocké au temps t dans le réservoir;
𝑖 𝑡 et 𝑞(𝑡) sont respectivement le flux rentrant et le flux sortant du réservoir, exprimés

dans la même unité;
𝐾: est la constante de temps du modèle, le seul paramètre du modèle à estimer.
Si le stock représente la hauteur d’eau stockée sur le bassin (c’est-à-dire la lame d’eau),
les un de ces différentes variables sont souvent en [𝑚𝑚]pour𝑆, en [𝑚𝑚/𝑕] pour les flux i
et𝑞 en [𝑕] pour 𝐾. Les flux ont alors la dimension d’une vitesse – ce sont exemple l’intensité
des précipitations pour i, le débit spécifique pour𝑞. Si le stock représente le volume d’eau
stocké au sein du bassin, les unités sont souvent en [𝑚3 ] pour𝑆, en [𝑚3 /𝑠]pour les flux, 𝑖 et
𝑞 et en [𝑠]pour 𝐾.
En combinant les équations de continuité et de vidange du réservoir, une équation

différentielle linéaire du premier ordre est obtenue. Elle constitue l'équation d'un réservoir
linéaire :( Hingray et al,2009)
𝑑𝑞 (II. 32)
𝐾. +𝑞 =𝑖
𝑞𝑡
53
II.5.3.2. Hydrogramme Unitaire
Proposée vers 1932 par J.K. Shermann (1932), la théorie de l'hydrogramme unitaire
permet de déterminer un hydrogramme de crue lorsque seules les précipitations sont
connues. C'est une contribution importante à l'étude du phénomène de l'écoulement de
surface car elle sert de base aux fonctions de transfert les plus couramment utilisées.
L'hydrogramme unitaire reste cependant un modèle globale applicable exclusivement pour
l'agrégation et la discrétisation en sous-bassins.(Musy, 1998)
Hypothèse
Dans la théorie de l’HU, le terme unitaire désigne la durée de référence 𝜏 de l’averse et

son volume normé (1 mm, 1 pouce etc…). Des variations de 𝜏 ou du volume conduisent à
un HU différent.
La théorie de l'HU repose sur les hypothèses suivantes:
 La pluie nette est uniforme sur le bassin versant et constante dans le temps;
 Les caractéristiques de l'HU sont invariantes dans le temps : la durée de

l'hydrogramme de ruissellement direct (temps de base𝑡𝑏 ) est constante et
indépendante de l'intensité de la pluie;
 La réponse du bassin versant aux pluies est linéaire: pour une même durée de
référence 𝜏, les ordonnées de l'HU sont proportionnelles au volume de la pluie nette
(Figure II-9);
 L’hydrogramme unitaire reflète l’ensemble des caractéristiques physiques du bassin

versant.
Figure ‎II-9:Hypothèse de linéarité de l'Hydrogramme Unitaire (Musy, 1998)
54
La linéarité et l’invariance dans le temps sont les deux principes fondamentaux de la théorie
de l’HU. Le principe de linéarité est illustré par la figure ci-dessus. L’Hydrogramme Unitaire
normé provient d’une averse de volume normé 1 et d’une durée de référence𝜏. Pour une
même durée de l'averse de référence, si le volume de pluie nette passe de 1 à 2, alors
l’hydrogramme de ruissellement voit ses ordonnées simplement multipliées par 2 (principe
de la linéarité). (Musy, 1998)
Dans le jargon hydrologique, les trois fonctions de réponse d'un bassin linéaire et invariant
sont les suivant: (Hingray, 2009)
 La réponse unitaire
 La réponse impulsionnelle
 La réponse indicielle
II.5.3.2.1. Hydrogramme Unitaire (réponse unitaire)
Définition
L'Hydrogramme Unitaire d'un bassin versant pour une durée de référence 𝜏 est la
réponse unitaire𝑢𝜏 𝑡 du bassin pour cette durée vis-à-vis de la transformation pluie nette-
débits. En d'autres termes, l'Hydrogramme Unitaire (HU) d'un bassin versant est
l'hydrogramme de ruissellement direct résultant d'une pluie nette de hauteur h, de durée 𝜏,
uniformément répartie dans l'espace et l'intensité constante sur la durée𝜏. (Hingrayet al,
2009)
Figure ‎II-10:Définitions des caractéristiques de l'Hydrogramme Unitaire. (Hingray et al, 2009)
55
𝑡𝑝 : Le temps de montée (le temps entre le début et la pointe de l'hydrogramme);

𝑡𝑑 : Le temps de décrue (le temps entre la pointeet la fin de l'hydrogramme);
𝑡𝑐 : Le temps de concentration (le temps entre la fin de la pluieet la fin de
l'hydrogramme);
𝑡𝑙 : Le temps de réponse (le temps entre le centre de gravité de l'averseet la pointe de
l'hydrogramme);
𝑡𝑏 : Le temps de base de l'hydrogramme 𝑡𝑏 = 𝑡𝑝 + 𝑡𝑑 = 𝜏 + tc (Figure II-11) (la
durée totale de l'hydrogramme). (Hingray; 2009 ; Musy, 1998)
Pour une durée de référence 𝜏 est une hauteur unitaire h fixées, les caractéristiques de l'HU
correspondant sont invariantes dans le temps.
Le terme unitaire désigne ici la durée de référence𝜏 de l'averse. Des variations de durée de
référence 𝜏 ou de hauteur h conduisent à un HU différent. En particulier, le volume de
ruissellement direct dépend de ces variables. (Hingray, 2009)
II.5.3.2.2. Méthodes d'identification de l'Hydrogramme Unitaire
Le cas général : inversion matricielle
Différents méthodes sont possibles pour l'identification de l'Hydrogramme Unitaire.

Elles font appel à la technique de convolution. Si la pluie nette et l'hydrogramme de
ruissellement sont décrits respectivement par M intensités moyennes 𝑖𝑚 et N valeurs de
débits 𝑄𝑛 N équations peuvent être alors formulées pour 𝑄𝑛 , 𝑛 = 1, … … , 𝑁.Elles font
apparaitre les N – M + 1 inconnues 𝑢𝑘 qui définissent l'HU (noter que comme𝑄𝑁+1 = 0
alors les termes 𝑢𝑁−𝑀+2 = 𝑢𝑁−𝑀+3 = 𝑢𝑁−𝑀+4 = 𝑢𝑁 = 𝑢𝑁+1 sont nuls). Ce système
d'équations peut s'écrire sous forme matricielle de façon suivante :
𝑖1 0 0 ⋯ 0 0 ⋯ 0 0 𝑄1
𝑖2 𝑖1 0 ⋯ 0 0 ⋯ 0 0 𝑄2
𝑖3 𝑖2 𝑖1 ⋯ 0 0 ⋯ 0 0 𝑢1 𝑄3
. . . ⋯ . . ⋯ . . 𝑢2 .
. ⋯ ⋯ 𝑈2 .
. ⋯ ⋯ . .
𝑖𝑀 𝑖𝑀−1 𝑖𝑀−2 ⋯ 𝑖1 0 ⋯ 0 0 ∙ . = 𝑄𝑀
𝑖1 𝑖𝑀 𝑖𝑀−1 ⋯ 𝑖1 𝑖1 ⋯ 0 0 . 𝑄𝑀+1
. ⋯ ⋯ . . .
. ⋯ ⋯ . . .
. ⋯ ⋯ . 𝑢𝑁−𝑀+1 .
0 0 0 ⋯ 0 0 ⋯ 𝑖𝑀 𝑖𝑀−1 𝑄𝑁+1
0 0 0 ⋯ 0 0 ⋯ 0 𝑖𝑀 𝑄𝑁
56
Ou encore sous forme synthétique :
𝑖 ∙ 𝑢 = 𝑄
𝑇
Avec 𝑢 = 𝑢1 , 𝑢2 , … 𝑢𝑁−𝑀+1 le vecteur des 𝑁– 𝑀 + 1 valeurs définissant
l'Hydrogramme Unitaire cherché.
En général, le système d'équations est surdéterminé car le nombre d'équations (𝑁)est

plus grand que le nombre d'inconnues (𝑁 – 𝑀 + 1).
Soit 𝑢 uneestimation de l'HU donnant l'estimation 𝑄 de l'hydrogramme observé 𝑄 par

convolution avec la pluie nette 𝑖 = 𝑖1 , 𝑖2 , … 𝑖𝑀 . L'identification de l'HU consiste alors à
trouver l'estimation 𝑢 qui minimise l'erreur 𝜀entre l'hydrogramme simulé 𝑄 et
l'hydrogramme observé 𝑄 . Le système étant généralement surdéterminé, il est pour cela
possible d'utiliser les méthodes classiques d'optimisation numérique. Ces méthodes
n'assurent cependant pas que toutes les ordonnées de l'HU cherché soient positives.
(Hingray et al, 2009)
II.5.3.2.3. L'Hydrogramme unitaire instantané (réponse impulsionnelle)
L'Hydrogramme Unitaire Instantané (HUI) est la réponse impulsionnelle 𝑢0 𝑡 du

bassin versant vis-à-vis de la transformation pluie nette-débits.
Figure ‎II-11: La réponse impulsionnelle.(Hingray, 2009)
En d'autres termes, l'HUI 𝑢0 𝑡 du bassin versant est l'HU pour une pluie unitaire de duré
infinitésimale. L'HUI 𝑢0 𝑡 est la limite de l'HU 𝑢𝜏 𝑡 lorsque la durée de l'averse unitaire 𝜏
tend vers 0. L'HUI est aussi la dérivée première de la courbe en S :
57
𝑑𝑠 𝑡
𝑢0 𝑡 = lim 𝑢𝜏 𝑡 = (II. 33)
𝜏→0 𝑑𝑡
II.5.3.2.4. La courbe en S (réponse indicielle)
L'hydrogramme en S est l'hydrogramme de ruissellement d'une pluie continue et d'une

durée supérieure à 𝑡𝑐 . Cet hydrogramme prend la forme d'une courbe en S qui devient
constante. En fait, au-delà de 𝑡𝑐 , la contribution du bassin versant devient maximale puisque
toute la surface contribue à l'écoulement. Ceci n'est vrai que dans le cas d'une averse
uniforme.
Pratiquement cet hydrogramme est obtenu par convolution d'une averse continue (Figure II-
12).
Figure ‎II-12:Construction de la courbe en S (Hingray et al, 2009)
La courbe en S permet de calculer un hydrogramme unitaire de norme 1 et de durée de

référence𝑇𝜏 , à partir de l'HU correspondant à 𝜏. (Musy, 1998)
II.5.3.3. Détermination de l'Hydrogramme Unitaire
II.5.3.3.1. Cascade de Nash
La résolution successive de l’équation du réservoir linéaire pour chacun des réservoirs

de la cascade de Nash permet d’obtenir l’expression de l’hydrogramme unitaire instantané du
modèle de Nash. Soit K le coefficient de stockage des réservoirs [𝑕] et n le nombre de
réservoirs. La sortie du premier réservoir pour une sollicitation unitaire instantanée est l’HUI
du réservoir linéaire donné par l’équation (34):
1 −𝜏
(II. 34)
𝑢1 𝜏 = ∙e 𝐾
𝐾
58
Lorsque cette fonction est utilisée comme entrée du second réservoir, le débit sortant
du second réservoir, obtenu par convolution de cette entrée et de l’HUI du second réservoir
(égal aussi à𝑢(𝑡)),a pour expression :
𝑡 𝑡
1 −1
1−𝜏 (II. 35)
𝑢2 𝜏 = 𝑖 𝜏 . 𝑢 𝑡 − 𝜏 . 𝑑𝜏 = 𝑢1 𝑡 ∙ ∙e 𝐾 ∙ 𝑑𝜏
𝐾
0 0
1 −𝑡
(II. 36)
𝑢2 𝑡 = ∙e 𝐾
𝐾2
En réitérant le procédé de calcul pour chacun des n réservoirs, on obtient

immédiatement l'expression du débit sortant du 𝑛è𝑚e réservoir, c'est-à-dire, l'HUI du
modèle de Nash :
𝑛−1
1 𝑡 −𝑡 (II. 37)
𝑢𝑛 𝑡 = ∙ ∙e 𝐾
𝐾. Ґ 𝑛 𝐾
Où Ґ 𝑛 est la fonction gamma définie par :
𝑛−1 ! ∀𝑛 ∈ 𝑅
Ґ 𝑛 =
∞ (II. 38)
𝑢𝑛−1 ∙ e−𝑢 ∙ 𝑑𝑢 ∀ 𝑛 ∈ 𝑁
0
II.5.4. L'Hydrogramme Unitaire Synthétique Régional de SCS

Le Soil Conservation Service américain (SCS) a effectué la synthèse d'un grand nombre
d'Hus moyens identifies sur des bassins versants de différentes régions des Etats-Unis.
Figure ‎II-13: Hydrogramme unitaire synthétique sans dimension proposé par le SCS, et son
approximation triangulaire.(Hingray et al, 2009)
59
Les HUS ont été décrits à cet effet dans un système de coordonnées adimensionnelles où les
abscisses sont exprimées par le rapport 𝑡/𝑡𝑝 et les ordonnées par 𝑞(𝑡)/𝑞𝑝 (avec 𝑡𝑝 temps de
montée et 𝑞𝑝 le débit de pointe de l’HU). L’HU synthétique régional résultant de l’analyse et
la courbe en 𝑆 correspondantes sont donnés à la (Figure ‎II-13).
Pour simplifier l’utilisation de l’HUS, le 𝑆𝐶𝑆 en a proposé une approximation triangulaire

qui dépend de 2 seuls paramètres :𝑞𝑝 et𝑡𝑝 . Le SCS propose pour un bassin donné de les
estimer par les équations régionales suivantes :
𝐷 (II. 39)
𝑡𝑝 = + 0,6 ∙ 𝑡𝑐
2
1,5𝑉𝑅 (II. 40)
𝑞𝑝 =
𝐷 + 1,2𝑡𝑐
Avec :
𝐷 : Ladurée de référence de la pluie unitaire à l’origine de l’hydrogramme
𝑉𝑅 : Le volume ruisselé
𝑡𝑐 : Le temps du concentration de bassin (défini ici comme le temps s’écoulant entre la

fin de l’averse et le point d’inflexion de la courbe de décrue).
Ces deux équations proviennent de la combinaison entre différentes équations issues de

l’expérience et diverses relations qui relient les variables décrivant l’HUS triangulaire :
𝑞𝑝 (II. 41)
Expression du volume de la crue : 𝑉𝑅 = 𝑡𝑝 +
2
Expression du temps de montée: 𝐷 (II. 42)
𝑡𝑝 = + 𝑡𝑙
2
Avec:
𝑡𝑟 :Durée de la décrue;
𝑡𝑙 : Temps entre le milieu de l’averse et la pointe de l’hydrogramme.
Le SCS estime qu’aux Etats-Unis vaut approximativement 0,6 ∙ 𝑡𝐶 et que37.5% de la

crue se situe avant la pointe de l’hydrogramme (par suite 𝑡𝑝 𝑡𝑏 = 0.375 et 𝑡𝑏 = 2.67 ∙
𝑡𝑝 𝑜ù 𝑡𝑝 le temps de base de l’hydrogramme). (Hingray et al, 2009)
60
II.5.4.1. Hydrogramme Géomorphologique
Suivons le cheminement d’une goutte de pluie dans le schéma de la figure 5.26, où un

cercle représente l’ensemble des cours d'eau de l’ordre i (ou un état i du système) selon la
classification de Horton-Strahler.
Le point de départ de la goutte de pluie peut être dans notre exemple un cours d’eau d’ordre
1,2 ou 3. En admettant que ce point de départ est aléatoire, on a les probabilités𝑝01 , 𝑝02 et
𝑝03 que la goutte d’eau commence son cheminement en 1, 2 ou 3. Ces probabilités
𝑝0𝑖 correspondent à l’état initial du système.
Suivons maintenant l’écoulement de la goutte d’eau dans le réseau hydrographique; par

exemple p12 est la probabilité qu’une goutte s’écoule d’un cours d’eau d’ordre 1 à un cours
d’eau d’ordre 2. On définit pij comme étant la probabilité de transition de l’état i à l’état j
avec𝑝𝑖𝑗 =0 lorsque j<i puisque l’eau ne remonte pas dans le réseau (j = i correspond au
maintien de la goutte d’eau dans un cours d’eau du même ordre). (Hingray et al, 2009)
Figure ‎II-14 : Schématisation d'un bassin versant du 3𝑒 ordre et probabilités de transition d'un
cours d'eau d'ordre 𝑖à un cours d'eau d'ordre 𝑗 > 𝑖 (Hingray, 2009)
II.6. Conclusion
Beaucoup de méthode des calculs hydrologiques utilisées en science de l'ingénieur se
basent sur des pluies de projet plus au moins complexes. Dans ce chapitre, nous avons
présenté une modélisation hydrologique de phénomène de l'inondation, en s’appuyant sur
les méthode statistiques les plus couramment utilisées pour l'estimation du débit des crues
dans un période de retour donné, ainsi que les différentes méthodes empiriques et
rationnelles.
61
III. Numéro du chapitre
— Albert Einstein
L e bassin versant est un objet complexe dont l'ensemble des

caractéristiques (géométrique, géologiques, hydrographique, etc)
joueront un rôle non seulement dans la réponse hydrologique du
bassin à une sollicitation des précipitations (régime des écoulements) mais
aussi, en amont et pour certaines d'entre elles ( altitude, exposition …),
directement dans le processus de formation de la pluie. La méthodologie
élaborée pour la présente étude est basée sur deux partie une présentation de
notre zone d'étude et celle d'une présentation de la morphologie du bassin
versant à fin de tiré les différents caractéristiques géométriques du bassin
versant ainsi que du cours d'eau.
III.1. Présentation de la zone d'étude

III.1.1. Présentation de la Wilaya
III.1.1.1. La situation
La région de Jijel fait partie du Sahel littoral de l'Algérie, elle se trouve sur la côte Est
algérienne, entre deux grands ports : celui de Bejaia (90 Km à l’ouest) et celui de Skikda (140
Km à l’Est). Elle est située au Nord-Est entre les latitudes 36° 10 et 36° 50 Nord et les
longitudes 5° 25 et 6° 30 Est. Les bassins versants de la région font partie du grand bassin
versant de l'Algérie (côtiers Constantinois), La wilaya de Jijel s'étant sur une superficie de
2.398,69 Km2 dont 82% sont des montagnes. (Fatih, 2009 ; Boudjoudjou, 2010).
62
Chapitre III — Zone d'étude et caractéristique morphométriques
Figure ‎III-1: Situation géographique de la région d’étude (source : google map)
III.1.1.2. Le Climat
Jijel est sous l’influence d’un climat maritime méditerranéen caractérisé par des étés secs
et tempérés, et des hivers doux et humides. Les pluies durent 06 mois totalisant une
pluviométrie de 1200 mm. La température est très douce en hiver (11° en moyenne au mois
de janvier) et la chaleur est tempérée par la brise de mer en été (26° en moyenne au mois
d’Août). (Ouarda et al, 2001)
Les formations végétales ont, pour la plupart, d'indiscutables rapports avec les climats, aussi
bien zonaux que régionaux. Leur étude n'est donc pas inutile. L'analyse du climat se fera sur
plusieurs paramètres, notamment les précipitations et les températures qui représentent les
paramètres les plus importants, de la part leur quantité annuelle, et leur répartition dans
l'espace et dans le temps. (Boudjoudjou, 2010)
III.1.1.2.1. Variations mensuelles et annuelles des précipitations
Les précipitations dans la région d'étude apparaissent d'une façon intense pendant la
période hivernale, qui s'étale de Décembre Février, où les stations de Jijel, Taher et Texenna
reçoivent une quantité non négligeable de pluies variant entre 541 mm et 1407.9 mm
annuellement. Les volumes d'eau reçue s'observent essentiellement en hiver
pourpratiquement les trois stations, au mois de décembre et Novembre, alors que pour
juillet et août, elles enregistrent une forte diminution(Figure ‎III-2). (Boudjoudjou, 2010)
63
Figure ‎III-2: Variations mensuelles des précipitations. Boudjoudjou, 2010)
64
Figure ‎III-3: Régime saisonnier des précipitations. (Boudjoudjou, 2010)

III.1.1.3. Réseau hydrographique
Le réseau hydrographique de la région est très dense. Il est dominé par une direction
Sud-Nord et des affluents de directions différentes favorisent l'écoulement des lames d'eau
précipitées qui déversent généralement dans la mer. Les plus importants oueds sont
(Boudjoudjou, 2010):
 Oued El-Kébir: qui prend naissance de la jonction d'Oued Rhumel et Oued Endja,
traverse El-Milia et El-Ancer et se jette à la mer dans la région de Beni-Belaid.
65
 Oued Djen-Djen qui prend sa source au Babors (Erraguene) est constituée de trois
étages (Boudjoudjou, 2010).
Figure ‎III-4: : Carte hydrographique de la région de Jijel au 1/200 000 (Boudjoudjou, 2010)
III.1.2. Présentation du bassin versant

III.1.2.1. Généralités sur l’oued Djendjen
Le bassin versant d’oued Djendjen se trouve dans la partie Nord-est de l’Algérie, il fait
partie de la wilaya de Jijel. Il étend d'Ouest Est, jusqu’à la station hydrométrique de Missa
(Tablout), puis s’oriente Sud – Nord, jusqu’à la mer. Le réseau hydrographique est dense,
traduisant un drainage élevé. Avant de s’ouvrir sur les plaines littorales, l’oued Djendjen est
très encaissé à l’amont, empruntant des ravins qui entaillent les reliefs de Texenna et Thar
Oussaf. Parmi ses affluents Principaux, citons : les oueds Agoug, Reha et Zatout. L’oued
Djendjen prend naissance dans le djebel Babor. La présence à l’amont de seuils rocheux
particulièrement résistants à l’érosion, accroît la force du ruissellement et induit des débits de
crues instantanées, élevés. La pente de l’oued Djendjen s’adoucit et s’estompe d’avantage en
pénétrant dans la plaine avale où elle devient plus faible. Ceci va provoquer, à ce niveau, un
ralentissement de la vitesse d’écoulement qui aura une incidence sur les relations des eaux
superficielles avec la nappe. (Fatih, 2009).
66
Figure ‎III-5: Bassin versant de l’Oued Djendjen (Fatih, 2009)
III.1.2.2. Situation géographique
Le bassin versant de l’oued Djendjen est situé au Nord-Est de la wilaya de Jijel, il

appartient au grand bassin versant côtiers Constantinois Ouest, Ce bassin versant s’étend du
5°30’ au 5°58’ de longitude et du 36°48’ au 36°22’ de latitude. Il est limité par :
 Au Nord, par la mer méditerranée ;
 Au Sud-Est par le grand bassin versant du Kébir-Rhumel ;
 Au Sud-Ouest par le bassin versant de La Soummam ;
 A l’Est par le bassin versant de Guebli ;
 A Nord-Ouest par le bassin versant des côtiers Algérois. (Derradji et Al, 2013)
III.1.2.3. Situation géographique de la plaine de l’oued Djendjen
La plaine de l’oued Djendjen est située dans la région de Jijel, au Nord Est de l’Algérie.
Elle est distante de 360 Km à l’Est d’Alger et bordée au Nord par la mer Méditerranée.
L’examen de la carte topographique montre que le bassin versant de l’oued Djendjen est
entouré par des monts dont l’altitude varie de 226 mètres à 1992 mètres. A l’arrière du bassin
versant, le massif des Babors, s’élève de façon abrupte, jusqu’à 2006 mètres d’altitude. Cette
chaîne montagneuse présente un relief très accidenté, compartimenté par des gorges
profondes, empruntées par les oueds.(Fatih, 2009)
67
III.1.2.4. La géologie du Bassin versant
La région de Jijel se trouve dans la zone dite : Branche Sud- méditerranéenne de

l’orogénie alpine, appelé "les Maghrébines". C’est la zone des massifs métamorphiques
Kabyles qui font partie des montagnes plissées du littoral méditerranéen (Derradji et Fenazi,
2013)
III.1.2.5. L’Hydroclimatologie du Bassin versant
Le bassin versant de l’oued Djendjen est caractérisé par un climat méditerranéen, avec
un hiver doux à très doux et un été sec.
 Il le reçoit une hauteur de précipitation d’environ 974 mm/an à la station de l’Achouat.

La grêle et la neige se produisent surtout sur les reliefs, au cours de la saison froide;
 Sur la partie littorale du bassin versant la température moyenne annuelle est de 18.74
°C
 L’apport liquide annuel de l’oued est important. Il est de l’ordre de 110 millions de m3.
(Derradji et al, 2013)
III.1.2.6. L’hydrogéologie du bassin versant
En raison de sa forme complexe, son étendue et sa grande variété lithologique, le bassin

versant présente des comportements et des caractéristiques hydrogéologiques complètement
différents à l’amont et à l’aval.(Derradji et al, 2013)
III.1.2.6.1. Systèmes aquifères amont
On peut l’appeler aquifères des roches à perméabilité en grand, ou réservoir à porosité

de fracturation. Il s’agit du réservoir qui constitue toute la partie amont du bassin versant. Il
est constitué par toutes les formations lithologiques (roches vertes, Schistes, calcaires,
dolomies; gypse…etc.).
III.1.2.6.2. Systèmes aquifères aval
On peut aussi l’appeler réservoir à porosité inter granulaire ou d’interstice, plaine ou

nappe alluviale tout simplement.
68
III.1.2.7. Le couvert végétal et l’occupation du sol
III.1.2.7.1. Le couvert végétal
Le couvert végétal (en particulier les forêts et les cultures) vient ajouter son influence à
celle de la nature géologique et topographique du bassin versant pour conditionner la
rapidité du ruissellement superficiel, le taux de l’évaporation et la capacité de rétention du
bassin. (Fatih, 2009)
III.1.2.7.2. Les différents types d’occupation du sol
En se basant sur les données statistiques disponibles au niveau de la direction des

statistiques agricoles et des enquêtes économiques de la willaya de Jijel, nous distinguons
trois domaines essentiels (Tableau III-1). (Fatih, 2009)
Figure ‎III-6 : Carte d'occupation de sol "bassin versant de Oued Djendjen"
Tableau ‎III-1: Le taux d’occupation du sol du bassin versant de l’oued Djendjen.
Domaines Espèces Superficie (ha) Taux (%)

Bois 13300 11.56
Domaines
Chênes lièges 43700 38
forestiers
Maquis 58000 50.43
Surface totale
99024 41
agricole
Domaines
Surface utile
agricoles 43649 44
agricole
Superficie irriguée 6210 14
Terre bâtis Terrains urbains 2110 9
69
III.2. Etude des caractéristiques morphologique du

bassin versant de l'oued DJENJEN
Les caractéristiques morphologiques sont des caractéristiques très importants à étudier,
lors d'une étude hydrologique, L’étude des caractéristiques de notre bassin versant va nous
permettre de souligner les paramètres géomorphologiques qui influence directement le
temps et la forme de réponse du bassin vis à vis d'un événement pluvieux donné.
Figure ‎III-7: Carte des limites de la zone d'étude.
La détermination de ces caractéristiques du bassin versant consiste à calculer, ou évaluer,

Les paramètres explicatifs, indispensables à l’évaluation des débits de crues.
III.2.1. Caractéristiques Géométriques

La superficie d'un bassin versant ainsi que son périmètre joue un rôle primordial dans
les études des bassins versants car elle représente l'aire de réception des précipitations et
d'alimentation des cours d'eau. De ce fait, leurs influences se manifestent dans la nature de la
relation entre le débit et le temps, donc un petit bassin versant réagit très vite à une averse
par rapport à un grand bassin. Pour notre bassin versant, il s’étale sur une superficie de
321Km2et un périmètre de 102 Km. On peut dire alors que la réponse à une averse sera
influencée par son étendue.
La longueur du thalweg principal mesurée est :(Touabia, 2004)
𝐿𝑝 = 44.45 km
70
III.2.2. Caractéristiques de forme

III.2.2.1. Indice de compacité de Gravelius :
L’indice admis par les hydrologues pour caractériser la forme d’un bassin versant d’une
part et de renseigner sur la rapidité de la concentration des écoulements est l’indice de
Gravelieus. (Laborde, 2009)
Le bassin versant est de forme allongée caractérisé par un indice de compacité de Gravillus.
P
K g = 0.28 III. 1
A
donc
Kg=1,59
III.2.2.2. Rectangle équivalent
C'est une transformation purement géométrique en un rectangle de dimensions L et

l'ayant la même surface que le bassin versant. Il permet de comparer les bassins versants
entre eux du point de vue de l’écoulement. Les courbes de niveau sont des droites parallèles
aux petits côtés du rectangle et exutoire est l'un de ces petits côtés (Touabia, 2004).
Le périmètre et la surface du rectangle sont respectivement :
𝑃 = 2(𝐿𝑟 + 𝑙) III. 2
et
𝐴 = 𝐿𝑟 ∙ 𝑙 III. 3
La longueur L et la largeur l en Km sont données par la résolution de P et S :
2
𝑆 1.128
𝐿𝑟 , 𝑙 = 𝐾𝑐 ∙ 1± 1− III. 4
1.128 𝐾𝑐
Les dimensions des rectangles équivalents du bassin versant sont :
𝐿𝑟 =43.65 Km
71
et
l = 7.35 Km
La figure ci-dessous montre la réalisation des dimensions de rectangle équivalent du bassin

versant:
Rectangle équivalent
L=43.65km
L=7,35km
13000m
1100m
1200m
1600m
1000m
1700m
1800m
900m
1500m
2000m
1900m
300m
1400m
Figure ‎III-8: Rectangle équivalant du bassin versant d'oued Djenjen.
III.2.3. Caractéristiques de relief

III.2.3.1.1. La courbe hypsométrique
Le relief est caractérisé par la courbe hypsométrique du bassin. Cette courbe donne la
répartition des altitudes en fonction des surfaces partielles du bassin versant. Les résultats
sont inscrits dans Tableau ‎III-2, qui donne la courbe hypsométrique représentée par la
Figure ‎III-9 (Touabia, 2004)
72
Tableau ‎III-2: La répartition des altitudes en fonction des surfaces partielles.

Altitudes Superficies Superficies Superficies
Altitudes [m]
moyennes [m] Partielles [km2] Partielles [%] Cumulées [%]
300-400 350 6,272 1,95% 1,95%
400-500 450 9,801 3,05% 5,01%
500-600 550 16,878 5,26% 10,27%
600-700 650 43,348 13,50% 23,77%
700-800 750 54,673 17,03% 40,80%
800-900 850 53,603 16,70% 57,50%
900-1000 950 45,399 14,14% 71,64%
1000-1100 1050 35,757 11,14% 82,78%
1100-1200 1150 24,028 7,49% 90,27%
1200-1300 1250 12,065 3,76% 94,03%
1300-1400 1350 9,869 3,07% 97,10%
1400-1500 1450 4,928 1,54% 98,64%
1500-1600 1550 1,826 0,57% 99,20%
1600-1700 1650 0,644 0,20% 99,41%
1700-1800 1750 0,814 0,25% 99,66%
1800-1900 1850 0,725 0,23% 99,88%
1900-2000 1950 0,370 0,12% 100,00%
Altitude(m)
H(m)
100 S(%)
2000
10
1900
1800
1700
1600
1500
1400
H5% 1300
1200
1100
1000
900
H50% 800
700
600
H95% 500
400
300 Surface(%)
100
0.12
2.90
0.34
17.22
5.97
9.73
42.50
59.20
98.05
28.36
89.73
76.23
Figure ‎III-9: Courbe hypsométrique du bassin versant
73
III.2.3.2. Les altitudes caractéristiques
Permet d'analyser les lois réglant les précipitations et le ruissellement superficiel.
Les déférentes altitudes caractéristiques le bassin versant d'oued Djendjen sont :
 L'altitude maximale 𝐻𝑚𝑎𝑥 représente le point le plus élevé du bassin est de 2006 m
tandis que l'altitude minimale considère le point le plus bas est de 237 m
 L’altitude médiane ou 𝐻50% est de 853,20 m qui correspond à l'abscisse 50 % de la

surface totale du bassin sur la courbe hypsométrique.
 L’altitude moyenne 𝐻𝑚𝑜𝑦 est de de 878.09 m
 La dénivelée est la différence d’altitudes (𝐻) comprise entre la surface cumulée, soit:
𝐷 = 836,80 𝑚
Figure ‎III-10: Carte caractéristique des altitudes du bassin versant d'Oued Djendjen
III.2.3.3. Les indices de pente
La pente ne cesse à se présenter comme facteur capital en relation avec le temps de

réponse d’un bassin, elle joue un rôle important dans le ruissellement.
Les différents indices de pentes permettent de caractériser les pentes d’un bassin afin
d’établir des comparaisons et des classifications. (Touabia, 2004)
74
III.2.3.3.1. Indice de pente globale I g
D
Ig = III. 5
Lr
D : Dénivelées spécifiques en m ;
Lr : Longueur du rectangle équivalent en Km.
Ig = 19,171 m/km
III.2.3.3.2. Pente moyenne de cours d'eau principal Ic
Elle exprime le rapport entre la dénivelée et longueur comprise entre 2 points

suffisamment éloignés sur le cours d'eau principal :
∆𝐻
𝐼𝑐 = III. 6
∆𝐿
H : Dénivelée entre 2 points suffisamment distants (m);
III.2.3.3.3. Dénivelée spécifique 𝐷𝑠
Elle dérive de la pente globale Ig en la corrigeant de l'effet de surface admis étant

inversement proportionnel à A
𝐷𝑠 = 𝐼𝑔 𝐴 III. 7
𝐷𝑠 = 343.47
La classification des reliefs du bassin versant selon la dénivelée spécifique montre que le
bassin a des reliefs très forts, le tableau suivant représente les différentes fourchettes de la
classification approuvée.
Tableau ‎III-3: Classification du relief par l'ORSTOM selon Ds
R1 Relief très faible Ds < 10 m

R2 Relief faible 10 m <Ds < 25 m
R3 Relief assez faible 25 m <Ds < 50 m
R4 Relief modéré 50 m <Ds < 100 m
R5 Relief assez fort 100 m <Ds < 250 m
R6 Relief fort 250 m <Ds < 500 m
R7 Relief très fort 500 m <Ds
75
III.2.3.4. Le Degré de développement du réseau
III.2.3.4.1. Densité de drainage 𝐷𝑑
Elle est définie comme étant Le rapport entre la longueur totale de tous les talwegs "Li’’
du bassin versant, à la surface totale "A". Elle reflète la dynamique du bassin, la stabilité du
chevelu hydrographique et le type de ruissellement de surface. Elle est exprimée en
Km/Km².
𝑛
𝑖=1 𝐿𝑖
𝐷𝑑 = III. 8
𝐴
Avec
𝑛
𝑖=1 𝐿𝑖 ∶ Somme des longueurs de tous les cours d’eau d’ordre i (Km)
𝐷𝑑 = 0.22 𝑘𝑚/𝑘𝑚2
L : distance entre ces 2 points (Km).
𝐼𝑐 = 2%
III.2.3.4.2. Temps de concentration
Le temps de concentration intervient dans la détermination du débit projet. Il est défini

comme étant le maximum de durée nécessaire à une goutte d’eau pour parcourir le chemin
hydrologique entre le point le plus éloigné du bassin versant et l’exutoire de ce dernier.
Théoriquement, on estime que le temps de concentration est la durée comprise entre la fin
de la pluie nette et la fin du ruissellement.
Mais pratiquement, il peut être déduit de mesures sur le terrain ou s’estimer à l’aide de
formules le plus souvent empiriques.
Le temps de concentration est influencé par diverses caractéristiques morphologiques dont

principalement (Khayra, 2015):
 La surface ;
 La forme (surface et longueur) ;
 Le relief du bassin (longueur et pente).
IL peut être notamment calculé par les formules morphométries suivantes :
76
Tableau ‎III-4: Différente formule du temps de concentration
Auteur Formule
4 𝐴 + 1.5𝐿
𝑇𝑐 =
0.8 𝐻𝑚𝑜𝑦 − 𝐻𝑚𝑖𝑛
𝑇𝑐 : Temps de concentration en heure
Giandotti (Touabia, 2004) A : surface du bassin versant en (Km²)
𝐈𝐈𝐈. 𝟗 L : longueur du cours d’eau principal (Km)
Hmin : altitude moyenne
Hmoy : altitude minimal
𝐴
𝑇𝑐 = 76.3 ∙
Ventura(Berbane, 2015) 𝐼𝑀
𝐈𝐈𝐈. 𝟏𝟎 𝑇𝑐 : Temps de concentration (min)
A : surface du bassin versant (km2)
𝐼𝑀 : pente moyenne du bassin (%)
3
Lr ∙ 𝐴
𝑇𝑐 = 64.8
𝐼𝑀
Passini(Iratni, 2014)
𝑇𝑐 : Temps de concentration (min)
𝐈𝐈𝐈. 𝟏𝟏 2
A : surface du bassin versant (km )
I : pente moyenne du bassin (%)
L : longueur du bassin (km).
𝐿𝑟 1.115
𝑇𝐶 = 0.067 ∙ 0.385
𝐻𝑚𝑎𝑥 − 𝐻𝑚𝑖𝑛
Basso (Berbane,2015)
𝑇𝑐 : temps de concentration en Heures ;
𝐈𝐈𝐈. 𝟏𝟐
𝐿𝑟 : longueur du bassin versant (Km) ;
𝐻𝑚𝑖𝑛 : altitude minimal (Km);
𝐻𝑚𝑎𝑥 : Altitude maximal (Km).
𝑇𝑐 = 0.6615 ∙ 𝐿0.77 𝐼𝑔 −0.385
Kirpich (Iratni,2014) L: longueur du chenal principal
𝐈𝐈𝐈. 𝟏𝟑 hydrographique en (Km);
Ig: indice de le la pente globale (%).
77
0.333 0.75
𝐴×𝐿
𝑇𝑐 = 4
Socolovsky ( Iratni,2014) 𝐼
𝐈𝐈𝐈. 𝟏𝟒 L: longueur du talweg principale en (Km) ;
A : surface du bassin versant en (Km2);
I: indice de pente en (m/km).
1
4
𝐴𝐵𝑉 × 𝐿
𝑇𝑐 = 3 3
Formule de L’ANRH : (SOGREAH 10𝐼𝐵𝑉 2
1986/1989) (Iratni,2014)
𝐼𝐵𝑉 : la pente du bassin versant en (%) ;
𝐈𝐈𝐈. 𝟏𝟓 L: longueur du chenal principal
hydrographique en (Km);
A: Surface du bassin versant en (Km2).
Les résultats de différentes méthodes sont affichés dans le Tableau ‎III-5 ci-dessous:
Tableau ‎III-5: Résultats de temps de concentration
Formule Giandotti Ventura Passini Basso Kirpich Socolvsky ANRH
Tc(h) 6,83 3,78 4,37 3,69 5,65 7,01 4,44
Dans ce cas, la valeur du temps de concentration prise est la moyenne entre les valeurs
calculées par les différentes formules, alors Tc est :
𝑇𝑐 = 5.10 𝑕eures
III.2.3.4.3. Le rapport de confluence
Sur la base de la classification des cours d'eau, Horton (1932) et Schumm (1956) ont
établi différentes lois présentées au tableau ci-après :
78
Tableau ‎III-6: Formule des rapports de confluence (Chow, 1988)
Loi Formule
𝑁𝑛
𝑅𝑏 = 𝑁
𝑛 +1
Rb : rapport de confluence des cours d'eau
Loi des nombres ("bifurcation ratio") ;

Nn : nombre des cours d'eau d'ordre n ;
Nn+1 : nombre des cours d'eau d'ordre
suivant.
𝑙 𝑛 +1
𝑅𝑙 = 𝑙𝑛
Rl : rapport des longueurs des cours d'eau ;

Loi des longueurs Ln : longueur moyenne des cours d'eau
d'ordre n ;
𝐴𝑛 +1
𝑅𝑎 = 𝐴𝑛
Loi des aires Ra : rapport des aires des cours d'eau ;

An: aire tributaire moyenne des cours d'eau
d'ordre n.
Avec :
n : ordre d'un cours d'eau n varie entre 1 et w (w est l'ordre du cours d'eau principal,
classification selon Strahler).
Les résultats obtenus sont présentés sur le tableau III.7.
Tableau ‎III-7: Résultats de rapport de confluence
Loi des nombre 𝑹𝒃 Loi des longueurs 𝑹𝒍 Loi des aires 𝑹𝒂

3.83 4.27 5.02
Le rapport de confluence est un nombre sans dimension exprimant le développement du

réseau de drainage. Il varie suivant l'ordre considéré. C'est un élément important à considérer
pour établir des corrélations d'une région à une autre. Selon Strahler (1964), le Rb varie de 3 à
5 pour une région où la géologie n'a aucune influence.
79
Figure ‎III-11: Les ordres de cours d'Oued Djendjen selon la classification de Horton.
Des informations récapitulatives concernant la morphométrie et l’hydrométrie du bassin

versant d'oued Djendjen sont résumées dans le tableau suivant :
Tableau ‎III-8: Caractéristiques géomorphologiques de bassin versant d'oued Djendjen
Caractéristiques morphologiques Symbole unité valeur

Surface A Km2 321
Périmètre P Km 102
Indice de compacité de Gravelius Kg / 1.59
Rectangle Longueur L Km 43.65
équivalent Largeur l km 7.39
Altitude maximale Hmax m 2006
Altitude moyenne Hmoy m 878.09
Altitude Altitude médiane H50% m 853.20

Altitude minimale Hmin m 237
Altitude de 5% H5% m 1331.60
Altitude de 95% H95% m 494.80
Densité spécifique Ds / 343.47
Indice de pente global Ig m/km 19.17
pente moyenne du BV Im % 4
pente moyenne de cours d'eau Ic % 36
80
Caractéristiques morphologiques Symbole unité valeur

La longueur principale du cours
L km 44.45
d'eau
Densité de drainage Dd km/km2 0.22
Temps de concentration Tc h 5.10
Rb 3.83
Rapports de confluence Rl / 4.27
Ra 5.02
III.3. Conclusion
L'étude des caractéristiques géomorphologiques du bassin versant de l'Oued Djendjen
nous a permis de comprendre le comportement hydrologique du bassin. Ce bassin versant
est de forme allongé avec une surface de 321 km2 et une pente moyenne de 4%. L'étude a
été menée sur un cours d'eau d'une longueur principale de 44.45 m et une pente moyenne de
36%.
81
IV. Numéro du chapitre
Nothing is softer or more flexible than water, yet nothing can resist it
— Lao Tzu
I l existe un grand nombre de méthodes hydrologiques qui conduisent à

des résultats différents pour l'estimation des crues de même probabilité
apparente. Le choix de la méthode doit être adapté aux caractéristiques
climatiques et régionales et aux donnés disponibles. Nous allons utiliser
quelques méthodes d'estimation des crues à savoir, celles basées seulement sur
les données d'observation des débits (les méthodes statistique "loi de Gumbel"),
celles basées sur les méthodes empiriques, celles basées sur les données
d'observation des pluies et des débits (méthode du Gradex) et en fin celles
basées sur les données d'observation des pluies (méthode de l'hydrogramme
synthétique).L’objectif principal de ce chapitre est de proposer et valider une
méthode d'estimation des débits de crue et des volumes de ruissellement, à
partir d’observations hydrométriques et météorologiques de bassin versant.
IV.1.Analyse et traitement des données

Les méthodes de calculs des débits se basent souvent sur les données disponibles dans
les zones où l’on souhaite effectuer l’étude. Le passage par la collecte de données est le
chemin le plus difficile pour chaque début d’études.
Nous avons travaillé sur une série des débits chronologiques observés dans le bassin
versant Djendjen pendant 21 ans (période de fonctionnement de 1972 à 1992 par la station
hydrométrique de Missa (code 03-04-08). Les donnée de base de cette étude sont recueillies
auprès du services de l'Agence National des Ressources Hydrauliques (A.N.R.H).
Cette série seras analysée pour extraire les crues observé sur 21 années. Nous avons pu
obtenir 370 crues sur le long de cette période.
82
Chapitre IV — Estimation du débit de pointe
Les données pluviométriques sont disponibles sur la période de 1914 jusqu’au2010.(

Annexe-A, tableau-1)
Pour mieux comprendre la méthodologie du travail nous avons l'incarner comme suit :
Données hydrométrique
station Missa
(code 03-04-08)
Séparation des épisodes

de crues(370 crues)
Calcul statistique
-Débit moyen horaire de
chaque crue.
-Débit de pointe de
chaque crue
Extraction des
débits
Débits de débits de pointe débits

pointe des crues maximum
horaire maximales journalier
Figure ‎IV-1: Diagramme explicative de la méthodologie d'extraction des données débit.
IV.2.Calcul de débit de pointe au droit de la station Missa

Le calcul des débits de pointe a été réalisé en appliquant la méthode probabiliste
d’ajustement statistique de Gumbel sur les débits maximums journaliers, les débits de
pointes maximales journalières et les débits de pointe horaire, ainsi que la méthode du
Gradex sur les débits maximums journaliers et les débits de pointes maximales journalières.
Le passage du débit journalier maximal extrapolé vers le débit de pointe correspondant a été
effectué par le biais d’un coefficient de pointe 𝐾𝑝 . Ce coefficient se définit comme étant le
83
rapport du débit de pointe maximale et du débit maximum journalier. Il traduit une relation
entre les deux débits.
débit de pointe maximale

𝑐𝑜é𝑓𝑓𝑖𝑐𝑖e𝑛𝑡𝑑e pointe =
débit maximum journalier
IV.2.1. Ajustement statistique

Pour une extrapolation directe des crues extrêmes à partir d’une série de 21 ans de
débits, on procède à un ajustement statistique selon la loi de Gumbel. Les graphes suivants
présentent l’ajustement des débits à la loi de Gumbel dans la station Missa.
A/ Débit maximum journalier
Figure ‎IV-2:Ajustement statistique de débit maximum journalier de la station Missa.
84
B/ Débit de pointe maximale journalier
Figure ‎IV-3: Ajustement statistique de débit de pointe maximale journalier de la station Missa.
C/ Débit de pointe horaire
Figure ‎IV-4:Ajustement statistique de débit de point horaire de la station Missa.
On observe que les débits paraissent alignés sur la droite d’ajustement statistique, Ce qui
confirme que la distribution des débits suit la loi de Gumbel. L'examen visuel de la
85
distribution des valeurs sur le graphique ne peut pas constituer une évidence sur la
conformité des données en question à la distribution par une loi de Gumbel. Il est
indispensable d'utiliser des tests d'adéquation adaptés. Nous avons utilisé le test d'adéquation
de khi-deux pour vérifier cette propriété.
IV.2.2. Teste d'adéquation par la loi de Gumbel

Le test du Khi-deux 𝜒 2 est un test non paramétrique qui permet de tester l'hypothèse
H0, selon laquelle les données observées sont engendrées par un modèle faisant intervenir
une loi de probabilité.
Hypothèses du teste Khi-deux:
H0: la distribution des valeurs suit une loi de Gumbel;

H1: la distribution des valeurs ne suit pas une loi de Gumbel
le tableau(IV-1) représente les résultat obtenus pour le test d'adéquation à la loi de Gumbel
pour les trois série de débit:
Tableau ‎IV-1:Teste khi2

Teste khi2
Débit de pointe Débit maximum Débit de pointe
Debit
horaire journalier maximale journalier
Degré de libeté 𝛾 3 3 3
N°-parametre de la
2 2 2
loi M
N°- classe K 6 6 6
χ2 𝐶𝑎𝑙𝑐𝑢𝑙 é 8,48 7,29 7,29
0,0633 0,0663
𝑃 χ2 𝐶𝑎𝑙𝑐𝑢𝑙 é 0,00379
Hypothése H1 Ho Ho
𝑑𝑜𝑛𝑐
𝑑𝑜𝑛𝑐
𝑑𝑜𝑛𝑐
𝑙 ′ 𝑎𝑗𝑢𝑠𝑡𝑒𝑚𝑒𝑛𝑡 𝑒𝑠𝑇 à 𝑐𝑜𝑛𝑠𝑖𝑑é𝑟𝑒
𝑙 ′ 𝑎𝑗𝑢𝑠𝑡𝑒𝑚𝑒𝑛𝑡 𝑒𝑠 à 𝑐𝑜𝑛𝑠𝑖𝑑é𝑟𝑒
𝑙 ′ 𝑎𝑗𝑢𝑠𝑡𝑒𝑚𝑒𝑛𝑡 𝑒𝑠 𝑟é𝑗𝑒𝑡é
𝑃 χ2 < 1%
𝑃 χ2 > 5%
𝑃 χ2 > 5%
Coclusion
86
D'après le tableau IV-1 la loi de Gumbel n’ajuste plus avec la série de débit de pointe horaire.
Les paramètres qu’on peut tirer à partir des séries des débits représentées dans l'annexe A.
(Tableau -2) Lorsqu’on applique la loi Gumbel, sont tabulés comme suite :
Tableau ‎IV-2: Les paramètres obtenus par la loi de Gumbel
Débits de pointe Débits maximaux

Débits de pointe horaire
maximale journalière journaliers
𝛼 𝑥0 𝛼 𝑥0 𝛼 𝑥0
104.503 180.941 142.297 187.828 164.253 203.849
Les quantiles des crues sont déterminés à partir de l’équation d’ajustement (II.9)
pour différentes périodes de retour, les tableaux suivant représentent les valeurs obtenues:
Tableau ‎IV-3: Les quantiles obtenus par l'ajustement des débits de pointe maximale journalière
Débits de pointe maximale journalière
Période de retour T Fréquence empirique Intervalle de
Qmax (m3 /s )
(ans) F(%) confiance
10 90 416.12 296 - 536
20 95 491.34 340 - 643
50 98 589.71 396 - 782
100 99 661.68 437 - 887
Tableau ‎IV-4: Les quantiles obtenus par l'ajustement des débits maximaux journaliers.
Débits maximaux journaliers
Période de retour T Fréquence empirique
Qmax ,j (m3 /s) Intervalle de
(ans) F(%) confiance
10 90 508.06 345 - 671
20 95 610.49 405 - 816
50 98 743.07 480 - 1010
100 99 842.42 536 - 1150
87
Tableau ‎IV-5:Les quantiles obtenus par l'ajustement des débits de pointe horaire
Débits de pointe horaire
Période de retour T Fréquence Intervalle de

Qmax ,h (m3 /s)
(ans) empirique F(%) confiance
10 90 573.48 385 -762
20 95 691.713 454 - 929
50 98 844.76 541 - 1150
100 99 959.44 606 - 1310
La Figure suivant représente les ajustement calculés par la loi de Gumbel pour des trois
données de débit:
1200,000
Qmax,j
1000,000
Qp(max,j)
800,000 Qmax,h
Q(M3/s)
600,000
400,000
200,000
0,000
0 20 40 60 80 100 120
T(h)
Figure ‎IV-5: Ajustement par la loi de Gumbel
La figure(IV-5) représente la déférence entre les quantiles calculés par l'ajustement des
trois débits : débit maximal journalier, débit de pointe maximale journalier et le débit de
pointe horaire par la loi de Gumbel, Cette différence de la pente 𝛼 (le Gradex) est
proportionnelle aux valeurs des débits calculé.
88
IV.2.3. La méthode Gradex

La méthode de Gradex est une approche hydro-pluviométrique probabiliste du calcul
des débits de crues. Elle s'applique dans les bassins jaugés pour calcule les quantiles de
période de retour à cent ans, ayant une station hydrométrique de référence.
Le schéma suivant représente les différentes étapes que nous avons appliquées pour calculer
les débits centennales.
Méthode du Gradex
Ajustement par la loi

de Gumbel
Ajustement pour Ajustement pour les

les séries séries de débit
pluviométrie
Précipitation Débit de point

maximale Débit maximal des crues
journalière. journalier maximales
Estimation du Gradex Débit décennal

des débits(𝛼) (débit pivot)
Calcul des quantiles
Figure ‎IV-6:Organigramme explicatif les différentes étapes appliqués de la méthode de Gradex
89
IV.2.3.1. Ajustement des précipitations maximales journalières par la loi de

Gumbel
L'ajustement des précipitations maximales journalières d'une série de 65 années de

mesure (observée à la station pluviométrique de Texanna) selon une distribution de Gumbel
permettre d'estimées paramètre qui souligne dans de tableau suivante:
Tableau ‎IV-6:les paramètres obtenus par la loi de Gumbel pour les Pmax,j
Paramètre de la loi de Gumbel pour les pluies maximaux journalier(𝑷𝒎𝒂𝒙,𝒋 )
𝛼 𝑥0
22.66 70.80
Les résultats des pluies obtenues pour différentes périodes de retour par l'ajustement
statistique par la loi de Gumbel sont tabulés comme suit :
Tableau ‎IV-7: les résultats d'ajustement des pluies maximales journalières par la loi de Gumbel
Période de Fréquence Intervalle de

𝐏𝐦𝐚𝐱,𝐣 (𝐦𝐦)
retour T (ans) empirique F(%) confiance
10 90 121.81 101-137
20 95 138.12 119-157
50 98 159.24 135-183
100 99 175.06 147-203
On représente l'ajustement statistique des pluies maximales journalières à la loi de Gumbel

dans la station Texenna par la figure ci-dessous.
90
Figure ‎IV-7: Ajustement statistique des pluies maximales journalières de la station Texenna
L’ajustement des précipitations selon la loi de Gumbel (Tableau IV-7) nous donne une
pente de 𝛼 = 22,66, qui est égale à la pente 𝛼 ′ estimé par la méthode Gradex.
La valeur du paramètre 𝑥0 est obtenue à partir des coordonnées du point de pivot qui
correspondant à la fréquence de non-dépassement de 2.25 pour les débits et pour 𝛼 ′ = 𝛼
On obtient les paramètres de la loi d’ajustement extrapolée des débits qui est parallèle à celle
des pluies. et en calculons les quantile pour la période de retour de 100 ans.
Tableau ‎IV-8: paramètres de la loi Gradex
Méthode du Gradex
Débits de pointe maximale journalière Débits maximaux journaliers
𝛼′ 𝑥0 𝛼′ 𝑥0
22.66 83.43 22.66 -17.540
IV.2.4. Calcul des quantiles pour les différentes périodes de

retour
Les résultats que nous avons trouvé d'après l'ajustement des débits selon la méthode de
Gradex pour des deux cas, débit maximaux journalier et le débit de pointe maximale
journalière sont représentés dans les tableaux suivant (tableau IV-9) :
91
Tableau ‎IV-9: Les quantiles obtenus par les débits maximaux journaliers.
Débits maximaux journaliers

période de
Fréquence Q Q Kp
retour T Q max (m3 /s )
empirique(%) (mm/j) (m3/s) (moyenne)
(ans)
10 90 33.46 124.32 508.06
20 95 49.78 184.94 755.75

4.09
50 98 70.89 263.39 1076.37
100 99 86.72 322.18 1316.69
le coefficient de pointe Kp est estimé à partir le rapport entre débit de pointe des crues
maximales et le débit maximaux journaliers et le. Coefficient se définit comme étant le
rapport du débit de pointe maximale et du débit maximum journalier. Il traduit une relation
entre les deux débits. son estimation peut être problématique pour le cas des bassins non
jaugés, donc des méthodes basées sur les caractéristiques du bassin versants peuvent être
employée pour son estimation.
Tableau ‎IV-10: Les quantiles obtenus par les débits de pointe maximaux journaliers.
période de Fréquence Q Q
retour T (ans) empirique(%) (mm/j) (m3/s)
10 90 112,00 416,12
20 95 128,316 476,73
50 98 149,4332 555,19
100 99 165,2576 613,98
Les graphes suivants montrent les droites obtenues en appliquant la méthode de Gradex
pour les débits maximaux journaliers et les débits de pointes.
92
200
Pluie (mm/j)
Débit 180
160
140
120
100
80 tg(α)
60 α
40
20
0
-2,000 -1,000 0,000 1,000 2,000 3,000 4,000 variable
5,000reduite(y)
6,000
Figure ‎IV-8Principe de méthode Gradex (les débits maximaux journaliers)
200
Pluie(mm/j)
180
Débit
160
140
120
100
tg(α)
80
60 α
40
Q 10
20
0
2.25
-2,000 -1,000 0,000 1,000 2,000 3,000 4,000 5,000 6,000
Figure ‎IV-9:Principe de méthode Gradex (les débits de pointe journaliers)
93
IV.2.5. Application des méthodes empiriques

Dans ce volet, on va appliquer la méthode de Fuller, Maillet et Gauthler et Mac Math et
celle de Giandotti pour l'estimation des débits des crues de différentes périodes de retour.
Les expressions analytiques de ces formules sont représentées au chapitre précédent
(Tableau IV-11).Le choix des méthodes de calcul se faits selon les données disponibles.
Alors on recueil les résultats obtenus dans le tableau ci-dessus :
Tableau ‎IV-11 : les résultats obtenus de débit de crues basé sur les méthodes empiriques
Maillet et
T(ans) Fuller Mac Math Giandotti
Gauthler
10 434.272 269.232 360.76 297.87
20 492.373 327.856 409.07 337.76
50 569.180 392.124 471.62 389.40
100 627.281 434.47 518.47 458.96
IV.2.6. Modélisation pluie-débit

En vue de reproduire le comportement du bassin versant face à des sollicitations
pluviométriques, la modélisation pluie-débit peut s’avérer un outil assez puissant et fiable.
Cette modélisation est basée sur des observations historiques ; quand ces observations ne
sont pas disponibles une approche géomorphologique peut être envisagée.
IV.2.6.1. Séparation des écoulements de base
L’écoulement de base est une composante de l’écoulement total qui constitue la

contribution du stock sous terrain par écoulement retardé, cette composante est souvent
isolée pour obtenir uniquement l’écoulement dit de surface.(Figure ‎IV-10)
94
Figure ‎IV-10: Les composantes de l’écoulement( Cour Musy,2005)
La méthode utilisée est celle du filtre digitale proposé par Arnold et Allen (1999), il s’agit
d’une méthode récursive donnée par l’équation suivant:
1 
qt    qt 1  (Qt  Qt 1 ) (IV-1)
2
Où:
𝑞𝑡 : Le débit filtré (réponse rapide) à l’instant t;
𝑄𝑡 :‎Le‎débit‎total‎à‎l’instant‎t;
β : le paramètre du filtre digitale qui varie de 0.925 à 1 et qui dépend des conditions
géologiques du bassin versant.
Dans la modélisation évènementielle au début de l’épisode le débit 𝑞0 est considéré nul.
IV.2.6.2. Procédure de calibration
Les évènements tabulés ont été sélectionnés et ils seront utilisés pour la calibration du
modèle, qui consiste à trouver les paramètres optimaux qui nous donnent des erreurs
minimales. Deux fonctions objectives ont été utilisées pour la calibration .
95
La procédure de calibration est résumé dans organigramme suivante:
Figure ‎IV-10:Organigramme explicatif du procédure de calibration.
96
Tableau ‎IV-12 : Fonction objective
Function objective Equation Optimum
( Robs  Rsim ) i
N 2
1
Bias (IV-2) N
 i 1 Robsi
0
 (Q obs Qsim ) 2
Nash- Sutcliffe 1 i 1
1
N
 (Q
efficiency (NSE)
obs Qobs ) 2
(IV-3) i 1
IV.2.6.2.1. Fonction de production
Pour la fonction de production SCS à 2 paramètres, Nous obtenons les valeurs

optimales suivantes :
𝐶𝑁 = 81
𝜆 = 0.005
Avec:
𝐵𝑖𝑎𝑠 = 33%
𝑁𝑆𝐸 = 70%
Tableau ‎IV-13 : calibration et validation de la fonction de production
ID P(mm) Robs(mm) S Ia Rsim(mm) %err

1 20.80 11.47 59.58 0.30 5.2 54.227
2 40.10 10.09 59.58 0.30 15.9 57.951
3 21.60 6.46 59.58 0.30 5.6 13.193
4 14.10 7.96 59.58 0.30 2.6 67.400
5 16.80 2.17 59.58 0.30 3.6 64.610
6 41.40 16.84 59.58 0.30 16.8 0.351
7 46.00 22.34 59.58 0.30 19.8 11.192
8 15.20 3.02 59.58 0.30 3.0 1.175
97
25,0
20,0
Ruissélement simulé (mm)
15,0
10,0
5,0
0,0
0,00 5,00 10,00 15,00 20,00 25,00
Ruissélement obsrvé (mm)
Figure ‎IV-11 :Ruissèlement simulé vers sue ruissèlement observé
Partant du même principe, on procède à la calibration de la fonction de transfert afin de

satisfaire les mêmes critères.
Le modèle de transfert adopté est le modèle de Nash à 2 paramètres qui seront calibrés à
travers les évènements disponibles.
E. Fonction de transfert
Une approche géomorphologique permet d’estimer à priori les paramètres de l’hydro-

gramme unitaire instantané qui sont dérivées des caractéristiques morphométriques du
bassin versant. La pente moyenne du bassin versant est utilisée pour l’estimation de la vitesse
moyenne de l’écoulement qui sera utilisée dans l’estimation des paramètres du
l’hydrogramme unitaire instantané géométrique.
Corradini, (2002) a proposé une formulation du temps d’atténuation (Lag-time) pour les
bassins versants méditerranéens
K    1.19 A0.33 (IV-3)
K : le temps d’atténuation en heure;
A : la surface du bassin versant en Km²;
98
η : est un paramètres à optimiser qui dépend des caractéristiques morphométries du

bassin versant (Brocca et al., 2011).
Une autre formulation a été proposé par SCS(1972) du coefficient
𝐾 = 0.6𝑇𝑐 (IV-4)
Tc [h] : le temps de concentration du bassin versants.
Ces formulations disponibles ont été prises en compte pour la calibration du modèle avec la
méthode (‘fmincon’ fonction dans MATLAB®) qui est une méthode basée sur un gradient
standard automatique pour l’optimisation non linéaire.(Brocca,2011)
Figure ‎IV-12: Calibration modèle cascade de NASH
Nous obtenons les coefficients K=3.9164 h et N=2.1 pour un critère NSE de l’ordre de
90% comme le montre la figure 2.
99
Figure ‎IV-13:Validation du modèle de NASH

Les paramètres K et N seront appliquée lors de l’estimation de l’hydrogramme de crue
généré par les crues de projets utilisées
IV.3.Pluie de projet
La pluie de projet est un évènement pluvieux fictif et est une représentation synthétique
de formes imposées (triangulaires, trapézoïdale,…). Elle est censée représenter l’ensemble
des informations pluviométriques disponibles (la pluie de projet est un évènement type) et
de produire à l’exutoire du bassin (le même effet que la pluie réelle (débit, volume, etc.…) à
laquelle il est possible d’associer une période de retour (Bourrier, 1997 ; Chocat, 1997 et
Valiron 1989).(Boutaghane et al; 1945)
L'élaboration des pluies de projet repose sur l'analyse statistique des divers éléments d'un
épisode pluvieux: durée de pluie, hanteur maximale précipitée sur intervalle de temps,
d'intensité de précipitation.
La répartition temporelle des précipitations ne peut pas être le même pour les différents
Avers Il est donc impossible de trouver deux pluie de projet avec le même schéma de
distribution.(chow,1988)
IV.3.1. Hyétogramme des pluie synthétiques

Le service de conservation des sols du Département de l'agriculture des États-Unis a
développé des hyétogrammes de pluie synthétiques à utiliser aux États-Unis pour les pluie de
100
6 et 24 heures. Ces hyétogrammes ont été dérivés des informations présentées par
Hershfield (1961) et Miller, Frederick et Tracey (1973) et d'autres données de crues.
Le tableau (3) dans l'annexe présente les hyétogrammes cumulatifs. Il y a quatre pluie d'une
durée de 24 heures, appelées types I, IA, II et III, respectivement. Ils devraient être
appliqués:
 Les types I et IA sont pour le climat maritime du Pacifique avec des hivers humides
et des étés secs.
 Le type III concerne le golfe du Mexique et les régions côtières de l'Atlantique, où

les tempêtes tropicales produisent de grandes quantités de pluie sur 24 heures.
 Le type II est pour le reste de la nation. (Chow,1988)
IV.3.2. Méthode de la pluie de projet double triangle

symétrique
Le principe de « la pluie de projet double triangle » a été développé en France par Des
bordes Chocat (1974) à partir de l’analyse statistique chronologique de pluies réelles. Elle est
construite à partir des courbes intensité-durée-fréquence et peut donc être adoptée à
différentes régions.
En partant du principe de « la pluie de projet double triangle » de Desbordes Chocat (1997),

propose qu’on puisse considérer cette pluie parfaitement symétrique lorsqu’on utilise le
modèle à réservoir linéaire, en supposant que la pointe du débit à l’exutoire apparaît pendant
la pluie intense (Figure IV-15).(Kwan , 2008)
101
Figure ‎IV-14:la pluie de projet double triangle (Kwan Tun lee,2008)
La période de retour sur la durée totale est prise égale à celle sur la durée de pluie
intense. L'hythogramme de la pluie de projet double triangle symétrique peut être
déterminée à partir les Courbe IDF comme suite:
𝑡
𝑑 − 𝑑𝑚 , 0 ≤ 𝑡 < 𝑡2
𝑡1 𝑡2
𝑡 − 𝑡1 𝑡𝑑
i t 4𝑑𝑚 , 𝑡2 ≤ 𝑡 < (IV-5)
𝑡𝑑 − 2𝑡1 2 2
𝑡𝑑 − 𝑡1 − 𝑡 𝑡𝑑
4𝑑𝑚 , ≤ 𝑡 < 𝑡𝑑 − 𝑡2
𝑡𝑑 − 2𝑡1 2 2
Avec:
1
𝑡1 = 𝑡 − 𝑡𝑚
2 𝑑 (IV-6)
𝑡𝑚 2
𝑡2 = 𝑡1 + 𝑑 − 𝑑𝑚
4𝑡1 𝑑𝑚 (IV-7)
Où:
102
𝑖(𝑡) : Intensité des précipitations à l'instant t;
𝑡𝑑 : La durée totale de hydrogramme;
𝑡𝑚 : Durée de hydrogramme triangulaire central;
𝑑 : Hanteur maximale précipitée pour la durée td;
𝑑𝑚 : Lame précipité correspondant à la durée 𝑡𝑚 ;
Pour 𝑡𝑑 et 𝑡𝑚 donnés, les valeurs correspondantes de 𝑑 et 𝑑𝑚 peuvent être obtenues à partir

d'une courbe IDF pour une période de retour spécifiée.
L'équation intensité-durée de la pluie utilisée dans cette étude peut être s'écrite par Montana
est:
𝐼 𝑡 = 𝑎𝑡 −𝑏 (IV-8)
Avec a et b sont des coefficient de Montana (Kwan Tun lee,2008)
IV.4.Application de la méthode SCS-CN et HUI du

Modèle de NASH sur les pluies de projet
Nous reprenons les résultats des paramètres qu'on a déterminés par calibration des
deux fonctions (production- transfert) dans le tableau ci-dessous:
Tableau ‎IV-14: Les paramètres du model pluie-débit
Fonction de transfert Fonction de production

Coefficient de Nombrer de
𝜆 Curve number (CN)
stockage (K) réservoir(N)
0.005 81 3.9164 2.1
Après la détermination des paramètres de chaque modèle (SCS et NASH). Nous

passerons à une autre étape pour appliquer les modèle pluie-débit sur les pluies de projet.
Nous avons appliqué la fonction de production du SCS-CN au pluies synthétiques avec les
la répartitions temporelle de type (I-IA-II-III) d'une part et une répartition aléatoire de pluie
de projet d'autre part et finalement sur la pluie de projet double triangle symétrique pour
calculer les pluie nettes distribués sur 24h.(Annexe-tableau 4)
103
La transformation de pluie de projet par les différentes parties du bassin versant, en

hydrogramme à l'exutoire du bassin est basée essentiellement à la fonction de transfert.
Après avoir exprimé les notions de base relative à la détermination d'un hydrogramme
unitaire de Nash (chapitre II), nous avons appliqué la technique de convolution pour calcule
le débit de point sur la période de retour de 100 ans. Le diagramme ci-dessous montre la
procédure de la modélisation pluie-débit.
Pluie de projet
Paramètres du Fonction de
production (SCS- Pluies nettes
modèle (CN-𝜆)
CN)
Fonction de
Paramètres du transfert
modèle (K-N) (Model de
NASH)
Débits de pointe(T)
Figure ‎IV-15: Organigramme de l'application des modèles P-Q
Les résultats obtenus de la modélisation Pluie-débit sont représentés dans le tableau et la

figure suivantes:
Tableau ‎IV-15: Les débits obtenu par la modélisation P-D pour T=100ans
Modèle Pluie-Débit
type double Répartition

type I type II type III type I-A
triangle aléatoire
684.93 824.69 824.77 581.05 766.71 862.62
104
Figure ‎IV-16: Hydrogrammes de crue obtenu par la modélisation P-D pour T=100ans
IV.5.Interprétation des résultats

Les résultats des débits centennaux obtenus par les différentes méthodes d'estimation
sont récapitulés dans le tableau suivant :
105
Tableau ‎IV-16: tableau récapitulatif des débits centennaux selon plusieurs méthodes
Débits maximaux journaliers (m3/s) 842.42 1

Gumbel 661.68 2
(m3/s)
Débits de pointe annuelle (m3/s) 959.44 3
Débits maximaux journaliers (m3/s) 1316.63 4

Gradex Débits de pointe maximale journalière
613.98 5
(m3/s)
Fuller (m3/s) 627.28 6
Méthodes Maillet et Gauthler(m3/s) 434.47 7
empiriques Mac Math(m3/s) 518.47 8
Giandotti(m3/s) 458.96 9
Type III(m3/s) 824.77 10
Type I(m3/s) 684.93 11
Modèles pluie- Type I-A(m3/s) 581.05 12

débit Type II(m3/s) 824.69 13
Type double Triangle(m3/s) 766.71 14
Type Aléatoire(m3/s) 862.62 15
Le tableau IV.16 montre qu’il y a une grande différence entre les résultats issus des
différentes méthodes. Cette différence est due d’un côté aux données de départ utilisées ainsi
qu’aux hypothèses et fondement de l’approche utilisée.
Chacune de ces méthodes comportent des avantages et des inconvénients d’où l’objet de
cette étude. Les approches utilisées peuvent être classée sous quatre grandes catégories ;
méthodes statistiques; méthodes empirique et méthodes basée sur la modélisation pluie-
débit.
Dans le cas des méthodes statistiques, l’ajustement par la méthode de Gumbel et la méthode
du Gradex ont été appliquée : pour la méthode Gumbel, afin d’estimer le débit de pointe,
l’ajustement a été appliqué aux débits extraits à partir d’une analyse statistique de la série
chronologique des débits enregistrés à la station Missa. Ces débits extraits sont les débits
106
moyens des crues maximales (Débits maximaux journaliers), les débits de pointe des crues
maximales(Débits de pointe maximale journalière) et les débits de pointes horaire , qui sont
des débits différents et qui donnent des résultats nettement différents.
Nous remarquons qu’il ya une différence dans l’estimation du débit centennale quand on
utilise différents débits de pointe pour la méthode de Gumbel (voir lignes 2 & 3 du tableau
IV.16).
Les estimations de débits de crues obtenus par la méthode de Gumbel pour une période de
retour de 100 ans sont relativement proches et donnent des résultats assez semblables
lorsque qu’on utilise les débits moyens des crues maximales et les débits de pointes. D'un
autre coté, l'utilisation des débits de pointe des crues maximales journalière sous-estime le
débit DE POINTE d'environ 30% (voir lignes 2 et 3 du tableau IV-15.
Cette déférence entre ces quantile (case 1-2-3) dépend de la pente α (le gradex)de la courbe
d'ajustement des trois débit: débit de pointe (𝛼 = 164,253); débit maximaux journalier
(𝛼 = 142,297)et le débit de pointe maximale journaliers (𝛼 = 104,503). Le débit calculé
est proportionnel à la valeur de la pente (voir figure IV-5 et tableau IV-2).
Pour la méthode Gradex, qui est une méthode intéressante dans le cas des bassins non
jaugés, elle a été appliquée pour estimer le débit de pointe en ayant pour point pivot le débit
calculé par ajustement statistique de la loi Gumbel pour une période de retour de 100 ans.
Si on n'a pas ou très peu d'enregistrement de crues du bassin étudié, il sera possible d'estimer
le débit de crue désiré par analogie avec d'autres bassins. Différentes méthodes régionales
permettent d'effectuer cette estimation. Les méthodes régionales les plus couramment
utilisées sont les méthodes empiriques.
L'estimation du débit de crue centennal par l'ajustement statistique de la loi Gumbel pour les
débits maximaux journaliers et les débits de pointe ( max) donne des résultats assez proche à
ceux obtenus par la méthode Gradex pour le débit de pointe (max) voir la case 4 .
Les résultats des débits centennaux obtenus après l'application des deux méthodes
statistiques sur la même série de donnés (débits de pointes maximaux journaliers) sont très
proches avec une différences de moins de 8%.
107
Les méthodes empiriques sont souvent utilisés, ont données des valeurs très faibles
comparant avec les méthodes statistique qui sont plus fiables car elle prennent en
considération les données historique ( débit pour gumbel) et ( pluies et débits décennale
pour gradex)
par contre Les formules empiriques sont utilisées lorsqu’on ne possède que peu ou pas de
données sur les débits des crues dans une région (bassin non jaugé) elle sont basées sur des
données régionales. Ces formules sont établies pour de nombreux cours d’eau et dans divers
pays, permettant d’estimer soit des débits maximums de crues soit des débits fréquentiels à
partir de certaines caractéristiques du bassin versant en les complétant parfois par certaines
données météorologiques en particulier la pluviométrie.
Les modèles pluies débits sont l’outil incontournable dans ces dernières décennies et ont
permis de faire une estimation correcte et qui se rapproche le plus des résultats statistiques,
mais c’est l’unique approche qui permet de construire un hydrogrammes de crue, qui, parfois
peut se révéler indispensable pour le dimensionnement de certains ouvrages et pour le calcul
de laminage.
Nous avons utilisé un procédé de calibration pour estimer les paramètres optimaux du
modèle pluie-débit. Ces paramètres sont 𝜆 et CN pour la fonction de transfert et K et N
pour la fonction de projection. Après la détermination de ces paramètres nous avons
appliqué la méthode SCS-CN et le model de NASH pour calculer les débits de crue
centennales et de construire les hydrogrammes de crues de la période de routeur de 100 ans .
L’inconvénient la modélisation pluie-débit c’est la calibration des modèles qui nécessitent

une base de données assez large mais pour des projets d’envergure elle permet de faire une
estimation plus fiable.
IV.6.Conclusion
L'estimation du débit des cures est une tâche très essentiel lors d'une étude hydrologie d'une
bassin versant. Dans Ce chapitre nous avons appliqué divers méthode pour estimé le
débit selon quatre approche hydrologique: les méthode statistique; les méthode empirique et
les méthodes basée sur la modélisation pluie-débit.qui donnent les résultat défirent. Les
méthodes basées sur la modélisation pluies-débits est celle qui donne le meilleur résultat par
apport au autre.
108
V. Numéro du chapitre
The process of scientific discovery is, in effect, a continual flight from wonder
—Albert Einstein
L a modélisation d’un phénomène aléatoire consiste à choisir un

modèle, connu pour refléter une image proche de la réalité et qui
permet de mieux comprendre le mécanisme et le fonctionnement du
phénomène. Cette étude représente une schématisation simplifiée d’un système
réel de l’écoulement à surface libre. Notre objectif est d’avoir les lames d’eaux
correspondantes au débits centennaux issus de l’étude hydrologique en
utilisant l’outil informatique HEC-RAS (Hydrologic Engineering Center –
River Analysis System) C’est un outil de modélisation hydraulique des
écoulements graduellement variés et à surface libre sous le modèle 1D de
Saint-Venant.
V.1. L’objectif de l’étude Hydraulique

L’étude hydraulique consiste à faire le diagnostic de la zone d’étude pour les
événements hydrauliques. Elle consiste à déterminé les hauteurs d’eau le longe du l'oued
étudier et les différentes caractéristiques hydraulique liées à ce écoulement . Pour cela, dans
cette partie du travail, nous avons utilisé deux logiciels, le premier logiciel de Système
d’Information géographique (SIG qui combiné avec HEC-GEORAS) et le deuxième logiciel
de modélisation de rivières (HEC-RAS) qui permettent de reproduire l’état naturel de la
zone à l’état numérique.
Le logiciel de modélisation hydraulique HEC-RAS permet de connaître le comportement

d’un cours d’eau et le cas échéant de prévoir les inondations et de pouvoir anticiper les
éventuels dégâts.
109
Chapitre V-Etude hydraulique
Une simulation hydraulique d’un modèle est un calcul, dont les résultats sont
caractéristiques du débit, de la géométrie du cours d’eau ainsi que des conditions aux limites
du modèle. Les données topographiques et les conditions aux limites sont les données de
base pour concevoir le modèle dans HEC-RAS. Il existe différents types de données
topographiques qui dépendent de l’extension géographique du modèle, de sa résolution et
surtout de la précision des valeurs. (Yahiaoui, 2012)
V.2. Présentation d'outilles d'information utilisées

Nous avons adopté une méthodologie basée sur le couplage entre ARCGIS et HEC-RAS à
l’aide de HEC-GEORAS. L’ensemble de ces outils constitue un dispositif informatique
intégré et cohérent
HEC-
ARCGIS HEC-RAS
GEORA
S
Figure ‎V-1: Représentation schématique du dispositif informatique utilisé
V.2.1. Présentation du modèle HEC-GEORAS

Le logiciel HEC-GeoRAS est un ensemble de procédures, d'outils et d'utilitaires pour le
traitement des données géo-spatiales dans ArcGIS en utilisant une interface utilisateur
graphique. Il permet d’une part, la préparation des données géométriques à exporter vers
HEC-RAS, et la récupération des résultats de simulation effectuée par HEC-RAS.
L’élaboration du fichier de la géométrie à exporter, est basée sur un modèle numérique de
terrain qui doit être de préférence sous forme de TIN.
110
Les données résultantes de la simulation par HEC-RAS, peuvent être exploitées par HEC-
GeoRAS sous ArcGIS pour analyser et cartographier les zones inondables, en plus d’autres
possibilités.(Benmia,2012)
V.2.1.1. Présentation du modèle HEC-RAS
Le modèle HEC-RAS (Hydrologic Engineering Center River Analysis System) est un

modèle unidimensionnel conçu par Hydrologic Engineering Centre de l’US Corp Engineers
pour modéliser l’écoulement à surface libre permanent et non permanent qui est conçue
pour exécuter les tâches suivantes: (Bachi, 2011)
Model
HEC-RAS
Utilisation
Calcule les Simulation de calcul la Exécution d'une

hauteur d'eau l'écoulement no capacité de analyse de
permanant transport solide qualité de l'eau
Figure ‎V-2: Organigramme récapitulative les tâche du logiciel HEC-RAS (US ACE, 2010)
Le HEC-RAS consiste en général à déterminer la hauteur d’eau dans la rivière à chaque

débit entrant et aussi quelques variables caractéristiques de l'écoulement , comme la charge,
la vitesse de l’eau, ou autre. Enfin un écoulement monodimensionnel donne une valeur
unique de chaque variable pour chaque profil en travers du modèle.
Les bases théoriques d’HEC-RAS sont publiées dans les brochures d’aide du logiciel. Malgré
cela, on va citer les principales concernant le régime permanent. L’HEC-RAS a la capacité de
détermination unidimensionnelle des lames d’eaux à partir du calcul de la variation de la
hauteur d’eau entre deux sections successives par la résolution, et d’une procédure itérative,
de l’équation d’énergie suivante :(Bachi, 2011).
111
𝛼2 𝑉1 2 𝛼2 𝑉2 2 𝑉. 1
𝑦2 + 𝑍2 + = 𝑦1 + 𝑍1 + + 𝑕𝑐
2𝑔 2𝑔
Tel que :
𝑦1 et 𝑦2 : Tirant d’eau dans chaque section;
𝑍1 , 𝑍2 : Côte de fond;
𝑉1 , 𝑉2 : Vitesse moyenne dans chaque section;
𝛼1 , 𝛼2 : Coefficients relatives à la répartition des vitesses dans une section;
𝑔 : Accélération de la pesanteur;
𝑕𝑐 : Perte de charge totale.(USACE,2011)
Figure ‎V-3: Représentation des termes de l’équation d’énergie(US ACE, 2002)
La perte de charge est composée de perte de charge linéaire et singulière.
𝛼2 𝑣22 𝛼1 𝑣12
𝑕𝑐 = 𝐿𝑆𝑓 + 𝑐 − 𝑉. 2
2𝑔 2𝑔
112
Avec :
2
𝑄
𝑆𝑓 = 𝑉. 3
𝐾
L : Distance entre deux sections,
C : Coefficient de perte par expansion ou contraction,
𝑆𝑓 : Pente représentative de frottement entre deux sections tel que :

𝑄: est le débit d'eau; 𝑄 = 𝑉 ∙ 𝐴 𝐸𝑞𝑢𝑎𝑡𝑖𝑜𝑛 𝑑𝑒 𝑐𝑜𝑛𝑡𝑖𝑛𝑢𝑖𝑡é
𝑄
d'ou: 𝑉= 𝑉. 4
𝐴
𝐸 et K : Coefficient de débit calculé comme suite :
1.486 2
𝐾= 𝐴𝑅𝑕 3 𝑉. 5
𝑛
Avec :
n : Coefficient de rugosité de Manning;

A : Section mouillé;
R : Rayon hydraulique.(Bachi,2011)
V.3. Etape de simulation PAR HEC-RAS

La simulation de l’écoulement d’eau par HEC-RAS passe par les étapes suivantes :
 Définition de MNT et la géométrie de la rivière sur Arc-Gis à l'aide de HEC-

GEORAS
 Importation du projet qui traité sur HEC-REORAS sur le logiciel HEC RAS,
 Choix des paramètres hydrométriques et conditions aux limites,
 Simulation hydraulique par HEC-RAS
V.4. Présentation générale

Un système intégré désigné à effectuer des calculs et des analyses hydrauliques, stocker
gérer les données et restituer les résultats sous forme de tableaux ou de graphiques.
113
Figure ‎V-4: L’interface du logiciel HEC-RAS 5.2.0
Projet File : file name.prj
Geometry Data Steady Flow Data

fileename.g01 file name. F01
Plan 1
file name.p01
fil
Run File 1
fil
File name.r01
Output File 1
file name. O01
Figure ‎V-5: Schéma illustrant les fichiers d'un projet(HEC-RAS users mnual)
Tout comme la plupart des logiciels (dont ArcGis d’ailleurs), HEC-RAS utilise le concept de
Projet afin de réaliser une simulation du comportement hydraulique d’un cours d’eau. Tel
que défini par HEC-RAS, un projet est un ensemble de fichiers permettant de simuler le
comportement hydraulique d'un cours d'eau sous diverses conditions. La figure ci-dessus
illustre la hiérarchie des fichiers d'un projet.
Il faut souligner l’existence de deux type de fichiers essentielles : les fichiers de données de
géométrie (Geometry Data files) et ceux de données d’écoulement (Flow Data files), qui
correspondent effectivement aux deux étapes primordiales de préparation de la simulation
114
par Hec-Ras, et qui seront détaillé ultérieurement. Les fichiers Plan, représentent des
scénarios programmés par l’utilisateur, c’est-à-dire, ce sont des combinaisons entre des
fichiers de géométrie et des fichiers d’écoulement. Les autres fichiers sont générés
automatiquement par le logiciel. .(Benmia,2012)
V.5. Elaboration du fichier de la géométrie

L’extension HEC-GeoRAS nous a permis la préparation des entités géométriques
nécessaires à la construction du modèle hydraulique. Après avoir mis en place les couches
utilisées comme fonds d’extraction de l’information (cartes topographiques, MNT, images
satellitaires, stations hydrologiques…), les fonctionnalités fournies par l’extension HEC-
GeoRAS nous ont permis d’une manière organisée, de numériser les entités suivantes :
 Le réseau hydrographique : à partir essentiellement du MNT en prenant en

considération le tracé des rivières sur la carte et l’image satellitaire.
 Les lignes berges (bank lines) : en se basant encore essentiellement sur la topographie
du MNT. Ces lignes permettent de délimiter le chenal d’écoulement.
 Les lignes d’écoulement (Flow path centerlines) : ils enveloppent le chenal

d’écoulement et les lignes berges, permettant au logiciel de calcul et des distances
entraînent chaque profil en travers et celui à son aval.
 Les lignes des profils en travers (Cross section cut lines) : c’est l’entité la plus
importante à laquelle il faut accorder le maximum d’attention puisque c’est le profil en
travers qui détermine l’altitude de la surface d’écoulement et la variation de sa pente,
ainsi que la topographie de la plaine d’inondation.
V.6. Simulation hydraulique

V.6.1. Le choix du coefficient de la rugosité (coeff. de Manning
Strickler)
Un des paramètres à définir représente la rugosité du lit. Ce coefficient n, appelé
coefficient de Manning, contribue à avoir une modélisation la plus proche de la réalité en
prenant en compte la vitesse d’écoulement de l’eau sur le fond de la rivière et sur les berges.
Il faut signaler en outre, que Hec-Ras permet à l’aide d’une fonctionnalité très avancée, de
prendre en considération l’atténuation de la rugosité à mesure que la profondeur de l’eau
augmente, ce qui rend le modèle encore plus réaliste.
115
1 2 1
𝑉𝑖𝑡𝑒𝑠𝑠𝑒 = ∙ 𝑅 𝑕 3 ∙ 𝑆2 𝑉. 6
𝑛
Avec:
Rh : Rayon hydraulique.
S :Surface
Le coefficient de Strickler (k) est également utilisé dans la littérature et correspond à l’inverse
du coefficient de Manning : Le coefficient de Strickler (k) est également utilisé dans la
littérature et correspond à l’inverse du coefficient de Manning :(Bachi,2011)
1
𝑘= 𝑉. 7
𝑛
Dans cette partie, nous avons proposé Le coefficient de Manning selon le tableau suivant:
Tableau ‎V-1: Coefficients de Manning

𝟏
Type de cour d'eau Valeurs usuelles de K(𝒎 𝟑 /𝒔)
Canal bétonné, très Liss 75 à 100

Canal en terre 50 à 75
Cours d'eau régulier, bien entretenu 30 à 50
Cours d'eau régulier, Bien entretenu 40 à 50
Cours d'eau ordinaire 30 à 40
Cour d'eau avec embâcles 20 à 30
On prend un coefficient de Manning n=0.0035
V.6.2. Donné hydrométrique et Les condition au limite

Après avoir introduit la géométrie du cours d’eau, l'étape suivante est de spécifier les
débits utilisés pour calculer les profils d'écoulement. Pour cela, on a engagé la simulation
hydraulique en régime d'écoulement mixte avec une valeurs de débits qui calculé par suivant
la loi Gradex et qui correspond au périodes de retours 100 ans. Les limites amont et /ou
aval de chaque tronçon doivent soumettre à des conditions qui peuvent être de différentes
natures. HEC-RAS propose une large gamme de conditions aussi bien pour l’amont que
pour l’aval.
116
Figure ‎V-6: L’introduction des conditions aux limites pour l’écoulement
V.6.3. Les profils en long

La ligne d'eau représente une vue en profil du canal d'étude avec le niveau et la surface
d'eau de l'Oued pour chaque station.
117
Figure ‎V-7: le profil en long du tronçon modélisé d’Oued Djendjen
118
V.6.4. Les profils en travers

Cette option d’affichage des résultats par profil transversal, permet d’acquérir le
maximum de détail à chaque point choisi du tracé du cours d’eau.
La figure suivent représente les profils en traverse quelque station le long du cour d'eau:
Station -2000 Station--6000

FINAAAAAAAAAAAAAAAALE Plan: Plan 01 6/21/2018 FINAAAAAAAAAAAAAAAALE Plan: Plan 01 6/21/2018
.035 .035
270 Legend 342 Legend
WS PF 1 341 WS PF 1
269
Ground Ground
Elevation (m)
Elevation (m)
Bank Sta 340 Bank Sta
268
339
267
338
266 337
0 20 40 60 80 100 120 140 0 20 40 60 80 100 120 140
Station (m) Station (m)
Station-10000 Station-1400
.035 .035
408 468.0 Legend
Legend
407 WS PF 1 467.5 WS PF 1
Ground
Elevation (m)
Ground
Elevation (m)
406 467.0
Bank Sta Bank Sta
405 466.5
404 466.0
403 465.5
402 465.0
0 20 40 60 80 100 0 50 100 150 200
Station-3600 Station-38000
.035 .035
740.5 792 Legend
Legend
740.0 791
WS PF 1 WS PF 1
739.5 Ground 790 Ground
Elevation (m)
Elevation (m)
739.0 Bank Sta 789 Bank Sta
738.5 788
738.0 787
737.5 786
737.0 785
0 20 40 60 80 100 120 140 160 0 20 40 60 80 100 120
Figure ‎V-8: Profile en traverse de quelque sections
119
V.6.5. Vitesse d'écoulement

La figure V-7, Représente la variation de la vitesse d'écoulement des eaux en fonction
de la distance à partir l'exutoire sur la bordure gauche, au centre et sur la bordure droite de la
rivière.
Plan 01 6/21/2018
12 Legend
Vel Chnl PF 1
10
Vel Left (m/s), Vel Chnl (m/s), Vel Right (m/s)
0
0 5000 10000 15000 20000 25000 30000 35000 40000
Main Channel Distance (m)
Figure ‎V-9: Profil du vitesse

V.6.6. la courbe de tarage
HEC-RAS a permis de fournir une représentation de la variance de la hauteur d’eau en
(m) en fonction de débit en m3/s
Plan 01 6/21/2018
791 Legend
W.S. Elev
790
789
W.S. Elev (m)
788
787
786
785
0 200 400 600 800 1000 1200 1400
Q Total (m3/s)
Figure ‎V-10: courbe de tarage
120
V.6.7. principaux paramètres hydrauliques

Le logiciel permet aussi de créer des tableaux récapitulatifs de l’ensemble des paramètres
hydrauliques pour un barrage ou autre ouvrage, pour une surface de stockage, pour chaque
station (profil en travers), ou pour l’ensemble des stations à la fois.
Les principaux paramètres hydraulique calculés pour chaque station sont représenté au
tableau ci dessous :
Tableau ‎V-2: les principaux paramètre hydraulique calculé pour chaque

station de La hanteur La Nombre
Vitesse Largeur au
Surface pente profondeur de Régime
la rivière d'eau (m) d'écouleme plafond
mouillé m2 I(m/m) hydraulique Froude- D'écoulement
nt (m/s) B
Dh Fr
2000 2.63 5.16 254.94 0.011 2.3 110.62 1.09 torrentielle
4000 2.26 4.33 304.4 0.017 1.24 245.04 1.24 torrentielle
6000 2.8 4.98 264.59 0.010 2.34 113.16 1.04 torrentielle
8000 2.09 5.55 237.36 0.021 1.54 154.51 1.43 torrentielle
10000 3.67 6.94 189.66 0.017 2.52 75.25 1.4 torrentielle
12000 4.79 4.94 266.35 0.009 2.48 107.55 1 critique
14000 1.48 6.27 209.84 0.033 1.32 158.97 1.74 torrentielle

Torrentielle
15400.11 3.89 5.52 238.44 0.009 3.07 77.66 1.01
torrentielle
18000 1.6 5.3 248.55 0.021 1.43 173.61 1.41
20264.2 4.7 6.39 205.99 0.012 3.05 67.64 1.17 torrentielle
22000 7.76 7.37 178.67 0.012 3.94 45.36 1.19 torrentielle
24000 4.1 5.02 262.31 0.009 2.56 102.32 1 critique
25470.09 2.95 6.27 210.12 0.018 2.14 98.33 1.37 torrentielle
28000 1.44 3.66 359.66 0.011 1.36 265.29 1 critique
30000 3.81 1.1 1193.72 0.0003 3.37 354.06 0.19 fluvial
32000 2.58 2.42 544.76 0.002 2.35 231.44 0.5 fluvial
33476.13 3.12 4.97 264.9 0.017 1.56 169.89 1.27 torrentielle
36000 2.12 5.29 248.72 0.017 1.7 146.52 1.3 torrentielle
38000 4.82 5.24 251.49 0.008732 2.78 90.59 1 critique
Nous avons observé une élévation importante de la hauteur d'eau dans la station(22000),
ainsi que augmentation de la vitesse d'écoulement, donc il y a un risque de déboulement du
lit de la rivière dans les zones qui situe dans cette station.(Tableau V-8)
121
Les zone qui situe dans les stations dans lequel se produit un changement brusque de la
pente, la vitesse et surélévation brusque de la ligne d'eau dans le cancel qui donné en
position fixé de longe du cancel ,subissent des débordement du niveau d'eau dans la rivière
(voire l'annexe A-tableau 10) qui contient avec les stations interpolées.
Ce de débordement peut être expliqué par un phénomène hydraulique comme il s'agit du

ressaut hydraulique lorsque le régime d'écoulement se transforme de régime super
critique(torrentiel) vers le régime subcritique (fluvial).
Le ressaut hydraulique est une surélévation brusque de la surface libre d‘un écoulement
permanant, qui se produit lors du passage du régime torrentielle au régime fluvial. Il est
accompagné d‘une agitation marquée et de grandes pertes d‘énergie. L‘écoulement torrentiel
à l‘amont du ressaut et l'effet de la gravité sur l'écoulement d'eau dans la rivière est
caractérisé par un coefficient cinétique (ou nombre de Froude Fr>1) . ce représente la
relation entre la vitesse de l‘écoulement et la célérité de propagation des onde de surface
libre.
Le phénomène de ressaut hydraulique est provoque un érosion du fond dans la rivières . les
particule solide provenant de l'amont peuvent se déplacer à l'avale où elles s'accumule et
changeant ainsi la morphologie du cours d’eau.
L'accumulation des particules solide en l'avale de la rivière produit une élévation du niveau
d'eau et elle est capable d'apparition un débordement en aval.
V.7. Conclusion
La réalisation et l’étude de la modélisation hydraulique du tronçon d’Djendjen qui coule
près du centre rural de Texenna permet de faire le diagnostic de cette zone d’étude pour
tout événement hydraulique. Cette simulation consiste à calculer les hauteurs d’eau et
déterminer les zones susceptibles d'être inondables en utilisant le logiciel HEC-GEORAS
qui couplé par le Système d’Information géographique « SIG » et le logiciel de modélisation
de rivières HEC-RAS qui permettent de reproduire l’état de la zone naturelle à l’état
numérique.
Les résultats de simulation hydraulique montre que le débordement du lit de cour dépond
par les deux phénomènes hydraulique qui sont le ressaut hydraulique et l'érosion qui
provoque le charriage des particules.
122
A travers ce mémoire, nous avons fait une analyse critique des différentes méthodes
d'estimation des crues dont dispose l'ingénieur pour estimer les crues de projet, qui ont
pour base les observations pluviométrique et/ou hydrométrique dans le bassin versant
d'Oued Djendjen.
L'étude hydrologique de la zone d'étude nous a permis d'extrait les différentes

caractéristiques géomorphologiques pour mieux comprendre le comportement
hydrologique du bassin versant.
L'estimation des crues de projet au bassin nous a montré une différence des résultats
obtenus en termes de convergence ou en termes de divergence. Cette différence est due
d’un côté aux données de départ utilisées ainsi qu’aux hypothèses et fondement de
l’approche utilisée.
Les méthodes statistiques (Gumbel et Gradex) sont des méthodes qui ont donnés une
fiabilité aux résultats car elles prennent en considération les données historiques
(données hydrométrique et donnés pluviométriques).
Compte tenu les différents problèmes existants dans la modélisation hydrologique et

indisponibilité des donnée sur le bassin étudié, les méthodes régionales et analogiques,
sont les méthodes adoptées dans ce cas. Les méthodes régionales les plus couramment
utilisées sont les méthodes empiriques.
Les différentes formules de l'approche empirique appliquée au bassin versant d'Oued

Djendjen ont données des valeurs très faibles comparant avec les méthodes statistiques
car elles sont basées sur des données régionales de la zone d'étude.
Les méthodes basées sur la modélisation pluies-débits nous ont également donné des
résultats assez fiables et qui se rapproche le plus des résultats statistiques.
Les résultats fournis par l'application de ces méthodes dépendent essentiellement des
données disponibles au départ, mais il faut insister sur la qualité (fiabilité) de ces données.
Parfois on est obligé d'appliquer des formules empiriques pour l'estimation des crues,
donc on recommande pour choisir les méthodes basées la modélisation pluie-débit ou les
123
méthodes statistique car elles prennent en considération la complexité du phénomène de
l’écoulement, en effet les méthodes statistiques sont des régressions qui synthétisent les
débits fréquentiel soit à partir des débits historiques, soit à partir des donnée historiques
de la pluviométrie. Cependant les modèles pluies débits quant à eux, permettent de
reproduire le comportement du bassin versant étudié face à une pluie donnée, ceci en
prenant en considération la géomorphologie du bassin versant, mais l’estimation des
paramètres nécessite la disponibilité de données certain nombre d’épisodes pluvieux. Ce
qui est intéressant est de souligner que pour les bassins versants non jaugés de la région,
une méthode d’analogie hydrologique peut être envisagée pour estimer les débits de crues
à partir du modèle hydrologique validé pour la région de Djendjen.
En fin de compte, nous avons terminé notre travail par la simulation hydraulique à l'aide
d'un outil de modélisation HEC-RAS en collaboration avec le Système d’Information
géographique « SIG » ARCGIS, afin d'extrait les différents paramètres hydrauliques
concernant Oued Djendjen, et d'indiquer au différentes zones susceptibles d'être
inondables qui se mis la région sous risque d'inondation.
124
ACREMAN, M., DUNBAR,M.,HANNAFORD, J.,MOUNTFORD,O., WOOD, P., HOLMES, N.et
KING, J. (2008), "Développions environnemental standards for abstractions frome UK
Rivers to implémente the eu water frame Works Directive/Développement de
standards environnementaux sur les prélèvements d'eau en rivière au royaume uni
pour la mise en œuvre de la directive cadre sur l'eau de l'Union Européenne ",
Hydrological Sciences Journal, 53(6), pp.1105-1120.
ANCTIL ,F., ROUSSELLE, J, AND NICOLAS LAUZON. Hydrologie: Cheminements de l'eau.

Presses inter Polytechnique, 2012
Arnold, J. G., & Allen, P. M. (1999). Automated methods for estimating baseflow and ground
water recharge from streamflow records. JAWRA Journal of the American Water
Resources Association, 35(2), 411-424
ASNOUNI, F. : "Etude du transport solide en suspension dans le bassin versant d'oued AL

ABD ".Thése de dactorat, 2014.
Bachi, M. : " PROBLEMATIQUE DU RISQUE INONDATION EN MILIEU URBAIN ; CAS DE

L’AGGLOMERATION DE SIDI BEL ABBES". Mémoire de Magister, Université
AboubakrBelkaïd-Tlemcen, 2011.
BENMIA ,K .:" Evaluation de la performance des barrages de protection contre les

inondations cas de la ville de ghazaouet". Mémoire magister université aboubakr
belkaïd , Tlemcen ,2012
BENNIS ,S.: "Hydraulique et hydrologie". troisième édition, Presse de l’université du Québec,

Canada, 2014.
BENNIS ,S.:"Hydraulique et hydrologie". Deuxième édition, Presse de l’université du

Québec, Canada, 2009.
BERDANE,KH. :"Crues et inondations en zone aride : cas du bassin de l’Oued Deffa dans le
Sud-ouest Algérien". Mémoire de Magister, Ecole Nationale Supérieure D’hydraulique
Arbaoui Abdallah, 2015.
BOIS,P.H.: "hydrologie générale". Ecoule nationale supérieure d'hydraulique et de

mécanique de Grenoble,2000.
BOUCENNA, F. :''Cartographie par les différentes méthodes de vulnérabilité a la pollution

d’une nappe côtière cas de la plaine alluviale de l’Oued Djedjen( Jijel, Nord-est
Algérien )".Mémoire de Magister, Université d'Annaba, 2009.
BOUDJEDJOU, L. : '' Etude de la flore adventice des cultures de la région de Jijel". Mémoire
de Magister, Université de Sétif, 2010.
125
BOULGHOBRA, N. : "Protection de la ville de Skikda contre l’inondation. essai de PPR".
Mémoire de Magister Aménagement Des Milieux Physiques, Université De Batna,
2006.
BOUTAGHANE ,H ., OUERDACHI ,L ., MANSOURI ,R .:" Contribution a l’étude du

dimensionnement des ouvrages d’assainissement pluvial : les méthodes des
hydrogramme". Revue scientifique semestrielle et technique, Le Journal de l’Eau et de
l’Environnement,1945
Brocca, L., Melone, F., Moramarco,T.:"Distributed rainfall-runoff modelling for flood

frequency estimation and flood forecasting". Hydrol, Process, 25, 2801–2813.
https://doi.org/10.1002/hyp.8042
CHOW, V. T., MAIDMENT, D. R, and MAYS, L. W. "Applied hydrology, 572 pp." Editions
McGraw-Hill, New York ,1988.
COMITE FRANÇAIS DES GRANDE BARRAGES .:"barrage et réservoir: méthode du gradex.

N°2,novembre,1994.
CORRADINI C, MORBIDELLI R, SALTALIPPI C, MELONE F. 2002.:" An adaptive model for flood

forecasting on medium size basins. In AppliedSimulation and Modelling", Ubertini L
(ed). IASTED Acta Press:Anaheim (CA); 555–559.
DAVIE, T.:"Fondametals of hydrologie". Deuxièmeédition,Rorutledge Fundamentals Of

Physical Geography,2002.
DEGOUTTE, G. : ''DiagnostiqueAménagement St Gestion Des Rivières''. hydrauliquement

morphologie fluviales appliquées, Deuxième édition, Lavoisier, Paris, 2012.
DUBREUIL,P.: "Initiation à l'analyse hydrologique".O.R.S.O.M, 24 rue bayard, paris VIIIe,

1974.
FENAZI, B. et DERRADJI, E. F. : ''La qualité chimique des eaux du bassin versant de l’oued
Djendjen (w. Jijel)''. 2016.
HINGRAY, B et Picouet.C et Musy.A.:" Hydrologie: Une science pour l'ingénieur ". Vol. 21,
PPUR presses polytechniques, 2009.
HORTON, R.E.:'' Drainage‐basin characteristics''. Eos, transactions American geophysical

union, 13(1), pp. 350-361, 1932.
IRATNI, N. : "Modélisation hydrologique de quelques sous bassins versant de la Tafna",

Mémoire de Magister, Université des Sciences et de la Technologie d’Oran-Mohamed
BOUDIAF, 2014.
LABORDE, J.P. : '' Elément d’hydrologie de surface''. Centre De Recherche Scientifique,

Université De Nice-Sophia, 2000.
LABORDE, J.P. : "Eléments d'hydrologie de surface". Ecole Polytechnique De L’universite De

Nice - Sophia Antipolis, 2009.
126
LAHLAH, S. : "Actes des Journées Techniques/ Risques Naturels : Inondation, Prévision,
Protection /Batna", 2004
LAHLAH, S. : " Introduction sur l’étude des inondations en Algérie". Agence Nationale des
Ressources Hydraulique (ANRH), Alger,2000.
Lang, C. : "Étiages et tarissements : vers quelles modélisations ? l’approche conceptuelle et

l’analyse statistique en réponse a la diversité spatiale des écoulements en étiage des
cours d’eau de l’est français. hydrologie". Université de Metz, 2007.
Lee, K . T, and Ho, J.Y. "Design hyetograph for typhoon rainstorms in Taiwan." Journal of
Hydrologic Engineering VOL 13, no 7,P .647-651,2008.
MCYLAN, P. , FAVOE, A.C , MUSY, A. : ''Hydrologie Fréquentielle : Une Science Prédictive''.

Presses Polytechniques et universitaires soudants, 2008.
MERABET, A. : "Etude de la protection de la ville de sidi bel abbés contre les inondations".
Mémoire de Magister, Université de Djilali Liabes-Sidi Bel Abbés, 2006.
MEYLAN,P et MUSY.A.: "hydrologie fréquentielle".H.G.A, bucaraest.1999.
MONTOROI, J.P."Rôle des sols sur la genèse des inondations." Publications du symposium
européen sur les problèmes actuels de la protection contre les inondations ,paris,2012
MUSY, A. et HIGY, CH. :"Hydrologie : une science de la nature". Presses Polytechniques et

Universitaires Romandes,2004.
MUSY,A." Hydrologie appliquée". 1998-001. Edition H.G.A, Bucarest ,1998.
OUARDA, T. B. M. J., GINGRAS, H., BOBEE, B.ET LEMONIER, M." Synthése des méthodes
simples de régionalisation" .Rapport de recherche No R-600, INRS-Eau, Québec,
Novembre 2001.
PÉGUY, CH-P. "Jeux et enjeux du climat ". Elsevier Masson, 1989.
QADEM, A. : "Quantification, modélisation et gestion de la ressources en eau dans le bassin

versant haut Sebou ".Thèse de Doctorat, Université De Sidi Mohamed Ben Abdallah
et Université De Lorraine, Maroc,2015.
RAGHUNATH ,H.M. Hydrology: principles, analysis and design. New Age International, 2006.
REMINI, B. : ''Quelques Applications en Hydraulique ''. Office de Publications Universitaires,

2002.
Roche, M. : "Hydrologie de surface". Gauthier-Villars Editeur, Paris,1963.
ROCHE, M. :"Les processus du raisonnement chez l’hydrologue", Bulletin De La Section

Géographique, Tome LXXX, 1969.
ROCHE, P.Q. , MIQUEL, J. , et GAUME, E. : ''Hydrologie Quantitative : Processus, Modèles et

Aides a la décision''. Spring-valag France, 2012.
127
SALOMON, J.N. "L’homme face aux crues et aux inondations". Presses Universitaire de
Bordeaux et Université Michel de Montaigne – Bordeaux, 1997.
SALOMON, J.N. : "L'Homme face aux crues et aux inondations". Université De Michel De
Montaigne, Bordeaux-3, 1997.
SATIN,M .,SELMI ,B.:"guide technique de l'assainissement". deuxième édition. presse

polytechniques et universitaires romandes,2008.
SCHUMM,S.A. : '' Evolution of drainage systems and slopes in badlands at Perth Amboy, New
Jersey''. Geological society of America bulletin, 67(5), pp.597-646 , 1956.
SCS,. SCS National Engineering Handbook. Section 4, Hydrology, in: National Engineering
Handbook,1972
SPPPI (2004) : " Prise En Compte Du Risque Inondation Dans Les Etudes De Dangers Des
Installations Classées ", Cahier technique N°8, DRIRE Haute-Normandie - 21, avenue
de la Porte des champs - 76037 ROUEN CEDEXT.
STRAHLER, A. N.: ''Quantitative geomorphology of drainage basin and channel

networks''. Handbook of applied hydrology, 1964.
TAHAR, S. : "Impact des inondations sur l'espace urbain ' le cas de la wilaya de Sidi Bel Abbes
'". Mémoire Magister, Université d'Oran, 2013.
TOUABIA, B. : "Manuel pratique d’hydraulique". Institut de l’ENSH, Blida, 2005.
US Army Corps of Engineers.: " HEC-RAS River Analysis System : User's Manual ". Manuel
d’utilisation de HEC-RAS, Hydrologic Engineering Center,2010
US Army Corps of Engineers.: "HEC-GeoRAS Extension for support of HEC-RAS using ArcView
:User's Manual ", Manuel d’Utilisation de HEC-GeoRAS, Hydrologic Engineering
Center, 2002
YAHIAOUI, A. :"Inondations torrentielles cartographie des zones vulnérables en algérie du

nord ( cas de l'Oued Mekerra, Zilaya de Sidi Bel Abbés)". Thèse de Doctorat, Ecole
Nationale Polytechnique,Université De Bechar, 2012.
128
Annexe A
A
Tableau 1: les données pluviométrique(1914-2010)
Ordre X(Km) Y(Km) Code Station Année P max,j Ordre X(Km) Y(Km) Code Station Année P max,j
1 776.2 377.85 30302 TEXENNA 1914 90.5 33 776.2 377.85 30302 TEXENNA 1947 111
2 776.2 377.85 30302 TEXENNA 1915 75.2 34 776.2 377.85 30302 TEXENNA 1948 70.4
3 776.2 377.85 30302 TEXENNA 1916 70.5 35 776.2 377.85 30302 TEXENNA 1949 47.7
4 776.2 377.85 30302 TEXENNA 1917 78.5 36 776.2 377.85 30302 TEXENNA 1950 76.5
5 776.2 377.85 30302 TEXENNA 1918 60 37 776.2 377.85 30302 TEXENNA 1951 63.7
6 776.2 377.85 30302 TEXENNA 1919 92 38 776.2 377.85 30302 TEXENNA 1952 82
7 776.2 377.85 30302 TEXENNA 1920 72.2 39 776.2 377.85 30302 TEXENNA 1954 70
8 776.2 377.85 30302 TEXENNA 1921 90 40 776.2 377.85 30302 TEXENNA 1955 56.3
9 776.2 377.85 30302 TEXENNA 1922 55.3 41 776.2 377.85 30302 TEXENNA 1968 28
10 776.2 377.85 30302 TEXENNA 1923 58.5 42 776.2 377.85 30302 TEXENNA 1972 135.3
11 776.2 377.85 30302 TEXENNA 1924 70.6 43 776.2 377.85 30302 TEXENNA 1973 83
12 776.2 377.85 30302 TEXENNA 1925 93.4 44 776.2 377.85 30302 TEXENNA 1975 93.2
13 776.2 377.85 30302 TEXENNA 1926 121.8 45 776.2 377.85 30302 TEXENNA 1977 58.5
14 776.2 377.85 30302 TEXENNA 1927 89.8 46 776.2 377.85 30302 TEXENNA 1978 99.8
15 776.2 377.85 30302 TEXENNA 1929 89 47 776.2 377.85 30302 TEXENNA 1979 45.3
16 776.2 377.85 30302 TEXENNA 1930 159.8 48 776.2 377.85 30302 TEXENNA 1980 80.6
17 776.2 377.85 30302 TEXENNA 1931 119 49 776.2 377.85 30302 TEXENNA 1981 48.2
18 776.2 377.85 30302 TEXENNA 1932 70.7 50 776.2 377.85 30302 TEXENNA 1985 64.4
19 776.2 377.85 30302 TEXENNA 1933 73.5 51 776.2 377.85 30302 TEXENNA 1986 58.3
20 776.2 377.85 30302 TEXENNA 1934 75.5 52 776.2 377.85 30302 TEXENNA 1990 152.8
21 776.2 377.85 30302 TEXENNA 1935 55 53 776.2 377.85 30302 TEXENNA 1998 79.5
22 776.2 377.85 30302 TEXENNA 1936 133 54 776.2 377.85 30302 TEXENNA 1999 80
23 776.2 377.85 30302 TEXENNA 1937 81.7 55 776.2 377.85 30302 TEXENNA 2000 96.4
24 776.2 377.85 30302 TEXENNA 1938 61 56 776.2 377.85 30302 TEXENNA 2001 60.8
25 776.2 377.85 30302 TEXENNA 1939 64.6 57 776.2 377.85 30302 TEXENNA 2002 117.7
26 776.2 377.85 30302 TEXENNA 1940 79.5 58 776.2 377.85 30302 TEXENNA 2003 100.7
27 776.2 377.85 30302 TEXENNA 1941 60 59 776.2 377.85 30302 TEXENNA 2004 103.2
28 776.2 377.85 30302 TEXENNA 1942 57 60 776.2 377.85 30302 TEXENNA 2005 84.4
29 776.2 377.85 30302 TEXENNA 1943 48 61 776.2 377.85 30302 TEXENNA 2006 98.7
30 776.2 377.85 30302 TEXENNA 1944 62.5 62 776.2 377.85 30302 TEXENNA 2007 142.6
31 776.2 377.85 30302 TEXENNA 1945 76 63 776.2 377.85 30302 TEXENNA 2008 165
32 776.2 377.85 30302 TEXENNA 1946 78 64 776.2 377.85 30302 TEXENNA 2009 140
65 776.2 377.85 30302 TEXENNA 2010 96
129
Tableau 2 : Donnée hydrométrique de 21 années (1972-19972)
Débits de Débits de
Débits
pointe pointe Coefficient de
Année maximaux
maximale maximale point
journalier
journalière journalière
1972 33.300 158.500 92.4 4.76
1973 75.356 206.500 107.8 2.74
1974 69.011 371.200 125.6 5.38
1975 84.428 395.800 141.8 4.69
1976 57.655 258.800 158.5 4.49
1977 19.388 125.600 160 6.48
1978 33.473 182.500 160 5.45
1979 49.859 209.800 170.5 4.21
1980 84.108 455.000 182.5 5.41
1981 51.519 257.200 190 4.99
1982 45.082 163.000 193.3 3.62
1983 60.148 141.800 209.8 2.36
1984 163.657 496.000 253.4 3.03
1985 198.515 552.500 257.2 2.78
1986 51.787 190.000 371.2 3.67
1987 43.845 88.000 395.8 2.01
1988 21.518 160.000 455 7.44
1989 110.570 160.000 517.5 1.447
1990 64.287 223.000 552.5 3.469
1991 27.180 78.000 776 2.870
1992 42.596 193.300 800.8 4.538
130
Tableau 3: pluie de projet
Double Répartition
T(h) Type I Type II Type IA Type III
triangle Aléatoire
0 0 0 0 0 0.00 0.00
1 1.7 1.1 2.2 1 0.52 0.85
2 1.8 1.2 2.9 1 1.05 2.22
3 2 1.2 3.2 1.2 1.57 2.22
4 2.1 1.3 3.3 1.1 2.09 2.22
5 2.3 1.6 4 1.4 2.62 1.79
6 2.6 1.6 4.8 1.5 3.14 0.51
7 3.1 2 6.4 1.7 3.67 0.51
8 3.8 2 15.7 2.6 4.19 0.51
9 6 2.7 9.5 3.3 4.71 0.51
10 26.1 3.4 5.7 4.1 5.24 1.51
11 10.9 5.5 4.6 6.1 11.87 2.58
12 5.8 42.7 4.1 35 19.53 2.58
13 4.6 11.3 3.7 15.1 11.87 2.58
14 3.8 4.9 3.5 6 5.24 4.94
15 3.3 3.1 3.3 3.7 4.71 7.83
16 3.1 2.5 3.1 3.8 4.19 1.69
17 2.7 2.2 3 1.8 3.67 14.49
18 2.5 1.9 2.8 1.8 3.14 12.56
19 2.3 1.6 2.6 1.7 2.62 11.17
20 2.1 1.5 2.4 1.8 2.09 7.52
21 2 1.2 3.8 1.1 1.57 12.08
22 1.9 1.2 1.9 1.1 1.05 0.00
23 1.8 1.2 1.8 1 0.52 7.15
24 0.017 1.1 1.7 1.1 0.00 0.00
131
Tableau 4: Application du fonction de production ( SCS-CN) au pluies synthétiques avec les la répartitions
temporelle de type I
Pluie de Pluie Ruissellement Ruissellement

temps type I%
projet cumulé cumulé (P nette)
0.000 0.000 0.000 0.000 0.000 0.000
1.000 1.700 2.976 2.976 0.115 0.115
2.000 1.800 3.151 6.127 0.519 0.404
3.000 2.000 3.501 9.628 1.263 0.744
4.000 2.100 3.676 13.305 2.331 1.067
5.000 2.300 4.026 17.331 3.787 1.456
6.000 2.600 4.552 21.883 5.740 1.953
7.000 3.100 5.427 27.309 8.426 2.686
8.000 3.800 6.652 33.962 12.154 3.728
9.000 6.000 10.504 44.465 18.803 6.649
10.000 26.100 45.691 90.156 54.032 35.229
11.000 10.900 19.082 109.237 70.424 16.392
12.000 5.800 10.153 119.391 79.380 8.956
13.000 4.600 8.053 127.444 86.576 7.196
14.000 3.800 6.652 134.096 92.575 5.998
15.000 3.300 5.777 139.873 97.819 5.244
16.000 3.100 5.427 145.300 102.773 4.954
17.000 2.700 4.727 150.026 107.108 4.335
18.000 2.500 4.377 154.403 111.137 4.029
19.000 2.300 4.026 158.429 114.856 3.719
20.000 2.100 3.676 162.106 118.262 3.406
21.000 2.000 3.501 165.607 121.513 3.252
22.000 1.900 3.326 168.933 124.609 3.096
23.000 1.800 3.151 172.084 127.549 2.939
24.000 1.700 2.976 175.060 130.330 2.781
132
temporelle de type I-A

temps type I-A%
0 0 0 0 0 0
1 2.200 3.851 3.851 0.200 0.200
2 2.900 5.077 8.928 1.092 0.892
3 3.200 5.602 14.530 2.744 1.652
4 3.300 5.777 20.307 5.030 2.286
5 4.000 7.002 27.309 8.426 3.396
6 4.800 8.403 35.712 13.203 4.777
7 6.400 11.204 46.916 20.464 7.261
8 15.700 27.484 74.401 41.076 20.612
9 9.500 16.631 91.031 54.769 13.693
10 5.700 9.978 101.010 63.277 8.508
11 4.600 8.053 109.062 70.271 6.993
12 4.100 7.177 116.240 76.586 6.315
13 3.700 6.477 122.717 82.343 5.757
14 3.500 6.127 128.844 87.835 5.492
15 3.300 5.777 134.621 93.050 5.215
16 3.100 5.427 140.048 97.979 4.928
17 3.000 5.252 145.300 102.773 4.795
18 2.800 4.902 150.201 107.269 4.496
19 2.600 4.552 154.753 111.460 4.191
20 2.400 4.201 158.954 115.342 3.882
21 3.800 6.652 165.607 121.513 6.171
22 1.900 3.326 168.933 124.609 3.096
23 1.800 3.151 172.084 127.549 2.939
24 1.700 2.976 175.060 130.330 2.781
133
temporelle de type II

temps type II%
0.000 0 0.000 0.000 0.000 0.000
1.000 1.100 1.926 1.926 0.043 0.043
2.000 1.200 2.101 4.026 0.220 0.176
3.000 1.200 2.101 6.127 0.519 0.300
4.000 1.300 2.276 8.403 0.971 0.451
5.000 1.600 2.801 11.204 1.687 0.717
6.000 1.600 2.801 14.005 2.564 0.876
7.000 2.000 3.501 17.506 3.856 1.293
8.000 2.000 3.501 21.007 5.342 1.485
9.000 2.700 4.727 25.734 7.610 2.269
10.000 3.400 5.952 31.686 10.830 3.220
11.000 5.500 9.628 41.314 16.724 5.893
12.000 42.700 74.751 116.065 76.431 59.708
13.000 11.300 19.782 135.847 94.161 17.729
14.000 4.900 8.578 144.425 101.972 7.812
15.000 3.100 5.427 149.851 106.947 4.975
16.000 2.500 4.377 154.228 110.976 4.029
17.000 2.200 3.851 158.079 114.532 3.557
18.000 1.900 3.326 161.405 117.612 3.080
19.000 1.600 2.801 164.206 120.212 2.599
20.000 1.500 2.626 166.832 122.653 2.441
21.000 1.200 2.101 168.933 124.609 1.956
22.000 1.200 2.101 171.034 126.568 1.959
23.000 1.200 2.101 173.134 128.530 1.961
24.000 1.100 1.926 175.060 130.330 1.800
134
temporelle de type III
Pluie synthétique Pluie de Pluie Ruissellement Ruissellement

Temps
type III% projet cumulé cumulé (P nette)
0 0 0 0 0 0
1 1 1.751 1.751 0.035 0.035
2 1 1.751 3.501 0.163 0.129
3 1 2.101 5.602 0.434 0.270
4 1 1.926 7.528 0.782 0.349
5 1 2.451 9.978 1.353 0.571
6 2 2.626 12.604 2.107 0.754
7 2 2.976 15.580 3.120 1.013
8 3 4.552 20.132 4.954 1.834
9 3 5.777 25.909 7.699 2.746
10 4 7.177 33.086 11.639 3.940
11 6 10.679 43.765 18.335 6.696
12 35 61.271 105.036 66.761 48.426
13 15 26.434 131.470 90.201 23.440
14 6 10.504 141.974 99.734 9.532
15 4 6.477 148.451 105.661 5.927
16 4 6.652 155.103 111.783 6.122
17 2 3.151 158.254 114.694 2.911
18 2 3.151 161.405 117.612 2.918
19 2 2.976 164.381 120.374 2.762
20 2 3.151 167.532 123.305 2.931
21 1 1.926 169.458 125.099 1.794
22 1 1.926 171.384 126.895 1.796
23 1 1.751 173.134 128.530 1.635
24 1 1.926 175.060 130.330 1.800
135
Tableau 8 : Application du fonction de production( SCS-CN) au pluies synthétiques avec les la
répartitions temporelle du double triangle.
Pluie
temps synthétique
type DT%
0 0.000 0.000 0.000 0.000 0.000
1 0.524 0.917 0.917 0.006 0.006
2 1.047 1.833 2.750 0.097 0.091
3 1.571 2.750 5.500 0.418 0.321
4 2.094 3.667 9.166 1.149 0.731
5 2.618 4.583 13.749 2.478 1.329
6 3.142 5.500 19.249 4.573 2.096
7 3.665 6.416 25.666 7.576 3.002
8 4.189 7.333 32.999 11.588 4.012
9 4.712 8.250 41.248 16.681 5.093
10 5.236 9.166 50.415 22.897 6.216
11 11.873 20.785 71.200 38.527 15.630
12 19.533 34.195 105.394 67.073 28.546
13 11.873 20.785 126.179 85.442 18.369
14 5.236 9.166 135.346 93.707 8.265
15 4.712 8.250 143.595 101.215 7.508
16 4.189 7.333 150.929 107.937 6.723
17 3.665 6.416 157.345 113.854 5.916
18 3.142 5.500 162.845 118.948 5.094
19 2.618 4.583 167.428 123.208 4.260
20 2.094 3.667 171.094 126.625 3.417
21 1.571 2.750 173.844 129.193 2.568
22 1.047 1.833 175.678 130.908 1.714
23 0.524 0.917 176.594 131.766 0.858
24 0.000 0.000 176.594 131.766 0.000
136
Tableau 9: Application du fonction de production ( SCS-CN) au pluies synthétiques avec une la répartition
temporelle Aléatoire du pluie.
répartition Pluie de Pluie Ruissellement Ruissellement

temps
Aléatoire% projet cumulé cumulé (P nette)
0 0.000 0.000 0.000 0.000 0.000
1 0.851 1.490 1.490 0.023 0.023
2 2.221 3.888 5.378 0.399 0.376
3 2.221 3.888 9.266 1.173 0.774
4 2.221 3.888 13.153 2.282 1.108
5 1.792 3.138 16.291 3.385 1.103
6 0.507 0.887 17.178 3.727 0.342
7 0.507 0.887 18.065 4.081 0.355
8 0.507 0.887 18.952 4.448 0.367
9 0.507 0.887 19.840 4.826 0.378
10 1.507 2.637 22.477 6.017 1.190
11 2.575 4.508 26.985 8.256 2.239
12 2.575 4.508 31.494 10.721 2.465
13 2.575 4.508 36.002 13.379 2.658
14 4.940 8.648 44.650 18.927 5.548
15 7.830 13.707 58.357 28.654 9.727
16 1.691 2.960 61.317 30.873 2.219
17 14.494 25.372 86.689 51.130 20.256
18 12.557 21.983 108.672 69.930 18.800
19 11.174 19.562 128.234 87.287 17.357
20 7.520 13.165 141.399 99.210 11.923
21 12.078 21.144 162.543 118.667 19.458
22 0.000 0.000 162.543 118.667 0.000
23 7.150 12.517 175.060 130.330 11.662
24 7.150 12.517 187.577 142.079 11.749
137
Tableau 10 : principaux paramètres hydrauliques pour un ensemble de stations
Vitesse d'écoule-
La hauteur d'eau
Surface mouillé
La profondeur
débit de point
D'ecoulement
pente I(m/m)
hydraulique
Nomber de
plafond B
Largeur au
Froude-Fr
River Sta
Régime
Profile
(M /s)
(m/s)
ment
(m)
DH
M
3
2
2000 PF 1 1316.69 2.63 5.16 254.94 0.011 2.3 110.62 1.09 supercritique
2500.00* PF 1 1316.69 2.83 4.81 273.77 0.010 2.16 126.55 1.04 supercritique
3500.00* PF 1 1316.69 2.84 3.78 348.56 0.009 1.58 220.23 0.96 supcritique
12000 PF 1 1316.69 4.79 4.94 266.35 0.009 2.48 107.55 1 Critique
12500.0* PF 1 1316.69 3.98 4.24 310.6 0.006 2.55 121.65 0.85 Supcritique
13000.0* PF 1 1316.69 2.93 4.61 285.49 0.010 2.13 134.1 1.01 Supercritique
13500.0* PF 1 1316.69 2.57 4.13 319.02 0.008 2.07 154.22 0.92 Supcritique
14000 PF 1 1316.69 1.48 6.27 209.84 0.033 1.32 158.97 1.74 Supercritique
15400.11 PF 1 1316.69 3.89 5.52 238.44 0.009 3.07 77.66 1.01 Supercritique
15833.4* PF 1 1316.69 4.11 4.05 325.47 0.004 3.29 98.81 0.71 subcritique
16266.7* PF 1 1316.69 3.2 4.5 292.52 0.007 2.61 111.96 0.89 subcritique
16700.0* PF 1 1316.69 3.08 3.95 333.08 0.006 2.55 130.56 0.79 subcritique
17133.3* PF 1 1316.69 2.63 4.04 326.21 0.007 2.27 143.62 0.86 subcritique
17566.6* PF 1 1316.69 2.47 3.74 352.19 0.006 2.17 161.95 0.81 subcritique
20264.2 PF 1 1316.69 4.7 6.39 205.99 0.012 3.05 67.64 1.17 supercritique
138
Suite du tableau 10
La
Vitesse profonde Largeur Nombe
débit de Surface pente
La hauteur d'écoule- ur au r de Régime
River Sta Profile point mouillé
d'eau (m) ment I(m/m) hydrauliq plafond Froude D'ecoulement
(M3/s) M2
(m/s) ue B -Fr
DH
24000 PF 1 1316.69 4.1 5.02 262.31 0.009 2.56 102.32 1 Critique
24490.0* PF 1 1316.69 4.33 4.96 265.34 0.008 2.65 100.22 0.97 Supcritique
24980.0* PF 1 1316.69 4.23 4.64 283.88 0.007 2.71 104.72 0.9 Supcritique
28000 PF 1 1316.69 1.44 3.66 359.66 0.011 1.36 265.29 1 Critique
28500.0* PF 1 1316.69 3.03 1.58 833.5 0.001 2.75 303.18 0.3 Supcritique
29000.0* PF 1 1316.69 3.39 1.34 981.87 0.001 2.98 329.93 0.25 Supcritique
29500.0* PF 1 1316.69 3.63 1.2 1100.49 0.0001 3.17 346.75 0.21 Supcritique
30000 PF 1 1316.69 3.81 1.1 1193.72 0.0003 3.37 354.06 0.19 Supcritique
30500.0* PF 1 1316.69 2.99 1.57 840.65 0.001 2.65 316.67 0.31 Supcritique
31000.0* PF 1 1316.69 2.51 2.07 637.34 0.002 2.24 284.55 0.44 Supcritique
31500.0* PF 1 1316.69 2.47 2.3 573.67 0.002 2.23 257.37 0.49 Supcritique
32000 PF 1 1316.69 2.58 2.42 544.76 0.002 2.35 231.44 0.5 Supcritique
32492.0* PF 1 1316.69 1.73 5.05 260.79 0.018 1.53 170 1.3 Supercritique
36000 PF 1 1316.69 2.12 5.29 248.72 0.017 1.7 146.52 1.3 Supercritique
38000 PF 1 1316.69 4.82 5.24 251.49 0.010 2.78 90.59 1 Critique
139

M HDR - Urb 2018 09

Transféré par

Droits d'auteur :

Formats disponibles

M HDR - Urb 2018 09

Transféré par

Informations du document

Titre original

Copyright

Formats disponibles

Partager ce document

Partager ou intégrer le document

Options de partage

Avez-vous trouvé ce document utile ?

Ce contenu est-il inapproprié ?

Droits d'auteur :

Formats disponibles

M HDR - Urb 2018 09

Transféré par

Droits d'auteur :

Formats disponibles

République Algérienne Démocratique et Populaire

Ministère de l’Enseignement Supérieur et de la Recherche Scientifique

Département de Génie Civil et Hydraulique

MÉMOIRE DE FIN D’ÉTUDES

ETUDE DE LA PROBLEMATIQUE DU RISQUE

Soutenu publiquement le 28/ 06 / 2018 devant le jury composé de :

‫دراسة مخاطر الفيضانات هذه باستخدام برنامجي ‪ ArcGis :‬و ‪.HEC-RAS‬‬

‫تدفق‪ ،SCS-CN ،Gradex ،Gumbel ،‬هيدروغرام وحدوي‪ ،ArcGis HEC-RAS ،‬فيضان‬

Nous remercions nos collègues de la première promotion de Master en hydraulique à l’université de

A tous mes frères et sœurs.

A ceux qui m’aime et m’ont aidé sans qu’ils le sachent.

« Tu m’as donné la vie, la tendresse et le courage pour réussir.

À ma chère amie Amel

« Notre passé se ressemble,

Notre présent nous assemble

Notre futur c'est ensemble. »

À mes frère Abd-El Basset et Abd-Allah.

À tout mes amis.

TABLE DES MATIERES .......................................................................................................................................... VII

TABLE DES FIGURES ............................................................................................................................................... X

TABLE DES TABLEAUX......................................................................................................................................... XIII

INTRODUCTION GENERALE ................................................................................................................................. XV

I. NUMERO DU CHAPITRE ................................................................................................................................ 1

PROBLEMATIQUE &GENERALITES SUR LES INONDATIONS .................................................................................... 1

I.1. CYCLE HYDROLOGIQUE ....................................................................................................................................... 1

II. NUMERO DU CHAPITRE ...............................................................................................................................27

II. ETAT DE L'ART SUR LA MODELISATION DES INONDATION ...........................................................................27

II.1. METHODE STATISTIQUE ................................................................................................................................... 27

III. NUMERO DU CHAPITRE ...............................................................................................................................62

CHAPITRE III .........................................................................................................................................................62

III. ZONE D'ETUDE ET CARACTERISTIQUE MORPHOMETRIQUES ........................................................................62

III.1. PRESENTATION DE LA ZONE D'ETUDE .................................................................................................................. 62

IV. NUMERO DU CHAPITRE ...........................................................................................................................82

IV. ESTIMATION DU DEBIT DE POINTE...........................................................................................................82

IV.1. ANALYSE ET TRAITEMENT DES DONNEES .............................................................................................................. 82

V. NUMERO DU CHAPITRE .............................................................................................................................109

V. ETUDE HYDRAULIQUE ...............................................................................................................................109

V.1. L’OBJECTIF DE L’ETUDE HYDRAULIQUE .............................................................................................................. 109

CONCLUSION GENERALE ....................................................................................................................................123

Figure ‎I-2: Cycle de l'eau(Laborde,2000) ................................................................................................................ 2

Figure ‎I-3: Lit mineur d'un cours d'eau (Bachi, 2011)............................................................................................. 5

Figure ‎I-4: Lit majeur d’un cours d’eau (Bachi, 2011)............................................................................................. 6

Figure ‎I-5: Lit moyen du cours d'eau(Bachi, 2011) ................................................................................................. 6

Figure ‎I-6: Représentation de termes géographique usuels (Degoutte, 2012) ...................................................... 7

Figure ‎I-7: Exemple de courbe de tarissement (Musy, 1998)................................................................................. 8

Figure ‎I-8:Hydrogramme de crue. (Bennis, 2009) .................................................................................................. 9

Figure ‎I-10: Inondation de Plaine. (Merabet ,2006) ............................................................................................. 15

Figure ‎I-12: Inondation par crues torrentielles. (Merabat, 2006) ........................................................................ 17

Figure ‎I-14: Types des inondations ....................................................................................................................... 19

Figure ‎I-18: Inondations de bab el Oued - 2011 ................................................................................................... 25

Figure ‎I-19: Inondations de Ghardaïa - 2008 ........................................................................................................ 25

Figure ‎I-20: Inondations de Skikda (Les allées) - 1984.......................................................................................... 25

Figure ‎I-21: Inondations de Jijel (Chekfa) ............................................................................................................. 25

Figure ‎II-1:Courbe QDF(Hingray et al ,2009) ........................................................................................................ 31

Figure ‎II-5: Distribution de la Pluit nette par la méthode de horton.(Musy,1998) .............................................. 43

Figure ‎II-6:Méthode de l'indice φ.(Musy,1998).................................................................................................... 44

Figure ‎II-7:Méthode du SCS.(Musy,1998) ............................................................................................................. 46

Figure ‎II-9:Hypothèse de linéarité de l'Hydrogramme Unitaire (Musy, 1998) ..................................................... 54

Figure ‎II-11: La réponse impulsionnelle.(Hingray, 2009) ...................................................................................... 57