Hadj Ben Abdelmoula Malika
Hadj Ben Abdelmoula Malika
Hadj Ben Abdelmoula Malika
Thème
Modélisation hydrologique des écoulements superficiels du
bassin versant du Haut Cheliff.
Soutenu le :………………………………………..
PRESENTE PAR :
Mlle : HADJ BEN ABDELMOULA MALIKA.
Avant tout je remercie ALLAH tous puissant, de m’avoir guidé toutes ces années d’étude et de
m’avoir donné la volanté la patiente, le courage pour terminer ce travail.
Ainsi mon profond respect va à mes enseignants classe Master 2 Pour leurs conseils,
orientations et aides.
Nos sincères remerciements vont à l’endroit de toutes les membres de ma famille pour le
soutien moral même lors des moments difficiles. Ainsi que tous les personnes qui nous ont
aidés de près ou de loin à élaborer ce travail.
Enfin, merci également le président de jury et les membres qui feront l’honneur d’examiner
consacrer votre temps à évaluer ce mémoire ; trouvent l’expression de ma parfaite
considération. Merci pour tout.
.
DEDICACE
J e tien à dédié ce modeste travaille en signe de respect à ma famille qui sont toujours été ma
source d’inspiration pour leurs sacrifices et leurs soutiens dont ils ont fait preuve pour que je
réussisse tous le long de ma vie et mes études et les mots me manquent pour exprimer toute la
reconnaissance, la fierté et le profond amour que je vous porte.
Résumé
Le principal objectif de cette étude et de chercher une relation pluie-débit permettant la
reconstitution des débits à l’exutoire de bassin versant du Haut Chéliff. Cette étude consiste à
évaluer les conditions en particulier géographiques, géologiques et hydro-climatiques qui
influent sur l’écoulement dans le bassin d’étude.
Cet objectif ne pourra être atteint qu’à travers une modélisation hydrologique,
d’observer le comportement du bassin et de contrôler les événements hydro-climatiques en
déterminant les paramètres optimaux par la simulation hydrologique.
D’après les résultats de simulations obtenus par le modèle SCS-LR sur le bassin
versant du Haut Chéliff, on constate que le modèle pluie-débit distribué parcimonieux fournit
des résultats très satisfaisantes, avec des paramètres conservant une forte signification
physique, et se prête à des applications à vocation opérationnelle dans le domaine de la
prédétermination et la simulation des écoulements.
Les valeurs du critère de Nash obtenus pour les 22 épisodes sont en moyenne
supérieure à 80 %, indiquant une signification assez suffisante des paramètres du modèle, ou
ce dernier peut être interprété en fonction des caractéristiques géographiques, géologiques et
hydro-climatiques du bassin d’étude, ce qui traduit l’efficacité et la robustesse du logiciel
ATHYS sur le bassin versant Haut-Chéliff.
Abstarct :
The main objective of this study and to seek a rainfall-flow relationship allowing
reconstitution of flows at the outlet of the Haut Chéliff watershed. This study consists of
evaluating the geographical, geological and hydro-climatic conditions that influence flow in
the study basin.
This objective can only be achieved through hydrological modeling, observation of the
basin's behavior and control of hydro-climatic events by determining the optimal parameters
by hydrological simulation.
Based on the results of the SCS-LR model on the Haut Cheliff catchment, it is found
that the distributed rainfall-parsimonious model provides very satisfactory results, with
parameters retaining high physical significance, To operational applications in the field of
predetermination and simulation of flows.
The values of the Nash criterion obtained for the 22 episodes are on average greater
than 80%, indicating a rather sufficient significance of the parameters of the model, or the
latter can be interpreted according to the geographical, geological and hydro - climatic
characteristics of the basin, Study, which reflects the efficiency and robustness of the ATHYS
software on the Haut-Chéliff watershed.
ملخص:
الهدف الرئيسي من هذه الدراسة هو البحث عن عالقة األمطار والسيول إلعادة هيكلة التدفقات عند
مخرج حوض الشلف االعلى .وتهتم هذه الدراسة الى تحديد األوضاع وبخاصة الجغرافية ،الجيولوجية
والمائية المناخية التي تؤثر على التدفق في منطقة الدراسة.
وال يمكن أن يتحقق ذلك إال بالنمذجة الهيدرولوجية ،ولمالحظة سلوك الحوض والسيطرة على
األحداث المائية المناخية في تحديد المعايير المثلى للمحاكاة الهيدرولوجية.
المحاكاة الهيدرولوجية لألحواض ،تتمثل في تحديد االمطار الى سيول من خالل استخدام النماذج
الهيدرولوجية .لهذا قمنا في دراستنا باختبار نموذج SCS-LRالمدمج في منصة برنامج. ATHYS
ووفقا لنتائج المحاكاة التي تم حصول عليها بواسطة نموذج SCS -LRنستخلص أن نموذج
األمطار والسيول لحوض الشلف االعلى يقدم نتائج مرضية للغاية ،مع الحفاظ على المعالم بمعنى المادي
القوي ،ويفسح مجاالت التطبيقات الموجهة نحو التشغيل في مجال التحديد المسبق وتدفق المحاكاة.
الكلمات المفتاحية:
حوض ،الشلف االعلى ،النمذجة الهيدرولوجية ،نموذج األمطار والسيول.ATHYS ،SCS-LR ،
Sommaire
Sommaire
Remerciements .............................................................................................................................
Résumé .........................................................................................................................................
Abstract ........................................................................................................................................
Introduction General…………………………………………………………………………-1-
I Aspect théorique sur modélisation hydrologique .......................................................... - 3 -
I.2 La modélisation hydrologique : .............................................................................. - 3 -
I.2.1 Les modèles pluies-débits ................................................................................. - 4 -
I.2.2 Le choix de modèle ........................................................................................... - 5 -
I.3 Littérature sur la modélisation pluie-débit. ............................................................ - 6 -
I.3.1 Les modèles hydrologiques ............................................................................... - 6 -
I.3.2 Classification des modèles hydrologiques ........................................................ - 7 -
I.3.3 Les modèles conceptuels spatialisés ............................................................... - 11 -
I.3.4 Fonction de production SCS ........................................................................... - 13 -
I.3.5 Fonction de transfert ........................................................................................ - 15 -
II. Présentation générale de la zone d’étude ................................................................... - 18 -
II.1 Domaine géographique de la zone d’étude : ........................................................ - 18 -
II.2 Contexte géologique et hydrogéologique ............................................................. - 19 -
II.2.1 Contexte géologique ........................................................................................ - 19 -
II.2.2 Contexte hydrogéologique .............................................................................. - 26 -
II.3 Contexte hydro-climatologique ............................................................................ - 28 -
II.3.1 Contexte hydrologique .................................................................................... - 28 -
II.4 Contexte climatologie .......................................................................................... - 34 -
II.4.1 Etude des précipitations .................................................................................. - 34 -
II.4.2 Températures ................................................................................................... - 37 -
II.4.3 Evaporation ..................................................................................................... - 38 -
II.4.4 Vents................................................................................................................ - 39 -
II.5 Les données hydrologiques .................................................................................. - 41 -
II.5.1 Réseau de mesures hydrologiques ................................................................... - 41 -
III. Modélisation des débits ............................................................................................... - 47 -
III.1 Mise en place des données pour les simulations .................................................. - 48 -
III.1.1 Données géographiques ................................................................................... - 48 -
III.1.2 Données hydro climatiques ............................................................................. - 49 -
III.2 Propriété et choix du modèle ................................................................................ - 51 -
III.2.1 Présentation du mode de calcul modèle .......................................................... - 52 -
III.2.2 Calage du modèle ........................................................................................... - 54 -
III.2.3 - Sensibilité du modèle aux paramètres ........................................................... - 55 -
Sommaire
F : Le paramètre d’ajustement.
Vm : La vitesse de transfert sur chaque maille,
m : Le mètre.
s : Le seconde.
Km : Le temps d’amortissement.
Tm : Le temps de propagation.
Liste des Symboles et Acronymes
pk : La pente de la maille k,
Ds : la dénivelée spécifique.
A : la superficie du bassin versant.
Ia/S : la proportion des pertes initiales par rapport à la capacité maximale du réservoir sol.
ds : Le coefficient de la vidange exponentielle de la lame infiltrée, simulant la vidange du
réservoir sol (évaporation, percolation),
α : l’exposant de la pente de la maille, adimensionnel.
INTRODUCTION GENERALE
En effet, les modèles pluie-débit sont des outils très utilisés, parce que les données
débimétriques sont moins accessibles que les données pluviométriques. L’intérêt des modèles
n'est pas seulement de reproduire des débits qui n'ont pas été observés sur un bassin donné
(reconstitution sur des longues périodes, simulation d’événements extrêmes), mais aussi de
simuler les effets de modifications du bassin ou du climat sur les écoulements.
De ce fait, nous proposons ici une approche qui permet l’étude et la quantification des
ressources en eau dans la région du haut Cheliff à travers une modélisation pluies débits.
Cependant, De nombreux modèles ont été appliqués à présent, mais restent souvent
inefficaces pour la simulation surtout pour des crues rapides. Dans la plupart des cas, ces
modèles comportent de trop nombreux paramètres, et sont susceptibles de dériver hors du
domaine d'observation par suite d'une mauvaise estimation de ces paramètres (équivalence ou
équifinalité des jeux de paramètres). Les modèles conceptuels distribués constituent un bon
moyen de réduire ce genre de problème. Ainsi, le modèle distribué SCS-LR, disponible dans
la plate-forme de modélisation MERCEDES du logiciel ATHYS a été utilisé dans le cadre de
cette étude.
Pour atteindre cet objectif notre étude est structurée en quatre chapitres :
-1-
Introduction générale
-2-
Chapitre I :
ASPECT Théorique SUR
LA Modélisation
HYDROLOGIQUE.
Chapitre I : Aspect théorique sur la modélisation hydrologique
Introduction
Le fonctionnement hydrologique des bassins versants reste encore assez mal connu et
incontournable dès lors que l’on s’intéresse à des problématiques relatives à la gestion des
ressources en eau et la maîtrise du risque hydrologique. La modélisation pluie-débit constitue
un substitut efficace à des mesures de débits fastidieuses est devenue un outil indispensable
pour divers domaines ; gestion des ressources en eau, prédétermination des crues et
aménagement de l’environnement ou de protection de l’environnement et de prévision de
catastrophes naturelles (simulations de crues à court terme, prévision d’étiages).
(Amiar,2014).
On distingue quatre types de variables qui existent en totalité ou en partie dans chaque
modèle :
- Variables d’entrées : Le modèle fait appel à ces variables qui dépendent du temps
et/ou de l’espace (pluie, ETP, caractéristiques et hydrodynamiques du milieu, …) ;
- Variables de sorties : Le modèle répond par un ensemble de variables (débits, flux
ou concentration en polluants, …) ;
-3-
Chapitre I : Aspect théorique sur la modélisation hydrologique
-4-
Chapitre I : Aspect théorique sur la modélisation hydrologique
✓ Simulation des débits pour des reconstitutions historiques ; (Prévision des crues) ;
✓ Détection de l’influence des changements d’occupation du sol ;
✓ Simulation de l’impact du changement climatique ;
✓ Pré détermination des débits de crue ou d’étiage (fréquence et durée) ; - (Prévision des
étiages).
Notre étude est portée sur le modèle ATHYS (Atelier hydrologique spatialisé), pour la
reconstitution des débits du bassin versant du haut Cheliff, et la quantification des ressources
en eau de ce bassin versant.
-5-
Chapitre I : Aspect théorique sur la modélisation hydrologique
Un modèle est une représentation d’un phénomène physique, dans le but de comprendre les
processus qui le régissent. Cette représentation peut être physique, analogique ou mathématique.
Dans le premier cas, le modèle est une maquette qui reproduit d’une manière plus au moins
adéquate la réalité. Les modèles analogiques se basent sur les similitudes entre le phénomène à
étudier et un autre phénomène physique. La modélisation mathématique est un outil essentiel
pour la connaissance des phénomènes naturels, elle essaye d’établir un lien entre les variables
d’entrée et de sortie par des relations mathématiques.
Au cours de ces dernières années, les efforts de la recherche pour la compréhension du cycle de
l’eau dans les milieux naturels (bassins versants, rivières, nappes, etc.), associés aux
développements de l’informatique, se sont concrétisés par l’apparition d’une multitude de
modèles mathématiques. On peut distinguer trois types d’utilisation des modèles mathématiques
en hydrologie :
- La modélisation comme outil de prévision : elle nous donne une anticipation des
évolutions futures du débit d’un cours d’eau. Il s’agit de l’utilisation opérationnelle la plus
courante des modèles hydrologiques.
Dans la plupart des cas cependant, les modèles développés sont basés sur des régressions
linéaires entre les variables indépendantes (pluie, débits amont) et les variables dépendantes
(débits aval), et font peu appel aux connaissances sur les processus hydrologiques.
-6-
Chapitre I : Aspect théorique sur la modélisation hydrologique
inondation. …). Cette modélisation doit décrire de manière fidèle et réaliste les différentes
étapes liées à la transformation de la pluie en débit. On l’utilise aussi pour obtenir des
informations intéressantes pour le dimensionnement des ouvrages hydrauliques.
Modèles Modèles
déterministes stochastiques
-7-
Chapitre I : Aspect théorique sur la modélisation hydrologique
-8-
Chapitre I : Aspect théorique sur la modélisation hydrologique
Ce type de modèle reproduit donc au mieux le comportement d’un système, plutôt qu’il n’avance
d’explications causales sur son comportement. Le modèle CEQUEAU est un parfait exemple de
modèle conceptuel que nous présenterons dans la section suivante.
Les modèles globaux offrent à l’utilisateur un choix très attractif, car il présente une structure
très simplifiée, il ne demande pas trop de données, faciles à utiliser et à calibrer. La
représentation du processus hydrologique est très simplifiée. Il peut souvent mener à des
résultats satisfaisants, et spécialement si l’objectif majeur est la prévision d’une crue.
-9-
Chapitre I : Aspect théorique sur la modélisation hydrologique
Les modèles conceptuels spatialisés ou semi-spatialisés représentent un grand progrès sur les
modèles globaux quand il s’agit d’analyser le fonctionnement interne d’un bassin. Le bassin
versant est discrétisé en unités spatiales (mailles ou sous-bassins) considérées comme
homogènes, qui se vident les unes dans les autres de l’amont en aval. Ainsi, on a la possibilité de
tenir compte de la répartition spatiale des facteurs et de suivre la genèse et la propagation des
débits à l’intérieur du bassin. C’est le cas du modèle CEQUEAU et du modèle HEC-HMS.
Les modèles à base physique spatialisés sont des modèles qui décrivent les mécanismes internes
d’un système (bassin versant) ayant comme base les lois de la mécanique, de la physique, de la
thermodynamique, etc. Du point de vue théorique, ces modèles sont indépendants de tout calage
parce que leurs paramètres sont mesurables. Ils permettent une description théorique unifiée de
la plupart des flux observés dans un bassin versant et servent à modéliser les principaux
processus hydrologiques comme :
▪ L’écoulement de surface à partir des équations de Saint-Venant ;
▪ L’écoulement en milieu saturé à partir des équations de Darcy ;
▪ L’écoulement en milieu non saturé à partir des équations de Richards ;
▪ L’évapotranspiration à partir des équations de conservation ou de relations entre flux.
A partir d’un découpage de l’espace en mailles de taille adaptée au problème à traiter, ils
simulent les diverses composantes du cycle de l’eau sur chaque maille (ruissellement,
infiltration, évapotranspiration) et convolent les transferts de maille à maille, jusqu’au réseau
hydrographique constitué. Ces modèles sont complexes à mettre en œuvre et exigent
d’importantes quantités de données. Ils sont bien adaptés à la simulation de la diversité d’un
bassin.
Cependant, ces modèles ne peuvent traiter les grands bassins en raison de leur grande
hétérogénéité morphologique et météorologique. Parmi ces modèles on peut citer le modèle
Mike SHE.
Pour dépasser les limites de chacune des approches précédentes (modèles conceptuels trop peu
réalistes, modèles à base physique trop complexes), il est intéressant d’essayer une modélisation
hydrologique qui peut être :
- 10 -
Chapitre I : Aspect théorique sur la modélisation hydrologique
La méthode de discrétisation spatiale varie d’un modèle à l’autre : mailles carrées, sous- bassins
versants, éléments de versant, plans versant et canal, unités hydrologiques ou aires contributives
(modèle TOPMODEL).
- 11 -
Chapitre I : Aspect théorique sur la modélisation hydrologique
I.2.3.1 ATHYS
L’Atelier Hydrologique Spatialisé, ATHYS, regroupe un ensemble de modèles hydrologiques
spatialisés associés à des traitements de données hydro-climatiques et géographiques.
Le débit généré pour un événement pluvieux (transformation pluie-débit) est calculé en trois
étapes :
- Pour chaque maille, le modèle de production permet d’estimer la quantité de pluie qui
va contribuer au ruissellement.
- Le modèle de transfert calcule l’hydrogramme produit par chaque maille à l’exutoire
du bassin versant. Ce calcul est réalisé à partir du résultat obtenu en appliquant la
fonction de production.
- 12 -
Chapitre I : Aspect théorique sur la modélisation hydrologique
- Les apports de chaque maille sont sommés pour obtenir le débit à l’exutoire.
Le modèle de production SCS relie le cumul de la pluie efficace Peau cumul de la pluie
brute Pb, par l’équation :
- 13 -
Chapitre I : Aspect théorique sur la modélisation hydrologique
Ce schéma a également été complété, dans MERCEDES, par un réservoir sol, dont le
niveau décrit l’évolution du stock hydrique en continu. Le réservoir est alimenté par l’eau
qui s’infiltre, et il est vidangé par la reprise évaporatoire, la percolation vers la nappe
profonde et l’eau exfiltrée, remise en écoulement gravitaire. La vidange est calculée
globalement à l’aide du modèle du réservoir linéaire.
Avec
- 14 -
Chapitre I : Aspect théorique sur la modélisation hydrologique
Pour la cohérence du schéma, la vidange appliquée au réservoir sol doit également être
appliquée au cumul des pluies, pour que le coefficient de ruissellement C(t) soit égal à 0
lorsque le réservoir sol est vide. Le cumul des précipitations P(t) doit donc être diminué en
relation avec la vidange du réservoir sol. Ceci est réalisé par l’introduction d’une vidange
Vid2 appliquée au cumul de pluie :
Avec
Enfin, on considère qu’une quantité exfiltrée, Exf(t), qui correspond à une fraction du
volume vidangé, est remise en écoulement gravitaire :
Le ruissellement total Rtot(t) [L.T-1] produit par une maille au temps t est donc égal à :
- 15 -
Chapitre I : Aspect théorique sur la modélisation hydrologique
𝐼𝑘 60
Tm = ∑ ⋅ 𝑑𝑡
𝑘 𝑉𝑘
Km, le temps d’amortissement, figurant les variations des vitesses dans le temps.
L’amortissement est simulé par un réservoir linéaire, dont la capacité Km.
𝑘𝑚 = 𝑘0 . Tm
Où lk est les longueurs des k-cellules entre la cellule m et l’exutoire (figure.2), V0la vitesse de
parcours, qui est censé être uniforme sur l’ensemble du bassin versant (paramètre de routage).
- 16 -
Chapitre I : Aspect théorique sur la modélisation hydrologique
Conclusion :
Dans ce chapitre nous avons présenté les différents modèles utiles à l’hydrologue pour
le traitement des problèmes liés à la simulation de la transformation de la pluie en débit.
La modélisation pluie-débit est considérée dans ce cas comme un des moyens les
mieux adaptés pour décrire une réalité complexe du fonctionnement d’un bassin versant qui
dépend des caractéristiques hydrologiques et pédologiques propre à ce bassin.
- 17 -
Chapitre II :
Présentation générale de
la zone d’étude
Chapitre II : Présentation générale de la zone d’étude
- 18 -
Chapitre II : Présentation générale de la zone d’étude
Le bassin versant du Haut Cheliff est situé dans le domaine du tell septentrional et
correspond à un sillon intra montagneux subsides, allongés grossièrement d’Est en Ouest, sur
plus de 60 Km de long et une moyenne de 25 Km de large, et localisé entre le massif de
Boumaâd et les contrefort de l’Ouarsenis d’après A. Perodon 1957, M.Mattaouer 1958. Du
point de vue litho stratigraphique, la dépression du Haut Cheliff est constituée dans son
ensemble par des terrains d’âges Mio-Plio-quaternaire.
II.2.1.1 Stratigraphie
Dans cette étude on s’est basé essentiellement sur les caractéristiques litho-
stratigraphiques, d’après l’analyse du log litho-stratigraphique nous pouvons ressortir les
différentes formations géologiques suivants (Figure II-2).
a) Le Primaire
Le Primaire est constitué d’une alternance des schistes noirs graphiteux et de bancs de
quartzites et argiles.
a) Le Trias
Le Trias est généralement constitué par des dolomies et des calcaires dolomitiques. Il affleure
dans le massif de Zaccar.
b) Le Jurassique
Le Jurassique est constitué, dans le massif de Zaccar par des calcaires massifs, compactes,
fracturés et karstifiés et surmontés par des schistes gréseux et des marno-calcaires.
L’ensemble de la série atteint une épaisseur de 700 m environ.
c) Le Crétacé
Le Crétacé affleure sur les bordures latérales de la plaine, il est représenté de la base au
sommet par :
- 19 -
Chapitre II : Présentation générale de la zone d’étude
d) Le Miocène
Le Miocène peu atteindre 1300m d’épaisseur et il est constitué de la base au sommet par :
Le Miocène inférieur repose sur des terrains plus anciens Éocènes et Crétacés et
débute par une formation conglomératique d’épaisseur environ 1000m.
Le Miocène supérieur est constitué par une série représentée de la base au sommet par des
niveaux conglomératiques épais de 1 à 2m à ciment sablo marneux, surmontés par un niveau
marneux de couleur bleu d’une épaisseur environ 50m.
Des marnes blanches alternant avec quelques bancs siliceux couronnés par un dépôt de
grès et de sable d’une épaisseur d’environ 100m.
e) Le Mio-Pliocène
Le Mio-Pliocène est constitué par des galets quartzeux, des conglomérats, des grés et
des argiles détritiques ainsi que le travertin ancien déposé au niveau des sources du Zaccar.
f) Le Pliocène
Il est représenté par une série marneuse qui débute par des niveaux détritiques gréseux
ou sableux et parfois conglomératiques, ensuite on a un passage progressif dans la partie
supérieure à des marnes sableuses et sables argileux formant tous les termes de passage aux
grés Astiens.
- 20 -
Chapitre II : Présentation générale de la zone d’étude
- L’Astien
Il est représenté par une série marine finement sableuse ou gréseuse de couleur jaune
fauve. Les grés sont souvent à ciment calcaire, localement passer à des calcaires gréseux.
Les dépôts sableux rouges à stratifications entre croisées qui succèdent aux
conglomérats précédents appartiennent au villafranchien.
Ces dépôts reposent le plus souvent sur les limons gris. Parfois les conglomérats du
Pliocène supérieur n’existent pas et les formations rouges du villafranchien se trouvent alors
en contact des grés astiens ; d’ailleurs, au Sud de Chélif sur la bordure de l’Ouarsenis, le
Pliocène disparait complètement. A ce niveau, le Quaternaire recouvre directement les
calcaires gréseux du Miocène.
g) Le Quaternaire
- Le Quaternaire ancien
Il est représenté essentiellement par des alluvions conglomératiques : Galets, graviers et des
sables. Les alluvions du Quaternaire ancien sont observées au pied du massif de Zaccar et en
particulier à l’Est de Sidi-Lakhder. Elles forment les vastes collines en bordure Sud de la
plaine d’El Khemis.
- Le Quaternaire récent
✓ Des limons bruns à rouges, les niveaux les plus anciens sont discordants sur le
Quaternaire ancien. Ces limons recouvrent la presque totalité de la dépression sous
forme d’alternance sablo gréseuse.
✓ Des Grés, visibles surtout sur les berges des oueds, forment des accumulations
de niveaux très fins dont les épaisseurs peuvent atteindre cinq mètres. Ce sont
d’anciennes vases asséchées.
- 21 -
Chapitre II : Présentation générale de la zone d’étude
- 22 -
Chapitre II : Présentation générale de la zone d’étude
Figure II-2 la litho stratigraphique de la zone d’étude (M. Mattauer, 1953, 1955,1958). .
- 23 -
Chapitre II : Présentation générale de la zone d’étude
a) La série autochtone
La série autochtone est constituée par des marnes et des calcaires formant de vaste
bombement recoupé par plusieurs failles des terrains essentiellement schisto-quartzitiques du
Crétacé et de calcaires du Jurassique néocomien.
b) Le complexe A
- 24 -
Chapitre II : Présentation générale de la zone d’étude
- 25 -
Chapitre II : Présentation générale de la zone d’étude
- 26 -
Chapitre II : Présentation générale de la zone d’étude
Toute fois notre étude dont l’objectif consiste à évaluer les réserves en eau souterraine est
axée essentiellement sur l’aquifère alluviale du haut Chéliff qui constitue l’aquifère le plus
important et le plus exploité dans toute la région de la plaine du haut Chéliff.
- 27 -
Chapitre II : Présentation générale de la zone d’étude
La quantification des apports en eau et leur distribution sont intimement liées aux
différents paramètres du climat et à ceux qui régissent les partages multiples et successifs des
eaux.
Le bassin versant du Haut Chéliff est composé de six (06) sous bassins versants
figure (II-1), dont Trois situés sur la rive droite et Trois sur la rive gauche.
Figure II-5 : Carte des sous bassin Versant de la région des Haut Cheliff.
La forme qui peut être traduite par l’indice de compacité de « Graveluis » Kc, étant
définie comme le rapport du périmètre du bassin versant par celui d’un cercle qui aurait la
même superficie.
P P
Soit : KC 0,28
2 A A
- 28 -
Chapitre II : Présentation générale de la zone d’étude
Plus (Kc) approche de l’unité plus le bassin serait ramassé et le temps de concentration
serait plus court. Un bassin très allongé ne réagira pas, toutes choses égales d’ailleurs, de la
même manière qu’un bassin de forme ramassée (M.Roche, 1963).
Tableau II-1 : Valeurs de l’indice de compacité sous bassins versants du haut Chéliff.
(MEHAIGUENE.2006)
Les valeurs de l’indice de compacité (Tableau II-1) montrent que les sous bassins versant ont
généralement une forme allongés avec des valeurs variant entre 1.33 à 1.78.
De ce fait les bassins versant de la région d’étude présentent, dans l’ensemble, des
formes assez allongées favorisant un écoulement lent des eaux de ruissellement. Mais
l’influence du relief est plus déterminante.
La répartition par tranches d’altitude de la superficie des bassins, plus significative que
l’altitude moyenne, ressort nettement des courbes hypsométriques ; elle est, par ailleurs,
visualisée par le diagramme des fréquences altimétriques établi pour les bassins de la région
d’étude.
Il s’avère que les secteurs, d’altitude supérieure à 1200 m, sont peu étendus dans les
bassins Aval Boughzoul; Elles correspondent à des massifs montagneux importants, seuls les
sommets des monts de l’Ouarsenis dépassant les 1600 m d’altitude.
Enfin, les secteurs, d’altitude inférieure à 400 m dans le bassin Aval Boughzoul,
représentant un peu plus du quart (29,45%) de la superficie totale du bassin. Elles concernent
essentiellement la moitie septentrionale du bassin (figure II-6).
- 29 -
Chapitre II : Présentation générale de la zone d’étude
Superficie (%)
0% 50% 100%
2000
1800
0, 01%
0, 04%
1600
0, 10%
0, 23%
1400
0,42%
Alt 1200 1, 04%
itu 3, 81%
de 7, 03% Altitude médiane
1000
(m 9, 01% courbe hypsométrique
) 9,61%
800
10,19%
10,40%
600
9,32%
9,33%
400 Altitude moyenne
9,63%
9,36%
200
6,86%
3,59%
0
0 2000 4000 6000 8000 10000 12000 14000 16000 18000 20000 22000 24000
Superficie(km 2)
Ce rectangle est défini pour avoir comparé facilement les bassins entre eux, et puisque
les valeurs de Kc sur tous les bassins de la région d’étude (Haut Cheliff) sont supérieures à
1,12 dans ce cas, ces bassins versants sont assimilables à un rectangle, ainsi les bassins sont
trois fois plus longs que larges.
Les valeurs de la longueur et de la largeur du bassin sont définies d’après les formules
suivantes :
- 30 -
Chapitre II : Présentation générale de la zone d’étude
Tableau II-2 : Longueur équivalente et largeur équivalente par sous bassin versant
(MEHAIGUENE.2006)
Le relief joue un rôle prépondérant dans l’aptitude au ruissellement des terrains. Son
appréhension peut être faite à l’aide de l’indice de pente global (Ig).
D
Ig
L
D : dénivelée simple ou bien la distance verticale, exprimée en mètres, sépare les altitudes
ayant 5% de la surface du bassin au dessous d’elles ; elle peut être assimilée à la distance
séparant les deux courbes de niveau maitresses extrêmes qui délimitent le bassin.
Pour permettre une comparaison plus poussée des bassins entre eux du point de vue
relief, on détermine la dénivelée spécifique (Ds) puisse déduit de l’indice de pente globale.
Ds Ig A
- 31 -
Chapitre II : Présentation générale de la zone d’étude
Tableau II-3 : Classes de relief au sens des hydrologues de l’O.R.S.T.O.M. (IRD actuellement)
Ds Classes de relief
Ds< 10 m Classe R1 – relief très faible
10 <Ds< 25 m Classe R2 – relief faible
25 <Ds< 50 m Classe R3 – relief assez faible
50 <Ds< 100 m Classe R4 – relief modéré
100 <Ds< 250 m Classe R5 – relief assez fort
250 <Ds< 500 m Classe R6 – relief fort
Ds> 500 m Classe R7 – relief très fort
(MEHAIGUENE.2006)
Tableau II-4 : Indice de pentes globales et classe de relief par sous bassin versant
Sous Ds Classe
Bassin Superficie Ig
bassin (m/km) de relief
0114 1383,69 17,67 657,11 R7
0115 767,29 15,01 415,86 R6
Haut 0116 857,64 19,41 568,51 R7
Chéliff 0117 743,97 23,87 650,97 R7
0118 680,55 23,73 619,02 R7
0119 898,54 21,38 640,85 R7
(MEHAIGUENE.2006)
Toutefois, d’après les valeurs de l’indice de pente global (Ig) et la dénivelée spécifique
(Ds) (Tableau II-3et II-4) les résultats permettent de distinguer deux catégories de classe de
relief : relief fort et relief très fort.
Les faibles valeurs de la dénivelée spécifiques (Ds), ou les reliefs forts, correspondent
au sous bassins Oued Chéliff Harbil (0115).
Les valeurs les plus fortes de l’indice de pente global (Ig) et la dénivelée (Ds), ou les
reliefs très forts, correspondent aux sous bassins : Cheliff Harraza ,AribEbda , Cheliff Ghrib,
AribChéliff et Deurder.
Dans ce même sens, et du fait que ce genre de classement comparatif s’applique mieux
à des bassins de petite taille, l’importance relative de la surface des bassins étudiés « égalisé »
en quelque sorte le relief, d’où les faibles différences entre les classes du relief trouvées.
Cette brève analyse a permis de relever les nuances morpho métriques existant entre
les unités de la région d’étude. Cependant le bassin du Haut Cheliff, montre des caractères
morpho métriques assez particuliers, favorisant, à priori, le ruissellement de surface.
- 32 -
Chapitre II : Présentation générale de la zone d’étude
Le bassin versant du Haut Cheliff est traversé par un cours d’eau principale qui
est l’oued Cheliff, il constitue un collecteur d’eaux de ruissellement des affluents de sa
rive droite et sa rive gauche.
Sur la rive droite, les principaux affluents sont l’oued Harbil et oued Arib. Sur
la rive gauche, les principaux affluents sont l’oued Deurdeur et oued Harreza
(figure II-6).
- 33 -
Chapitre II : Présentation générale de la zone d’étude
Le bassin versant du Haut Chéliff est équipé par dix-neuf stations pluviométriques
(figure II-5). Ces stations font partie du réseau hydro climatologique national géré par
l’Agence Nationale des Ressources Hydriques (ANRH).
- 34 -
Chapitre II : Présentation générale de la zone d’étude
Les caractéristiques des dix-neuf stations sont mentionnées dans le tableau II-5.
Coordonnée
N° Code Nom de station
x Y Z
1 011405 GHRIB BGE 487250 318400 460
2 011501 TAMEZGUIDA 497000 323900 490
3 011510 DOMAINE MOSBAH 469800 321550 313
4 011512 AIN SULTAN 465500 325950 285
5 011514 DJENANE BEN OUADAH 475650 325200 320
6 011601 MARABOUT BLANC 460175 315600 368
7 011603 BORDJ EL AMIR AEK 461100 285300 1080
8 011605 THENEIT EL HAD 439800 285650 1150
9 011606 SIDI MOKREFI 463200 304250 425
10 011609 DEURDEUR BGE 458800 300050 575
11 011702 ARIB CHELIF 439500 332550 230
12 011703 B. EL AMIR KHALED 455850 313800 370
13 011711 SIDI LAKHDER 452000 329650 250
14 011715 EL ABABSA 443900 318500 313
15 011717 EL KHEMIS ANRH 456200 329000 300
16 011718 HARREZA BARRAGE 445350 321400 312
17 011801 ARIB EBDA 439650 335550 280
18 011803 SIDI MEDJAHED 452150 337050 850
L’analyse des précipitations s’est basée sur les données disponibles relatives aux observations
des stations pluviométriques El Khemis ANRH (011717), Harreza Barrage (011718), Arib
Cheliff (011702)et El Ababsa(011715) pour la période (1985-2011) (Annexes II-1 et II-2).
La figure ci-dessous expose d’une part la variation des cumuls annuels pour chaque
station de mesure et d’autre part la variation des moyennes de ces cumuls annuels.
On constate tout d’abord de fortes disparités en fonction des années, par exemple il ya
deux fois plus en 1995 qu’en 1993 (à Arib Cheliff602 mm en 1995 et 194 mm en 1993). La
variabilité temporelle de la pluviométrie est un facteur important à prendre en compte dans le
cadre de cette étude : les années sont indépendamment sèches ou humides.
- 35 -
Chapitre II : Présentation générale de la zone d’étude
500
400
300
200
100
0
1985
2008
1986
1987
1988
1989
1990
1991
1992
1993
1994
1995
1996
1997
1998
1999
2000
2001
2002
2003
2004
2005
2006
2007
2009
2010
2011
Figure II-9 : Précipitations moyennes annuelles (période 1985-2011)
L’analyse des pluies annuelles pour la période 1985 à 2011 (figure II-9), montre
que les années 1986, 1995,2008 et 2010 sont les plus pluvieuses sur l’ensemble des
stations pluviométriques, avec un maximum de 613 mm enregistré en 2008 à la station de
KHEMISANRH. Par contre les années les moins pluvieuses sont celles de 1993 et
1999avec un minimum 176 mm enregistré au niveau de la station de Harreza.
Code Nom de station Sépt. Oct. Nov. Déc. Jan. Fév. Mar. Avr. Mai Juin Juil. Aout
011715 EL ABABSA 20.6 28.2 40.6 46.0 51.4 49.5 39.6 36.1 21.8 6.6 2.3 5.9
011702 ARIB CHELIFF 19.6 28.0 42.5 44.3 45.6 46.3 40.0 40.6 23.1 6.6 1.5 6.1
011717 EL KHEMIS ANRH 20.3 28.4 41.9 45.2 46.1 47.7 39.8 40.8 22.6 5.7 1.6 6.3
011718 HARREZA 24.3 28.0 43.9 46.9 48.6 51.2 42.0 33.1 23.4 9.8 4.1 7.1
- 36 -
Chapitre II : Présentation générale de la zone d’étude
On constate aussi, que les mois de Mai, Novembre et Aout sont les mois les plus pluvieux
pour les quatre stations pluviométriques, avec un maximum de 51.4 mm en mois Janvier pour
la station HARRAZA et AL ABABSA.
60
11715
EL
50 ABABS
Précipitations (mm)
A
40
30
20
10
0
Sépt Oct Nov Déc Jan Fév Mar Avr Mai Juin Juill Aout
II.4.2 Températures
La température, second facteur constitutif de climat, influe sur le développement des
végétations. Ce sont les températures extrêmes plus que les moyennes qui ont influence sur la
végétation sauf si elles sont exceptionnelles et de courte durée (Greco, 1966).
Les températures moyennes atteignent une valeur maximale de 29,5 °C en Juillet et Aout et
9,4 °C en mois de Janvier.
- 37 -
Chapitre II : Présentation générale de la zone d’étude
30
25
20
15
10
5
0
2003
1989
1990
1991
1992
1993
1994
1995
1996
1997
1998
1999
2000
2001
2002
2004
2005
2006
2007
2008
2009
2010
2011
2012
2013
2014
II.4.3 Evaporation
L’évaporation est la quantité d’eau restituée à l’atmosphère sous forme de vapeur, son
étude est essentielle pour connaitre le potentielle hydrique d’une région ou d’un bassin
versant.
- 38 -
Chapitre II : Présentation générale de la zone d’étude
350,0
300,0
Evaporation (mm)
250,0
200,0
150,0
100,0
50,0
0,0
Sept Oct Nov Déc Jan Fév Mars Avr Mai Juin Juil Aout
- Une période de forte évaporation allant du mois Mai au Mois d’Octobre dont elle
atteint une valeur maximale de 293mm.
- Une période de faible évaporation allant du mois de Novembre au mois d’Avril ou
elle atteint une valeur minimale de 44,36 mm.
A l’échelle annuelle, on enregistre une faible évaporation durant la période allant de 2007 à
2014 d’où une valeur minimale de 57,07 mm enregistrée durant l’année 2009. L’évaporation
la plus forte à été enregistrée durant l’année 1989avec une valeur de 190,6 mm
II.4.4 Vents
Dans la vallée du Cheliff, les pressions atmosphériques sont plus élevées le matin (à 7
h) que sur le littoral, la situation s’inverse à 13 heures car l’abaissement absolu est très
important dans le Chéliff est assez faible sur le littoral (mai et juin). Ce phénomène explique
l’importance de certain vents locaux vent d’Est à Ain Defla du Nord-Est à Khemis -Miliana)
(Remaoun M-2007).
D’une manière générale, la région d’étude se caractérise par des vents légers, généralement au
moins 10 km.
- 39 -
Chapitre II : Présentation générale de la zone d’étude
4,5
4
3,5
3
Vents (Km/h)
2,5
2
1,5
1
0,5
0
Sept Oct Nov Déc Jan Fév Mars Avr Mai Juin Juil Aout
- 40 -
Chapitre II : Présentation générale de la zone d’étude
- 41 -
Chapitre II : Présentation générale de la zone d’étude
D’après les données hydrométriques disponibles, les débits maximums mensuels sont
observés pendant les mois de Janvier, Février et Mars pour les stations Arib Cheliff et Djenan
Ben Ouadah Situé sur le cours d’eau principale du Cheliff.
Par contre, pour les stations situées sur la rive gauche du Cheliff (Marabou Blanc et
BirOuled Tahar) les débits maximums sont enregistrés durant les mois de Septembre, Octobre
et Novembre (tableau II-9).
350
300
250
200
150
100
50
0
Sept. oct. nov. déc. janv. fév. Mars Avril Mai Juin juil. Aout
Les débits minimums mensuels ne dépassent guerre 1 m3/s pour l’ensemble du bassin du Haut
Cheliff (figure II-15). Le tableau ci-dessous présente la variation des débits minimums
moyens mensuels des cinq stations pour la période (1991-2013).
- 42 -
Chapitre II : Présentation générale de la zone d’étude
Tableau II-10 : Variation des débits minimums moyens mensuels pour la période (1991-2013).
Code Nom station Sept. oct. nov. déc. janv. fév. Mars Avril Mai Juin juil. Aout
011601 MARABOUT B. 0,15 0,01 0,01 0,03 0,04 0,05 0,06 0,06 0,04 0,32 0,17 0,07
011715 EL ABABSA 0,00 0,00 0,00 0,00 0,00 0,01 0,02 0,01 0,00 0,00 0,00 0,00
011514 DJ. B. OUADAH 0,16 0,13 0,29 0,34 0,37 0,47 0,48 0,35 0,56 1,08 1,30 0,64
011905 B. O. TAHAR 0,00 0,00 0,00 0,02 0,04 0,15 0,13 0,12 0,08 0,02 0,00 0,00
011702 ARIB CHELIF 0,01 0,02 0,09 0,40 0,61 0,59 0,35 0,15 0,02 0,02 0,02 0,01
1,4
MARABOUT BLANC
EL ABABSA
1,2
DJENANE BEN OUADAH
BIR OULED TAHAR
1,0
Débits Min (m3/s)
ARIB CHELIF
0,8
0,6
0,4
0,2
0,0
sept oct nov dec janv fev mars avril mai juin juil aout
Le tableau suivant présente la variation des débits moyens mensuels des cinq stations pour le
période (1991-2013).
Code Nom station Sept. oct. nov. déc. janv. fév. Mars Avril Mai Juin juil. Aout
011601 MARABOUT B. 0,96 0,69 0,38 0,20 0,76 0,60 0,75 0,70 0,74 0,99 0,78 0,93
011715 EL ABABSA 0,05 0,06 0,07 0,05 0,11 0,22 0,15 0,09 0,05 0,01 0,00 0,01
011514 DJ. B. OUADAH 1,44 2,82 1,22 1,18 3,02 3,45 2,67 2,21 2,66 2,39 2,66 1,91
011905 B. O. TAHAR 0,14 0,13 0,23 0,22 0,63 0,98 0,97 0,62 0,39 0,09 0 0,02
011702 ARIB CHELIF 1,87 3,59 2,53 4,50 6,91 9,74 7,96 4,70 2,19 1,59 1,93 1,78
- 43 -
Chapitre II : Présentation générale de la zone d’étude
12 MARABOUT BLANC
EL ABABSA
10 DJENANE BEN OUADAH
Débits moyens (m3/s)
BIR OULED TAHAR
8 ARIB CHELIF
0
sept oct nov dec janv fev mars avril mai juin juil aout
Les écoulements annuels dans le bassin du Haut Cheliff sont analysés à travers l’analyse des
débits enregistrés au niveau des stations Arib Cheliff, Djenan Ben Ouadah , Marabout Blanc,
El Ababsa et Bir Ouled Tahar durant la période (1991-2013).
Les plus forts débits écoulés sont enregistrés durant l’année (1994-1995),par contre les plus
faibles débits écoulés sont enregistrés durant l’année (1993-1994).
350 B. O. TAHAR
300 DJ. B. OUADAH
250
200
150
100
50
0
2007
1991
1992
1993
1994
1995
1996
1997
1998
1999
2000
2001
2002
2003
2004
2005
2006
2008
2009
2010
2011
2012
- 44 -
Chapitre II : Présentation générale de la zone d’étude
L’analyse de la variation des débits minimums moyens annuels (figure II-18) des cinq
stations étudiées pour le période (1991-2013) ; montre que les débits minimaux sont au-
dessous d’une valeur seuil de 1 m3/s sur l’ensemble des stations hydrométriques du bassin du
Haut Cheliff.
ARIB CHELIF
0,80
0,60
0,40
0,20
0,00
1994
2001
2008
1991
1992
1993
1995
1996
1997
1998
1999
2000
2002
2003
2004
2005
2006
2007
2009
2010
2011
2012
Figure II-20 : Variation des débits minimums moyens annuels (1991-2013)
Le tableau suivant présente la variation des débits moyens annuels des cinq stations étudiées
pour le période (1991-2013).
8,0
MARABOUT BLANC
7,0 EL ABABSA
DJENANE BEN OUADAH
Débits moyens (m3/s)
6,0
BIR OULED TAHAR
5,0 ARIB CHELIF
4,0
3,0
2,0
1,0
0,0
- 45 -
Chapitre II : Présentation générale de la zone d’étude
Durant la période 1991-2013 et sur l’ensemble des stations de la région d’étude les débits
écoulés les plus forts sont enregistrés durant l’année (1991-1992) surtout à la station Arib
Cheliff.Les plus faibles débits écoulés sont enregistrés durant l’année (1993-1994) marquant
dans la station Marabout Blanc.
Le tableau suivant présente la variation des apports moyens annuels des cinq stations étudiées
pour le période (1991-2013).
B. O. TAHAR
DJ. B. OUADAH
150
100
50
Conclusion :
- 46 -
Chapitre III :
Modélisation des débits
Chapitre III : Modélisation des débits
La fonction SCS est conçue pour simuler le ressuyage des sols entre deux averses, et
comporte trois paramètres :
- Vo la vitesse de transfert et
- Ko un coefficient de diffusion.
✓ Pour chaque maille, le modèle de production permet d’estimer la quantité de pluie qui
va contribuer au ruissellement.
✓ Le modèle de transfert calcule l’hydrogramme produit par chaque maille au niveau du
point de référence (exutoire du sous-bassin). Ce calcul est réalisé à partir du résultat
obtenu en appliquant la fonction de production.
✓ Les apports de chaque maille sont sommés pour obtenir le débit au point de référence.
- 47 -
Chapitre III : Modélisation des débits
Le MNT utilisé dans notre cas est un extrait du MNT du Nord Algérie, Provenant de la
base de données USGS/NASA SRTM Data, crées en 2007. C’est un MNT de maillage
régulier de 10m x10m.
Le maillage à été réduit ensuite à 100m x100m pour diminuer la taille des fichiers utilisés.
L’ensemble des fichiers ont le même format (Ascii raster) avec un maillage de 100m par
100m :
- 48 -
Chapitre III : Modélisation des débits
Tableau III-1 : Fichier entête des données hydro climatiques, la période (1985 à 2007)
- 49 -
Chapitre III : Modélisation des débits
Les données hydro climatiques sont enregistrées dans un fichier au format (.txt). Ce
fichier peut être converti sous le module VISHYR. Ce module propose différentes opérations
de correction, calcul, gestion et visualisation de données hydro-climatiques locales (figure III-
4).
Après avoir traité et analyser les données disponibles nous avons pu constituer notre base de
données correspondant à :
- 50 -
Chapitre III : Modélisation des débits
Le modèle continu simulera en complément les processus qui s’étalent sur une
temporalité plus longue : redistribution de l’eau dans les sols, évapotranspiration…
L’avantage des modèles continus est qu’ils décrivent les variables d’état du système à
tout moment, notamment au début des épisodes pluvieux intenses. A l’inverse, le modèle
événementiel doit être initialisé en fonction d’indicateurs d’humidité initiale du bassin, qui
peuvent être le débit de base, un indice des pluies antérieures, un indicateur saisonnier ou
empirique.
De ce fait notre choix est porté sur des modèles événementiels et spécialement le
modèle distribué dérivé du SCS. Ce modèle est réputé être souple, susceptible de s’adapter à
tout type de bassin sans préjuger de la nature des processus de formation de crues, et enfin, ce
- 51 -
Chapitre III : Modélisation des débits
modèle a déjà été utilisé sur des bassins comparables (Gaume, 2002 ; Marchandise, 2003 ;
Lequien,2003).
Le modèle choisi est un modèle distribué, reposant sur la discrétisation du bassin en mailles
carrées régulières. Le modèle est construit en 5 étapes :
Pour un cumul de pluie P inférieur à 0.2S, le coefficient C n’est pas défini. En réalité, il est
pris égal à 0. Cela correspond à une imbibition « forfaitaire » des premières pluies par le
bassin versant.
- 52 -
Chapitre III : Modélisation des débits
- Lorsque les pluies cessent, il n y a plus de génération de ruissellement avec cette loi.
- Le coefficient de ruissellement est toujours croissant dans le temps, et reste stable
pendant les interruptions de la pluie.
Cette dernière remarque nous a conduits à introduire une vidange de réservoir. Celle-ci
permet de gérer le ressuyage des sols et de générer des écoulements par exfiltration. On
définit le réservoir stoc(t) qui contient l’eau effectivement infiltrée dans le sol à l’instant t du
fait de la loi SCS. Ce réservoir est vidangé de manière linéaire sur le pas de temps Dt d’une
quantité vid(Dt) telle que :
Où ds est un paramètre global du modèle à caler. Cette vidange permet de simuler la décrue
par une reprise partielle de la vidange de taux w(0<w<1). w est le deuxième paramètre de
vidange à caler (Figure 33)
- 53 -
Chapitre III : Modélisation des débits
Ko est un paramètre global du schéma à caler. Pour chaque maille du MNT, la direction de
drainage est orientée vers l’une des 8 voisines, celle qui a l’altitude la plus basse. Les
directions de drainage permettent de définir les trajectoires d’écoulement de chaque maille
jusqu’à l’exutoire du bassin versant.
L’hydrogramme de crue est la somme sur la maille exutoire, dans le temps et dans
l’espace, des hydro grammes élémentaires des mailles appartenant au bassin versant.
- 54 -
Chapitre III : Modélisation des débits
Les paramètres du modèle sont calés à partir d’une base de données de 22 événements
enregistrés à Arib Cheliff (011702). L’ensemble de ces événements ont été choisis dans la
base de données initiale des 239 événements suivant les critères suivants :
Le paramètre Vo : agit sur le calage dans le temps de pic de crues et c’est le seul
paramètre dans ce modèle qui influe sur le calage dans le temps.
Conclusion
Dans ce chapitre on a essayé d’exposer le contexte général de logiciel hydrologique
ATHYS et les modèles à utilisées, la nécessité d’une base de données élaboré selon les
exigences du ce logiciel, la préparation et l’introduction des données géographiques et
hydro-pluviométriques nécessaires à notre modélisation.
- 55 -
Chapitre V :
Résultat et discussion
Chapitre IV : Résultats et discussions
- Le calage avec optimisation de type séparée d’où on optimise d’abord les paramètres de
production à l’aide d’une fonction critère calculée sur les volumes écoulés, puis les
paramètres de transfert à l’aide d’une fonction critère calculée sur le débit.
- Le calage sur les épisodes sur un groupe d’épisodes les paramètres optimisés sont
commun à un groupe d’évènement.
Pour le modèle SCS-LR, on admettra que les paramètres à calibrer sont S, la/S, ω,
ds,V0,α,β et K. Les deux paramètres S et V0, ont des effets relativement bien différenciés : S
joue sur le volume et ne modifie pas le calage dans le temps, Vo joue sur le temps de montée
et la position du maximum, mais ne modifie pas les volumes, l’observation la plus frappante
de ces tests est la faible sensibilité de ds, il est donc difficile d’en conclure un sens physique
(pour ce paramètre uniquement).
L’optimisation des paramètres est conditionnée par l’optimisation des critères d’erreurs
de simulation, à savoir le critère EAM (Ecart Arithmétique Moyen), EQM (Ecart Quadratique
Moyen) et le critère de Nash. Ce dernier est considéré comme critère pour tester l’efficacité
du modèle.
- 56 -
Chapitre IV : Résultats et discussions
Lorsque le critère de Nash est égal à 1, la modélisation est optimale : les valeurs
calculées sont égales aux valeurs mesurées. Un Nash inférieur à 0 traduit une mauvaise
simulation.
- 57 -
Chapitre IV : Résultats et discussions
Epi 09 : 29/01/1986 : 23h00 à 18/02/1986 :22h00, Epi 19 : 05/02/1987 : 06h00 à 09/02/1987 :22h00,
Nash=0.64 Nash=0.85
Epi 27 : 14/11/1987 : 21h00 à 16/11/1987 :23h00, Epi 81 : 26/01/1992 : 05h00 à 31/01/1992 :23h00,
Nash=0.88 Nash=0.92
Epi 121 : 26/02/2003 :23h00 à 03/03/2003 :23h00, Epi 149 : 01/03/1999 : 11h00 à 05/03/1999 :21h00,
Nash=0.8 Nash=0.89
Débit observé
Débit calculé
- 58 -
Chapitre IV : Résultats et discussions
D’autre part, en enregistre une vitesse moyenne d’écoulement de 3 m/s sur l’ensemble
des événements analysés. Cela se traduit par les faibles pentes caractérisant une grande partie
de la surface du bassin du Haut Cheliff (plaine du Haut Cheliff).
D’une manière générale, on remarque certaines différences entre les débits observés et
ceux calculés surtout dans la courbe de décrue de chaque hydro gramme (Figure IV-1),
malgré que les valeurs du critère de Nash, sont assez acceptables. Ces différences peuvent être
dues soit à des pertes par évaporation soit à la mauvaise connaissance de l’état initiales du sol
qui à été mal paramétrée lors de la simulation.
- 59 -
Chapitre IV : Résultats et discussions
Conclusion
La modélisation distribuée, via la méthode des mailles interactives, est la méthode qui a
donné un Nash plus important pour tous les évènements. Il a atteint 0,96 pour le cas de
l’épisode 11 et 0,91 pour l’évènement 230 ; on remarque que le transfert des débits d’une
maille à une maille adjacente, en incluant les pertes donne de bons résultats
D’après les résultats de calage, C’est difficile d’interpréter les résultats de calage surtout
les paramètres qui varient d’un évènement à l’autre, les valeurs de S obtenus sont presque
invariables entres les événements care on est dans la même région ou les conditions physique
sont fixé.
Sur l’ensemble des épisodes analysés (annexe IV-1)on remarque que les pics des crues
ainsi que les décrues sont bien simulés surtout pour les évènements de forts débits que ceux
de faibles débits ou les valeurs de critère NASH sont très importants, le NASH moyenne pour
les faibles crues(débit<100m3/s) égale à 0,77, par contre pour des fortes crues Q>100m3/s le
Nash moyenne égale à 0.90. Ainsi, et sur la base des résultats obtenus on constate que le
modèle SCS-LR s’ajuste bien sur les fortes crues que sur les faibles crues dans le bassin
versant Haut Chéliff.
- 60 -
Conclusion générale
CONCLUSION GENERALE
Nous pouvons conclure, à travers cette étude qui a pour objectif de modélisé les
écoulements superficiels de bassin versant haut Chéliff ont abouti à la conclusion suivante :
Cette étude a permis d’abord de faire une synthèse bibliographique sur la modélisation
hydrologique, la variabilité et la description des processus hydrologiques selon le type du
modèle utilisé. Ainsi notre choix, s’est porté sur le modèle utilisé SCS-LR du logiciel
ATHYS. En tenant compte de l’objectifs de notre étude hydrologique, disponibilité de
données et la robustesse et la simplicité du modèle ce qui permettre de négliger ou simplifier
certain paramètre dans la modélisation.
Le bassin versant du haut Chéliff constitue un des principaux bassins du grand bassin
hydrographiques du Cheliff-Zahrez, avec une superficie de 5331 Km2. En effet, le bassin du
Haut Cheliff, montre des caractères morpho métriques assez particuliers, favorisant, à priori,
le ruissellement de surface.
Une forme allongée avec une altitude moyenne de l’ordre de 640 m dans le bassin Haut
Cheliff, les secteurs, d’altitude supérieure à 1200 m, sont peu étendus. Elles correspondent à
des massifs montagneux importants, seuls les sommets des monts de l’Ouarsenis dépassant
les 1600 m d’altitude. Ainsi le relief est fort à très fort.
Du point de vue climatique le bassin du Haut Cheliff, est caractérisé par un climat semi -
aride avec un hiver froid et pluviaux et un été chaud et sec. La température moyenne annuelle
est de l’ordre de 19°C, les précipitations moyennes annuelles enregistré au niveau de la plaine
est de400mm.
Dans ce contexte, la modélisation des débits par le logiciel ATHYS qui s’appuie sur les
données géographiques et hydro climatiques, peut rendre en compte le comportement
hydrologique de bassin du Haut Cheliff.
Les résultats de simulation obtenus par l’application du modèle SCS-LR en mode non
continu sur le bassin d’étude sont assez suffisants :
Les paramètres appliqués dans la modélisation des différents épisodes, sont assez
cohérents avec l’interprétation physique proposée, ce qui permet d’envisager généralisation
sur d’autres épisodes. Le calage du modèle reste cependant nécessaire, mais peut se limiter
aux épisodes sélectionnés.
L’application et, surtout, l’interprétation du modèle supposent que la formation des crues
soit effectivement régie par des processus liés à la saturation des formations superficielles.
- 61 -
Conclusion générale
Dans d’autres cas (ruissellement hortonien), les indicateurs auxquels sont reliés les paramètres
du modèle et les conditions initiales du système peuvent être très différents.
Cependant, effet la valeur moyenne du critère de Nash de 0, 83, est assez acceptable.
Alors on peut conclue que le modèle utilisé SCS-LR de logiciel ATHYS est représentatif pour
le bassin du Haut Cheliff. Ainsi, on peut donc vérifier l’objectifs recherché relatif à la
reconstitution correcte des débits à l’exutoire du bassin en tenant compte les processus
hydrologiques ainsi les paramètres de calage qui pouvant avoir lieu.
- 62 -
Références bibliographiques
REFERENCES BIBLIOGRAPHIQUE
ALI TARIK.B, 2006, Modélisation pluie-débit mensuelle et journalière par les modèles
conceptuels et les systèmes neuro-flous (Application aux bassin algériens), Doctorat, INA,
Alger,30p-36p.
EDIJATNO, 1991. Mise au point d’un modèle élémentaire pluie-débit au pas de temps
journalier. Thèse de Doctorat, Université Louis Pasteur/ENGEES, Strasbourg, 242 p.
GAUMEE.,2002 Eléments d’analyse sur les crues éclair. Thèse de doctorat, ENGREF,
CEREVE.
GAUME E., LIVET M., DESBORDES M., VILLENEUVE J.P., 2004. Hydrological
analysis of the river Aude, France, flash flood on 12 and 13 November 1999. J. Hydrol., 286
(2004) 135–154.
- 63 -
Références bibliographiques
HSU, K.I…, GUPTA, H.V, SOROOSHIAN S., 1995.Artificiel neural network modeling of
the rainfaill-runoff process. Water Resour.Res.31(10), pp.2517-2530.
MAKHLOUF, Z.ET MICHEL, C., 1994. A two-parameter monthly water balance model
for French watersheds. Journal of Hydrology 162, 299-318.
MICHEL, C.1989. Hydrologie appliquée aux petits bassins versants ruraux, Cemagref,
Antony.
MINNS, A.W.AND HALL M.J., 1996. Artificial neural networks as rainfall-runoff models
hydrological Sciences journal 41(3): pp.399-418.
NASCIMENTO N.O. 1995. Appréciation à l’aide d’un modèle empirique des effets d’action
anthropiques sur la relation pluie-débit à l’échelle du bassin versant. Thèse de Doctorat,
CERGRENE/ENPC, Paris, 550 p.
NASCIMENTO N.O.De Oliveira. ’Application à L’aide D’un Modèle Empirique Des Effets
D’actions Anthropiques Sur La Relation Pluie-Débit à L’échelle D’un Bassin Versant ‘ Ecole
Nationale des Ponts et Chaussés, 1995.
STEENHUIS T.S.., WINCHELL M., ROSSING J., ZOLLWEG J.A., WALTER M.F.,
1995.SCS Runoff Equation Revisited for Variable-Source Runoff Area.j.Irrig.and
Drain.Engrg. 121,234(1995)
- 64 -
Références bibliographiques
THOMAS Z., MOLÉNAT J., CAUBEL V. AND MÉROT P., Modelling investigation of
the influence of a hedge surrounding bottomland on soil-water fluxes, Hydrological
Processes, (sous presse).
- 65 -
ANNEXES
Annexe N°II-1 : Les précipitations moyennes annuelles pour chaque station.
ARIB AL KHEMIS-
HARREZA
Année CHELIFF ABABSA ANRH
1985 456,9 436,3 470,60 450,10
1986 528,3 484,6 520,00 433,70
1987 297,1 285,8 311,50 275,00
1988 433 343,8 429,10 317,20
1989 306 278,9 286,60 278,00
1990 374,1 375 422,90 321,10
1991 487,4 280,84 477,60 363,30
1992 254,5 238,9 288,20 261,30
1993 194,4 225,2 192,10 253,11
1994 467,7 394,8 443,50 441,00
1995 602,5 459,3 590,50 298,28
1996 314,3 275,7 312,90 284,70
1997 406,5 333,5 433,40 315,90
1998 379,07 300,6 415,20 291,30
1999 229,1 185,4 244,10 176,60
2000 358,7 440,9 417,90 403,00
2001 335,37 237,9 317,00 251,60
2002 394,5 389,7 381,80 416,80
2003 307,8 315,8 378,20 318,60
2004 366,7 322,9 381,90 332,00
2005 445,2 356,3 421,40 374,50
2006 387,9 317,8 379,00 319,00
2007 416,2 362 443,10 395,90
2008 529,4 495,5 612,90 543,90
2009 464,6 441,7 496,00 508,70
2010 593,2 428,7 516,60 386,40
2011 503,5 434,4 445,02 389,58
Annexe N°II-2 : les débits maximums annuels des cinq stations étudiées pour le période (1991-
2013).
Numéro d’événement 09