REPUBLIQUE TUNISIENNE

MINISTERE DE L’AGRICULTURE ET DES RESSOURCES HYDRAULIQUES

INSTITUTION DE LA RECHERCHE ET DE
L’ENSEIGNEMENT SUPERIEUR AGRICOLES

INSTITUT NATIONAL
AGRONOMIQUE DE TUNISIE

Protection Contre les Inondations

SLIMANI Mohamed,
Professeur à l’INAT
Aouissi Jalel
Maitre-assistant à l’INAT

1
 Version Sept2020

Protection Contre les Inondations

Chap I Les Inondations dans le Monde ..................................................................................... 4
Introduction ............................................................................................................................ 4
HISTOIRE .............................................................................................................................. 4
Les inondations dans le Monde : ............................................................................................ 5
Les inondations en Tunisie ..................................................................................................... 6
INONDATION 2003.............................................................................................................. 9
Description de quelques inondations en Tunisie .................................................................. 11
CH II RAPPELS STATISTIQUES ......................................................................................... 15
Variables aléatoires : ............................................................................................................ 15
Les paramètres statistiques descriptives :............................................................................. 15
Les lois statistiques : ajustement : ........................................................................................ 16
Liaisons entre deux variables : ............................................................................................. 17
CH III Détermination du débit décennal sur un ..................................................................... 19
basin versant non jaugé Méthode Socose ................................................................................. 19
Introduction : ........................................................................................................................ 19
Les paramètres utilisés : ....................................................................................................... 19
Procédure pratique : ............................................................................................................. 22
Calcul des débits de période de retour > 10ans : .................................................................. 23
Application : ......................................................................................................................... 24
CH : IV Les formules empiriques ............................................................................................ 26
Estimation de l’apport moyen annuel : ................................................................................ 26
Estimation des apports fréquentiels : ................................................................................... 28
Estimation des débits fréquentiels : ...................................................................................... 29
Estimation du temps de concentration : ............................................................................... 30
CH V : Calcul des apports moyens et fréquentiels et des ....................................................... 32
Débits fréquentiels sur quelques BV du Sud tunisien. .............................................................. 32
(Médenine et Tataouine) .......................................................................................................... 32
Les formules à appliquer ...................................................................................................... 32
Les paramètres utilisés : ....................................................................................................... 32
Les bassins versants : ........................................................................................................... 32
CH : VI Les Ajustements statistiques et méthode de Gradex .................................................. 35
Introduction : ........................................................................................................................ 35
Ajustement statistique : (série longue n ≥ 20ans) : ............................................................. 36
Année ..................................................................................................................................... 36
Méthode du Gradex : ............................................................................................................ 38
Application: .......................................................................................................................... 39
CH : VII Dimensionnement des ouvrages................................................................................ 41
(lac collinaire) .......................................................................................................................... 41

2
 Etudes préliminaires : ........................................................................................................... 41
Le volume utile de la retenue Vu (m3) : ............................................................................... 42
Dimensionnement grossier de la digue : .............................................................................. 42
Etude hydrologique détaillée:............................................................................................... 43
Conception et dimensionnement de la digue : ..................................................................... 44
Volume de la digue (Vd) : .................................................................................................... 45
Ouvrages annexes : (évacuateur de crue) : ........................................................................... 45
Projet ..................................................................................................................................... 49
Chapitre VIII : Calcul des débits et apports fréquentiels sur les bassins versants du Sud ..... 49
Introduction .......................................................................................................................... 49
Les Apports moyens et fréquentiels : ................................................................................... 54
Les débits fréquentiels :........................................................................................................ 63
PROJET DE CONSTRUCTION D’UN OUVRAGE DE PROTECTION CONTRE LES
INNONDATIONS : CAS D’UN LAC COLLINAIRE ............................................................ 69
Les Etapes ............................................................................................................................ 69
Chapitre I : Présentation générale ........................................................................................ 70
chapitre II: Etude morphologique du bassin versant ............................................................ 75
Chapitre III : Etude hydrologique du bassin versant ............................................................ 80
Étude de la pluviométrie moyenne ....................................................................................... 80
temps de concentration : ....................................................................................................... 82
L’intensité des pluies ............................................................................................................ 83
Étude des apports liquides .................................................................................................... 84
Chapitre IV : Etude topographique ...................................................................................... 89
Chapitre V : Etude géologique et géotechnique ................................................................... 95
Chapitre VI : Etude de la digue .......................................................................................... 103
Chapitre VII : Etude des ouvrages annexes ....................................................................... 115

3
 Chap I Les Inondations dans le Monde
Introduction
Depuis le début de ce siècle les inondations ont causé la mort de 2, 250,000 personnes dans le
monde.
En 1931, 4 000,000 morts en Chine
En 1975 : 30,000 morts au Bengladesh.
En 1972 : plusieurs centaines de mort au philippine.
En 1982 : 5,000 morts au Pérou.
En 1987 : 1000 morts à l’Equateur, 500 morts au Colombie et 1200 au Venezuela.
En 1988 plus que 500 morts à Rio de Jeneiro.
En 1979 : 1500 morts en chine et la destruction de 300 villages.

Le risque inondation est la conséquence de deux composantes :

- la présence de l'eau, qui s'écoule habituellement dans son lit mineur, mais qui peut aussi en
sortir occasionnellement pour recouvrir une partie ou la totalité du fond de la vallée (lit
majeur).

- la présence de l'homme, qui s'installe dans l'espace alluvial qui a été progressivement
façonné par le cours d'eau, pour y implanter toutes sortes de constructions, d'équipements ou
d'activités. Cette occupation humaine joue un double rôle : d'une part, elle constitue le risque
en exposant des personnes et des biens aux inondations, d'autre part, elle l'aggrave, en amont
comme en aval, en modifiant les conditions d'écoulement de l'eau. Une inondation peut se
traduire par :

 un débordement du cours d'eau, une remontée de la nappe phréatique, une stagnation

des eaux pluviales : inondations de plaine,
 des crues torrentielles (Vaison-la-Romaine),
 un ruissellement en secteur urbain (Nîmes),

L'ampleur de l'inondation est fonction de :

 l'intensité et la durée des précipitations,

 la surface et la pente du bassin versant,
 la couverture végétale et la capacité d'absorption du sol,
 la présence d'obstacles à la circulation des eaux.

HISTOIRE

Il existe de nombreuses histoires de déluge, de l’Inde à la Grèce. Une des légendes les plus
connues est celle de Noé et de son arche, racontée dans le premier livre de la Bible. Toute la
terre est inondée après 40 jours de pluies diluviennes, seul le vieil homme, sa famille et des
douzaines d’animaux différents sont sauvés en se réfugiant sur un bateau en bois. Des
archéologues découvrirent les preuves d’une inondation terrible du fleuve Euphrate en
Mésopotamie (vers 2000 à 3000 ans av. J.C.), certains historiens pensent qu’elle peut être à
l’origine de l’histoire de Noé. Dans la littérature grecque, la légende prétend que Zeus, roi des

4
dieux grecs, fut tant dégoûté par les hommes qu’il créa l’inondation. Deucalion, père
d’Hellên, ancêtre de la race hellénique, avait construit un navire à bord duquel il monta avec
son épouse pendant le temps de l’inondation. La première inondation enregistrée à Florence(
Italie ), à cause du fleuve Arno, remonte au 2ème siècle ap. J.C. Léonard de Vinci, grand
ingénieur, architecte et artiste du 16ème siècle, établit des plans pour pouvoir contrôler les
crues du fleuve Arno, mais tous ses plans furent rejetés en blocs.

Les inondations dans le Monde :

- Egypte : La crue du Nil

Le Nil est le plus long fleuve du monde, de ses 6 741 km, seuls 1 508 coulent en territoire
égyptien. Chaque année, à la suite des pluies torrentielles qui s’abattent sur les montagnes
d’Ethiopie et la région des lacs équatoriaux, le Nil sort de son lit et inonde en quelques mois
toute la vallée. La crue atteint l’Egypte fin mai-début juin jusqu’en octobre. Au milieu de
l'été, le fleuve se remplit progressivement d'une quantité si énorme d'eau et de limon en
suspension que ses rives ne peuvent plus le contenir. Il inonde lentement. Pendant la longue
période où les eaux restent au même endroit au fond de la vallée, le limon se dépose sur le sol,
de sorte qu'après le retrait des eaux, le pays apparaît recouvert d'une couche de sédiments
riches et fertiles que les anciens Egyptiens appelait "kemet" ou Terre Noire. Le Nil a permis la
naissance d’une des plus grandes civilisations de l’antiquité. Hérodote (484 - 420 av JC),
historien grec et grand voyageur, a écrit : L’ EGYPTE EST UN DON DU NIL. Depuis
l’époque des pharaons, on trouve dans toute l’Egypte, des nilomètres.

Ce sont en général des escaliers qui descendaient dans les eaux du Nil et qui permettaient de
mesurer la crue du Nil et de calculer ainsi les impôts demandés aux paysans. Depuis les
années 1960, la quantité de sédiments a considérablement baissée du fait de la construction du
Haut barrage d’Assouan.

- La mousson indienne

L'Asie des moussons, qui s'étend au sud et à l'est, est tropicale, chaude toute l'année, sèche
l'hiver et très arrosée l'été. Dans les régions de Mousson, les crues des fleuves peuvent - être
catastrophiques. Chaque année, la mousson d’été dure environ 3 mois. Des semaines entières,
de nombreuses cités sont inondées. Des barrages et des digues sont construits pour protéger
les champs et les habitations. En été l'échauffement des parties méridionales du continent
asiatique entraîne la disparition des anticyclones et leurs remplacement par une série de
dépression, qui s'alignent suivant un axe allant du nord - ouest de l'Inde. Ces pluies sont
apportées par la mousson, nom donné à un vent qui souffle sur l'Asie et change de direction
avec les saisons. De juin à octobre, il est tellement chargé d'humidité que les passants qui ne
peuvent s'abriter ont l'impression de recevoir des seaux d'eau sur la tête. Pendant la mousson,
il peut tomber jusqu'à 3 cm d'eau en 5 minutes. C’est dans ce pays que l’on relève les taux de
chutes de pluie les plus élevées au monde, plus de 600 cm de pluie par an dans certaines
villes. Cette mousson est néanmoins attendue dans certains endroits, car elle permet de faire
des stock d'eau pour l'année, et l'irrigation et la fertilisation de certaines plaines

- Madagascar :
La grande île de Madagascar a été traversée par la tempête tropicale Gloria les 4 et 5 mars
2000, du nord-est au sud-ouest. Suite aux pluies diluviennes, les prairies ont été inondées, les

5
huttes de terre effondrées, les maisons en dur crevassées et de vastes étendues de cultures sont
devenues spongieuses.

- Mozambique :

Après le passage du cyclone Eline, début février 2000, des précipitations diluviennes se sont
abattues sur le Mozambique qui connaît ainsi les pires inondations de son histoire. Devant la
montée des eaux, les arbres servent de derniers [Link] mars, plus de 240 000 personnes
avaient du fuir leurs villes ou leurs villages totalement inondés.

- Inondations au Venezuela (1999 ) :

Cette catastrophe est l'une des plus importante du siècle en l'Amérique latine, depuis celle
survenue au Pérou en 1970. Près de 50 000 morts, 200 000 personnes évacuées du littoral et
de très nombreux dégâts. 20 000 personnes bloquées au principal aéroport du pays. 6 000
militaires sont chargés d'assurer la sécurité

- Inondations au Mexique (1999 ) :

Des pluies diluviennes et des inondations ont également affecté le Mexique et sa capitale,
provoquant la mort de centaines de personnes en quelques heures. Ces pluies se sont abattues
sur l'Amérique centrale pendant deux semaines, et ont provoqué dans cette région les plus
importants dégâts depuis le passage de l'ouragan Mitch qui avait fait 26000 morts et disparus
fin octobre 1998, entraîné des dégâts estimés à plusieurs dizaines de millions de dollars et
l'évacuation de 50 000 mille personnes.

Les inondations en Tunisie

Depuis plusieurs siècles la Tunisie a été le théâtre de beaucoup d’inondations.

Celles qui sont les plus importants dans l’histoire sont :

Année Région Année Région

861 Centre Kairouan 1931 Sahel et Nord.
921 Rakkada (Kairouan). 1932 Centre.
953 Kairouan. 1958 Les côtes orientales.
1020 Plusieurs régions. 1959 Sahel, Gabès, Kairoman.
1325 Gabès. 1962 Tout le Sud.
1456 Centre. 1964 Kairouan.
1612 Plusieurs régions. 1965 Zarzis.
1700 Nord. 1969 Centre et Sud.
1726 Sud Ouest. 1973 Centre et Sud.
1761 Nord Ouest. 1982 Sfax.
1902 Tout le pays. 1990 Centre Sud.

Comme la plupart des pays du globe, la Tunisie n’a pas été épargnée par les inondations,
plusieurs régions de notre pays ont été touchées par cette catastrophe naturelle comme le
montre le tableau suivant :

6
Tableau 2 : Historique des crues en Tunisie

7
Année Dégâts et observations
861 Une inondation remplit le pont de Kairouan

921 Les pluies violentes qui tombèrent sur RAKKADA et Kairouan ruinèrent les constructions.

953 Il y a eu à Kairouan une longue période de vent violents et froids acoompagnés de pluie, puis il tombât une grande quantité de
neige.
988 La rivière pénétra à Mansoura et ruina les maisons
1020 Il y a eu à Friguia un violent orage acoompagné de pierres telles qu’au paravent on n’en avait jamais vues d’aussi grosses et
d’aussi abondantes, elles causèrent des accidents mortels.
1325 Le célèbre voyageur Ibn Batouta dut rester bloqué pendant 10 jours à Gabes à cause des pluies incessantes.

1456 On signale un ouragan accompagné de pluies avec des grêlons gros comme un œuf de poule et même d’avantage.
1612 Des pluies abondantes à la suite desquelles les maisons de Djerba s’écoulèrent. En certains points les pluies durèrent 5 mois.

1700 Des pluies excessives, dans les plaines les récoltes pourrirent.

1726 La pluie tombât pendant 5 mois d’une façon continue. Les pertes furent énormes. A Tozeur une petite pluie de 2 heures avait
démolis des maisons, si elle s’était prolongée plus longtemps, de la ville il ne serait resté qu’un tas de boue.

1931 Les pluies furent particulièrement denses en Novembre dans la région entre Sousse et Sfax et sur une partie de la dorsale. Le
lac de Kelibia débordait et coupait la GP1 à la hauteur de Sidi Bou Ali.
1958 Les régions de Tunisie qui furent affectées par le désastre furent Sfax avec 280 familles sinistrées, Moknine avec 830 familles,
Mahdia avec 124 familles, El Jam et Ksour Essef avec chacune 50 familles.

1959 Le désastre frappa la région du Sahel.

1959oct Les pluies diluviennes qui tombèrent au début du mois d’Octobre , sur les régions de Kairouan et Gabes, enflèrent
démesurément les oueds, qui ne pouvant se déverser dans la mer par suite de marée exceptionnel, inondèrent Gabes, en y
faisant 4000 personnes sans abris dont la moitié étaient des enfants. De nombreuses maisons ainsi que des magasins et dépôts
furent détruit en grand nombre. Il y eu pour quelques 2 millions de dinars de dégâts matériels. D’autres inondations de
moindres importances furent enregistrées à Sfax et Nefta.

1962 A la fin du mois de Novembre 1962 des pluies diluviennes s’abattirent sur le sud de la Tunisie, inondant une vaste région
située prés de la ville de Gabes ou les eaux envahirent toutes la région. Bilan de la catastrophe : une cinquantaine de morts,
5000 sans abris, des dégât matériels considérables….
1969 Toutes les régions de Tunisie ont subies des dégâts, les régions les plus touchées sont : La région de Kairouan (elle a reçue la
pluie au de habituelle de 4 ans en 36 jours), Sousse, Sfax, Gafsa, puis le Kef et Kasserine. Les oueds ont monté de 11 mètres au
dessus de la normale et des villages sont disparus, notons que le débits max enregistré dans l’oued Zeroud été de 18000 m3/s.
32 ponts principaux ont été détruis et la carte de la Tunisie a été changée par la création de nouveaux oueds. La force de l’eau
torrentielle a arraché des blocs de béton encré de 100 tonnes. L’infrastructure de l’irrigation estimé à un coût de 7 millions de
dollars a été dévastée en 6 heures 542 personnes ont perdues la vie, 303974 personnes se sont trouvés sans logés, 70504
maisons détruites. Les récoltes agricoles, les routes, les chemins de fer et les ponts été endommagés. Le total des dommages
est estimé à 100 millions de dollars dont 10% dans le domaine agricole. 200Km de routes ont été détruites et gravement
endommagés en plus de plus de 32 ponts détruits. La destruction de la voie ferré porté un grand préjudice à l’économie du
pays. L’exportation de la matière première a été stoppée. Le phosphate et le fer de Gafsa ne sont pas exploités et les ouvriers
sont en chômage.120 Km de voie ferré ont été détruits et le trafic a été stoppé pendant 2 mois. 15 à 20.000 ouvriers se sont
trouvés sans abris. Plusieurs pays ont porté leurs aides à la Tunisie comme l’Algérie, Australie, Canada, Finlande,
Allemagne……. Des chiffres très intéressants sont à retenir de cette inondation mettant en évidence l’ampleur des
inondations : 18000m3 débit de pointe de Zeroud ; 4 milliards de m3 sont écoulés vers la mer …..
1973déc Les pluies torrentielles ont commencé à tomber pendant la première quinzaine de décembre, le centre et le sud du pays sont les
plus affectés (gouvernorat de Gafsa, Kasserine, Kairouan, Sousse, Sfax et Gabes) .Après c’est le nord du pays qui a affecté,
c’est le cas de BOU SALEM ou plusieurs maisons ont été envahis par les eaux.
Les oueds ont quitté leurs lits. Des villages et des fermes agricoles ont été inondés .Les routes furent coupés, les ponts détruits
et le téléphone et le télétexte interrompus. Il y avait aussi des pertes en vie humaine énormes prés de 100 morts.

1989avr Une forte averse s’est abattue sur le golf de Tunis à l’après midi du mardi 11 Avril 1989. Cette averse a débuté vers 13h et s’est
poursuivie jusqu'à 21h avec une accalmie entre 18 :30 et 19 :15. L’ensemble des observations recueillis sur la zone en question
montre que l’épicentre de l’averse s’est développé sur la partie médiane et l’amont des Oueds Roriche et Gareb qui drainent
toute la région. Le corps principale de l’averse se situe entre 15h et 17 :45 et représente 75% de l’averse totale, les intensités
des pluies enregistrées sont comprises entre 45mm/h sur la périphérie de la région et une centaine de mm dans la zone de
Menzeh6. Les périodes de retour de ces intensités ne sont pas exceptionnelles et ne dépassent pas les 5 années.
Les intensités max en 2 heures sont par contre très élevées et ont atteint un niveau de l’épicentre des valeurs exceptionnelles de
40mm/h ce qui correspond à une reccurence sup. à 100 ans. Pour L’ensemble de l’averse on peut estimer cette intensité
raisonnablement entre 15 et 20 mm/h ce qui représente une intensité de reccurence de 30 années environ.
Les coefficients de ruissellement ont certainement atteint de fortes valeurs en raison de la forte urbanisation, en plusieurs
endroits, ils ont à coup sur frôlé les 100%. Même la partie non urbanisée des bassins versants présente des cœfficients élevés en
raison de leurs reliefs pentus et de la force des précipitations. Toutes ces conditions ont été à l’origine des forts écoulements et
des inondations qui se sont produites sur cette région.

8
INONDATION 2003
Au début de l’année hydrologique 2003/2004, des précipitations remarquables se sont
abattues sur le grand Tunis et ont provoqué des inondations et des dégâts considérables
surtout pour les nouveaux noyaux urbains qui ne respectent pas rigoureusement les plans
d’aménagement préétablis par le ministère de l’équipement.
Sur le plan géographique, la région couverte par les événements pluviométriques
exceptionnels est formée par les gouvernorats de l’Ariana, de Tunis et de Ben Arous.
Sur le plan hydrologique, les oueds les plus importants de la zone affectée par l’événement
sont :
-Bassin versant de l’Oued Guiriana limité à Sebkhet Sijoumi
-Bassin des Oueds Guereb Rourich
-Bassin de l’oued Ennkhilet qui se déverse à Sebkhet Ariana et draine un bassin versant de
17,1 km2 et d’autre bassins d’une superficie égale à 73,6 Km2.
Cette inondation est caractérisée par trois épisodes pluvieux différents :
-Episode des jours 16 et 17 Septembre,
-Episode du 24 Septembre,
-Episode des jours 11 et 12 décembre 2003.

a- Episode du 16 et 17 décembre 2003

Les isohyètes du cumul pluviométrique des 2 jours successifs 16 et 17 Septembre a montré

que le secteur ayant reçu plus que 200 mm est localisé au niveau de l’Ariana, El Menzah,
Ettahrir, Mutuel ville et Tunis Carthage, tandis que l’isohyète 140 a englobé en plus le secteur
de Mannouba et de Jbal Ammar. L’isohyète 100 a couvert pratiquement tout le grand Tunis
sauf le secteur formé par Jdeida, Tebourba et Borj El Amri.
b- L’épisode pluvieux du 24 Septembre 2003 :

9
Cet événement a eu lieu une semaine après l’événement du 17 Septembre, les états des sols
sont très favorables au ruissellement. La durée de cet événement a été très courte et n’a pas
dépassé les 2 heures. Ces pluies torrentielles sont caractérisées par des intensités horaires très
élevées variant de 16,2 mm/h à Raoued à 79 mm/h à Manouba.

c-Episode pluvieux des jours 11 et 12 décembre 2003 :

Les deux épicentres de cet événement sont localisés au niveau de Ben Arous et du Lac
Sijoumi. Au niveau de ce dernier le niveau de l’eau a largement dépassé les 10 mètres ce qui a
envahi les quartiers situés sur la rive Ouest tel que Sidi Hcine et El Mghira.

Tableau n° 3: Pluie de Septembre 2003 (station de Saida Mennoubia)

10
 ARIANA VILLE TUNIS VILLE BEN AROUS MANNOUBA
VILLE
Moyenne
mensuelle de 39 39.2 29 37
Septembre (mm)
Septembre 2003
447.7 213.32 281.2 357.5
(mm)

Description de quelques inondations en Tunisie

Les inondations restent toujours un problème d'actualité surtout au Centre et au Sud

du pays. Les pluies y sont très irrégulières, les cours d'eau assez larges et peu profonds et les
sols sont nus et à relief accidenté. Les inondations restent gravées dans la mémoire des
indigènes et ils les utilisent comme des repères historiques : 1969, 1973, 1982, 1990 et 2003
pour ne citer que les plus récentes et les plus importantes. En effet on se souvient d'un
événement quand il a apporté un changement important dans la vie : c'est le cas des
inondations par les dommages humains et matériels souvent assez lourds. C'est pour cela
qu'on reprend ici sous forme d'une étude sommaire les inondations de1969, les inondations
de 1973, les inondations de 1990 et celles de 2003.

Inondations de 1969
Les dix derniers jours du Septembre 1969 ont été très pluvieux et presque la totalité de pluie
de ce mois a été tombé durant cette période. En fait, le tableau suivant nous donne quelques
chiffres de la pluviométrie tombée sur différents régions de la Tunisie :

Tableau n°1: Comparaison de la pluie du Septembre 1969 avec la normale du mois de

Septembre et la moyenne annuelle

station Pluie Pluie Normale du (2)/(3) (2)/(1)

moyenne moyennedu mois de
annuelle mois de Septembre
Septembre (mm)(3)
(1) (mm)(2)
Nord Tunis 415 34.2 49.1 1.44 0.12
de la Beja 626 38.2 49.9 1.31 0.08
Tunisie
Centre Kairouan 286 36.7 135 3.68 0.47

11
de la Sousse 327 41.4 130.3 3.17 0.4
Tunisie
Sud de Gabès 175 14.2 106.5 7.5 0.61
la Gafsa 152 13.5 98.7 7.31 0.65
Tunisie

Aussi pendant le mois d’octobre trois périodes pluvieuses ont été observées :

*Du 03 au 08 octobre : la pluviométrie est comprise entre100 et 200mm au niveau du sud et

du centre tunisien notamment à Sfax, Gabès et Djerba.

* Du 11 au 17 octobre : c’est une période moins pluvieuse

*Du 21 au 31 octobre : comparable avec celle de septembre avec des pointes de 200 à 300mm

Les régions qui ont enregistrées des inondations sont

 Au centre et sud

 La région de kairouan crue d’oueds Zéroud-Merguellil-Nebhana

 La région de Meknassy (crue d’oued Leben)
 Région de SidiBouzid (crue d’oued Fekka)

 Au nord tunisien

 Dépression de Boussalem (crue d’oued Tessa, Mellègue,

Medjerda)
 Région sud de Tunis (crue d’oued Meliane)
 Région nord ouest de Tunis (crue de Medjerda)

Plusieurs dégâts sont enregistrés lors de ces crues :

Au niveau humain, à Sidi Bouzid, les inondations ont causé la mort de plus de 250 personnes
avec une pluie maximale journalière de 200 mm.

Au niveau du réseau routier, 32 ponts sone endommagés et plus de 100 km de routes sous des
dépots inondés ou démolis. Et aussi au niveau ferroviaire, en effet sur la ligne tunis-Sfax le
pont de 40m sur l’oued Essed a été détruit ainsi que les voies de la gare de Sidi Bou Ali.

12
Sur la ligne goulette, Jerissa, Jerda, 3km de voie été détruite. (PIAS, JSTUCMANN,J. ,1970)

Les inondations de novembre et décembre 1973

La Tunisie a connue trois périodes pluvieuses importants qui sont :
-les pluies du 19, 20 et 21 novembre ;
-les pluies du 3, 4 et 5 Décembre ;
- les pluies du 11,12 et 13 Décembre.
Ces pluies sont tombées essentiellement sur le Cap-Bon, bassin versant du Meliane, du
Merguellil, du Zeroud, du Baiech et surtout le littoral tunisien.
Au cours des deux premières périodes pluvieux quelques oueds sont ruisselés surtout au
niveau du Cap-Bon et la région de Gabès. Les dégâts ne sont pas élevés à part quelques oueds
qui ont occasionnés des coupures des routes lors de leurs débordements. Par contre les pluies
de 11,12 et13 Décembre ont été marquées par leur continuité et leur forte intensité surtout sur
le Sahel et le Sud du pays (région de Gafsa et Djerid). En effet, des crues dévastatrices ont été
produites sur les oueds :
-Miliane et ses affluents ;
-Tous les oueds cotières de Sahel de Sousse et de Sfax ;
-Le Merguellil et Zeroud ;
-Oued Baiech dans la région de Gafsa et ses deux principaux affluents (Sidi
Aich et oued Kébir) ;
-les oueds de la région de Gabès.

Les dégâts ont été catastrophiques dans la région de Gafsa ou les crues ont
provoqués (Kallel M.R, 1973):
-le décès de dizaines de personnes dans des véhicules et des engins emportés ;
-la destruction de plusieurs ouvrages ;
-Coupures des routes desservant la ville de Gafsa et par l’isolement complet du gouvernorat
du reste du pays.

Aussi à la suite des crues des oueds côtières entre Sousse et Gabès, la GP1 a été coupé en
plusieurs points et la voie de chemin de fer Sfax-Gafsa a été endommagée en plusieurs
endroits.

13
 Les inondations de Janvier 1990
Les pluies, pendant les quatre jours (21, 22, 23 et 24 Janvier 1990), ont dépassé la normale du
mois de Janvier de 26 fois à Chaal(365mm), 16 fois à Maknassy(308 mm), 18 fois à Gafsa, 13
fois à Tozeur... etc (N. SAKISS, 1990).
Ces inondations ont causé la mort de 24 personnes et la disparition de 32, l'effondrement total
ou partiel de trente neuf milles logements, la disparition de quatorze mille ovins et la
fermeture à la circulation de vingt huit routes à grande et à moyenne circulation. Le bilan des
pertes des terres agricoles, des arbres, ... est encore plus lourd. (Slimani,2003)

Tableau n° 2 : Comparaison de la pluie du 21-22-23 janvier 1990 avec la normale du

mois et la moyenne annuelle

station Pluie moyenne Pluie du 21-22-23 Normale du mois (2)/(3) (2)/(1)

annuelle janvier 1990 due janvier
(mm)(2) (mm)(3)
(1)
Sidi Bouzid 245 160.8 18.1 8.9 0.7
Meknassy 205 308 18.7 16.5 1.5
Chaal 205 365 14 26 1.8
Gafsa 170 352.1 14 18 1.5
Tozeur 95 134.1 10 13.4 1.4
Gabès 185 127.3 19 6.7 0.7

14
 CH II RAPPELS STATISTIQUES

Variables aléatoires :

Une variable aléatoire est définie comme une variable qui décrit un domaine défini sur [0,+∞[
ou ] -∞, +∞[ d’une façon aléatoire. Pour les variables hydropluviométriques, le domaine de
définition est souvent [0, +∞[ exemple : Pluie débit, ruissellement, infiltration…etc.
On s’intéresse à la pluie et au débit. Les différentes variables aléatoires utilisées sont :
-La pluie : Pluie annuelle, pluie mensuelle, pluie journalière, pluie maximale journalière,
intensité…etc.
-Les débits : débit instantané, débit moyen, débit maximal, débit de point…etc.
Les périodes utilisées pour ces variables peuvent être : heure, jour, mois, année,…Les
données de base de la pluviométrique proviennent du pluviogramme, alors que celles du débit
proviennent de l’hydrogramme.

Les paramètres statistiques descriptives :

Soit une variable aléatoire à n observations : X : x1, x2, …, xi, …, xn alors on peut définir :

a- la moyenne : notée X ou m
1 n
X   xi
n i 1
b- L’écart type : notée Sx ou T

n
1
On définit la variance V X 
n
 (x
i 1
i  X ) 2 on a Sx = √Vx

1 n 1 n 2 n n
On peut écrire VX =  i
n i 1
( x  X ) 2
  i 
n  i 1
x  2
i 1
x i X  
i 1
X 2]
n n n
1 1 1
= 
n i 1
xi2  2 X  xi   X 2
n i 1 n i 1

1 n 2 n 2 1 n 2
=  xi  2 X  X   xi  X 2
2

n i 1 n n i 1

1 n 2
V X   xi  X 2
n i 1

c- Moment d’ordre r : noté Mr ou mr

1 n
Mr =  ( xi  X ) r
n i 1
1 n
M3 n 
( xi  X ) 3
d- Coefficient d’asymétrie : Cas = 3  i 1
SX ( Sx)3

15
 1 n
M4 n 
( xi  X ) 4
e- Coefficient d’aplatissement : Cap = 4
 i 1
Sx (S X ) 3

Pour une variable X quelconque, la représentation du Cas est :

F(x) est densité de probabilité ou fréquence calculée avec f(xi) = ni/n avec ni : le nbre
d’apparition de xi.

asy = 0→ répartition symétrique des valeurs autour de la moyenne.

asy <0→la majorité des valeurs sont situées au dessous de la moyenne.

Ex : X : 235, 425, 238, 325, 206, 408, 1080, 230, 980 :X = 600
<X <X <X <X <X <X <X <X <X
asy>0→La majorité des valeurs sont situées au dessus de la moyenne.
X: 235, 425, 238, 325, 206, 408, 45, 230, 25 X = 200

>X >X >X >X >X >X <X >X <X

Pour une variable x quelconque la représentation du Cap est
Cap→ disposition des valeurs autour de la moyenne.
X = 250 Cap = 3→ 230, 210, 215, 245, 255, 260, 280, 270, 275
X = 250 Cap > 3→ 170, 180, 195, 205, 220, 230, 245, 250, 255, 265, 270, 280, 275, 290
X = 250 Cap < 3→ 235, 245, 240, 2 5, 265, 250, 266, 270
Cap et Cas sont 2 paramètres qui indiquent le type de la loi statistique qu’on cherche à
ajuster :

Loi normale Cas→ 0

Cap→ 3
Loi weibull Cas>0
Cap<0
Loi de Gamma Cas<0
Cap<3 ou Cap>3 selon le paramètre .

Les lois statistiques : ajustement :

a-Les lois statistiques :

loi normale, loi log-normale, loi Gamma, loi poisson, loi Gumbel, loi Frech… (cours
hydrologie statistique).
b- Méthode d’ajustement :
Méthode des moments.
Méthode du maximum de vraisemblance.
(cours hydrologie statistique)

c-Pratique de l’ajustements :
-On calcule les paramètres statistiques de la variable : X, 5X, Cas, Cap.
-On choisi la loi en fonction des données et des paramètres statistiques :

*données pluviométriques moyennes → loi normale.

*données pluviométriques tronquées (au delà d’un seuil) → loi log-normale.
*données pluviométriques minimal → loi de weibull.

16
*données pluviométriques maximal → loi de Gumbell, loi Gamma.
*débits moyens, débits instantanées → loi normale.
*débits au delà d’un seuil → loi log-normale.
*débits minimum → loi de weibull.
*débits maximal → loi Gumbell, Gamma.

-On calcule les paramètres de la loi choisie par la méthode des moments en max de
vraisemblance.
-On classe la série par ordre croissant.
-On calcule les fréquences de non dépassement. La valeur x i dans la série classée aura
i  0,5
f ( xi )  avec i : le rang de xi.
n
On transforme les valeurs de F(xi) en Ui = Log (F(xi)) → loi normale → loi log-normale.
Ou Ui = -Log(-Log(F(xi)) → loi Gamma → loi Gumbell, Frêchet.
On représente sur le même graphique et en utilisant les échelles des abscisses transformées ,xi
en fonction de Ui et la loi théorique.

-On juge «à l’euil » la qualité de l’ajustement.

Si la loi est acceptée alors on cherche les quantiles x10, x20, x50, x100 pour des périodes de
retour 10, 20, 50, 100, si non on essaye une autre loi.
 l’unité de la période de retour T =10, 20, 50, 100,…etc est l’année où la série des x i est
une série annuelle (c.à.d une valeur/an) .
1
T est donnée par T 
1 F
1
Si F = 0,90 → T   10ans
1  0,90
1
Si F = 0,95 → T   20ans
1  0,95
1
Si F = 0,98 → T   50ans
1  0,98
1
Si F = 0,99 → T   100ans
1  0,99

1
Si F = 0,999 → T   1000ans
1  0,999

Liaisons entre deux variables :

Soient deux variables à n observations chacune :

X : x1, x2,…xi,…xn
Y : y1, y2,…yi,…yn
1 n
On peut définir la covariance entre X et Y avec cov (X,Y) =  xi yi  XY
n i 1
On peut définir la Corrélation r / -1 ≤ r ≤ 1 entre X et Y
Cov ( X , Y )
r
SX SY

17
 1 n
Cov (X,Y) =  xi yi  XY
n i 1
1 n 2 1 n 2
SX= 
n i 1
xi  X 2 et SY = 
n i 1
yi  Y 2

1 n 1 n
X=  xi
n i 1
; Y=  yi
n i 1

On peut définir si r  0,50 Une régression linéaire simple entre X et Y et on écrit :

y i  axi  b  ei
On cherche a et b par la méthode des moindres carrés (cours hydrologie statistique).
Cov ( X , Y ) 1 n
On a a= avec Var( X )   xi2  X 2
Var ( X ) n i 1
b  Y  aX
Cov( X , Y ). SY Cov( X , Y ) SY
a  .
SX SX . SY SX SY SX
SY
ar
SX

18
 CH III Détermination du débit décennal sur un
basin versant non jaugé Méthode Socose

Introduction :

La méthode est utilisée sur des bassins versants de superficie allant de 2 à 2000 km2, dans une
zone rurale et sans observation de débits.
L’information utilisée est : les caractéristiques du bassin versant et les paramètres
pluviométriques du bassin versant ou des zones avoisinantes.

Les paramètres utilisés :

On détermine le débit de crue en utilisant les 2 paramètres suivants :

-D : la durée pendant laquelle le débit a dépassé Qmax/2.
-J : la rétention potentielle du bassin versant (BV).
On prend Qmax = Qd = débit décennal.
Le principe de la méthode :
On part d’une pluie de projet quelconque : un enregistrement pluviométrique d’une averse
donnée :
Le calcule de ruissellement à partir d’une pluie donnée se fait selon la formule empirique
suivante :

( P(t )  0,2 J ) 2
R(t ) 
P(t )  0,8 J

J : interception potentielle du BV en t.
P(t) : la pluie cumulée à l’instant t.
R(t) : ruissellement cumulé à l’instant t.

Ex:

th 0 1 2 3 4 5 6 7 8 9 10 11 12
P(st) (mm) 0 10 15 33 45 36 31 26 22 15 12 10 5
P(t) 0 10 25 58 103 139 170 196 218 233 245 255 260
R(t) (mm) 0 0 1,2 1,5 15,1 32,8 51,2 68,35 83,82 94,8 103,8 111,4 115,2

19
 Pluie horaire d'une averse (mm)

50
Intensité (mm/h)

40
30
20
10
0
1 2 3 4 5 6 7 8 9 10 11 12 13
temps(h)

J=200m

Représenter graphiquement R(t) en fonction de t.

La quantité d’eau totale ruisselée = 115,2mm.
La pluie totale tombée est 260mm.
115,2
Le coefficient de ruissellement global Cr   44%
260
Cr est fort → la formule donnée ci-dessus (établie en France) ne peut pas être appliquée
sur des bassins versants tunisien du centre et du sud (où Cr est très faible =10 à 15%). Elle ne
peut appliquée que sur des BV de l’extrême Nord.
Ce ruissellement va être transféré dans le cours d’eau à l’aide d’une fonction de transfert.
fonction fonction
de production de transfert
Pluie brute → ruissellement → débit
( P(t )  0,2 J ) 2
R(t )  -hydrogramme unitaire.
P(t )  0,8 J
- hydrogramme de Nash.
Ou autres formules.
- hydrogramme H2U.

t
( )4
1 D
On choisie une fonction de transfert empirique du type : q (t ) 
1,53D t
1  ( )8
D
T : temps
Qd
D : durée en heures pendant laquelle le débit a dépassé
2

Remarque : J et D sont donnés par les formules empiriques ci-dessous.

On va prendre à titre d’exemple D = 3heures puis on va calculer les coefficients de la
fonction de transfert.

T 0 1 2 3 4 5
q(t) 0 2,610-3 0,94 0,32 0,026 7,710-3

20
 fonction de transfert

1
0,9
0,8
0,7
0,6
u(t)

0,5
0,4
0,3
0,2
0,1
0
1 2 3 4 5 6
t(h)

La réponse du BV est très rapide comme la montre la fonction de transfert.

On fait la convolution entre q(t) et r(t) (ruissellement horaire) pour avoir Q(t) dédit en mm/h.
Q(t )  S10 3 Q(t )  S 3
Q(t) (mm/h) →  m /s
3600 3,6
EX : surface du BV= 760km2

On cherche r (t)

T 0 1 2 3 4 5 6 7 8 9 10 11 12
R(t) 0 0 1,21 0,37 13,59 17,7 18,4 17,1 15,5 11 9 7,5 3,8

Temps de réponse du BV = 4 heures

Convolution : Q(t )  i  j q(i )r (i  j )

Q(1) = q(1) r(1)
Q(2) = q(2) r(1)+q(1) r(2)
Q(3) = q(3) r(1)+q(2) r(2)+q(1) r(3).
Q(4) = q(4) r(1)+q(3) r(2)+q(2) r(3)+q(1) r(4)
Q(5) = q(5) r(1)+q(4) r(2)+q(3) r(3)+q(2) r(4)+q(1) r(5)
.
Q(13) =…
Q(0) = 0
Q(1) = 0 r(0) = r(1) = r(2) = r(3) car le ruissellement est nul.
Q(2) = 0 q(1) = 2,6 10-3
Q(3) = 0
Q(4) = 35,3310-3mm/h = q(4) r(1)+…+ q(1) r(4)
Q(5) = 12,83mm/h = 0 +2,610-3 . 13,6 = 35,3310-3mm/h
= q(5) r(1)+… + q(2) r(4) + q(1) r(5)
= 0,94 . 13,6 + 2,610-3 . 17,7
Q(6) = 0 + 0 + 0 + q(3) r(4) + q(2) r(5) + q(1) r(6)
0,32 . 13,59 + 0,94 . 17,7 + 2,610-3 . 18,4 = 21,21mm/h
Q(7) = q(7) r(1) + q(6) r(2) + 0 + q(4) r(4) + q(3) r(5) + q(2) r(6) + q(1) r(7)
=0 + 0 + 0 + 0,026 . 13,59 + 0,32 . 17,7 + 0,94 . 18,4 +

21
 2,610-3 . 17,1
Q(7) = 23,5 mm/h.
Q(8) = 22,56mm/h.
Q(9) =20,70mm/h.
Q(10) =15,9mm/h.
Q(11) = 12,53mm/h.
Q(12) =10,34mm/h.
Q(13) =q(5) r(9) + q(4) r(10) + q(3) r(11) + q(2) r(12)
= 7,710-3 . 12,83 + 0,026 . 9 + 0,32 . 7,5 + 0,94 . 3,8
= 6,3mm/h.
Q(14) = q(5) r(10) + q(4) r(11) + q(3) r(12)
7,710-3 . 9 + 0,026 . 7,5 + 0,32 . 3,8 = 1,48mm/h.
Q(15) = q(5) r(11) + q(4) r(12) = 7,710-3 . 7,5 + 0,026 . 3,8
Q(16) = q(5) r(12) = 7,710-3 . 3,8 = 0,03mm/h.

T 1 2 3 4 5 6 7 8 9 10 11 12 13
Qmm/h 0,035 12,80 21,21 23,35 22,56 20,7 15,9 12,53 10,34 6,3 1,48 0,15 0,03
Qm3/s 7,38 2707 4477 4926 4760 4367 3356 2643 2181 1329 312 31 6

hydrogramme de la crue ( en mm )

25

20

15
Q(t) mm

10

0
1 2 3 4 5 6 7 8 9 10 11 12 13
temps

L’hydrogramme a été établi en prenant des valeurs approximatives de J et de D.

Procédure pratique :

1-Ressemblement des données :

S : superficie du BV en km2
L : longueur du talweg le plus long.
L est mesurée sur une carte depuis l’extrémité sur la périphérie jusqu’à l’exutoire (km).
P : pluie décennale journalière du poste pluviométrique le plus proche(mm).
b : l’exposant dans la formule I-D-F (intensité – durée – fréquence).
i = at-b

22
 Pa: pluie moyenne annuelle du poste pluviométrique le plus proche(mm).
ta : température moyenne annuelle en °c .
2 – Calcule des grandeurs D et J :
Pa 1
D est donnée en heure par : Ln (D) = -0,69 + 0,32 Ln (S) + 2,2 
P ta
J est exprimée en mm :
S Pa
J = 260 + 21 . Ln ( )  54
L P
3 – Calcul intermédiaires :
On calcule K :
24 b P
K=
S
21(1  3 )
30 D
0,2 J
Puis ρ = 1-
K (1,25D)1b
Résultats : les débits :
Ks 2
on a Qd : débit décinnal : Qd(m3/s) = Q  si   0
(1,25D) b 15  12 
Qd = 0 si ρ ≤0
Avec : Q : coefficient correcteur proche de 1 en fonction de ρ et de b.

Calcul des débits de période de retour > 10ans :

On suppose que J et D restant constants.

J : Interception potentielle ne dépend pas de T.
Q max
D : Durée pendant laquelle QT est supérieur à reste constante.
2
24 b P'
On calcule la valeur de K en utilisant la forme suivante : K
s
21(1  )
303 D
Avec : P’ = P* - P
P* : pluie de période de retour T = 50 ou 100 (par exemple).
Si on suppose en plus que les pluies maximales suivant la loi de Gumbel (hydrologie
statistique première A) alors on écrit :
║P’ =
1
4,6P  2,25P * ( P *  P) ln( ln( F ))
2,35
0,98 ou 0,99
50ans 100ans

Ks 2
On calcule en suite ρ : si ρ >0 alors QT = Q
(1,25D) b 15  12 
Si ρ ≤ 0 alors QT = 0

23
Application :

Calculer Q10, Q50 et Q100 pour le bassin versant dont les caractéristiques sont :

S = 760km2 P10 = 160mm =P P100 =250mm Pa =400mm

L = 60km P50 = 190mm b =0,55 ta = 22°c
Etapes :
1) calcul de D et J
2) calcul de K
3) calcul de ρ
4) calcul de Q10
5) calcul de P’(en remplaçant P* par P50 )
6) calcul de K, puis de ρ, puis de Q50.
7) calcul deP’ ( en remplaçant P* par P100)
8) calcul de K, puis ρ, puis Q100.
Solution :
400 1
1) ln (D) = -0,69 + 0,32 ln(760) + 2,2 
160 22
D =9 heures.
760 400
J = 260  21ln ( )  54  J  227,93
60 160
(24) 0,55 .160
2) K=  K  30,3
760
21(1  3 )
30 9
0,2 . 27,93
3)   1    0,5
30,3(1,25.9)10,55
30,3.760 (0,5) 2
4) Q10 = 1. .  Q10  168,9m 3 / s
(1,25.9) 0, 55
15  12.0,5

5) P’ =
1
4,6.160  2,25.190  (190  160) ln( ln(0,98))  P'  181mm
2,35
(24) 0,55 .181
6) K=  K  34,33
760
21(1  3 )
30 9
0,2 . 227,93
  1    0,55
34,33 (1,25.9)10,55
34,33. 760 (0,55) 2
Q50  1. .  Q50  248,2m 3 / s
(1,25 . 9) 0 , 55
15  12 . 0,55

7) P' 
1
4,6 .160  2,25. 250  (250  160) ln( ln(0,99))  P'  250mm
2,35

24
 (24) 0,55 . 250
8) K  K  47,41
760
21(1  )
30 3 9
0,2 . 227,93
  1    0,67
47,41(1,25 . 9)10,55
47,41. 760 (0,67) 2
Q100  1. .  Q100  613,89m 3 / s
(1,25. 9) 0, 55
15  12 . 0,67

25
 CH : IV Les formules empiriques

Estimation de l’apport moyen annuel :

Soit un bassin versant avec un cours d’eau sans observations de ruissellement, alors on utilise
les formules empiriques.
a-Formule de TIXERONT :
P3
R avec : P : pluie moyenne annuelle en m.
3E 2
E : évaporation annuelle en m.
R : Ruissellement en m.
En fonction des régions :
P3
-Le Nord on a E  1m donc R 
3
P3
-Le Centre on a E  1,15m donc R 
4
3
P
-Le Sud on a E  1,3m donc R 
5
b- Formule algérienne :
R  s.P.(1  10  k P )
2

R : ruissellement en m.
S : superficie du BV en km2.
P : précipitation moyenne annuelle en m.
k : a-0,01 log(s) avec a = 0,185
s: superficie en km2.
c- Méthode de Chaumont:
R = σ (P- μ)3/2.
R : ruissellement en m.
P : précipitation moyenne annuelle en m.
σ : coefficient d’écoulement.
μ : facteur de rétention du sol.
Le coefficient d’écoulement σ dépend de la pluie, on a :

Pmm 200 300 400 450 500 550 600 650 700 750 800 850 900
Σ 0,03 0,06 0,095 0,12 0,16 0,21 0,25 0,30 0,34 0,38 0,42 0,45 0,47

Cette méthode n’est applicable qu’au Nord de la Tunisie et qu’au centre.

μ : il dépend de la nature du sol et de la végétation.

La formule de Chaumont est applicable dans le Nord et dans le Centre.

Μ(mm) Nature du BV
0 Massifs calcaires et zones karstiques
100 Marne et argile sans boisement
150 Marnes et argiles cultivés ou boisés
150 à 200 Marnes et argiles cultivés ou boisés avec des grés en couverture
200 à 300 Grés avec cultures ou forêts

26
d- Formule de TURC :
Le ruissellement est donné par de déficit d’écoulement(D) :
1) D = P-Q Q : ruissellement.
1 0,9 1
2) 2
 2  2
D P T
3) T = 300 + 25T + 0,05T3.
4)
Q : lame d’eau écoulée ou ruisselée en mm.
P : pluie moyenne annuelle en mm.
D : moyenne annuelle du déficit en mm.
T : moyenne annuelle de température en °c.

On estime D par (2) puis Q par (1).

Si le déficit D < P → Q = P-D.
Si D > P → Q = 0.
Formule Turc 2 : (pour les pluies faibles).
5 1
P Zm 3
Q = ( )3 ( )
30 100
Q : lame ruisselée en mm.
Zm : altitude médiane en m.
P : pluie moyenne en mm.
Zm est donnée par la courbe hypsométrique du BV.

e- Méthode de Fersi : (applicable pour le Sud tunisien)

R  KP IG
K ≈ 16,3910-3 pour le Sud.
P : pluie moyenne annuelle en mm.
R : ruissellement en mm.
IG : indice de pente globale en m/km. (voir cours hydrologie générale).

f- Méthode de Ghorbel :
Le = a  b( s  so)x  xo   c
Le : lame ruisselée en mm.
S : superficie du BV en km2.
So : superficie limite de S.
a, b, c : constantes qui dépend de x.
x = P 3 H méd  H ex
P : pluie moyenne annuelle en m.
Hméd : altitude médiane en m.
Hex : altitude de l’exutoire en m.
xo : valeur limite de x en fonction de s et x.
xo et so sont mal définies. C’est une formule qui contient beaucoup de paramètres donc
beaucoup d’incertitudes.

g- Méthode de D.G.R.E :

27
Pour les BV de l’extrême Nord et les affluents gauches de la Medjerdah on peut définir 3
groupes :
Groupe 1 : BV peu ruisslnts (perméables) :
Le = 0,26 P-80
Groupe 2 : BV moyennement ruisselants :
Le = 0,26P
Groupe 3 : BV très ruisselants :
Le = 0,26P + 150
Le : lame ruisselée en mm.
P : pluie moyenne annuelle en mm.
Ce son des formules trop simplifiées et ne sont pas applicable pour certaines conditions au
limites.
Ex : pour le groupe 3 : pou P = 100mm → Le = 175 mm
Le > P (pas logique)
h- formule de FRIGUI :

 BV de l’extrême Nord et de l’Ichekeul :

domaine d’utilisation : s = 30 à 1100km2.
Alt = 100 à 400m.
Pluie entre 500 et 1300mm.
Le = 8,18 (P-500)0,58.
Le et P en mm.
 BV de la Medjerdah : s =18 à 21,000km2.
P = 300 à 1300mm.
Le = 7,92 10-7 P2,84 ;Le et P en mm.

Estimation des apports fréquentiels :

Les formules sont utilisées lorsqu(on ne dispose pas d’analyse statistiques sur les séries
d’apport observés. Pour les BV tunisien on a souvent utilisé les deux formules suivantes :

a-Pour le Sud tunisien :

Formule de Fersi :Le (T) = 0,21 Lo T0,89.
Lo : apport moyen annuel calculé par l’une des formules de III .1 en mm.
T : Période de retour en années.
Le (T) : lame écoulée ou apport fréquentiel pour la période de retour T.

b-Formule de Ghorbel :
A(T )
On a :  RTV
Amoyen
A(T) : apport fréquentiel pour une période de retour T.
Amoy : apport moyen annuel calculé par les formules de III . 1.
RTV : rapport donnée selon les régions.

Zone I (le Nord) : RTV = 0,92 log (T) + 0,66. (log néperien)
Zone II (le Centre) : RTV = 1,27 log (T) + 0,45.
Zone III (le Sud) : RTV = 3,65 log (T) – 3,60.

28
Estimation des débits fréquentiels :

a - Méthode rationnelle :

QT = A I(T) . C/3.6
QT: débit de pointe pour une période de retour T en m3/s.
C : coefficient de ruissellement.
I(T) : intensité de l’averse estimée à partir des courbes I-D-F pour une période de retour T en
mm/h.
A : superficie du BV en km2.
Condition d’application : A < 200km2.
Tc : le temps de concentration du BV considéré.

b- Formule de Burkli-Ziegler :
QT  A . C . I (T ) . S / A
QT: débit pour une période de retour T en m3/s.
A : Superficie du B.V en km2.
I(T) : intensité à partir des courbes I-D-F.
S : pente du BV (en %).
C : coefficient de ruissellement.

c- Méthode de Ghorbel :
QT
 RTQ
Qmoy
[Link] 1
avec : Qmoy  S 0,8 (1,075( . )  0,232) Nord
L IC
P : pluie moyenne annuelle.
DH : Hmax-Hex
L : longueur du talweg.
IC : superficie en km2.
Qmoy = 85 Log (s) Sud.
S : superficie en km2.
Q en m3/S.

Les formules de RTQ :

zone I (Le Nord) : RTQ = 1,33 Log (T) + 0,46.
zone II (Le Centre) : RTQ = 1,074 T0,4 – 0,71.
zone III (Le Sud) : RTQ = 1,47 T0,4 – 1,35.

d- Formule de Fersi : (Sud) :

QT  2,410 3 Po . IG .T 0,89 . S .
avec : QT : débit fréquentiel en m3/s.
Po : pluie moyenne annuelle en mm.
IG : indice global de pente.
S : superficie du BV en km2.
T : période de retour en années.

e- formule de Kallel :

29
 T 0, 41
QT = qr .S
S
Avec : qr = 5,5 Nord
Qr = 12,5 Sud
S : superficie du BV en km2.
QT : débit fréquentiel en m3/s.
T : période de retour en années.

Estimation du temps de concentration :

a-Formule U.S.C:
Tc = 0,945 (L1,55/H0,385)
Tc: temps de concentration en heures.
L : longueur du taweg en km.
H : dénivellée entre les extrémités du taweg en m.
Condition d’application : dans les zones arides ou semi-arides et pour des petits
BV accidentés et sans végétation.

b-Méthode de Passini :
Tc  B A L / I
Tc : temps de concentration en heures.
B : coefficient régional (≈1).
L : longueur du talweg (km).
I : pente du BV en %.
A : superficie du BV en km2.
Condition d’application : BV de superficie 25<A<150km2.

c-Méthode de Ventura :
Tc  1,27 A / I
Tc : temps de concentration en heures.
A : superficie du B.V (km2).
I : pente du BV en %.
Condition d’application :
Superficie 1<A<25 km2.

d- Méthode de Giandotti :
Tc  (4 A  1,5L) /( 0,8 H )
Tc : temps de concentration en heures.
H : dénivellée en m.
L : longueur du talweg en km.
S : superficie du BV en km2.
Condition d’application :
10 < A < 50km2.

e- Méthode de Kent’s :
Tc  0,1203 ( L0,8 / I 0,5 )(1000 / CN 1, 69  9)1, 67
Tc : temps de concentration en heures.
L : longueur de talweg en km.

30
I : pente du BV en %.
25400
CN 
254  s
Condition d’application : A > 10km2.
↓ superficie.

31
 CH V : Calcul des apports moyens et fréquentiels et des
Débits fréquentiels sur quelques BV du Sud tunisien.
(Médenine et Tataouine)

Les formules à appliquer

apport moyen → Formule de Fersi.

apport fréquentiels → Fersi, Ghorbel.

débits fréquentiels → Fersi, Kallel, Ghorbel.

Les paramètres utilisés :

Apport moyen → Formule de Fersi → P ;

Apports fréquentiels → Fersi → Lo
Ghorbel → Lo
Débits fréquentiels → Fersi → P ; IG
Ghorbel → S
Kallel → S
P: pluviomètrie.
IG: indice de pente m/km.
S : surface en km2.

Les bassins versants :

Apport moyen → Fersi R = KP IG

K = 16,39 10-3

BV Superficie km2 IG m/km P mm Apport moyen 106 m3

EZ-ZEUSS:
-Bouramli 83 12,5 180 0,86
-Medjerdah 46 13,8 180 0,504
-Ez-zeuss 256 6,9 180 1,98
OUM-ZESSAR :
-Lahimar 26 18,3 150 0,27
-Negueb 1 27 19,6 160 0,31
-Nsim 24 15,3 150 0,23
-Mogar 76 16,6 150 0,76
-Negueb 2 142 15,5 160 1,46
-Abid 10 13,8 150 0,09
-Hallouf 1 50 14,3 140 0,43
-Hallouf 2 218 10,2 140 1,81
-Moussa 17 10,5 150 0,13
-Machchana 8 8,9 170 0,006
-Koutine 10,5 6,5 180 0,079
-Oum zessar 368 6,8 170 2,67
BOU HAMED:
-Egrager 26 11,5 80 0,112
-Biroucha 56 10,5 90 0,267
-El oussif 36 8,9 95 0,167
-El Ghordhab 13 9,8 120 0,098
-Ben Abaîra 54 10,5 130 0,372
-El Himame 42 10,5 140 0,312

32
-Bou hamed 257 8,9 170 2,13
-Esseder 76 7,5 180 0,614
-Bou hamed 463 8 160 3,43
FESSI:
-Ghomrassen 84 15,2 80 0,429
-Tataouine 1 496 12,1 80 2,26
-Affluent 1 164 11,5 90 0,82
-Affluent 2 104 10,5 90 0,497
-Tataouine 2 922 12,6 80 4,29
-Ghabou 106 11,5 90 0,53
-El Gouafel 36 10,5 90 0,172
-El Ghar 122 14,8 100 0,769
-Affluent 3 44 12,5 110 0,28
-Affluent 4 42 11,5 110 0,256
-Fessi 1 1112 10,5 120 7,08
-Affluent 5 24 9,8 120 0,147
-Affluent 6 36 8,7 130 0,220
-Fessi 2 1230 9 130 7,86
-El Gdah 1 68 8,6 150 0,490
-El Gdah 2 260 8,5 150 1,863
-El Gdah 468 8,3 150 3,31
-Saheg 220 9,8 160 1,806
Fessi 2230 8,4 170 18

Tableau des apports fréquentiels :

BV Apport moy 106 V5(*) 106 V10 106 V20 106 V50 106 V100 106
m3 m 3
m 3
m 3
m 3
m 3

1 0,86 0,97 1,4 2,6 5,87 10,88

2 0,504 0,57 0,82 1,52 3,44 6,38
3 1,98 2,25 3,22 5,98 13,52 25,05
4 0,27 0,3 0,44 0,81 1,84 3,42
5 0,31 0,35 0,5 0,94 2,12 3,9
6 0,23 0,26 0,37 0,69 1,57 2,9
7 0,76 0,86 1,24 2,3 5,19 9,61
8 1,46 1,66 2,38 4,41 9,97 18,47
9 0,09 0,1 0,15 0,27 0,61 1,14
10 0,43 0,49 0,7 1,3 2,93 5,44
11 1,81 2,06 2,95 5,5 12,36 22,9
12 0,13 0,15 0,21 0,4 0,88 1,64
13 0,006 0,007 0,01 0,02 0,04 0,076
14 0,079 0,09 0,13 0,24 0,54 0,99
15 2,67 3,04 4,35 8,06 18,23 33,78
16 0,115 0,13 0,19 0,35 0,78 1,45
17 0,267 0,30 0,43 0,8 1,82 3,38
18 0,167 0,19 2,72 5,04 11,4 21,13
19 0,098 0,11 0,16 0,3 0,67 1,24
20 0,372 0,42 0,6 1,12 2,54 4,7
21 0,312 0,35 0,51 0,94 2,13 3,95
22 2,13 2,42 3,47 6,43 14,54 26,95
23 0,614 0,7 1 1,85 4,2 7,77
24 3,43 3,9 5,6 10,4 23,42 43,4
25 0,429 0,48 0,7 1,3 2,93 5,43
26 2,26 2,57 3,7 6,82 15,43 28,6
27 0,82 0,93 1,34 2,5 5,6 10,37
28 0,497 0,56 0,81 1,5 3,39 6,29
29 4,29 4,88 7 12,9 29,3 54,28

33
30 0,53 0,6 0,86 1,6 3,62 6,7
31 0,172 0,19 0,28 0,52 1,17 2,17
32 0,769 0,87 1,25 2,32 5,25 9,73
33 0,28 0,32 0,46 0,84 1,9 3,54
34 0,256 0,29 0,41 0,77 1,75 3,24
35 7,08 8,06 11,54 21,4 48,38 89,66
36 0,147 0,17 0,24 0,44 1 1,86
37 0,220 0,25 0,36 0,68 1,54 2,86
38 7,86 8,94 12,81 23,7 53,67 99,46
39 0,490 0,56 0,8 1,48 3,34 6,2
40 1,863 2,12 3,04 5,63 12,72 23,57
41 3,31 3,76 5,4 10 22,6 41,88
42 1,806 2,05 2,94 5,45 12,33 22,85
43 18 20,5 29,34 54,4 122,9 227,76

(1)Lo = 0,86.106 m3.

Fersi : V5 = 0,21 Lo.T0,89 = 0,21.0,86 106.50,89 = 0,75 106 m3.

V10 = 0,21.0,86 106.100,89 = 1,4 106 m3.
V20 = 0,21.0,86 106.200,89 = 2,59 106 m3.
V50 = 0,21.0,86 106.500,89 = 5,87 106 m3.
V100 = 0,21.0,86 106.1000,89 = 10,88 106 m3.

Ghorbel :
RTV Sud = 3,65 log(T) –3,6 (log néperien)
T = 5 → RTV = 2,27 → V5 = 1,95 106 m3.
T = 10 → RTV = 2,27 → V10 = 4,13 106 m3.
T = 5→ RTV = 2,27 → V20 = 6,3 106 m3.
T = 5 → RTV = 2,27 → V50 = 9,18 106 m3.
T = 5 → RTV = 2,27 → V100 = 11,35 106 m3.
T = 5 on calcul RTV.

A(T ))
On a  RTV  A(T )  RTV . Lo
Lo
(1) Pour T = 5 , T = 10 , T20 → Ghorbel = 2. Fersi.
T = 50 Ghorbel = 1,8. Fersi.
T = 100 Ghorbel = Fersi.
(2) Formule de fersi V5 < Vmoyen (pas logique)

Formule de Fersi établie pour les régions de Sud donc on essaye de garder cette formule tout
en corrigeant son défaut pour T = 5ans.
(*) La formule de Fersi a été modifiée :
V5 = 0,21 Lo T1,05
Pour avoir V5 > Lo.

Débits fréquentiels : formule de Ghorbel

Qmoy = 85 Log (s)

zone III (Le Sud) : RTQ = 1,47 T0,4 – 1,35.

34
 B.V S Q5 Q10 Q20 Q50 Q100
1 83 220,3 292,7 388,9 566,2 752,4
2 46 164 217,9 289,5 421,5 560
3 256 386,9 514 683 994,5 1321
4 26 123,3 163,8 217,6 316,9 421
5 27 125,6 166,9 221,8 322,9 429
6 24 118,4 157,4 209 304,5 404,5
7 76 210,8 280 372 541,8 719,9
8 142 288 382,8 508,7 740,6 984
9 10 76,4 101,6 134,9 196,5 261
10 50 170,9 227 301,8 439,5 583,9
11 218 357 474,3 630,3 917,7 1219,3
12 17 99,7 132,4 176 256,2 340,5
13 8 68,3 90,8 120,7 175,8 233,5
14 10,5 78,3 104,1 138,3 201,4 267,6
15 368 463,8 616,3 818,9 1192,3 1584,2
16 26 123,3 163,8 217,6 316,9 421
17 56 180,9 240,4 319,4 465 618
18 36 145 192 ,7 256 372,9 495,5
19 16 96,7 128,5 170,7 248,6 330
20 54 177,6 236 313,7 456,7 606,8
21 42 156,7 208,2 276,6 402,8 535,2
22 257 387,6 515 684,3 996,4 1323,9
23 76 210,8 280 372 541,8 719,9
24 463 520,3 691,3 918,5 1337,4 1777
25 84 221,6 294,4 391,2 569,6 756,9
26 496 538,5 715,5 950,7 1384,2 1839,2
27 164 309,6 411,4 546,7 795,9 1057,6
28 104 246,6 327,6 435,3 633,8 842,2
29 922 734,2 975,6 1296,2 1887,3 2507,6
30 106 248,9 330,7 439,5 639,9 850,2
31 36 145 192,7 256 372,9 495,5
32 122 267 354,8 471,5 686,5 912
33 44 160,4 213 283 412,3 547,8
34 42 156,7 208,2 276,6 402,8 535,2
35 1112 806,3 1071,4 1423,5 2072,7 2753,9
36 24 118,4 157,4 209 304,5 404,5
37 36 145 192,7 256 372,9 495,5
38 1230 848 1126,8 1497,2 2179,9 2896,4
39 68 199,4 264,9 352 512,5 681
40 260 389,9 518 688,3 1002,2 1331,6
41 468 523,1 685 923,5 1344,6 1786,6
42 220 358,6 476,5 633,9 921,9 1224,9
43 2230 1141,9 1517,2 2015,9 2935,2 3899,9

CH : VI Les Ajustements statistiques et méthode de Gradex

Introduction :

Soit une série d’observation de débits ou d’apport maximum.

Ex : les apports maximum mensuels.
Les apports maximum annuels.
Les débits maximum mensuels.

35
 Les débits maximum annuels.
Si la série est suffisamment longue on peut effectuer directement un ajustement statistique sur
cette série puis on cherche les valeurs de période de retour élevées (50ans, 100, voir 1000ans).
Si la série n’est pas assez longue alors on utilise une méthode dite méthode de Gradex. Elle
consiste à combiner l’information hydrométrique (débits) et l’information pluviométrique.
Certaines hypothèses sont à supposer pour utiliser cette méthode.

Ajustement statistique : (série longue n ≥ 20ans) :

a / Données : débits maximum annuels (station Bou Salem sur Medjerda).

Année 1926 27 28 29 30 36 37 38 39 40 47 48 49 50 51 59 60 61
Débits 425 431 1220 1760 58 150 1420 310 566 1780 911 1700 851 383 191 1140 254 337
m3/s

62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82
603 672 587 449 685 119 167 118 1490 400 221 3180 129 620 210 743 172 303 208 41 208

83 84 85 86 87 88 89
319 580 917 80 776 147 272

31 32 33 34 35 41 42 43 44 45 46 52 53 54 55 56 57 58
578 2060 496 307 894 231 943 X 351 196 743 651 904 470 322 465 255 X

b / Lois utilisées:
Pour les séries de valeur maximum on utilise les lois des extrêmes : loi de Gumbel, loi de
Gamma, loi de Pearson…
La plus simple est la loi de Gumbel.

X  Xo
F ( X )  exp( exp( ))
a
avec a = 0,78 Sx
Xo = X - 0,577 . a.
Les étapes de l’ajustement:
1) calculer X et Sx puis a et Xo.
2) Classer la série par ordre croissant.
3) Faire le tableau des fréquences annulés
i  0,5
F(Xi ) 
n
et effectuer une transformation d’échelle U i   Log ( Log ( F ( X i )))

Rag valeur F(Xi) Ui

1 X1 1  0,5 1  0,5
 Log ( Log ( ))
n n
2 X2 2  0,5
n

4) Tracer sur le même graphique :

La droite de Gumbel Qi = aUi + Xo

36
 Les points expérimentaux Qi = f(Ui)
(Q = X)
5) Si l’ajustement est satisfaisant :
On cherche les débits de période de retour élevée.
Q50 pour T = 50ans.
1
On a la formule T 
1 F
Avec F : fréquence de non dépassement.
Donc si T = 50 → F = 0,98
U = -Log(-Log(0,98)) = 3,9
Q100 → T = 100
1
F = 1
100
U = -Log(-Log(0,99)) = 4,6
1
Q1000 → T = 1000 ans → F = 1   0,999
1000
U = -Log(-Log(0,999)) = 6,90

n = 62
X = 615,6
Sx = 568,3
a = 0,78 . Sx = 443,3
Xo = X – 0,577 a = 359,8
X  Xo
F ( X )  exp( exp( ))
a
rang valeur F(Xi) Ui
1 41 8,06 10-3 -1,57
2 58 2,41 10-2 -1,31
3 80 4,03 10-2 -1,16
4 118 5,64 10-2 -1,05
5 119 7,25 10-2 -0,96
6 129 8,87 10-2 -0,88
7 147 1,04 10-1 -0,81
8 150 1,2 10-1 -0,74
9 167 1,37 10-1 -0,68
10 172 1,53 10-1 -0,62
11 191 1,69 10-1 -0,57
12 196 1,85 10-1 -0,52
13 208 2,01 10-1 -0,47
14 208 0,21 -0,42
15 210 0,23 -0,37
16 221 0,25 -0,32
17 231 0,26 -0,28
18 254 0,28 -0,23
19 255 0,29 -0,19
20 272 0,31 -0,14
21 303 0,33 -0,1
22 307 0,34 -0,05
23 310 0,36 -0,01
24 319 0,38 0,03
25 322 0,39 0,07
26 337 0,41 0,11
27 351 0,42 0,16

37
 28 383 0,44 0,2
29 400 0,46 0,25
30 425 0,47 0,29
31 431 0,49 0,34
32 449 0,5 0,39
33 465 0,52 0,43
34 470 0,54 0,48
35 496 0,55 0,53
36 566 0,57 0,58
37 578 0,58 0,63
38 580 0,6 0,68
39 587 0,62 0,74
40 603 0,63 0,8
41 620 0,65 0,85
42 651 0,67 0,91
43 672 0,68 0,97
44 685 0,70 1,03
45 743 0,71 1,1
46 743 0,73 1,17
47 776 0,75 1,24
48 851 0,76 1,32
49 894 0,78 1,4
50 904 0,8 1,49
51 911 0,81 1,58
52 917 0,83 1,68
53 943 0,84 1,79
54 1140 0,86 1,91
55 1220 0,88 2,04
56 1420 0,89 2,2
57 1490 0,91 2,37
58 1700 0,92 2,58
59 1760 0,94 2,84
60 1780 0,96 3,19
61 2060 0,97 3,7
62 3180 0,99 4,81

Qi = 443,3 Ui + 359,8
T = 50 ans → F = 0,98 = 1 – 1/T
U = -Log(-Log(0,98)) = 3,9
T = 100 ans → F = 0,99
U = 4,6
T = 1000 ans → F = 0,998
U = 6,9

Méthode du Gradex :

Elle est utilisée si la série des débits est insuffisante (n < 20) donc on la complète par
l’information pluviométrique.
Les hypothèses de la méthode :
1) Le sol se sature à partir d’une crue de période de retour supérieure ou égale à 10 ans
donc tout se qui pleut se transforme directement en ruissellement. Donc il y a
parallélisme entre l’ajustement des pluies et celui des débits au delà de la crue
décennale (F = 0,9).
2) On suppose que les pluies suivent une loi à décroissance exponentielle.
3) Le BV doit être relativement perméable et le terrain n’est pas Karstique.

38
Pour avoir l’ajustement sur les débits de pointe on utilise un coefficient R = 1,50 et on aura
l’ajustement des débits de pointe en affectant une multiplication des débits maximum par le
rapport R = 1,50.
Puis on cherche Q50, Q100 et Q1000 en extrapolant la droite des débits de pointe.
Les étapes :
1) calculer X et SX de la série des pluies.
2) Calculer a et Xo des pluies.
3) Classer par ordre croissant les pluies et calculer F(Xi) et Ui.
4) Classer par ordre croissant les débits transformés en mm et calculer les F(Xi) et Ui.
5) Représenter sur le même graphique :
- La distribution Xi et Ui des pluies.
- La droite de Gumbel sur les pluies.
- La distribution Xi et Ui des débits.
6) tracer une droite à partir des coordonnées (F = 0,9 et Q10) parallèlement à celle des
pluies.
7) Effectuer une multiplication de cette droite par R = 1,5.
8) Calculer Q50, Q10 et Q100 en utilisant cette droite.

Q m 3 / s . t ( h ) . 3600
Q m /s 
3
2 6
.10 3
S ( Km ) .10
Q m 3 / s . t ( h ) . 3600
=
S ( Km 2 ) .10 3
Pour le cas de Bou Salem on a :
T(h) : durée de la crue = 120 heures.
S(Km2) : 16483 Km2.

Application:
On suppose qu’on connait les 20 premières années de mesure des débits mais on connait la
série des pluies de durée t = 120 heures.
La pluie à utilisée doit avoir une durée égale au temps de concentration du B.V et on prend la
pluie maximale annuelle de durée t.
Série des débits mesurés :

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19
452 431 1220 1760 58 578 2060 496 307 894 150 1420 310 566 1780 231 943 351 196

20
743

Série des pluies :

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25
62 54 124 218 46 110 86 72 66 120 34 92 60 162 122 76 150 54 42 56 84 62 134 66 38

26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50
128 136 126 36 88 60 130 88 52 54 50 62 76 50 36 46 28 148 58 58 200 40 46 42 50

51 52 53 54 55 56 57 58 59 60 61 62
412 356 501 532 516 62 72 690 22 106 24 18

39
1) X et Sx de la série des pluies :
X = 119,14 Sx = 135,81
2)a = 0,78 SX = 105,93
Xo = X – 0,577 a = 58,02
3) Qi = aUi + Xo
Qi = 105,93 . Ui + 58,02

Par ajustement direct :

Q50 = 2000 m3/s
Q100 = 2400 m3/s
Ajustement par méthode de gradex :
Q50 = 9700 m3/s
Q100 = 12500 m3/5
La méthode de gradex donne des valeurs surestimés car on a pris des valeurs de pluie de
durée t = 120 heures assez élevés (P = 516, 536, 690…).
Si on enlève ces valeurs et on refait l’ajustement on trouvera des valeurs des débits qui
s’approchent à celles estimées par la méthode d’ajustement direct.

40
 CH : VII Dimensionnement des ouvrages
(lac collinaire)
Etudes préliminaires :

On a généralement trois études à effectuer :

- Etude topographique.
- Etude géologique.
- Etude hydrologique.

En fonction de ces trois études on peut classer les sites :

- sites favorables → les études sont bonnes : c’est un site à retenir.
- sites satisfaisants → les caractéristiques des études sont moyennes à bonnes : c’est un site à
retenir avec réserves.
- sites défavorables → les études sont mauvaises : c’est un site à registrer.
Etude topographique :
On cherche pour les sites retenus le rendement topographique en calculant le coefficient
d’efficacité (Ef).
Un site favorable doit avoir Ef > 3.

Volume utile de la retenue

Ef 
Volume de la digue
donc il faut un levé topographique dans la cuvette. Ce levé va permettre d’avoir les différentes
surface de plan d’eau en fonction de la hauteur.

H i 1
Le volume d’eau comprit entre Si et Si + 1 est Vi = . ( S i 1  S i )
2
L’étude topographique permet :
- le choix de l’emplacement de la digue.
- La détermination de la zone inondable par la retenue.
- L’évaluation du remblai de la digue.
- La délimitation des zones d’emprunt.

Etude géologique :
Un site géologiquement favorable :
- permet la réalisation d’une retenue imperméable.
- Assure une bonne assise pour la digue.
- Fournir les matériaux nécessaires à la réalisation.
Il permet aussi d’étudier :
- L’étanchéité de la cuvette et la stabilité des talus.
- La stabilité et l’étanchéité des appuis et des fondations de la
digue.
- Les caractéristique matériaux utilisés.

Etude hydrologique :
on doit tout d’abord effectuer une étude hydrologique sommaire qui permettra de déterminer :
- La superficie du BV (km2) : S.

41
 - Pluviométrie moyenne annuelle (mm) : P.
- L ‘apport moyen annuel (m3/an) : A.
- L’évaporation annuelle (mm) : E.
- La volume utile de la retenue (m3) : Vu.
- Le débit de crue (m3/s) : Q.
- L’envasement moyen annuel (m3/an) : En.
La superficie du BV sera déterminée par planimétrage, la pluie moyenne annuelle par
observation au niveau d’un poste sur le BV ou le plus proche de ce bassin. L’apport moyen
annuel sera déterminé par des formules empiriques (Ghorbel, Fersi, TIXERON, formule
algérienne, TURC…).
L’évaporation sera déterminée à partir des observations sur le BV ou sur le poste le plus
proche.

Le volume utile de la retenue Vu (m3) :

il est déterminé en fonction de l’apport moyen annuel et de la pluviométrie annuelle sur la
région selon le tableau suivant :

P(mm) < 550 550-650 650-750 750-1000 >1000

Vu/A 1,3 1,2 1,0 0,80 0,50

Le débit de crue utilisé pour dimensionner l’évacuateur de crue pour une période de retour T
(souvent 50 ans pour les lacs) est déterminé par les formules empiriques (rationnelle, Ghorbel,
Kallel,…).
L’envasement moyen annuel (En) :
En = 1,4 . Es . S
S : surface en Km2.
Es : érosion spécifique en m3/an/Km2.

Dimensionnement grossier de la digue :

- Hauteur utile de la digue (Hu en m) : Hu correspond à la hauteur de la cote normale de la
retenue.
1 3V
On a Vu  H u . S  H u  u
3 S
- hauteur maximale de la digue (H en m) :
on peut prendre pour première approximation H = Hu + 2.
- Volume de la digue (Vd en m3) : Calculé par la formule :
(L . H )  H 
Vd  l  ( Fm  Fv)
2  3 
L : longueur de la digue en crête (m).
l : largeur de la digue en crête (m).
H : hauteur maximale de la digue (m).
1
Fm : fruit amont ( ) ; Pm : pente du talus amont.
Pm
1
Fv : fruit aval ( ) ; Pv : pente du talus avale.
Pv

Cas d’une vallée triangulaire.

Si on a une vallée trapézoidale alors on a :

42
 H H  H H 
Vd  ( L  L ' )l  ( Fm  Fv )  L' l  ( Fm  Fv )
2  3  2  2 
L’ : largeur du lit de l’oued (m).
Pour les valeurs de L, Fm et Fv on peut prendre dans une première approximation :
L = 3m
1
Fm = 2,5 → Pm 
2,5
1
Fv = 2  Pv 
2

Etude hydrologique détaillée:

KC A  1,12 2 
Description du B.V : l 1  1  ( ) 
1,12  XC 
- La superficie (Km2) S
- Le périmètre (Km) P
XC A  1,12 2 
- La longueur du Talweg (Km) L L 1  1  ( ) 
1,12  KC 
- Le coefficient de compacité Kc = 0,28 P / S
- Les dimensions du rectangle équivalent l, L.
- La hauteur maximale sur le BV (m) : Hmax.
- La hauteur minimal sur le BV (m) : Hmin.
- La hauteur moyenne (m) : Hm = (Hmax + Hmin)/2.
- La pente moyenne du B.V : I = (Hmax – Hmin)/√S.
- La hauteur à 95% , à 5% et à 50% : H95 , H5 , H50.
- L’indice global de pente :
Ig = (H95 – H5)/L
- La dénivelée Spécifique:
DS = Ig √S (m).
- On doit aussi classer les reliefs:

DS (m) Classe Relief

< 10 R1 Très faible
10-25 R2 Faible
25-50 R3 Assez faible

> 500 R7 Très fort

- On doit aussi tracer la courbe hypsométrique :

Altitude Superficie en % % cumul

705-650 6,6 6,6
650-600 14,7 21,3
600-550 20 43,3

450-410 100

43
- On doit aussi étudier la pluviométrie sur le B.V (la pluie mensuelle).
- On estime les crues de projet et les apports fréquentiels en utilisant les formules empiriques
(Fersi, rationnelle, Ghorbel, Kallel,…).
-On estime aussi le transport solide sur le B.V.

Conception et dimensionnement de la digue :

1-hauteur de la digue : (H) donnée par

H = RN + h + R
RN : hauteur de la retenue normale.
H : charge maximale au dessus du déversoir (côte des plus hautes eaux).
R : revanche.
a- hauteur normale (RN) :
c’est la hauteur qui correspond à V tel que : V = VU +VS.
VU : volume utile.
VS : volume des sédiments qui correspond à une durée égale à la durée de vie du lac (50
ans).
VU est déterminé en fonction de l’apport et la pluie et VS déterminé par formule empirique
(voir VI –1).
La hauteur de la retenue normale est déterminée au niveau de l’étude des plans d’eau
successifs et de l’envasement. On trace la courbe hauteur-volume.

Hi (m) Dénévelée H (m) Volume partiel V Volume cumulés

H0 ∆H = 0 Vo =0 0
H1 ∆H =H1 –H0 ∆V = (S1+ S0) ∆H/2 V = V1
H2 ∆H =H2 – H1 ∆V = (S1 + S2) ∆H/2 V = V1 + ∆V
Hi ∆H =Hi – Hi-1 ∆V = (Si+Si-1) ∆H/2 V = Vi-1 + ∆V

Pour chaque volume Vi de la retenue correspond une cote Ci tel que Hi = Ci = Co.
Si on prend ∆H = 1m : constant.
On détermine par planimétrage chaque surface du plan d’eau Si.

h : charge maximale → donnée par le dimensionnement de l’évacuateur de crue.

h : est généralement comprise entre 0,60m et 1,5m.
R : la revanche : c’est la tranche comprise entre la cote des plus hautes eaux et la ligne de
crête de la digue.

On a R = 1 + 0,3 √L (m)
L: longueur de la retenue en Km.
2- largeur en crête (l) :
Pour les lacs collinaires on adopte une valeur de l = 3m.
3- Pente des talus :
Elles dépendent des caractéristiques des matériaux constituant la digue.
Pour les lacs collinaires on adopte les valeurs suivantes.

H(m) Type de la digue Pente des Talus

Amont Avale
< 5m -homogène ½,5 ½
zonée(noyau d’argile) ½ ½

44
5 à 10m -homogène
*Granulométrie étendue ½ ½
*à fort % d’argile ½ ,5 ½
-zonée 1/2 ½,5

Volume de la digue (Vd) :

on adopte le volume calculé en VI-1 et on lui ajoute 15% pour le tassement.
Vd : = Vd . 1,15
Volume d’encrage :

  
Vr  L  l  ( )
 P 
L : longueur de la tranchée.
L : largeur de la base de la tranchée.
∆ : profondeur de la tranchée.
P : pente des berges de la tranchée.

Ouvrages annexes : (évacuateur de crue) :

Les différents éléments : - déversoir ou seuil.

- chenal.
- Coursier.
- Bassin de dissipation.
 Le déversoir :
Il est situé dans la partie supérieure de l’évacuateur de crue à la même cote que la retenue
normale. Le débit (Q) est fonction de la charge d’eau au dessus de seuil de déversoir (h) et du
profil de ce seuil.
Pour un déversoir de longueur L le débit q est donné par Q   . L . h . 2 gh

  0,38 déversoir à large seuil.

  0,48 déversoir à seuil déversant.
h : charge d’eau au dessus de seuil (m).
L : longueur du déversoir (m).
g en m/s2 Q en m3/s
Q est donnée par l’étude hydrologique = QT avec T: durée de vie de l’ouvrage (100, 1000,
même 10 000 ans).
Pour les lacs collinaires, on adopte T = 50 ans.
- Si on fixe L on aura h.
- Si on fixe h (pour les lacs h < 1,5 m) on détermine L.
QT   . L . h . 2 gh
QT
L
 . h . 2 gh
QT
- Si h = 1,5 m alors L 
0,48 .1,5 2.9,81.1,5
  0,48

45
 QT
L
3,9
Pour un lac collinaire QT  100 m 3 / s  L  25 m
 Le chenal d’écoulement :
Il se situe juste à l’avale du déversoir et il doit avoir une pente suffisamment faible (I < IC)
pour assurer un écoulement fluvial. IC : pente critique.
On adopte généralement une forme rectangulaire.
 Le coursier :
Le passage de l’écoulement fluvial se fait par un changement de pente et rétrécissement.
La forme de coursier est rectangulaire. Le coursier peut être réaliser : - en béton, en
maçonnerie, en terre revêtu d’herbes.
La profondeur normal de l’eau dans le coursier est donnée par Manning-Strickler :
2 1
1
Q  R3 I 2 S
n
n : coefficient de manning.
R : rayon hydraulique.
S : section mouillée.
I : pente du coursier.
Le rayon hydraulique et la surface sont fonction de h (la profondeur normale). On fixe la
forme géométrique et la pente.
S R 2 / 3  n . QT / I
n, QT, I sont connus donc n QT / I est connu.
On calcule plusieurs valeurs de S R2/3 pour différentes valeurs de h jusqu’à l’égalité entre
SR 2 / 3 et n QT / I
h trouvée est alors la hauteur normale.
1
Ex : n  , QT  100m 3 / s; i  0,01
60
Formule rectangulaire :
S  L.h  2 / 3 L.h 2 / 3
R  ( )
P  L  2h  L  2h
L.h 2 / 3
S R 2 / 3  L.h ( )
L  2h
L = 20m.

46
 ( L.h) 5 / 3 (20h) 5 / 3
S R2/3    n QT / I  16,66
( L  2h) 2 / 3 (20  2h) 2 / 3
(20.0,5) 5 / 3
h  0,5  6,09
(20  2.0,5) 2 / 3
(20.1) 5 / 3
h  1m  18,76
(20  2.1) 2 / 3
(20.0,9) 5 / 3
h  0,9m  15,44
(20  2.0,9) 2 / 3
(20.0,95) 5 / 3
h  0,95  17,28
(20  2.0,95) 2 / 3
h  0,92  16,41
Donc h = 0,92m
 Le bassin de dissipation :
A la sortie du coursier, il faut prévoir un bassin pour dissiper l’énergie de l’eau.
On utilise un bassin de dissipation qui par création de ressaut hydraulique va transformer
l’écoulement torrentiel au niveau du coursier en un écoulement lent pour le restituer dans le lit
de l’oued.

Dimensionnement :
L : longueur du bassin de dissipation.
L : largeur du bassin.
h1 : tirant d’eau amont du bassin.
h2 : tirant d’eau aval du bassin.

La vitesse d’eau dans le bassin est : V1  0,8 2 gh(h2  h1)

Le débit Q qui passe à travers le bassin est Q = V1 . h1 l.
On assimile h2 à une constante = H = la hauteur d’eau à partir du fond du bassin.
V1  0,8 2 g ( H  h1 )
Pour chercher h1 on donne des valeurs successives à h1 et on calcule à chaque fois V1 et Q et
on s’arrête quand on aura Q = QT = Q50 on prend l = 2 x longueur du coursier.
h
Puis on calcule h2  1 l  8F 2  1
2
V1
F : nombre de froud =
gh1
Et on prend H : hauteur du l’eau à partir du fond du bassin.

La longueur L du bassin est prise égale à L = 6(h2 – h1).

Résumé :

V1  0,8 2 g ( H  h1 )
Q  V1 . l . h1
différentes valeurs de h1 jusqu’à avoir Q = QT → h1.

47
 h1
On calcule h2  l  8F 2  1
2
L = 2 x longueur de coursier.
F = V1 / gh1
L  6(h2  h1 )

48
 Projet

Calcul des apports et des débits fréquentiels sur les cours d’eau du gouvernorat de Kébili

Documents utilisés :

- Cartes d’ état-Major au 1/100.000

- Données pluviométriques
- Courbes d’ ajustements et courbes IDF

Chapitre VIII : Calcul des débits et apports fréquentiels sur les bassins versants du Sud

Introduction

Le gouvernorat de Kébili est situé dans le sud Est du pays . Il est limité au Nord par le
gouvernorat de Gafsa , au Nord Ouest par le Gouvernorat de Tozeur , au Sud Ouest par les
frontières Tuniso-Algériennes , à l'Est par les gouvernorat de Gabès et Mednine et au Sud
par le gouvernorat de Tataouine. La plus grande partie est constitué par des plaines et des
dépressions fermées . Les montagnes sont situées au Nord et au Nord-Est . Il est situé dans
la région tunisienne la plus sèche et la plus chaude , le climat est caractérisé par des pluies
rares avec une très forte variabilité temporelle . Le régime thermique est très contrasté avec
des amplitudes thermiques diurnes , mensuelles et annuelles très élevées . La moyenne des
minima est égale à 3.1 °C et celle des maxima est de 42.2 °C L ‘ évaporation dépasse 1900
mm/an . Les unités hydrologiques les plus importantes dans la région sont Segui , Chereb ,
El Bahaier , Tebaga et Dahar et enfin Bouflidja . On s’intéressera aussi au bassin versant de
Hallouf ( gouvernorat de Gabès et Mednine )dont l’exutoire se trouve dans le gouvernorat
de Kébili . Le reste du Gouvernorat est occupé par des dépressions ( Chott Jerid, Chott
Fejjej, quelques sebkhas et garaas isolés) et des étendus de dunes de sables dans la région d’
El Faouar et Régim Maatoug .Par conséquent le réseau hydrologique dans ces plaines de
sables est diffus et sans aucune importance quant à l’apport hydrique ou l’érosion . On les
évoquera quand on s’intéressera à l’érosion éolienne et les mouvements des dunes de
sables.
Le climat de la région est saharien supérieur ( quotient d’Emberge = 7.6 ) caractérisé par une
aridité quasi-permanente , une continentalité affirmée et une omniprésence des vents . Ils
sont fréquents de novembre à avril de secteurs nord , nord-ouest et sud-ouest , accompagnés
généralement de tempêtes de sable. Les vents du type sirocco soufflent en moyenne 35 jours
par an

Les Principaux bassins versants :

- Les bassins versants de Segui :

a- Présentation Générale

Il est situé au Nord Ouest de Kébili à la limite avec le gouvernorat de Gafsa et à l’Est à la
limite du gouvernorat de Gabès. Les principales montagnes dans la région sont Djebel Berda
au Nord et Djebel Sif Ellaham et Oum Ali au Sud . La région de Ségui constitue alors un

49
couloir entre ces deux chaînes, de longueur 24 km( Est-Ouest) et de largeur variable entre 3 et
6 km (Nord-Sud) . La région est traversée du Sud au Nord par la route de Kébili – Gafsa .
C’est plaine vaste de superficie totale 22400 ha ( y compris les montagnes ). Elle est classée
dans l’étage bioclimatique saharien supérieur. Elle profite de l’effet positif du relief qui
favorise les précipitations locales . La moyenne pluviométrique annuelle y est de l’ordre de
150 mm. La région du ségui est subdivisée en terres collectives . Une plaine dite El Griaâ se
présente sous forme d’une enclave à l’intérieur du gouvernorat de Gafsa fait partie de la
région du Segui . Elle a une superficie de 1225 ha et draine le flanc nord de Djebel Asker.
Elle est située à 30 km à l’ouest de la route de Kébili-Gafsa De même un bassin versant (oued
Battoum) à l’extrême est de la région de Ségui , de superficie 1000 ha est occupée par la
collectivités de Beni Zid d’ El Hamma . Elle fait partie de l’ Unité hydrologique Segui.

Tableau 1 : les caractéristiques physiques des bassins versants de la région de Segui

N° Cours S kc Hmx Hmn(m H5(m) H95 (m) H50 (m) Hmy( D(m) Ig(m/km)
d’eau (km2) (m) ) m)
1 Lettouchi 16.5 1.31 772 190 610 200 300 481 410 57
2 Enf. saad 9.5 1.82 772 200 740 215 410 486 525 58.66
3 Regaieg 5.8 1.57 475 210 455 215 290 342 240 42.10
4 El Gsir 13.5 1.74 926 210 860 230 485 568 630 61.16
5 Nir Karma 5.2 1.83 600 215 542 220 320 381 322 48.79
6 Jedari 7.0 2.0 926 215 860 230 485 570 630 71.60
7 Houinit 4.3 1.35 636 250 515 262 320 443 253 64.87
8 El Blida 7.5 1.84 926 230 835 250 510 578 605 75.62
9 Retam 7.6 2.2 850 220 805 235 365 535 570 55.62
10 [Link] 9.95 1.33 455 200 397 215 480 327 182 31.38
11 Oum Ali 4.6 1.44 416 175 410 187 270 295 223 49.34
12 El Araguib 11.6 1.40 416 150 400 160 250 283 240 35.30
13 Jerouala 5.0 1.25 515 150 430 160 205 332 270 73.93
14 Sania 2.7 1.36 515 150 462 160 185 332 302 96.80
15 Bouzaiene 4.4 1.33 515 150 462 155 230 332 307 79.74
16 [Link] 7.6 1.32 502 150 500 162 312 326 338 67.600
17 Kh. 17.9 1.32 539 190 532 217 347 370 315 45.65
Besbass
18 Serag 248.0 1.4 926 125 375 160 295 525 215 6.82

b- Caractéristiques des bassins versants :

On dénombre 17 bassins versants dans la région . Ils ont tous une pente assez forte et logés en
grande partie dans les chaînes montagneuses de Berda , Halfaya , Sif El Laham , Hachichina
et Batoum ( avec un altitude maximale variant entre 400 m et plus que 900 m ) , avec des lits
encaissés surtout dans la partie amont . Tous les cours d’eau de la région ruissellent du Nord
vers le Sud ( Djebel Berda ) et du Sud vers le Nord ( Djebel Sif Laham et Oum Ali ). L’ oued
Sareg qui joue le rôle d’un collecteur principal est situé dans le couloir séparant les deux
chaînes , donc il draine les cours prenant naissance sur les crêtes des chaînes citées ci-dessus.
Il coule de l’ Ouest vers l’ Est et se jette dans Sebkhat Sidi Mansour.
La région du Ségui est subdivisé en deux partie :
-Segui Ouest : La région située à l’ouest de la route de Kébili-Gafsa
-Segui Est : La région située à l’Est de la route de Kébili-Gafsa

- Les bassins versants de Chareb :

50
Le bassin versant de Chareb constitue la partie qui draine les flancs sud d’une série de
montagnes s’étendant de l’ouest vers l’est : Sefra, Zitouna, El Asker, Halfaya, Sif El Laham ,
Oum Ali , Hchichina, Battoum . Ils se jettent tous dans Chott El Fejjej . Le chareb constitue la
partie des terres fertiles située entre ces chaînes de montagnes et chott El Fejjej . Tous les
cours d’eau sont de direction Nord-Sud. Cette région est subdivisée en deux parties :
- à l’ Ouest de la route de Kébili-Gafsa : Le Cherb Ouest
- à l’ Est de la route de Kébili-Gafsa : Le Cherb Est

 Chareb Ouest :

Elle est située à l’Ouest de la route Kébili-Gafsa et s’étant jusqu’aux limites administratives
des gouvernorats de Kébili et Tozeur. Elle est bordée par la chaîne montagneuse de Djebel El
Asker au Nord et Chott El Fejjej au Sud et s’étale sur une longueur de 45 km de l’Est à l’
Ouest . Sa superficie totale est de 23000 ha . Elle draine les flancs Sud des djebels Meghatta ,
El Hamra , El Asker , Ezzitouna , Sefra et Taferma . La zone est classée dans l’étage
bioclimatique saharien supérieur à hiver tempéré . La pluviométrie moyenne annuelle dans la
région est de l’ordre de 120 mm . Les bassins versants présentent tous une pente assez forte à
forte. Ils sont logés en grande partie dans la chaîne montagneuse El Asker.

Tableau 2 : les caractéristiques physiques des bassins versants de la région de Chereb Ouest

N° Cours d’eau Sup kc Hmax Hmin H5(m) H95 H50 Hmoy D(m) Ig(m/km)
(km2) (m) (m) (m) (m) (m)
1 ElMeghatta 36.0 1.21 508 70 387 102 180 289 285 31.9
2 Sinaouni 11.0 1.35 539 75 325 95 150 307 212 34.2
3 Tine 16.0 1.33 540 75 400 90 180 307 310 42.5
4 Dhiebi 9.3 1.19 582 75 502 95 260 328 407 94.6
5 ElHamra Est 8.6 1.34 412 35 360 42 132 223 318 58.9
6 Elhamra Ouest 14.8 1.45 582 75 464 88 256 328 376 45.8
7 SidiAhmed 8.3 1.67 462 50 406 62 125 256 344 90.5
8 Mahaïa 8.6 1.53 608 75 500 112 352 341 388 57.1
9 El Asker 12.8 1.17 625 50 500 72 288 337 428 87.3
10 El Mahser 9.2 1.15 645 75 540 132 312 360 408 99.5
11 EZZitouna 22.1 1.43 645 50 490 92 260 346 398 42.8
12 Taferma 6.9 1.38 593 50 460 62 174 321 398 76.5
13 El Hamil 8.1 1.18 593 30 424 39 134 311 385 98.7

 Chareb Est :

La zone concernée est située entre la route de Kébili-Gafsa à l’Ouest et la limite

administrative des gouvernorats de Kébili et Gabès à l’ Est. C’est une plaine aux piedmonts
des montagnes de Hachichina, Naïmia , Oum El Gloub et Lefaïa au Nord et à l’Est et Le
Djebel Barrani au Sud. Sa superficie totale est de 11000 ha .Elle est limitée au Sud par Chott
El fejjej qui est l’exutoire naturel de tous les cours d’eau drainant ces montagnes. La
pluviométrie annuelle moyenne est de 120 mm. La zone de superficie 3000 ha et située entre
chott El Fejjej au Sud et Djebel Barrani au Nord est appelée Chareb barrani , alors que le
reste est appelé Chareb dakhlani de superficie 8624 ha . Les pentes des bassins versants
sont assez fortes et ils sont logés dans la chaîne montagneuse de Hachichina , Naïmia , Oum
El Gloub et Lefaïa . Les lits de ces oueds sont encaissés dans la partie amont. Elle est
constituée par 5 bassins versants. La zone est classée dans l’étage bioclimatique saharien
supérieur à hiver tempéré.

51
Tableau 3 : les caractéristiques physiques des bassins versants de la région de Chereb
Est

N° Cours d’eau Sup kc Hmax Hmin H5(m) H95 H50 Hmoy D(m) Ig(m/km)
(km2) (m) (m) (m) (m) (m)
1 Bazina 30.0 1.12 455 70 282 84 142 262 198 30.5
2 Hachichina 19.4 1.59 455 120 440 140 272 287 300 27.7
3 Oum ElGloub 24.5 1.13 384 80 250 100 140 232 150 25.0
4 Naïmia 10.5 1.64 521 122 474 134 270 321 340 41.0
5 El Kseub 25.8 1.43 596 100 474 122 254 348 352 33.8

- Les bassins versants d’ El Bahaier :

C ‘est une région située au nord de la chaîne de Tebaga , limitée à l’ est par le gouvernorat de
Gabès. Elle est formée par des petits bassins versants qui drainent le flanc nord de Tebaga .
Les cours d’eau se présentent sous forme d’un chevelu d’affluents plus ou moins parallèles et
qui confluent par deux ou trois aux piedmonts de la chaîne. Ils traversent la route Kébili-
Gabés dans le sens sud-nord. La région d’El Bahayer est une petite partie de la région d’El
Bahira ( sur la carte d’état major) qui s’étend sur tous les versants nord de la chaîne Tebaga
entre les deux gouvernorats Kébili et Gabès. Nous citons ( tableau 6 ) les principaux cours
d’eau intéressant le gouvernorat de Kébili (El Bahayer ).

Tableau 4 : les caractéristiques physiques des bassins versants de la région d’ El Bahayer

N° Cours d’eau Sup kc Hmax Hmin H5(m) H95 H50 Hmoy D(m) Ig(m/km)
(km2) (m) (m) (m) (m) (m)
1 Abdellatif 15.6 1.34 273 50 203 75 119 161 128 17.5
2 Zigzaou 2.8 1.34 225 100 183 110 136 162 73 23.5
3 Radhouane 6.0 1.37 226 125 213 131 164 175 82 17.4
4 Saïdane 7.2 1.15 300 100 247 108 136 200 139 40.9

- Les bassins versants de Tebaga et Dahar :

Les bassins versants de Tebaga sont situés au Nord de Kébili. Les oueds prennent leurs
sources dans les versants Sud de cette chaîne montagneuse. Ils constituent la zone de recharge
de la nappe du complexe terminal. Donc les travaux entrepris ou à entreprendre sur ces
bassins visent surtout cet objectif. En effet , cette nappe est la plus sollicitée du gouvernorat .
Le long de cette chaîne , il y a des affleurements de la roche aquifère de la nappe du complexe
terminal. En effet , les principaux cours d’eau qui drainent le ruissellement et contribuent à l’
alimentation de la nappe du complexe terminal descendent des crêtes de Tebaga et des
Matmata et le flanc Ouest du Dahar. On cite particuliéremet Oued Brimba à Mansoura , Oued
El Hachana à Fej errih , Oued El Melah –Nord-Kébili , Oued Oum Goullem à Bazma-Kébili
et Oued El Melah-Jemna( Dahar). L’exutoire naturel de ces cours d’eau est Chott de Kébili ,
prolongement de chott Jerid au sud de Kébili . Certains cours d’eau se déversent directement
dans chott Jerid , tel que l’oued Brimba, d’ autres se déversent dans la plaine à l’ Est de Kébili
et provoquent la coupure du trafic routier entre kebili , Douz et Jemna. Par conséquent , la
protection de la ville de Kébili , des villages voisins ainsi que de l’ infrastructure routière de la
région est aussi l’un des objectifs des travaux de CES. Les bassins versants sont souvent peu
accidentés . La plaine dite de Bougarfa est à la limite entre le gouvernorat de Gabès et de
Kébili, localisée entre kebili et Douz . Elle est jalonnée par Djebel Daouaya à l’ouest , Djebel
Kheneg Bougarfa au Nord et Djebel El Gueha à l’ est . Elle draine les bassins versants de

52
Oued El Guelta , Oued El Guesseïet et Oued Errmel. C’est une plaine agricole subdivisée en
terres collectives . Les bassins versants de El Guelta et El Guesseïet présentent une pente
modérée avec des lits encaissés surtout dans la partie amont. Les cours d’eau concernés
ruissellent du Nord vers le Sud-Est . Ils traversent la piste Bazma-Dahar pour s’épandre dans
la plaine Bougarfa avant de se déverser dans l’Oued Errmal, au sud ouest de Bir Negoua.
Oued Errmal , recevant les eaux de El Guelta , El Guesseïet , Chebib , Medaoua et plusieurs
autres affluents secondaires , travese la route Douz-Matmata puis s’ étend dans la dépression
dite behiret ben Aïssa ,avant la confluence avec Oued El Hallouf à l’ouest de Margueb Dhieb.
Certains autres cours d’eau , tels que Oued Dzou et Sedret El M’bet ne déversent pas dans ces
chotts , mais dans des plaines et les eaux peuvent être utilisées pour le développement des
parcours et la recharge des nappes phréatiques . Le principal cours d’eau des versants sud de
Tebaga est l’ Oued El Maleh . Il est constitué par 3 branches : oued El Maleh , Oued
Mansoura et Oued Kerker. L’Oued El Melah-Jemna ( Dahar ) quant à lui il est constitué par
deux affluents principaux : l’oued Dzou (appelé Oued Guelib El Anz dans sa partie amont ,
avec un affluent principal: oued El Baghla) et l’oued El Hassane (appelé Oued Nesrif dans sa
partie avale ,avec les affluents suivants : Ben Gharbia , El Asnam , El Majen , Guedah Triki ,
Mziraâ , Tolba , Haj Mhamed) . Le bassin versant de cet Oued est particulièrement accidenté
dans sa partie haute .De l’ autre coté de la chaîne de Tebaga , on a les bassins versants d’ El
Bahayer cités ci-dessus.
Le principal Oued de la région et de tout le gouvernorat de Kébili est El Hallouf. Il prend ses
sources à Beni Khdeche ( gouvernorat de Mednine ) est constitué de plusieurs affluents
principaux : Oued Zmertene au sud ouest de Toujéne , Oued Oum Lebbas qui prend ses
sources à l’ouest de Techine , Oued Leghene qui prend naissance à Techine et Matmata et
enfin Oued Ahmadi ( Oued Essedd en aval ) prenant naissance à Tamezret. Par conséquent
oued El Hallouf et ses affluents prennent leurs sources dans les gouvernorats de Gabès et
Mednine, et il joue le rôle d’un récepteur de plusieurs oueds qui drainent les monts ouest de
Djeffara de Tamezret jusqu’à Béni Khdeche.

Tableau 5 : les caractéristiques physiques des bassins versants de la région de Tebaga et dahar
N° Cours d’eau Sup kc Hmax Hmin H5(m) H95 H50 Hmoy D(m) Ig(m/km)
(km2) (m) (m) (m) (m) (m)
1 ElHachana 6.8 1.07 165 30 125 32 46 57 93 29
2 Bimba 20.7 1.30 155 25 69 23 57 90 46 5.4
3 Elmaleh ( 65.0 1.32 257 25 155 50 88 141 105 7.2
Teb)
4 Elhallouf 1800.0 1.15 713 155 495 142 300 434 352 6.67
5 Oum Goulam 31.5 1.84 302 47 217 67 136 147 150 9.1
6 Dzou 53.7 1.68 394 120 345 145 225 257 200 11.4
7 Sedret M’ bet 22.0 1.19 252 90 221 105 145 171 116 17.8
8 BenGharbia 13.7 1.28 327 100 316 121 153 213 195 31.0
9 Guelib El Anz 7.5 1.53 355 100 338 126 170 227 212 33.6
10 ElGhiraine 8.5 1.15 300 125 273 131 178 212 142 38.4
11 Thmed 9.5 1.36 350 125 284 134 183 237 150 25.4
12 ElMerbah 8.6 1.34 394 150 357 165 211 272 192 36.2
13 Med Snoussi 10.2 1.40 394 150 373 157 209 272 216 33.8
14 Recifa 6.8 1.93 375 140 358 158 215 257 200 24.7
15 ElHassane 27.0 1.62 429 140 410 177 291 284 233 18.3
16 Haj M’hamed 14.5 1.36 450 175 405 195 250 312 210 27.8
17 ElMajene 16.5 1.93 280 120 254 133 163 200 121 9.53
18 ElMaleh (dah) 427.5 1.17 429 70 284 80 145 249 204 7.08
19 ElGuelta 25.7 1.60 445 160 384 176 262 302 208 16.9
20 ElGuessïet 17.0 1.63 425 150 345 160 238 287 185 17.8

53
- Les bassins versants de la région de Bouflidja :
La région concernée ( formée des bassins versants ) a une superficie totale de 2294 km2 , avec
plus de 90% des affleurements durs , donc favorables au ruissellement . La région est
caractérisée par des crues violentes ( le débit de pointe de l’Oued Belkacheb a atteint 1940
m3/s le 5 Mars 1979 ) . Tous les oueds se jettent dans la dépression endoréique dite garaet
Bouflija ( Ksar Ghilaine ) , située à la limite du Grand erg oriental. La zone inondable a une
superficie de 13 km2 avec une oasis de 100 ha et une zone touristique importante. Les bassins
versants présentent une pente modérée et ils sont logés dans la chaîne montagneuse du flanc
Ouest du Dahar . Les lits des bassins versants sont souvent larges, ceci favorise l’ infiltration
dans les roches aquifères. Les cours d’eau ruissellent de l’ Est vers l’ Ouest . Ils traversent
rapidement le vaste plateau du Dahar pour se déverser enfin dans le bassin endoréique de
Garaet Bouflija. Les oueds drainent le flanc ouest de la série des Djebel compris entre Djebel
Bazina au Sud et Djebel El Fjeiej au Nord. Les principaux cours d’eau sont Belkacheb amont
confluence , oued Larej , Oued Belkacheb aval confluence ou Magran , Oued Rechada et
oued Mahabés. L’ Oued Belkacheb amont est issu des Djebels El Fjeiej et Djebel Adroum . Il
collecte un grand nombre d ‘affluents . l’ Oued Larej a un bassin adjacent à celui de l’oued
Belkacheb. Le bassin est peu compact . L’ oued prend ses sources dans les flancs des
montagnes situées entre le Djebel Adroum au Nord et le Djebel Flouassen . Il rejoint oued
BelKacheb à 1 km aval de la piste Trapsa . Il devient alors Oued Belkacheb aval ou Magran ,
et se jette dans Garaet Bouflija à environ 4 km de la confluence. Le devient large après la
confluence, et barré sur sa partie avale par les dunes de sables en mouvement. L’ oued perd
une bonne partie de ses apports avant d’arriver à l’exutoire. L’ oued Rechada prend naissance
dans les djebels limitant au nord le bassin d’ oued Belkacheb et au sud celui d’oued Mehabés
. Il jouit très peu des apports des zones montagneuses . Il a un bassin allongé et peu compact.
Enfin l’ oued Mehabés est issu des Djebels entre le Djebel Timasjal au Nord et Djebel Bazina
au Sud . Il se jette directement dans Garaet Bouflija . Le bassin versant est très peu compact et
allongé.

Tableau 6 : les caractéristiques physiques des bassins versants de Bouflidja

N° Cours d’eau Sup kc Hmax Hmin H5(m) H95 H50 Hmoy D(m) Ig(m/km)
(km2) (m) (m) (m) (m) (m)
1 Belkacheb 824 1.29 689 230 - - - 459 137 3.5
2 Larej 670 1.34 599 240 - - - 419 136 5.27
3 Magran 1540 1.36 689 215 - - - 402 137 3.5
4 Rechada 108 1.37 427 248 - - - 337 57 5.5
5 Mehabés 576 1.69 605 215 - - - 410 97 4.0

Les Apports moyens et fréquentiels :

- Les formules utilisées :

-Formules des apports moyens :

On a plusieurs formules pour estimer l’ apport moyen , les plus connues sont :
3
P
1-Formule de Tixeront : R 2 ; avec P : Pluviométrie moyenne annuelle en m ; E :
3E
évaporation annuelle en m et R : ruissellement en m .
En fonction des région , l ‘ évaporation E est de 1 m au Nord et peut atteindre 2.5 m au Sud.

54
 2
kp
2-Formule Algérienne : RSP(110 ) ; avec R : ruissellement en m , S :
superficie du bassin versant en km2 , P : pluie moyenne annuelle en m , k = a-0.01log(S) ,
avec a = paramètre régional ( a~0.185 )
3
3-Formule de Chaumont : R (P) ; avec R : ruissellement en m , P : pluie
2

moyenne annuelle en m ,  : coefficient d’écoulement ,  : facteur de rétention du sol .

Cette méthode n’est applicable qu’ au Nord et au Centre de la Tunisie

4-Formule de Turc1 : Le ruissellement R est donné indirectement par le déficit

1  0.9  1
d’écoulement (D) en utilisant 3 formules : D = P - R et 2 et
D P 
2 2

 300 25  0.05

2
,
avec R : lame ruisselée en mm , P : pluie moyenne annuelle en mm , D : moyenne annuelle
du déficit en mm et  : moyenne annuelle de la température en °C.
Si on D < P alors R = P- D sinon R = 0

5-Formule de Turc2 : Elle est utilisée surtout pour les pluies faibles , la formule s ‘écrit
5 1
3 Zm 3
sous la forme R( P ) ( ) , avec R : lame ruisselée en m , Zm : Altitude
30 100
moyenne en m et P : pluie moyenne annuelle en m.

6-Formule de Fersi : Elle a été établie sur des bassins versants du sud Tunisien . Elle s’écrit
sous la forme R KP Ig , avec P : pluie moyenne annuelle en mm , R : ruissellement en
mm , Ig : indice de pente global ( en m/km ) et K : paramètre égal à 0.01639 pour le Sud
Tunisien . Cependant l’ application de cette formule sur des bassins versants de l ‘ Oued
Leben ( Gouvernorat de Sidi Bouzid ) a montré que la valeur K = 0.0116 est plus adaptée à
cette région ( Mammou Halima , 1992 )

Beaucoup d’ autres formules sont établies par Ghorbel , Frigui , DGRE , …etc mais elle
sont surtout applicable au Nord ( Ichkeul , Medjerda ) et au Centre ( Zeroud , Marguellil ,
Nebhana) .

-Formules des apports Fréquentiels :

Toutes les formules que nous proposons ici , utilisent l’ apport moyen comme information
de base . Par conséquent ce dernier doit être bien estimé. On citera les deux principales
formules utilisées : Formule de Fersi et Formule de Ghorbel .

1-Formule de Fersi : R(T)0.21 RT

0.89
avec R : apport moyen en mm calculé par l’une
des formules ci-dessus , T : Période de retour en années , R(T) : ruissellement fréquentiel
pour la période de retour T .
Cette formule présente une anomalie : en effet si on fait le rapport entre R et R(5) on a
R(5)/R = 0.88 . Par conséquent l’ apport moyen annuel est supérieur à l’ apport de période de
retour 5 ans …. !!! . On se proposera d’améliorer cette formule ou ne l ‘utiliser qu’ au delà d’
une certaine période de retour…

55
 R(T)
2-Formule de Ghorbel : RTV  , le rapport RTV dépend de la région . On a alors :
R
Zone I ( Le Nord ) : RTV = 0.92 Log( T ) + 0.66
Zone II ( Le Centre ) : RTV = 1.27 Log( T ) + 0.45
Zone III ( Le Sud ) : RTV = 3.65 Log ( T ) – 3.60

- Calcul des lames ruisselées :

Notation : Ral : apport moyen par la formule algérienne
Rtu : apport moyen par la formule de Turc2
Rfer1 : apport moyen par la formule de Fersi ( k =0.01694)
Rfer2 : apport moyen par la formule de Fersi ( k =0.0116)

Tableau 7 : les lames ruisselées par différentes formules dans les bassins versants de Segui.

N° Cours d’eau Sup (km2) Hmoy(m) pluie(mm) Ig(m/km) Ral(mm) Rtu(mm) Rfer1(mm) Rfer2(mm)
1 Lettouchi 16,50 481,00 120 57,00 10,278 17,014 14,849 10,509
2 Enf. ali saad 9,50 486,00 120 58,66 6,125 17,073 15,064 10,661
3 Regaieg 5,80 342,00 120 42,10 3,853 15,186 12,761 9,032
4 El Gsir 13,50 568,00 120 61,16 8,517 17,984 15,381 10,886
5 Nir El Karma 5,20 381,00 120 48,79 3,477 15,743 13,738 9,723
6 Jedari 7,00 570,00 120 71,60 4,598 18,005 16,642 11,779
7 Houinit 4,30 443,00 120 64,87 2,907 16,554 15,841 11,211
8 El Blida 7,50 578,00 120 75,62 4,906 18,089 17,103 12,105
9 Retam 7,60 535,00 120 55,62 4,967 17,629 14,668 10,381
10 Khanguet Gherib 9,95 327,00 120 31,38 6,397 14,961 11,018 7,798

11 Oum Ali 4,60 295,00 120 49,34 3,098 14,456 13,815 9,778
12 El Araguib 11,60 283,00 120 35,30 7,388 14,257 11,686 8,270
13 Jerouala 5,00 332,00 120 73,93 3,351 15,036 16,911 11,969
14 Sania 2,70 332,00 120 96,80 1,875 15,036 19,351 13,695
15 Bouzaiene 4,40 332,00 120 79,74 2,971 15,036 17,563 12,430

16 Baten El Karma 7,60 326,00 120 67,60 4,967 14,945 16,171 11,445
17 Khanguet El Besbass 17,90 370,00 120 45,65 11,093 15,590 13,289 9,405
18 Serag 248,00 525,00 120 6,82 127,871 17,518 5,136 3,635

Tableau 8 : les lames ruisselées par différentes formules dans les bassins de Chareb ouest.

N° Cours d’eau S (km2) Hmoy(m) Pluie(mm) Ig(m/km) Ral(mm) Rtur(mm) Rfer1(mm) Rfer2(mm)
1 ElMeghatta 36 289 120 31,9 21,313 14,357 11,109 7,862
2 Sinaouni 11 307 120 34,2 7,029 14,649 11,502 8,141
3 Tine 16 307 120 42,5 9,986 14,649 12,822 9,075
4 Dhiebi 9,3 328 120 94,6 6,004 14,976 19,130 13,539
5 ElHamra Est 8,6 223 120 58,9 5,579 13,168 15,094 10,683

56
 6 Elhamra Ouest 14,8 328 120 45,8 9,283 14,976 13,310 9,420
7 SidiAhmed 8,3 256 120 90,5 5,396 13,788 18,710 13,242
8 Mahaïa 8,6 341 120 57,1 5,579 15,171 14,862 10,519
9 El Asker 12,8 337 120 87,3 8,102 15,112 18,377 13,006
10 El Mahser 9,2 360 120 99,5 5,944 15,448 19,619 13,885
11 EZZitouna 22,1 346 120 42,8 13,511 15,245 12,867 9,107
12 Taferma 6,9 321 120 76,5 4,536 14,869 17,202 12,175
13 El Hamil 8,1 311 120 98,7 5,274 14,713 19,540 13,829

Tableau 9 : les lames ruisselées par différentes formules dans les bassins versants de Chareb
est.

N° Cours d’eau S (km2) Hmoy(m) Pluie(mm) Ig(m/km) Ral(mm) Rtur(mm) Rfer1(mm) Rfer2(mm)
1 Bazina 30 262 120 30,5 17,978 13,895 10,862 7,688
2 Hachichina 19,4 287 120 27,7 11,960 14,324 10,351 7,326
3 Oum ElGloub 24,5 232 120 25 14,878 13,343 9,834 6,960
4 Naïmia 10,5 321 120 41 6,729 14,869 12,594 8,913
5 El Kseub 25,8 348 120 33,8 15,615 15,274 11,435 8,093

Tableau 10 : les lames ruisselées par différentes formules dans les bassins versants d’El
Bahayeir

N° Cours d’eau Sup(km2) Hmoy(m) pluie(mm) Ig(m/km) Ral(mm) Rtur(mm) Rfer1(mm) Rfer2(mm)

1 Abdellatif 15,6 161 100 17,5 5,648 8,718 6,856 4,853

2 Zigzaou 2,8 162 100 23,5 1,124 8,736 7,945 5,623
3 Radhouane 6 175 100 17,4 2,304 8,963 6,837 4,839

4 Saïdane 7,2 200 100 40,9 2,735 9,371 10,482 7,419

Tableau 11 : les lames ruisselées par différentes formules dans les bassins de Tebaga-dahar

N° Cours d’eau Sup(km2) Hmoy(m) pluie(mm) Ig(m/km) Ral(mm) Rtur(mm) Rfer1(mm) Rfer2(mm)

1 ElHachana 6,8 57 100 29 3,729 6,167 8,826 6,247

2 Bimba 20,7 90 100 5,4 10,591 7,181 3,809 2,696
3 Elmaleh ( Tebaga) 65 141 100 7,2 30,802 8,341 4,398 3,113
4 Elhallouf 1800 434 100 6,67 655,582 12,133 4,233 2,996
5 Oum Goulam 31,5 147 100 9,1 15,680 8,457 4,944 3,499
6 Dzou 53,7 257 100 11,4 25,785 10,188 5,534 3,917
7 Sedret El M’ bet 22 171 100 17,8 11,212 8,895 6,915 4,894
8 BenGharbia 13,7 213 100 31 7,196 9,570 9,126 6,459

57
 9 Guelib El Anz 7,5 227 100 33,6 4,088 9,776 9,501 6,724
10 ElGhiraine 8,5 212 100 38,4 4,598 9,555 10,157 7,188
11 Thmed 9,5 237 100 25,4 5,104 9,917 8,260 5,846
12 ElMerbah 8,6 272 100 36,2 4,649 10,383 9,861 6,979
13 Med Snoussi 10,2 272 100 33,8 5,457 10,383 9,529 6,744
14 Recifa 6,8 257 100 24,7 3,729 10,188 8,146 5,765
15 ElHassane 27 284 100 18,3 13,577 10,533 7,011 4,962
16 El Haj M’hamed 14,5 312 100 27,8 7,589 10,869 8,642 6,116
17 ElMajene 16,5 200 100 9,53 8,565 9,371 5,060 3,581
18 ElMaleh (dahar) 427,5 249 100 7,08 175,999 10,082 4,361 3,087

19 ElGuelta 25,7 302 100 16,9 12,965 10,751 6,738 4,769

20 ElGuessïet 17 287 100 17,8 8,808 10,570 6,915 4,894

Tableau 12 : les lames ruisselées par différentes formules dans les bassins versants de Bouflija

N° Cours d’eau S (km2) Hmy(m) pl(mm) Ig(m/km) Ral(mm) Rtur(mm) Rfer1(mm) Rfer2(mm)

1 Belkacheb 824 459 50 3,5 27,943 3,894 1,533 1,085

2 Larej 670 419 50 5,27 23,119 3,777 1,881 1,331
3 Magran 1540 402 50 3,5 49,453 3,725 1,533 1,085
4 Rechada 108 337 50 5,5 4,294 3,513 1,922 1,360
5 Mehabés 576 410 50 4 20,126 3,750 1,639 1,160

- Calcul des coefficients de ruissellement et des volumes moyens :

La formule algérienne donne des estimations aberrantes pour les grands bassins versants.
Celle de Fersi 2 donne les estimations les plus faibles , alors que la formule de Turc surestime
légèrement , pour les petits bassins versants , la lame ruisselée par rapports aux valeurs
données par Fersi1 et Fersi2. Pour les grands bassins versants , la formule de Turc donne
même plus que le double de celles de Fersi . En examinant de prés cette formule , on
remarque que les paramètres qui interviennent sont la pluie et l ‘ altitude moyenne. Elle
ressemble à celle de Fersi qui fait intervenir la pluie et la pente ( expression indirecte de la
pente ) . Fersi a validé cette formule sur des bassins versants du Sud . On préfère donc cette
formule avec k =0.01694 . Par la suite on calculera les coefficients de ruissellements et les
volumes moyens écoulés en m3.

Tableau 13 : les coefficients de ruissellement et les apports moyens des bassins du Segui

Cours d’eau Sup (km2) pluie (mm) Rfer1(mm) Cr V(m3)

Lettouchi 16,50 120 14,849 0,12 245008,74
Enfidhet ali saad 9,50 120 15,064 0,13 143105,01
Regaieg 5,80 120 12,761 0,11 74016,62
El Gsir 13,50 120 15,381 0,13 207647,98
Nir El Karma 5,20 120 13,738 0,11 71437,94

58
Jedari 7,00 120 16,642 0,14 116497,00
Houinit 4,30 120 15,841 0,13 68116,23
El Blida 7,50 120 17,103 0,14 128274,34
Retam 7,60 120 14,668 0,12 111478,03
Khanguet Gherib 9,95 120 11,018 0,09 109625,04

Oum Ali 4,60 120 13,815 0,12 63550,30

El Araguib 11,60 120 11,686 0,10 135551,88
Jerouala 5,00 120 16,911 0,14 84555,25
Sania 2,70 120 19,351 0,16 52247,03
Bouzaiene 4,40 120 17,563 0,15 77277,13
Baten El Karma 7,60 120 16,171 0,13 122898,66

Khanguet El Besbass 17,90 120 13,289 0,11 237866,69

Serag 248,00 120 5,136 0,04 1273808,28

Soit un Apport moyen annuel de Vségui = 3322962.14 m3

Tableau 14 : les coefficients de ruissellement et les apports moyens des bassins du Cherb
Ouest

Cours d’eau Sup (km2) Pluie (mm) Rfer1(mm) cr V(m3)

ElMeghatta 36 120 11,109 0,09 399906,11
Sinaouni 11 120 11,502 0,10 126521,97
Tine 16 120 12,822 0,11 205151,48
Dhiebi 9,3 120 19,130 0,16 177905,23
ElHamra Est 8,6 120 15,094 0,13 129812,44
Elhamra Ouest 14,8 120 13,310 0,11 196994,75
SidiAhmed 8,3 120 18,710 0,16 155296,83
Mahaïa 8,6 120 14,862 0,12 127813,49
El Asker 12,8 120 18,377 0,15 235221,65
El Mahser 9,2 120 19,619 0,16 180492,67
EZZitouna 22,1 120 12,867 0,11 284363,83
Taferma 6,9 120 17,202 0,14 118697,07
El Hamil 8,1 120 19,540 0,16 158271,89

Soit un Apport moyen annuel de Vchareb-ouest = 2496449.41 m3

Tableau 15 : les coefficients de ruissellement et les apports moyens des bassins du Cherb est

Cours d’eau Sup (km2) Pluie(mm) Rfer1(mm) cr V(m3)

Bazina 30 120 10,862 0,09 325860,24
Hachichina 19,4 120 10,351 0,09 200817,62

59
Oum ElGloub 24,5 120 9,834 0,08 240933,00
Naïmia 10,5 120 12,594 0,10 132233,48
El Kseub 25,8 120 11,435 0,10 295011,03
Soit un Apport moyen annuel de Vchareb-est = 1194855.37 m3

Tableau 16 : les coefficients de ruissellement et les apports moyens des bassins du Tebaga-
Dahar

Cours d’eau Sup (km2) pluie(mm) Rfer1(mm) cr V(m3)

ElHachana 6,8 100 8,826 0,09 60018,74

Bimba 20,7 100 3,809 0,04 78839,92
Elmaleh ( Tebaga) 65 100 4,398 0,04 285863,40
Elhallouf 1800 100 4,233 0,04 7619287,75
Oum Goulam 31,5 100 4,944 0,05 155743,60
Dzou 53,7 100 5,534 0,06 297170,48
Sedret El M’ bet 22 100 6,915 0,07 152128,87
BenGharbia 13,7 100 9,126 0,09 125020,25
Guelib El Anz 7,5 100 9,501 0,10 71254,10
ElGhiraine 8,5 100 10,157 0,10 86330,35
Thmed 9,5 100 8,260 0,08 78472,85
ElMerbah 8,6 100 9,861 0,10 84807,00
Med Snoussi 10,2 100 9,529 0,10 97193,56

Recifa 6,8 100 8,146 0,08 55390,63

ElHassane 27 100 7,011 0,07 189307,69
El Haj M’hamed 14,5 100 8,642 0,09 125305,27
ElMajene 16,5 100 5,060 0,05 83485,17
ElMaleh (dahar) 427,5 100 4,361 0,04 1864368,26
ElGuelta 25,7 100 6,738 0,07 173163,13
ElGuessïet 17 100 6,915 0,07 117554,13

Soit un Apport moyen annuel de Vtebaga-dahar = 11800705.16 m3

Tableau 17: les coefficients de ruissellement et les apports moyens des bassins du El Bahayeir

Cours d’eau Sup (km2) pluie(mm) Rfer1(mm) cr V(m3)

Abdellatif 15,6 100 6,856 0,07 106960,29

Zigzaou 2,8 100 7,945 0,08 22246,97
Radhouane 6 100 6,837 0,07 41020,87
Saïdane 7,2 100 10,482 0,10 75469,78

60
Soit un Apport moyen annuel de Velbahayier= 245697.91 m3

Tableau 18: les coefficients de ruissellement et les apports moyens des bassins de Bouflija

Cours d’eau Sup (km2) pluie(mm) Rfer1(mm) cr V(m3)

Belkacheb 824 50 1,533 0,03 1263310,75

Larej 670 50 1,881 0,04 1260460,00
Magran 1540 50 1,533 0,03 2361041,94
Rechada 108 50 1,922 0,04 207564,97
Mehabés 576 50 1,639 0,03 944064,00

Soit un Apport moyen annuel de Vbouflija = 6036441.65 m3

- Calcul des volumes fréquentiels :

On va utiliser deux formules : La formule de Fersi et celle de Ghorbel pour la zone sud.
On estimera les apports fréquentiels pour les périodes de retour suivantes : 5 , 10 , 20 , 30 , 50
,100 et 1000 ans. Pour cela on donnera les rapports entre l’ apport moyen et l’ apport
fréquentiel pour chaque pas de temps et pour chaque formule .

5 10 20 30 50 100 1000
Fersi 0.87 1.63 3.02 4.33 6.83 12.65 98.22
Ghorbel 2.27 4.80 7.33 8.81 10.67 13.20 21.61

Si la formule de Fersi a donné une relative satisfaction dans l’estimation des apports moyens ,
elle présente cependant deux anomalies majeures dans l’estimation des apports fréquentiels .
-L’apport fréquentiel de période de retour 5 ans est inférieur à l’apport moyen ( rapport < 1) :
ceci est une aberration ….. !!!.
-L’apport fréquentiel de période de retour 1000 ans est presque égal à 100 fois l’apport
moyen : ceci semble être une forte surestimation
On va par conséquent utiliser la formule de Ghorbel.

Tableau 19: les apports moyen en m3 et fréquentiels en millions de m3 de Segui

N° Cours d’eau V(m3) V5 V10 V20 V30 V50 V100 V1000

1 Lettouchi 245008,74 0,556170 1,176042 1,893918 2,158527 2,614243 3,234115 5,294639
2 [Link] 143105,01 0,324848 0,686904 1,106202 1,260755 1,526930 1,888986 3,092499
3 Regaieg 74016,62 0,168018 0,355280 0,572148 0,652086 0,789757 0,977019 1,599499
4 El Gsir 207647,98 0,471361 0,996710 1,605119 1,829379 2,215604 2,740953 4,487273
5 Nir Karma 71437,94 0,162164 0,342902 0,552215 0,629368 0,762243 0,942981 1,543774
6 Jedari 116497 0,264448 0,559186 0,900522 1,026339 1,243023 1,537760 2,517500
7 Houinit 68116,23 0,154624 0,326958 0,526538 0,600104 0,726800 0,899134 1,471992
8 El Blida 128274,34 0,291183 0,615717 0,991561 1,130097 1,368687 1,693221 2,772008

61
 9 Retam 111478,03 0,253055 0,535095 0,861725 0,982121 1,189471 1,471510 2,409040
10 Khanguet Gherib 109625,04 0,248849 0,526200 0,847402 0,965797 1,169699 1,447051 2,368997
11 Oum Ali 63550,3 0,144259 0,305041 0,491244 0,559878 0,678082 0,838864 1,373322
12 El Araguib 135551,88 0,307703 0,650649 1,047816 1,194212 1,446339 1,789285 2,929276
13 Jerouala 84555,25 0,191940 0,405865 0,653612 0,744932 0,902205 1,116129 1,827239
14 Sania 52247,03 0,118601 0,250786 0,403870 0,460296 0,557476 0,689661 1,129058
15 Bouzaiene 77277,13 0,175419 0,370930 0,597352 0,680812 0,824547 1,020058 1,669959
16 Baten El Karma 122898,66 0,278980 0,589914 0,950007 1,082737 1,311329 1,622262 2,655840
17 Khanguet El Besbass 237866,69 0,539957 1,141760 1,838710 2,095606 2,538038 3,139840 5,140299
18 Serag 1273808,28 2,891545 6,114280 9,846538 11,222251 13,591534 16,814269 27,526997

Tableau 20: les apports moyen en m3 et fréquentiels en millions de m3 de Chareb Ouest

N° Cours d’eau V(m3) V5 V10 V20 V30 V50 V100 V1000

1 ElMeghatta 399906,11 0,907787 1,919549 2,931312 3,523173 4,266998 5,278761 8,641971

2 Sinaouni 126521,97 0,287205 0,607305 0,927406 1,114659 1,349989 1,670090 2,734140
3 Tine 205151,48 0,465694 0,984727 1,503760 1,807385 2,188966 2,708000 4,433323
4 Dhiebi 177905,23 0,403845 0,853945 1,304045 1,567345 1,898249 2,348349 3,844532
5 ElHamra Est 129812,44 0,294674 0,623100 0,951525 1,143648 1,385099 1,713524 2,805247
6 Elhamra Ouest 196994,75 0,447178 0,945575 1,443972 1,735524 2,101934 2,600331 4,257057
7 SidiAhmed 155296,83 0,352524 0,745425 1,138326 1,368165 1,657017 2,049918 3,355964
8 Mahaïa 127813,49 0,290137 0,613505 0,936873 1,126037 1,363770 1,687138 2,762050
9 El Asker 235221,65 0,533953 1,129064 1,724175 2,072303 2,509815 3,104926 5,083140
10 El Mahser 180492,67 0,409718 0,866365 1,323011 1,590140 1,925857 2,382503 3,900447
11 EZZitouna 284363,83 0,645506 1,364946 2,084387 2,505245 3,034162 3,753603 6,145102
12 Taferma 118697,07 0,269442 0,569746 0,870050 1,045721 1,266498 1,566801 2,565044
13 El Hamil 158271,89 0,359277 0,759705 1,160133 1,394375 1,688761 2,089189 3,420256

Tableau 21 : les apports moyen en m3 et fréquentiels en millions de m3 de Chareb Est

N° Cours d’eau V(m3) V5 V10 V20 V30 V50 V100 V1000

1 Bazina 325860,24 0,739703 1,564129 2,388556 2,870829 3,476929 4,301355 6,911496

2 Hachichina 200817,62 0,455856 0,963925 1,471993 1,769203 2,142724 2,650793 4,259342

3 Oum ElGloub 240933 0,546918 1,156478 1,766039 2,122620 2,570755 3,180316 5,110189
4 Naïmia 132233,48 0,300170 0,634721 0,969271 1,164977 1,410931 1,745482 2,804672
5 El Kseub 295011,03 0,669675 1,416053 2,162431 2,599047 3,147768 3,894146 6,257184

Tableau 22 : les apports moyen en m3 et fréquentiels en millions de m3 de Tebaga-Dahar

Cours d’eau V(m3) V5 V10 V20 V30 V50 V100 V1000

N°
1 ElHachana 60018,74 0,136243 0,288090 0,439937 0,528765 0,640400 0,792247 1,297005

62
 2 Bimba 78839,92 0,178967 0,378432 0,577897 0,694580 0,841222 1,040687 1,703731
3 Elmaleh ( Tebaga) 285863,4 0,648910 1,372144 2,095379 2,518457 3,050162 3,773397 6,177508

4 Elhallouf 7619287,8 17,295783 36,572581 55,849379 67,125925 81,297800 100,574598 164,652808

5 Oum Goulam 155743,6 0,353538 0,747569 1,141601 1,372101 1,661784 2,055816 3,365619
6 Dzou 297170,48 0,674577 1,426418 2,178260 2,618072 3,170809 3,922650 6,421854
7 Sedret El M’ bet 152128,87 0,345333 0,730219 1,115105 1,340255 1,623215 2,008101 3,287505

8 BenGharbia 125020,25 0,283796 0,600097 0,916398 1,101428 1,333966 1,650267 2,701688

9 Guelib El Anz 71254,1 0,161747 0,342020 0,522293 0,627749 0,760281 0,940554 1,539801
10 ElGhiraine 86330,35 0,195970 0,414386 0,632801 0,760570 0,921145 1,139561 1,865599
11 Thmed 78472,85 0,178133 0,376670 0,575206 0,691346 0,837305 1,035842 1,695798
12 ElMerbah 84807 0,192512 0,407074 0,621635 0,747150 0,904891 1,119452 1,832679
13 Med Snoussi 97193,56 0,220629 0,466529 0,712429 0,856275 1,037055 1,282955 2,100353
14 Recifa 55390,63 0,125737 0,265875 0,406013 0,487991 0,591018 0,731156 1,196992
15 ElHassane 189307,69 0,429728 0,908677 1,387625 1,667801 2,019913 2,498862 4,090939
16 El Haj M’hamed 125305,27 0,284443 0,601465 0,918488 1,103939 1,337007 1,654030 2,707847
17 ElMajene 83485,17 0,189511 0,400729 0,611946 0,735504 0,890787 1,102004 1,804115
18 ElMaleh (dahar) 1864368,3 4,232116 8,948968 13,665819 16,425084 19,892809 24,609661 40,288998

19 ElGuelta 173163,13 0,393080 0,831183 1,269286 1,525567 1,847651 2,285753 3,742055

20 ElGuessïet 117554,13 0,266848 0,564260 0,861672 1,035652 1,254303 1,551715 2,540345

Tableau 23 : les apports moyen en m3 et fréquentiels en millions de m3 de El Bahayeir

N° Cours d’eau V V5 V10 V20 V30 V50 V100 V1000

1 Abdellatif 106960,29 0,242800 0,513409 0,784019 0,942320 1,141266 1,411876 2,311412

2 Zigzaou 22246,97 0,050501 0,106785 0,163070 0,195996 0,237375 0,293660 0,480757
3 Radhouane 41020,87 0,093117 0,196900 0,300683 0,361394 0,437693 0,541475 0,886461
4 Saïdane 75469,78 0,171316 0,362255 0,553193 0,664889 0,805263 0,996201 1,630902

Tableau 24 : les apports moyen en m3 et fréquentiels en millions de m3 de Bouflija

N° Cours d’eau V V5 V10 V20 V30 V50 V100 V1000

1 Belkacheb 1263310,8 2,867715 6,063892 9,260068 11,129768 13,479526 16,675702 27,300145
2 Larej 1260460 2,861244 6,050208 9,239172 11,104653 13,449108 16,638072 27,238541
3 Magran 2361041,9 5,359565 11,333001 17,306437 20,800779 25,192317 31,165754 51,022116
4 Rechada 207564,97 0,471172 0,996312 1,521451 1,828647 2,214718 2,739858 4,485479
5 Mehabés 944064 2,143025 4,531507 6,919989 8,317204 10,073163 12,461645 20,401223

Les débits fréquentiels :

63
 - Les formules utilisées :

S i(T)C
1-Formule rationnelle : Q(T) avec Q( T) : débit de pointe en m3/s pour
3.6
une période de retour T en années , C : Coefficient de ruissellement , i(T) : intensité en mm/h
de l’averse maximale estimée à partir des courbes I-D-F pour une période de retour T et une
durée D égale au temps de concentration du bassin versant , S : Superficie du bassin versant
en km2.

ip
2-Formule de Burkli-Ziegler : Q(T)SCi(T)
S , avec Q( T) : débit de pointe en
m3/s pour une période de retour T en années , C : Coefficient de ruissellement , i(T) :
intensité en mm/h de l’averse maximale estimée à partir des courbes I-D-F pour une période
de retour T et une durée D égale au temps de concentration du bassin versant , S : Superficie
du bassin versant en km2 , ip : pente du bassin versant en %.

Q(T)
3-Formule de Ghorbel : RTQ , avec RTQ : un rapport régional dont la
(Qmax)moy
formulation est la suivante :
Zone I( Nord ) : RTQ = 1.33 Log( T ) + 0.46
Zone II ( Centre ) : RTQ = 1.074 T0.4 – 0.71
Zone III ( Sud ) : RTQ = 1.47 T0.4 – 1.35

Q(T) : est le débit de pointe en m3/s de période de retour T en années , T : Période de retour
en années .
(Qmax)moy : est le débit maximum moyen annuel : c’est la moyenne des débits de crues
observés . Il peut être donné par les formules empiriques . Pour le Sud cette formule est
(Qmax)moy = 85 Log(S) , S : Superficie du bassin versant en km2.

3
4-Formule de Fersi : Q(T)2.4 10 [Link].T .S , avec S : Superficie du bassin versant
0.89

en km2. , Q(T) : débit de pointe en m3/s pour une période de retour T en années , P : Pluie
moyenne annuelle , Ig : indice de pente global en m/km , T : période de retour en années

0.41
5-Formule de Kallel : Q(T)qr T S , avec S : Superficie du bassin versant en km2, Q(T) :
S
débit de pointe en m3/s pour une période de retour T en années , T : Période de retour en
années . qr est un paramètre régional , il est égal à 5.5 pour le Nord et 12.5 pour le Sud .

- Calcul des débits fréquentiels :

On se propose d’utiliser les 3 formules établies en Tunisie : Ghorbel , Fersi , Kallel pour
quelques cours d’eau . On comparera les résultats entre eux et avec ceux calculés par
[Link] ( 1998 ). On prend de chaque région 3 cours d’eau.

Tableau 25 : Tableau de comparaison entre les 3 formules

fersi ghorbel kallel

oued Q5 q10 q20 q5 q10 q20 q5 q10 q20

64
segui Lettouchi 113,458 210,257 389,644 149,886 242,413 364,504 98,227 130,512 173,410
Enfidhet ali saad 67,227 124,583 230,874 120,369 194,675 292,722 74,533 99,031 131,581
Khanguet El 98,576 182,678 338,535 154,241 249,456 375,093 102,309 135,937 180,617
Besbass
Serag 204,038 378,119 700,721 294,785 476,760 716,877 380,815 505,983 672,292
chareb ElMeghatta 138,538 256,735 475,776 191,599 309,876 465,943 145,090 192,780 256,144

Sinaouni 45,383 84,103 155,858 128,207 207,352 311,784 80,202 106,563 141,589
Tine 82,032 152,020 281,720 148,241 239,753 360,503 96,727 128,520 170,762
Dhiebi 106,133 196,683 364,488 119,231 192,835 289,955 73,744 97,983 130,189
Hachichina 64,827 120,136 222,634 158,543 256,414 385,556 106,510 141,518 188,033

Brimba 11,237 20,825 38,592 162,011 262,023 393,990 110,020 146,182 194,231
Elmaleh ( Tebaga) 47,048 87,188 161,575 223,191 360,970 542,770 194,959 259,040 344,183

Elhallouf 1206,960 2236,709 4145,017 400,762 648,158 974,599 1025,945 1363,158 1811,209
ElMaleh (dahar) 304,273 563,872 1044,954 323,899 523,846 787,679 499,984 664,321 882,674
ElGuelta 43,663 80,915 149,950 173,579 280,732 422,121 122,590 162,883 216,421
bahayeir Abdellatif 27,445 50,860 94,252 146,887 237,563 357,211 95,510 126,903 168,614
Zigzaou 6,615 12,258 22,717 55,050 89,034 133,875 40,464 53,764 71,435
Bouflij Larej 177,480 328,902 609,514 347,922 562,700 846,101 625,929 831,663 1105,019
Magran 270,928 502,077 930,436 392,421 634,668 954,315 948,960 1260,870 1675,300
Rechada 29,857 55,331 102,538 250,338 404,875 608,789 251,304 333,904 443,654

Si la formule de Fersi s’applique pour les ruissellements dans le Sud , elle donne ici des
valeurs trop élevées . Par contre la méthode de Kallel souestime ,mais légèrement , les débits .
On a pris comme référence , les valeurs calculées par [Link] , sur ces cours d’eau (
[Link] , 1998). Cependant , la formule de Kallel , donne des valeurs plutôt supérieures à
celles calculées par [Link] , pour les grands bassins versants ( Hallouf, Larej , Magran
,…etc ). Par conséquent , le paramètre Surface , prend une importance capitale ici . La
formule de Kallel semble être très sensible à ce paramètre . D’ailleurs , il intervient avec la
racine carrée . Or , comme on l’ a constaté lors des visites sur terrains , ce sont surtout les
petits cours d’eau qui causent le plus de problèmes d’érosion . Ces oueds sont d’ une pente
forte et descendent des montagnes . Par conséquent , on ne tolère pas cette souestimation. Les
grands cours d’eau revêtent aussi une importance primordiale surtout pour la protection
contre les inondations des villes et des infrastructures hotelliéres ( Ksar Ghilaine par exemple
) . Il faut signaler que [Link] a utilisé la formule de Ghorbel qu’ on a utilisé aussi , mais il
n’ a pas utilisé les mêmes rapports . En effet , pour le Sud les rapports en fonction des
périodes de retour sont , avec
Zone III ( Sud ) : RTQ = 1.47 T0.4 – 1.35 :

T 5 10 20 50
RTQ 1.44 2.34 3.52 5.67
RTQ(Touil) 1.0 2.2 3.7 6.7

On a trouvé des valeurs assez proches de celles estimées par [Link] .

65
La formule de Fersi , semble tenir compte de la pluviométrie , la pente , la superficie et de la
période de retour. Mais , il semble que l’ exposant 0.89 de la période de retour est exagérée.
En plus , la formule de Fersi fait intervenir l’ indice de pente globale Ig . Ce dernier est très
fort pour les petits cours d’eau dans les montagnes et il est faible pour les grands cours d’eau
dans les plaines. Lorsque la pente est trop forte ou la surface est très importante , la formule
de Fersi donne des valeurs trop exagérées . Donc il faut intervenir au niveau de l’exposant , de
la pente et de la superficie. On réduit l’exposant , et on cherche à assurer une meilleure
compensation entre la superficie et la pente. Déjà une première remarque s’impose , le
coefficient 0.0024 semble être élevé( dans la littérature) , en prenant 0.00024 on divise toutes
les valeurs par 10 et on aboutit à des valeurs beaucoup plus raisonnables.. pour les présenter
dans ce rapport.
On utilisera la méthode de Ghorbel ( adoptée par [Link]) et on calculera les débits
fréquentiels pour tous les bassins versants. On ne peut pas utiliser la formule de Fersi
améliorée , parce qu’elle donne des valeurs faibles pour certains bassins et forts pour
d’autres…De même , la formule de Kallel est très sensible à la superficie du bassin versant
et elle tend à surestimer les débits pour les grands cours d’eau.

Tableau 26 Les débits fréquentiels ( Méthode de Ghorbel ) pour la région du Ségui

N° Cours d’eau Sup (km2) qmax Q5 Q10 Q20 Q30 Q50 Q100 Q1000
1,44 2,34 3,52 4,38 5,68 7,92 21,95
1 Lettouchi 16,50 103,486 149,020 242,158 364,271 453,269 587,801 819,610 2271,521
2 Enf. ali saad 9,50 83,107 119,673 194,469 292,535 364,006 472,045 658,204 1824,188
3 Regaieg 5,80 64,891 93,444 151,846 228,418 284,224 368,583 513,940 1424,366
4 El Gsir 13,50 96,078 138,353 224,823 338,196 420,823 545,725 760,941 2108,920
5 Nir El Karma 5,20 60,860 87,639 142,413 214,228 266,568 345,686 482,013 1335,883
6 Jedari 7,00 71,833 103,440 168,090 252,853 314,630 408,013 568,920 1576,742
7 Houinit 4,30 53,845 77,537 125,997 189,534 235,840 305,839 426,451 1181,894
8 El Blida 7,50 74,380 107,107 174,050 261,818 325,785 422,480 589,091 1632,646
9 Retam 7,60 74,869 107,812 175,194 263,539 327,927 425,257 592,964 1643,378
10 Khanguet Gherib 9,95 84,815 122,134 198,467 298,549 371,490 481,749 671,734 1861,688
11 Oum Ali 4,60 56,334 81,122 131,823 198,297 246,745 319,979 446,169 1236,540
12 El Araguib 11,60 90,479 130,290 211,721 318,486 396,298 513,920 716,593 1986,012
13 Jerouala 5,00 59,412 85,554 139,025 209,132 260,227 337,463 470,547 1304,103
14 Sania 2,70 36,666 52,799 85,798 129,064 160,597 208,262 290,394 804,817
15 Bouzaiene 4,40 54,693 78,759 127,983 192,521 239,557 310,659 433,172 1200,522
16 Baten El Karma 7,60 74,869 107,812 175,194 263,539 327,927 425,257 592,964 1643,378
17 Khanguet El Besbass 17,90 106,493 153,349 249,192 374,854 466,437 604,877 843,421 2337,511
18 Serag 248,00 203,528 293,081 476,256 716,420 891,454 1156,041 1611,945 4467,448

Tableau 27 :Les débits fréquentiels ( Méthode de Ghorbel ) pour la région du Chareb Ouest

N° Cours d’eau Sup (km2) qmax Q5 Q10 Q20 Q30 Q50 Q100 Q1000
1,44 2,34 3,52 4,38 5,68 7,92 21,95
1 ElMeghatta 36 132,286 190,491 309,549 465,646 579,411 751,383 1047,703 2903,671

66
 2 Sinaouni 11 88,518 127,466 207,133 311,585 387,710 502,784 701,066 1942,978
3 Tine 16 102,350 147,384 239,499 360,273 448,294 581,349 810,614 2246,587
4 Dhiebi 9,3 82,321 118,542 192,631 289,770 360,566 467,584 651,983 1806,947
5 ElHamra Est 8,6 79,432 114,383 185,872 279,602 347,914 451,176 629,104 1743,540

6 Elhamra Ouest 14,8 99,472 143,240 232,765 350,142 435,688 565,002 787,820 2183,416
7 SidiAhmed 8,3 78,122 112,495 182,805 274,988 342,173 443,731 618,723 1714,770
8 Mahaïa 8,6 79,432 114,383 185,872 279,602 347,914 451,176 629,104 1743,540
9 El Asker 12,8 94,113 135,523 220,224 331,277 412,214 534,561 745,374 2065,777
10 El Mahser 9,2 81,922 117,968 191,697 288,365 358,818 465,317 648,822 1798,187
11 EZZitouna 22,1 114,273 164,554 267,400 402,242 500,517 649,073 905,045 2508,300
12 Taferma 6,9 71,302 102,675 166,847 250,984 312,304 404,996 564,713 1565,083
13 El Hamil 8,1 77,221 111,199 180,698 271,819 338,229 438,617 611,592 1695,006

Tableau 28 :Les débits fréquentiels ( Méthode de Ghorbel ) pour la région du Chareb Est

N° Cours d’eau Sup (km2) qmax Q5 Q10 Q20 Q30 Q50 Q100 Q1000
1,44 2,34 3,52 4,38 5,68 7,92 21,95
1 Bazina 30 125,555 180,800 293,799 441,955 549,932 713,154 994,398 2755,939
2 Hachichina 19,4 109,463 157,627 256,144 385,310 479,449 621,751 866,948 2402,716
3 Oum ElGloub 24,5 118,079 170,034 276,305 415,638 517,187 670,689 935,187 2591,837
4 Naïmia 10,5 86,801 124,994 203,115 305,540 380,189 493,030 687,465 1905,284
5 El Kseub 25,8 119,988 172,782 280,771 422,357 525,546 681,530 950,302 2633,729

Tableau 29 :Les débits fréquentiels ( Méthode de Ghorbel ) pour la région du Tebaga-Dahar

N° Cours d’eau Sup (km2) qmax Q5 Q10 Q20 Q30 Q50 Q100 Q1000
1,44 2,34 3,52 4,38 5,68 7,92 21,95
1 ElHachana 6,8 70,763 101,899 165,586 249,087 309,943 401,935 560,445 1553,254
2 Bimba 20,7 111,857 161,075 261,747 393,738 489,936 635,350 885,911 2455,272
3 Elmaleh ( Tebaga) 65 154,098 221,901 360,588 542,424 674,948 875,275 1220,453 3382,443

4 Elhallouf 1800 276,698 398,445 647,474 973,978 1211,938 1571,646 2191,449 6073,525
5 Oum Goulam 31,5 127,356 183,393 298,014 448,295 557,821 723,384 1008,663 2795,473
6 Dzou 53,7 147,048 211,749 344,092 517,608 644,069 835,232 1164,619 3227,700
7 Sedret El M’ bet 22 114,106 164,313 267,008 401,653 499,784 648,122 903,719 2504,625
8 BenGharbia 13,7 96,621 139,135 226,094 340,107 423,201 548,809 765,240 2120,836

9 Guelib El Anz 7,5 74,380 107,107 174,050 261,818 325,785 422,480 589,091 1632,646

10 ElGhiraine 8,5 79,001 113,761 184,861 278,082 346,023 448,723 625,685 1734,063
11 Thmed 9,5 83,107 119,673 194,469 292,535 364,006 472,045 658,204 1824,188

67
12 ElMerbah 8,6 79,432 114,383 185,872 279,602 347,914 451,176 629,104 1743,540
13 Med Snoussi 10,2 85,731 123,453 200,611 301,773 375,502 486,952 678,990 1881,796

14 Recifa 6,8 70,763 101,899 165,586 249,087 309,943 401,935 560,445 1553,254
15 ElHassane 27 121,666 175,199 284,698 428,264 532,897 691,062 963,594 2670,567
16 El Haj M’hamed 14,5 98,716 142,151 230,996 347,481 432,377 560,708 781,833 2166,822
17 ElMajene 16,5 103,486 149,020 242,158 364,271 453,269 587,801 819,610 2271,521
18 ElMaleh (dahar) 427,5 223,630 322,027 523,293 787,176 979,498 1270,216 1771,146 4908,669
19 ElGuelta 25,7 119,844 172,576 280,436 421,852 524,918 680,716 949,167 2630,583
20 ElGuessïet 17 104,588 150,607 244,736 368,150 458,096 594,061 828,338 2295,710
Tableau 30 :Les débits fréquentiels ( Méthode de Ghorbel ) pour la région de El Bahayeir

N° Cours d’eau Sup (km2) qmax Q5 Q10 Q20 Q30 Q50 Q100 Q1000
1,44 2,34 3,52 4,38 5,68 7,92 21,95
1 Abdellatif 15,6 101,416 146,038 237,312 356,983 444,200 576,041 803,211 2226,072
2 Zigzaou 2,8 38,008 54,732 88,940 133,790 166,477 215,888 301,027 834,285
3 Radhouane 6 66,143 95,246 154,774 232,823 289,706 375,691 523,851 1451,836
4 Saïdane 7,2 72,873 104,937 170,523 256,514 319,185 413,920 577,156 1599,568

Tableau 31 :Les débits fréquentiels ( Méthode de Ghorbel ) pour la région de Bouflija

N° Cours d’eau Sup (km2) qmax Q5 Q10 Q20 Q30 Q50 Q100 Q1000
1,44 2,34 3,52 4,38 5,68 7,92 21,95
1 Belkacheb 824 247,854 356,909 579,978 872,445 1085,600 1407,810 1963,002 5440,391
2 Larej 670 240,216 345,912 562,106 845,562 1052,148 1364,429 1902,514 5272,749
3 Magran 1540 270,939 390,153 633,998 953,706 1186,714 1538,935 2145,839 5947,117
4 Rechada 108 172,841 248,891 404,448 608,400 757,044 981,737 1368,901 3793,860
5 Mehabés 576 234,636 337,876 549,048 825,918 1027,705 1332,732 1858,316 5150,258

68
PROJET DE CONSTRUCTION D’UN OUVRAGE DE PROTECTION CONTRE LES
INNONDATIONS : CAS D’UN LAC COLLINAIRE

Les Etapes

1- Présentation générale
 localisation
 objectifs
 Caractéristiques climatiques
 Occupation du sol
2- Etude morphologique du bassin versant
 Caractéristiques de forme du bassin versant
 Relief du bassin versant
3- Etude hydrologique du bassin versant
 Pluviométrie moyenne
 Temps de concentration
 Intensités
 Apports liquides moyens
 Répartition mensuelle des apports annuels
 Apports fréquentiels
 Débits fréquentiels de pointe
 Hydrogramme des crues
 Les apports solides
 Qualité des eaux
4- Etude topographique
 travaux et relevés topographiques
 topographie du site
 topographie de la cuvette
 volume de la digue
5- Etude géologique et géotechnique
 géologie régionale
 géologie locale
- bassin versant
- retenue
 Etude géotechnique
- Travaux de reconnaissance
- Lithologie et affleurements
- Essais au laboratoire
- Résultats et interprétations
6- Etude de la digue
 généralités
 conception générale
- dimensionnement
- optimisation de la hauteur
- hauteur des plus hautes eaux
- revanche
- calage de la crête
- largeur en crête
- tranchée d’ancrage

69
  Etude de l’infiltration
- positionnement de la ligne de saturation
- dispositifs de drainage
 débit de fuite
 protection des talus
 tassement
 stabilité de la digue
- méthode de Fellinus
- résultats
7- Les ouvrages annexes
 Evacuateur de crue
- Définition et conception
- Emplacement
- Dimensionnement
 Ouvrage de prise et de vidange
- conduite de vidange
- tour de prise
- chambre des vannes

Chapitre I : Présentation générale

Localisation et accès au site du lac
Le site du lac est situé dans la partie centrale de la Tunisie ; dans le gouvernorat du
Kairouan appartenant à la zone nord de la région naturelle des basses steppes. Le site proposé
est rattaché administrativement au secteur Ouled Khalfallah, délégation de Haffouz, localisé
sur la carte d’état majeur Paviller a échelle 1/50.000 aux coordonnées Lambert : X=487.450
,Y =242.250
Objectifs de la réalisation du lac
L’objectif principal de la réalisation du lac est le développement agricole par l’apport
des eaux de surface en quantité stockée. On peut citer aussi d’autres objectifs qui sont visés
par l’exécution du lac :
 La mobilisation des eaux de ruissellement afin de promouvoir le développement
économique et social suivant une utilisation rationnelle des ressources et la mise en
valeur de l’intensification du système de l’exploitation.
 L’approvisionnement en eau du cheptel de la région et la création d’un microclimat.
 La maîtrise des eaux de ruissellement qui constituent l’agent principal de la
dégradation de l’environnement de la région.
 L’amélioration des conditions de vie de la population de la zone d’influence par les
infrastructures dérivées de la réalisation du lac collinaire et la création d’emploi.

70
 V e r s K a ir o u a n N

P i s t e d 'a c c é s
v e r s H a je b E l A y o u n

$ S i t e d u la c

3 0 3 K il o m e t e r s

$ E x u to ir e 1 . s h p
L i m it e d u b a s s i n v e r s a n t . s h p

Figure 1- Carte de localisation du site de lac collinaire Oued Oummilha

71
 Caractéristiques climatiques
Sur le plan climatique, la région du bassin versant de l’oued Oummilha appartient à
l’étage bioclimatique méditerranéen aride supérieur avec un été chaud et sec et un hiver
tempéré. Le régime pluviométrique présente une irrégularité annuelle et interannuelle avec
des précipitations à caractère orageux. La température est très élevée en été avec plusieurs
jours de sirocco par an.
Le bassin versant de l’oued Oummilha est un petit bassin dépourvu de station climatique
d’où le recours aux données enregistrées par les stations climatiques les plus proches de notre
site qui sont :

Tableau 1- Coordonnées géologiques des stations climatiques (INM, 2005)

Station Altitudes (m) Latitude (Gr) Longitude (Gr)
Kairouan 66 39.64 8.63
Sbiba 627 39.49 7.49
Les températures
La répartition moyenne, maximale et minimale des températures est basée sur les
enregistrements des stations de Kairouan et de Sbiba présentée dans les tableaux suivants :

Tableau 2- Répartition mensuelle des T° à la station de Kairouan (INM, 2005)

Mois
T °C S O N D J F M A M J J A [Link]
T°moy 25.5 21 15.9 12 10.7 12 14 16.9 20.8 25.3 28.3 28.6 19.3
T°max 32.6 27.2 22 17.8 16.5 18.3 20.7 24 28.4 33.6 37 36.9 26.3
T°min 18.3 14.7 9.6 6.5 4.9 5.8 7.3 9.6 13.1 17 19.4 19.4 12.1

Tableau 3- Répartition mensuelle des T° à la station de Sbiba (INM, 2005)

Mois [Link]
T °C S O N D J F M A M J J A
T°moy 22.1 17.1 12.3 9.9 8.9 9.2 10.8 13.4 18.1 23.3 26.1 26.1 16.4
T°max 28.8 23 17.8 15 13.8 14.6 16.7 19.3 25 31.1 34.4 34.2 22.8
T°min 15.4 11.2 6.7 4.8 3.7 4 5.3 7.4 11.2 15.5 17.5 17.7 10
En analysant les données fournies par les deux tableaux précédents on peut conclure
que les températures moyennes annuelles se situent entre 16.4°C et 19.3°C. On observe les
plus basses températures moyennes mensuelles en janvier avec une moyenne de l’ordre de
8.9 °C à 10.7°C et les plus élevées au mois de juillet avec une moyenne supérieure à 26°C.
Les étés sont très chauds avec des maxima moyens supérieurs à 30°C, les hivers sont froids
avec des minima moyens supérieurs à 4°C.

72
 Les températures minimales absolues peuvent atteindre des valeurs inférieures à 0°C à
partir du mois de décembre jusqu’au mois de mars, pendant cette période, il y a un risque
permanent de gelées au sol.
Les précipitations
C’est une zone peu pluvieuse recevant moins de 300 mm de pluviométrie annuelle. En
se référant aux données climatiques fournit au niveau de la station de Sidi Saad jaugeage, la
station la plus proche, les précipitations annuelles sont irrégulières variant entre 112.9 mm et
551.5 mm pour une période d’observation allant de 1973 à 2003. La moyenne annuelle au
niveau de cette station est de l’ordre de 230 mm pour la même période d’observation. Les
précipitation mensuelles sont également très irrégulières, mis à part les mois de Mars et
Octobre, en effet tous les mois peuvent être à secs d’une année à l’autre. Pour la région, la
caractéristique essentielle est la prédominance des totaux de pluies d’automne avec un
deuxième maximum au printemps, ces deux périodes de maximum sont séparées par une
période plus sèche en hiver. Les répartitions mensuelle et saisonnière de la pluviométrie sont
données dans le tableau suivant :

Tableau 4- Répartition mensuelle et saisonnière des pluies

Mois S O N D J F M A M J J A Total annuel
PMOY. Mensuelle (mm) 28.4 24.6 24.5 23.7 22.7 17.0 25.5 22.6 20.1 5.9 3.0 12.2 230.2
PMOY. Mensuelle (%) 12.3 10.7 10.6 10.3 9.9 7.4 11.1 9.8 8.7 2.6 1.3 5.3 100
PMOY. Saisonnière (%) 33.6 27.6 29.6 9.1 100
Les figures ci-dessous représentent la répartition mensuelle et saisonnière de la pluviométrie :

30
25
Pluie mensuelle
moyenne(mm)

20
15
10
5
0
S O N D J F M A M J J A
Mois

Pluie mensuelle

Figure 2 - La répartition mensuelle de la pluviométrie

L’évaporation
Les valeurs de l’évaporation relevées sur bac pour ces deux stations sont reportées dans
le tableau suivant :

73
Tableau 5- Répartition mensuelle de l’évaporation aux stations de Kairouan et Sbiba ( INM, 2005)
Mois Moyenne
Station S O N D J F M A M J J A annuelle
Kairouan 174.6 137.6 113.7 107 105.4 110.7 136.8 150.7 211.8 262.8 293.5 266.7 2071.3
Sbiba 129.9 98.5 82 82 77.9 62.6 85.9 87.4 118 158.8 177.7 175.6 1336.3

Elle dépend essentiellement de la température. Elle est maximale en été pour les stations
en plaine, et minimale pour les stations en altitude et en hiver à cause des températures plus
fraîches. Environ 50% des quantités s’évaporent entre le mois de mai et de août, et seulement
15% s’évaporent en hiver.
Les vents
La direction et la vitesse du vent moyen dépendent de la situation topographique de la
station. C’est la station de Kairouan, que les vitesses du vent sont les plus élevées. Sur les
montagnes, les vents doivent être plus violents.

Tableau 6- Répartition mensuelle de la vitesse des vents à la station de Kairouan ( INM, 2005)
Mois S O N D J F M A M J J A

Vitesse(m/s) 2.7 3.3 2.6 2.3 3.1 2.7 3.9 3 2.7 3 2.8 2.9

En hiver, les directions prédominantes du vent sont celles du nord-ouest et nord/nord-

ouest. En été les vents dominats sont ceux de nord-est et de sud. Pour le tiers des observations
on enregistre un temps calme, en effet, 70% des vents observés ont une vitesse inférieure à 3
m /s. Le sirocco souffle du sud et sud-ouest à partir du mois d’avril jusqu’au mois de
septembre, il est accompagné souvent de vents de sables. Le nombre de jours de sirocco
peuvent atteindre 30 à 45 j /an en moyenne.
L’occupation du sol
Le bassin versant de l’oued Oummilha se caractérise par une mise en valeur très limitée.
La partie nord-est caractérisée par des pentes faibles et une profondeur des sols importantes
constitue la zone de mise en valeur agricole. La partie sud du bassin versant est constitué par
le massif de montagne dont l’occupation est très réduite.
Les cultures annuelles
Les cultures annuelles, conduites exclusivement en sec, représentée par les céréales se
situent au nord du futur site du lac. La superficie des cultures annuelles est assez limitée et est
de l’ordre de 25.54 ha soit 3.3 % de la superficie totale du bassin versant.

74
 L’arboriculture fruitière
L’arboriculture fruitière correspond essentiellement à l’olivier et l’amandier. On les
rencontre au niveau de plusieurs endroits du bassin versant et essentiellement dans la partie
nord-est. La superficie est de 127.71 ha soit 16.5 % de la superficie totale.
La plantation pastorale
La zone pastorale est constituée essentiellement par des jeunes plantations d’Acacia.
Elle se situe au niveau de la zone montagneuse. Les plantations pastorales couvrent 219.04 ha,
soit 28.3% de la superficie totale du bassin versant.
Les parcours
La superficie des parcours est limitée à 16.25 ha soit 2.1% de la superficie totale du
bassin versant. Ils sont localisés essentiellement au niveau des berges des affluents et d’autres
endroits dispersés dans le bassin versant.
La zone inculte
La zone inculte correspond à une superficie de 385.45 ha soit 49.8% de la superficie
totale. Cette zone est constituée par des terrains en pentes caractérisées par une érosion
importante. Le tableau ci-dessous récapitule la répartition spatiale du couvert végétal du
bassin versant de l’oued Oummilha :

Tableau 7- La répartition de couvert végétal du bassin versant Oummilha

Occupation du sol Superficie

(Ha) (%)
- cultures annuelles 25.54 3.3
- arboricultures fruitières 127.71 16.5
- plantation pastorale 219.04 28.3
- parcours 16.25 2.1
- zone inculte 385.45 49.8
La figure ci-dessous explique la répartition spatiale des différentes composantes du couvert
végétal :

chapitre II: Etude morphologique du bassin versant

Les caractéristiques physiographiques d'un bassin versant influencent fortement sa
réponse hydrologique, et notamment le régime des écoulements en période de crue ou
d'étiage. Le temps de concentration qui caractérise en partie la vitesse et l'intensité de la
réaction du bassin versant à une sollicitation des précipitations, est influencé par diverses
caractéristiques morphologiques : en premier lieu, la taille du bassin (sa surface), sa forme,

75
son élévation, sa pente et son orientation. A ces facteurs s'ajoutent encore le type de sol, le
couvert végétal et les caractéristiques du réseau hydrographique. Ces facteurs, d'ordre
purement géométrique ou physique, s'estiment aisément à partir de cartes adéquates. (Musy,
2003)
Planimétrage du bassin versant
Limité au site du lac collinaire et en se basant sur l’analyse de la carte d’état major de
Paviller et de Djebel Trozza à l’échelle 1/50000 ème
, le bassin versant de l’oued Oummilha
couvre une superficie de 7.74 km2 et un périmètre de 13.91 km.
Les prospections du terrain ont permit de caractériser le bassin versant par :
 Une forte pente sur la grande partie.
 Une présence locale de broussaille.
 Une érosion locale relativement faible.
 Une absence de retenue de régularisation en amont du site du lac.
La forme du bassin versant
La forme d'un bassin versant influence l'allure de l'hydrogramme à l'exutoire du bassin
versant. Par exemple, une forme allongée favorise, pour une même pluie, les faibles débits de
pointe de crue, ceci en raison des temps d'acheminement de l'eau à l'exutoire qui seront plus
importants. (Musy, 2003)
1.1 L’indice de Gravélus ou de forme

Cet indice se détermine à partir d'une carte topographique en mesurant le périmètre du bassin
versant et sa surface. Il est proche de 1 pour un bassin versant de forme quasiment circulaire
et supérieur à 1 lorsque le bassin est de forme allongée. (Musy, 2003)

Pour le bassin versant de l’oued Oummilha, S=7.74 km2 et P=13.91 km, ainsi KG vaut 1.40 >1
indiquant que le bassin versant a une forme allongée.

1.2 Le rectangle équivalent

Pour le bassin versant de l’oued Oummilha les dimensions du rectangle équivalent sont les
suivants :
La longueur équivalente : L eq = 6.9 km
La largeur équivalente : l eq = 1.12 km

76
 Le relief

L'influence du relief sur l'écoulement se conçoit aisément, car de nombreux paramètres

hydrométéorologiques varient avec l'altitude (précipitations, températures, etc.) et la
morphologie du bassin. En outre, la pente influe sur la vitesse d'écoulement. Le relief se
détermine lui aussi au moyen d'indices ou de caractéristiques suivants :

1.3 La courbe hypsométrique

100%
90%
80%
70%
Superficie (%)

60%
50%
40%
30%
20%
10%
0%
230 245 260 275 290 305 320 335 350 365 380 395 410 425 440 455 470 485 500

Altitude (m )

Figure 1- Courbe hypsométrique du bassin versant de l'oued Oummilha

1.4 Les altitudes caractéristiques

 Les altitudes maximale et minimale

L’altitude maximale : H max = 504 m

L’altitude minimale : H min = 220 m

 L’altitude médiane

L’altitude médiane correspond à l’altitude lue au point d’abscisse 50% de la surface

totale du bassin, sur la courbe hypsométrique. L’altitude médiane : H med = 303 m

1.5 L'indice global de pente

L indice globale de pente : Ig = 41 m/km

1.6 La dénivelée spécifique

77
 La dénivelée spécifique est donnée par la formule suivante : Ds = Ig* S

Ds : La dénivelée spécifique en m,

Ig : L’indice global de pente en %,

S : La superficie du bassin versant

La dénivelée spécifique : Ds = 114 m ; donc le relief est assez fort.

Le réseau hydrographique
Le réseau hydrographique se définit comme l'ensemble des cours d'eau, permanents ou
temporaires, qui participent à l'écoulement. Le réseau hydrographique est sans doute une des
caractéristiques les plus importantes du bassin versant. La différenciation du réseau
hydrographique d'un bassin est due à quatre facteurs principaux qui sont la géologie, le climat,
la pente du terrain et la présence humaine. Afin de caractériser le réseau hydrographique, on
reporte son tracé en plan sur une carte à une échelle adéquate (Musy, 2005). L’élaboration du
réseau hydrographique est réalisée par le logiciel ArcView, la figure ci-dessous le présente :

$ Ex u toire 1 .sh p
Lim ite d u b a s s in v e rsa n t.s h p
Le ré s e au hy d ro g ra p h iq u e .s h p

3 0 3 K ilo m e te rs

Figure 2- Le réseau hydrographique de bassin versant de l'oued Oummilha

78
 1.7 Le profil longitudinal du cours d'eau
La structure du réseau hydrographique est représentée par son profil en long. C’est la
variation altimétrique du fonds du cours d’eau en fonction de la distance à l’exutoire (Muzy,
2003).

550
500
450
400
Altitude (m)

350
300
250
200
150
100
50
0
250 550 900 1900 3800 4600 5315 5320
Longueur de l'oued (m)

Figure 3- Profil en long de l’oued Oummilha

Divers paramètres descriptifs sont utilisés pour définir le réseau hydrographique :

L’altitude maximale du cours d’eau principal : H0max = 500 m

L’altitude minimale du cours d’eau principal : H0min = 220 m

La longueur du cours d’eau principal : L o =3.8 km

La pente moyenne du cours d’eau principal : I moy =

H 0
max  H min
0

L0

I moy = 21 m/km

Le tableau suivant récapitule les caractéristiques physiographiques du bassin versant de

l’oued Oummilha :

Caractéristiques physiographiques du bassin versant de l’oued Oummilha

Désignation Paramètre Unité Valeurs

Surface de bassin versant S km2 7.74
Périmètre de bassin versant P km 13.91

Longueur de Talweg L km 5.32

Altitude maximale de bassin versant Hmax m 504
Altitude minimale de bassin versant Hmin m 220

79
 Altitude médiane de bassin versant H50% m 303
Altitude à 5% de la superficie H5% m 493
Altitude à 95% de la superficie H95% m 266
Indice global de pente Ig m/Km 41
Indice de compacité de Graveluis KG Adimensionnel 1.4
Dénivelée spécifique DS m 114
La pente moyenne du cours d’eau principal I moy m/Km 21

Chapitre III : Etude hydrologique du bassin versant

L’étude hydrologique consiste à estimer tous les paramètres hydrologiques nécessaires

pour le dimensionnement définitif des ouvrages. Ces paramètres sont l’apport moyen et
fréquentiel, leur distribution, le débit maximal instantané, le volume et l’hydrogramme de
crue, l’apport solide et la qualité des eaux.

La réaction hydrologique d'un bassin versant à une sollicitation particulière est

caractérisée par sa vitesse et son intensité. Ces deux caractéristiques sont fonction du type et
de l'intensité de la précipitation qui le sollicite mais aussi d'une variable caractérisant l'état du
bassin versant : le temps de concentration des eaux sur le bassin. (Muzy, 2003).

Étude de la pluviométrie moyenne

Parmi les méthodes généralement proposées pour calculer la moyenne des pluies à partir
de l'ensemble des mesures ponctuelles obtenues à plusieurs stations pluviométriques sur le
bassin ou à proximité, on distingue la méthode de la moyenne arithmétique, la méthode du
coefficient de pondération ou l'utilisation d'isohyètes. Le choix de la méthode dépendra
notamment de la longueur de la série de données dont on dispose, la densité du réseau de
mesure, et la variation du champ pluviométrique.

1.1 Calcul de la moyenne arithmétique

Caractéristiques des stations pluviométriques pour une période d'observation de 30 ans -1973/2003-
(DRE,2005)

Station Altitude Latitude Longitude Moyenne annuelle observée

(m) (Gr) (Gr) (mm)
Code Nom
66588 Sidi Saad 238 39.3260 8.1770 230.16
Jaugeage
63038 Haffouz DRE 280 39.59806 8.1550 313.50
76066 Sidi Amor 117 39.3280 8.5785 278.50
Bouhajla

80
 Cette méthode donne une pluviométrie moyenne annuelle de 274 mm. Cependant elle
est souvent peu recommandée car elle est peu représentative. Il faut lui préférer des méthodes
graphiques (tracé d'isohyètes) ou statistiques qui permettent de donner un poids différent à
chacun des points de mesures (moyennes pondérée).

1.2 La méthode du coefficient de pondération

1
di2
La formule suivante permet le calcul du coefficient de pondération« Ci » : Ci 
1
 i di2

Ci : le coefficient de pondération pour la station « i » ;

di : la distance entre la station et le centre du bassin versant en Km.
La pluviométrie moyenne du bassin versant s’écrit comme suit :
Pmoy   Pi moy   Ci  Pi
i i

P moy : la pluviométrie moyenne annuelle en mm ;

P i : la pluviométrie moyenne de la station « i » en mm.
Le tableau suivant récapitule les résultats de cette méthode :

Pluviométrie moyenne annuelle par la méthode du coefficient de pondération

Station Pi (mm) di (km) Ci Pi moy

Sidi Saad 230.16 12 0.6 137

Jaugeage
Haffouz DRE 313.5 15.9 0.3 106
Sidi Amor 278.5 36 0.1 18
Bouhajla
La sommation des pluviométries moyenne au niveau de chaque station implique une
pluviométrie moyenne de 262 mm.
1.3 La méthode des isohyètes

C’est la méthode la plus rigoureuse qui demeure moins lourde en utilisant les moyens
actuels tel que le logiciel SURFER, elle est fondée sur l'utilisation des isohyètes. Les
isohyètes sont des lignes de même pluviosité. Grâce aux valeurs pluviométriques acquises au
niveau des stations avoisinantes de bassin versant de l’oued Oummilha, on peut tracer le
réseau d'isohyètes. Le tracé des isohyètes doit être dessiné avec le maximum de
vraisemblance compte tenu de la région, du réseau et de la qualité de la mesure. (Muzy, 2005)

81
 En traçant les courbes isohyètes du gouvernorat de Kairouan, on trouve que notre site
est situé au niveau d’une courbe isohyète ayant comme pluviométrie moyenne 262 mm/an.

Résultats de calcul de la pluviométrie moyenne annuelle par les différentes méthodes

Méthode utilisée Pluviométrie moyenne
calculée (mm)
- Méthode de la moyenne 274
arithmétique
- Méthode du coefficient de 262
pondération
- Méthode des isohyètes
262

Le résultat de la méthode des isohyètes confire celle de coefficient de pondération donc

on admet une pluviométrie moyenne annuelle de 262 mm pour le bassin versant de l’oued
Oummilha.

Isohyètes interannuelles gouvernorat de Kairouan

(1973-74 au 2002-03)

SIDI NAJI
NEBHANA Bge
PONT ROMAIN

SBIKHA DELEGATION
Essbikha
EL ALAM OFFICE
OUSLATIA FORET

El Ouslattia

OUSLATIA INRAT

KAIROUAN SM
KAIROUAN
HAFFOUZ
Haffouz DRE
HAFFOUZ SM
Chebika
El Alla

LAC OUMMILHA

HAJEB LAAYOUNSIDI SAAD JAUGEAGE SIDI AMOR BOUHAJLA

Hajeb Laayoun
Bouhajla
Nasrallah

Lac Oummilha
Délégations
Postes pluviométriques Cherarda

(262/230.2) pluie du site sur la pluie de la station la plus proche sidi saad Jaugeage ( corrigé)

Figure 3- Carte des isohyètes du gouvernorat de Kairouan (DRE, 2005)

temps de concentration :
Résultats de calcul du temps de concentration par différentes formules

Formules Temps de concentration tc (h)

82
 Formule de Giandotti 3.05
Formule de Kent’s 2.36
Formule de Kripich 0.43
Formule de Ventura 2.44
Formule de Turazza 1.20
Estimation du projet national de PCI 1.31

On note que la valeur donnée par la formule de Giandotti (t c = 3.05 h) , Kent’s (2.36
h)et Ventura (t c = 2.44 h) surestiment le temps de concentration, cellesdonnées par Kripich (t
c = 0.43 h) le sous estime, alors que les deux formules de Turazza (t c= 1.2 h) présente une
valeur proche de celle estimée lors du projet de protection contre les inondations, donc on
adopte tc = 1.3 comme un temps de concentration pour notre site.

L’intensité des pluies

Puisqu’on ne dispose pas, au bassin versant de l’oued Oummilha, d’une station

pluviométrique équipé de pluviographe pour nous permettre d’analyser la répartition
pluviométrique. On évalue l’intensité de la pluie au niveau du bassin versant par la courbe
Intensité-Durée-Fréquence de la station de Kairouan.

Figure 4- la courbe IDF de la station de Kairouan (Projet PNUD/FAO TUN/86/020, 1995)

Les résultats des intensités retrouvées pour différentes périodes de retour, mentionnée à
la courbe IDF, sont mentionnés au tableau suivant :

83
 Tableau 5- L’intensité de pluie pour différentes périodes de retour
t c = 79 min
Période de retour T (an) 1 5 10 20 50 100
Intensité en (mm/h) 12.5 19 23 29 35 43

Étude des apports liquides

Apport moyen annuel

Résultats de calcul de lame ruisselée par les différentes formules
Formule Lame d’eau ruisselée (mm) Apport moyen annuel (m3)
Tixeront 4.53 35 062
Fersi 27.5 212 850
Algérienne 7.14 55 263
Sammie 19.7 152 478
Turc 1 43.5 336 700
Mont Marin 32 247 680
Marguellil 23.5 181 900
DG/RE 9.3 72 000
Estimation Chikhaoui 15-20 116 100 – 154 800

Les formules de Tixeront, algérienne et de DG/RE donnent des valeurs très faibles. Les
formules de Turc 1 et de Mont Marin donnent des valeurs
élevées. Alors que les autres formules donnent des valeurs proches. En effectuant la moyenne
entre les résultats des formule de Marguellil, Fersi et Sammie on obtient une lame ruisselée de
l’ordre de 23.56 mm. Si en se réfère à la carte de ruissellement établie pour le bassin versant
de Zéroud et Marguellil (Chikhaoui 1986), on remarque que la lame ruisselé est comprise
entre 15 et 20 mm. Et étant donné que le site du lac Oummilha est très proche de l’oued
Marguellil, on peut juger que la valeur obtenue par cette formule représente mieux le contexte
du centre tunisien, donc on peut adopté la valeur donnée par cette formule ; la lame ruisselée
est égale à 23.5 mm.
D’où l’apport moyen annuel sur le bassin versant de l’oued Oummilha est de l’ordre de 182
000 m3.
La répartition mensuelle des apports
La connaissance des apports mensuels permet le dimensionnement optimal du lac
collinaire ainsi que sa gestion rationnelle. Du fait que le bassin versant de l’oued Oummilha
fait partie de grand bassin versant de Zeroud, les observations hydromètriques faites sur oued
Zeroud ont permis de déterminer la répartition mensuelle d’oued Oummilha.

84
 Répartition mensuelle des apports et de la pluviométrie du bassin versant de
Zéroud(Guide hydraulogique,Ghorbel 1991)
Désignation S O N D J F M A M J J A
Répartition des apports (%) 19 13 6 9 3 3 8 5 7 15 3 9
Répartition mensuelle de la pluviométrie 15 11.4 8.7 11.5 9.6 6.6 11 9.8 5.7 2.3 0.9 7.5
(%)

En tenant compte de la pluviométrie mensuelle, cette répartition doit être corrigée pour
les mois de : Janvier, Février, Mars, Avril, Juin et Juillet.
 Janvier, Février, Mars et Avril : la Répartition des apports de Zeroud sous-estime les
apports de ces quatre mois, ce qui ont été révisés en tenant compte du taux mensuel de la
pluviométrie.
 Juin et Juillet : du fait que la pluviométrie est faible à absente durant ces deux mois, le
ruissellement a été considéré comme nul.
En tenant compte de ces paramètres, la répartition mensuelle des apports, obtenu par la
multiplication de l’apport moyen annuel par le pourcentage correspondant à chaque mois est
donnée par le tableau suivant :

- Répartition mensuelle des apports du bassin versant de l’oued Oummilha

Désignation

Total
Mars

Août
Sept

Juin
Mai
Nov

Avr

Juil
Déc

Jan
Oct

Fév

19 13 6 9 9.6 6.6 11 9.8 7 0 0 9 100

apports (%)
corrigée des
Répartition

35 24 11 16 17 12 20 18 13 0 0 16 182
adoptée des
Répartition

(103 m3)
apports

Apport annuel fréquentiel

La connaissance de la répartition dans le temps des apports pour différentes périodes de
retour est primordiale. Ainsi, la gestion rationnelle, l’optimisation et le dimensionnement des
ouvrages hydrauliques sont fondés sur l’exploitation des apports au niveau de leurs
variabilités annuelles, saisonnières ou mensuelles.

85
Le calcul des apports fréquentiels pour le bassin versant du lac collinaire de l’oued Oummilha
s’effectue selon les méthodes suivantes :
 La méthode des paramètres régionaux
 La méthode d’ajustement statistique.

Résultats du calcul des apports annuels fréquentiel calculés par les différentes méthodes
Période Méthode des paramètres Méthode d’ajustement Apport
de régionaux statistique fréquentiel
retour adopté
T (ans) A T (103 m3)

2 151.060 142.669 146.864

5 242.060 223.616 232.828

10 313.040 277.209 295.123

20 382.200 328.617 355.408

50 467.740 395.159 431.449

100 544.180 445.024 494.602

Les valeurs données par les deux méthodes sont proches, on adopte la moyenne afin
d’obtenir les apports fréquentiels.
Étude des crues
l’étude des crues consiste à déterminer le débit de la crue de projet qui est le débit
maximum devant être évacué par l’évacuateur de crue après laminage tout en évitant les
risques de débordement et de submersion de la digue. La période de retour de la crue de projet
dépend d’une part de l’importance de l’ouvrage et d’autre part des menaces de la sécurité
publique et dommage causés par l’éventuelle rupture de la digue. Pour les lacs collinaires on
adopte généralement une période de retour de T=50 ans tout en vérifiant la capacité
d’évacuation des eaux pour T= 100 ans.

- Résultats du calcul des débits de crue par les différentes méthodes

Méthode utilisée Q2(m3/s) Q5(m3/s) Q10(m3/s) Q20(m3/s) Q50 (m3/s) Q100 (m3/s)
graphique - - - - 73.53 -
rationnelle - 22.7 37.62 48.37 69.44 87
Kallel - - 87.95 116.86 341.44 453.67
Ghorbel 8.67 21.31 34.4 51.74 83.50 116.50
Frigui - - 9.75 14.18 169.20 228.70
Fersi - - - 61.8 89.61 122.05
Froncou-Rodier - - 23.86 54.10 122.69 236.18
QT adopté (m3/s) 8.67 22 36 48.37 69.44 87
Les débits de crue obtenus par les formule de Kallel, Frigui et Froncou-Rodier sont
assez élevés. En effet, ces formules ont étés calées sur des bassins versants de grandes
superficies. Les méthodes rationnelle et Ghorbel donnent des valeurs relativement proches

86
pour les débits de crues T=5 et 10 , pour cela les débits adoptés seront la moyenne entre ces
deux formules. À partir de T=20 ans la formule de Ghorbel aboutit à des valeurs plus
importantes que celles de la méthode rationnelle,qui donne des dédits approximatifs à la
région d’étude, donc on adopte les débits vingtennal, cinquantennal et centennale estimés par
cette dernière.

L’hydrogramme des crues

L’hydrogramme des crues est utilisé pour la conception et le dimensionnement des
ouvrages ainsi que pour la gestion du stock d’eau dans la retenue, il sert aussi à l’estimation
de différentes caractéristiques de la crue tel que : la forme, le volume, le temps de montée et
de base (Projet PNUD/FAO TUN/86/020, 1995). Afin de tracer l’hydrogramme des crues, on
adopte la méthode de « Sokolovsky » qui assimile l’hydrogramme à deux équations
paraboliques, l’une pour la montée de la crue et l’autre pour la décrue. Les équations
s’écrivent de la manière suivante :

2
 t 
Q t = Q max (moy)  
 t 
 pour t  tm
 m 

3
 t t 
Q t = Q max (moy)   d  pour t  tm
 td 
 
Avec :
Q t : débit à l’instant t en m3/s ;
Q max (moy) : débit de pointe de période de retour T en m3/s ;
t m : temps de montée en heure ; tm = 1.3 h ;
t d : temps de décrue en heure ; td=1.95 h.
Pour simplifier l’étude, on adopte à l’hydrogramme de crue un temps de montée égale
au temps de concentration (tm=tc) et un temps de décrue égale à un coefficient multiplicatif au
temps de montée (td= c  tm) où « c » est un coefficient qui varie avec la taille du cours d’eau,
la perméabilité et le taux du boisement du bassin versant ( projet RAB/80/0.11,1987) et il est
généralement pris égale à 3/2 par les bureaux d’étude en Tunisie (Touihri, 2004).
La figure 14 présente les hydrogrammes des crues :

87
 100
90
80
Débit (m3/s) 70
60
50
40
30
20
10
0
0,0 0,3 0,5 0,8 1,0 1,3 1,5 1,8 2,0 2,3 2,5 2,7 3,0 3,2
Temps (h)

Q100 Q50 Q20 Q10 Q5 Q2

Figure 5- Hydrogrammes des crues

Evaluation des volumes des apports de crue

Il est primordial de transformer les débits de pointe estimés par les formules empiriques
en volume, ceci peut être effectué grâce aux formules suivantes :

QT
 582  I g  S 0.556
0.222

Hr
Et VT =Hr  S
QT : débit de crue de période de retour T en m3/s ;
Hr : volume de la crue ramenée à une hauteur moyenne sur le bassin versant en m ;
Ig : indice global de pente en m/km ;
S : superficie de bassin versant en m2 ;
VT : volume de la crue de période de retour T en m3.

Tableau 6- Evaluation des volumes des apports des crues

Période de retour T (ans) QT (m3/s) H r (m) V T (103)
50 69.44 0.016 123.84
100 87 0.021 162.61

Evaluation des apports solides

Le transport solide est par définition la quantité de sédiment transportée par un cours
d’eau. C’est la conséquence de l’érosion qui est fonction de plusieurs agents en particulier les

88
pluies, les ruissellements qui en découlent et les vents. Ainsi que des facteurs qui
conditionnent les quantités de particules arrachées principalement l’érosivité des pluies, la
nature des sols, les caractéristiques de la végétation, la topographie et l’action entropique. La
durée de vie des ouvrages hydrauliques en général et les lacs collinaires en particulier, est
directement liée à l’importance de phénomène de l’érosion. En raison des difficultés de
mesure du transport solide en absence des données d’observation, on a recours aux formules
empiriques qui nous permet d’évaluer l’apport solide.

Résultats de calcul de l’apport solide moyen annuel par les différentes méthodes
Formule Apport solide moyen annuel Apport solide moyen annuel
(t/km2/an) (m3/km2/an)
La formule de CES et FAO 1300 867.80
Formule de Tixeront 568.41 378.94
Formule de Fersi 2739.70 3932
Formule de DGRE 5898.50 1826.46

Pour le bassin versant de Merguillil et Zéroud, l’apport solide est de 1000 à 2000
t/km2/an. Par ailleurs, le suivi de l’envasement de barrage Nebhna sur une période de 25 ans
montre un taux d’érosion spécifique de 1414 t/km2/an et le suivi de barrage Sidi Saad,
construit sur oued Zéroud et contrôlant une superficie de 8950 Km 2, montre une érosion
spécifique de 800 t/km2/an (Bouzaiane, 1986). Donc la valeur d’apport solide moyen annuel
obtenu par la formule de CES et FAO est la plus acceptable. Ainsi on adopte 1300 t/km2/an
comme apport solide moyen.
Qualité des eaux
L’analyse des conditions géologiques du bassin versant de l’oued Oummilha (absence
des formations salifères) et l’enquête menée auprès des habitants peuvent signifier que la
salinité des eaux de l’oued ne dépasse pas 1.5 à 2 g/l. La salinité des eaux de l’oued Marguillil
à Haffouz est de 1 g/l. Autres valeurs enregistrés sur des cours d’eau dans la région d’étude
sont comprise entre 1.9 g/l pour les débits de crue et 2.7 g/l pour les débits d’étiage. Vue la
non disponibilité des analyses d’eau dans le bassin versant et par analogie aux autres bassins
versants, la valeur adoptée pour la salinité des eaux du futur lac collinaire Oummilha est la
moyenne observée dans la région, soit 2 g/l.

Chapitre IV : Etude topographique

L’étude topographique d’un lac collinaire permet de choisir la future retenue collinaire
pour déterminer avec précision l’implantation des côtes des principaux ouvrages, elle permet

89
également d’évaluer les dimensions de la retenue et de la digue ainsi que la caractérisation du
site et de la cuvette.
L’étude topographique d’un lac collinaire permet :
 Le choix de l’emplacement de l’axe de la digue.
 La détermination de la zone inondable par la retenue, d’où sa capacité en
fonction de la hauteur de la digue.
 L’évaluation du volume de remblai de la digue en fonction de sa hauteur.
 L’emplacement de l’évacuateur des crues.
 La délimitation des zones d’emprunt.
Travaux topographiques
L’étude topographique détaillée a été effectuée sur la base d’un plan côté au 1/500 avec
des courbes de niveau équidistantes de 1 m, établi à partir d’un levé topographique réalisé au
niveau du site et de la cuvette avec une densité moyenne de 150 points/ha.
Cependant le plan côté, en plus des deux points repères matérialisant l’axe de la
digue,comprend d’autres points repères situés en dehors du site, constituant des références
pour l’implantation des ouvrages et permettant de reconstituer les repères perdus par les
travaux.
Le plan côté à l’échelle 1/500 nous a permis de :
 Définir l’emplacement optimal de l’axe de la digue et de l’évacuateur de crues, compte
tenu de la hauteur approximative de l’ouvrage, et des particularités topographiques du site.
 Etablir le profil en travers de la vallée au droit de l’axe de la digue.
 Etablir les courbes hauteur- surface et hauteur- volume de la retenue.
 Déterminer le volume de la digue pour différentes hauteurs.
D’autre part, ce plan côté sert comme fond topographique pour :
 Le tracé en plan de l’emprise de la digue et de l’évacuateur de crues.
 La délimitation du niveau de la retenue normale.
Par ailleurs, le plan côté permet aussi l’élaboration de:
 Profils en long et en travers de la digue permettant de déterminer avec précision le volume
de la digue.
 Profils en long et en travers de l’évacuateur de crues, permettant d’évaluer le volume des
travaux de terrassement.

90
 254 Rive gauche Rive droite
252
250
248
Cote (mNGT)

246
244
242
240
238
236
0 20 40 60 80 100 120 140 160
Distance cumulé (m)

Figure 6- Profil en long de la digue

Topographie du site
La topographie du site a été déterminée d’une part sur la base des reconnaissances de
terrain et d’autre part à partir du plan côté et du profil en travers de la vallée au droit de l’axe
de la digue. Pour le lac collinaire de l’oued Oummilha, le site se trouve au niveau d’une vallée
en forme de V, pas trop large mais évasée en haut. En effet elle présente des talus de 40% puis
de 10% sur la rive droite et de 19% puis de 10% sur la rive gauche, ce qui indique que le
volume de la digue devient assez important au-delà d’une certaine hauteur. D’autre part, le lit
de l’oued présente une forme trapézoïdale est peu large de 2 m et peu profond de 2 à 3 m.
Topographie de la cuvette
La topographie de la cuvette a été étudiée d’une part d’après les reconnaissances de
terrain, et d’autre part à partir du plan côté.
Pour le futur lac collinaire de l’oued Oummilha la cuvette est assez profonde et peu
large présente de talus assez fort 20 à 30%, en plus le lit de l’oued présente une forte pente, ce
qui lui donne une retenue peu étendue et assez courte et ayant une capacité de stockage assez
limitée, et par conséquent la digue doit avoir une hauteur élevée pour pouvoir stocker la
grande partie de l’apport du bassin versant.
Pour déterminer la capacité de stockage de la cuvette nous avons procédé par étape :
 Dans une première étape on détermine la surface du plan d’eau correspondant à chaque
courbe de niveau délimitée par l’axe de la digue, ce qui permis de déterminer la variation de
la surface inondée en fonction de la hauteur d’eau, et par conséquent le traçage de la courbe
hauteur- surface de la retenue.

91
 Dans une deuxième étape, on calcule le volume de la retenue correspondant à chaque
hauteur d’eau dans la cuvette, ce qui permis de déterminer la variation de la capacité de
stockage de la cuvette en fonction de la hauteur de la retenue, et par conséquent le traçage de
la courbe hauteur -volume de la retenue.
En effet, à chaque surface du plan d’eau correspond un volume partiel Vi compris entre
deux courbes de niveau successives et équidistantes de 1m, est donné par la relation suivante :

Si 1  Si
Vi   h
2
Vi : volume d’eau entre les courbes de niveau i-1 et i pour une hauteur hi;
Si-1 : surface du plan d’eau correspondant à la courbe i-1 pour une hauteur hi-1;
Si : surface du plan d’eau correspondant à la courbe i pour une hauteur hi ;
∆ h : dénivelée entre les deux courbes de niveau i-1 et i, ∆ h =1m .
Le volume V de la retenue correspondant à une hauteur hi est donné par la relation
suivante :

V  Vi
i

Ainsi, en appliquant ces relations pour chaque hauteur du plan d’eau, on détermine le
volume de la retenue en fonction de sa hauteur. Le tableau et la figure suivants donnent le
volume de la retenue et la surface du plan d’eau en fonction de la hauteur.

Tableau 1- Variation du volume de la retenue et sa superficie en fonction de sa hauteur

Côte (m NGT) Hauteur (m) Volume cumulé (103 Surface du plan d’eau
m3) (ha)
236 0.00 0.000 0.000
237 1.00 0.275 0.055
238 2.00 1.075 0.105
239 3.00 2.750 0.230
240 4.00 5.650 0.350
241 5.00 9.675 0.455
242 6.00 14.700 0.550
243 7.00 20.800 0.670
244 8.00 28.200 0.810
245 9.00 36.850 0.920

246 10.00 46.700 1.050

247 11.00 58.050 1.220
248 12.00 71.000 1.370

92
 249 13.00 86.150 1.590
250 14.00 103.800 1.910
251 15.00 121.950 1.990
252 16.00 143.000 2.200

Lors de la description de la topographie de site on a remarqué qu’il se trouve au niveau

d’une vallée en forme de V, pas trop large mais évasée en haut. En effet elle présente des talus
de 40% puis de 10% sur la rive droite et de 19% puis de 10% sur la rive gauche, ce qui veut
dire que le volume de la digue devient assez important au-delà de 16 m de hauteur.

160 3,2

140 2,8
Volume (1000 m3)

120 2,4

Surface (ha)
100 2

80 1,6

60 1,2

40 0,8

20 0,4

0 0
0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16

Hauteur (m)

Volume Surface

Figure 7- Variation de la surface inondée et du volume stocké en fonction de la hauteur

Volume de la digue
Sur la base du plan côté et des profils en travers de la vallée au droit de la digue au
niveau de l’axe, qui sont des coupes verticales effectuées parallèlement à l’axe de la digue et
équidistantes de 10 m, ainsi on peut estimer le volume de la digue par la méthode des
tranches verticale qui exige les conditions suivantes :
- Une hauteur de décapage de 1 m ;
- Une pente de talus amont de 1/3 ;
- Une pente de talus aval de 1/ 2.5 ;
- Une largeur en crête qui varie entre 3 et 4 m selon la hauteur de la digue :
L=3 m pour H d <10 m ;
L= H d /3 pour H d >10 m.
Chaque tranche verticale présente une surface trapézoïdale Si :

93
 Avec :

Si=
Ei  Li   hi
2
Où E i est l’empattement donné par la formule suivante :
Ei = Li + (Pam+P av)  hi
Li : la largeur en crête de la digue en m ;
Pam : la pente de talus amont en % ;
P av : la pente de talus aval en % ;
hi : la hauteur de la digue en m.
Sachant la surface Si, on calcule le volume Vi : Vi=Si  ei
ei : l’équidistance entre deux verticales consécutives qui vaut 10 m.
Ainsi on peut déduire le volume de la digue pour une cote donnée par la sommation des
volumes partiels Vi. VT= Vi
i

Tableau 2- La la variation de la hauteur en fonction du volume de la digue

Cote (m NGT) Hauteur (m) Volume de la digue (m3)

244 8 7200
245 9 9880
246 10 46500
247 11 58000
248 12 71000
249 13 86000
250 14 104000
251 15 122000
252 16 143000
La figure suivante présente la variation de la hauteur en fonction du volume de la digue :

94
 18

16

14

12
Hauteur (m)

10

0
0 10000 20000 30000 40000 50000 60000
Volume (m3)

Hauteur-Volume de la digue Puissance (Hauteur-Volume de la digue)

Figure 8- Courbe Hauteur-Volume de la digue

Chapitre V : Etude géologique et géotechnique
L’étude géologique se fait à partir de la carte géologique de la région étudiée et sa notice
explicative. Cette étude a permis de faire :
 La description géologique de la cuvette.
 L’appréciation de l’étanchéité de site, de la cuvette et de la stabilité des appuis.
 La description de zones d’emprunt ainsi que la qualité des matériaux de construction
disponibles et l’élaboration des caractéristiques géotechniques du site.
Géologie régionale
Le site de l’oued Oummilha se situe dans la Tunisie centrale à 75 Km de la ville de
Kairouan. La région est caractérisée par deux directions de plissement :
 Une direction Nord/Sud matérialisée par les Djebels Cherahil et Siouf et dont la
terminaison péri-anticlinale butte au nord contre l’oued Zeroud.
 Une direction Nord-Ouest matérialisée par les Djebels Hallouf, Grijima et El Houereb ( au
nord du site).
La région a subi au cours de son histoire géologique des mouvements tectoniques
intenses. Les formations géologiques rencontrées dans la région s’étalent des dolomies gris
bleus du jurassique (Djebel El Houereb) aux ensembles conglomératiques de Dra’a Affane du
quaternaire.
Géologie locale
1.1 Bassin versant

95
 Le bassin versant du site Oummilha est traversé, en allant de l’amont vers l’aval (du plus
récent au plus ancien), par les formations géologiques suivantes :
 Formation Beglia constituée de sables grossiers à stratification oblique et argile ;
 Formation Ain Grab constituée de calcaire gréseux bioclastiques et d’argiles bioclastiques ;
 Formation Messiouta constituée par des alternances de sables, argiles et grés ;
 Formation Cherahil formée par des marnes et calcaires lumachelliques et dolomitiques ;
 Formation El Haria formée par des argiles gypseuses à horizons de calcaires
Le site se situe se situe du point de vue géologique au niveau du flanc Est du synclinal
Sidi Abdelhamid, Oued Faid et dont la base a été affectée par un accident tectonique majeur.

1.2 Cuvette de la retenue

En parcourant la cuvette de la retenue de l’amont vers l’aval, nous pouvons distinguer
les alternances des grès fins tendres jaunâtres débités en bancs métriques avec des lits
d’argiles grises à verdâtres dont l’épaisseur varie de quelques centimètres à une dizaine de
mètres. Quelques bancs de grès dolomitiques roux très durs font ressaut dans la topographie et
individualisent des ensembles gréseux plus au moins perméables et des ensembles argilo-
gréseux étanches pouvant constituer des écrans d’étanchéité aidés en cela par le pendage
amont de 45°. La queue de la future retenue se situe au niveau de la barre dolomitique grise
d’age Yprésien qui repose sur des argiles de la formation El Haria dont l’épaisseur varie entre
8 et 10 m.
1.2.1 Etanchéité de la cuvette
La nature lithologique des formations rencontrées dans les fondations de la future
cuvette de retenue assure une étanchéité acceptable. Les alternances de grès fins jaunâtres et
les lits d’argiles gris à gris verdâtres plaident en faveur d’une étanchéité de la cuvette de la
retenue. Les lits d’argiles, assez nombreux et intercalés entre les bancs de grès assez
perméables, forment des écrans d’étanchéité naturelle à tente circulation amont-aval.
De plus le pendage amont de 30° à 45° des couches rencontrées perpendiculaire au sens
normal de l’écoulement des eaux, présente un avantage certain pour l’étanchéité aussi bien de
la cuvette que de l’axe du site.
1.2.2 Géologie du site
Les affleurements visibles au niveau du site de l’oued Oummilha sont d’age Hauterivien
Valanginien et appartiennent à la formation Melloussi. Il s’agit de grés tendres et fins
jaunâtres à blancs grisâtres parfois nettement individualisés et parfois alternant avec des lits
d’argiles grises à gris verdâtre.

96
 L’appui rive gauche est recouvert par des éboulis de pentes constitués de silts argilo-
sableux, de cailloutis et de blocs rocheux dont le diamètre peut atteindre 40 cm., l’épaisseur
de ces éboulis peut atteindre 70 cm.
Au niveau du lit mineur, les alluvions sont constituées de sables fins à moyens, de galets
et de blocs de dolomies grises et de grès roux très durs parfois métriques.
Sur le berge de la rive droite, au niveau de l’axe du site, nous pouvons observer un
conglomérat à élément de grés dolomitiques arrondis de 0.2 m à 1 m de diamètre dans un
ciment de silt argileux à calcaire induré reposant en discordance sur des grès fins jaunâtres à
blancs grisâtres de la formation Melloussi, l’épaisseur peut atteindre 1.5 m. Au niveau de
l’appui rive droite, les affleurements ne sont pas visibles.
1.2.3 Zones d’emprunts
Les zones d’emprunts sont définies comme suit :
 Rive gauche : les matériaux identifiés sont de sables argileux à silteux avec des éléments
graveleux.
 Rive droite : les matériaux identifié sont formés d’argile grise à gris verdâtre. Ces
matériaux peuvent convenir moyennant un décapage de 0.5 à 1 m de colluvions.
Etude géotechnique
Les travaux de reconnaissance géotechnique sont élaborés dans le but de :
 Déterminer l’épaisseur des terrains de couverture (colluvions, alluvions……)
 Déterminer la nature lithologique exacte des différentes formations rencontrées.
 Prélever des échantillons remaniés et intacts au niveau des fondations du site et de les
identifier qualitativement et quantitativement.
Les travaux de reconnaissances géotechniques
Les travaux de reconnaissances géotechniques ont permis de :
 Etablir les profils et les coupes géologiques.
 Prélever les échantillons représentatifs des matériaux.
 Evaluer les quantités de matériaux disponibles.
 Déterminer la profondeur d’ancrage de la tranchée parafouille d’étanchéité.
Au total, huit puits de reconnaissances d’une profondeur de 2.3 m, ont été exécutés
manuellement au niveau de l’axe du site et au niveau de 5 zones d’emprunts :
 Trois puits de reconnaissances au niveau de l’axe du site ont été réalisés jusqu’à une
profondeur de 2.3m.
- Axe de la digue au niveau de la rive gauche (PA1)
- Axe de la digue au niveau du lit de l’oued (PA2)

97
 - Axe de la digue au niveau de la rive droite (PA3)
 Cinq puits de reconnaissances au niveau des zones d’emprunts ont été réalisés jusqu’à une
profondeur de 2.4m :
- zone d’emprunt amont rive gauche (PE2)
- zone d’emprunt amont rive droite (PE3)
- zone d’emprunt aval rive gauche (PE1)
- zone d’emprunt Med. Ben Ahmed Nsibi (PE4)
- zone d’emprunt Ammar Abidi (PE5)
Les reconnaissances géotechniques de la fondation de la digue et des zones d’emprunts
ont intéressé :
 l’examen visuel et la description des différentes couches rencontrées dans chaque puit.
 La réalisation des essais de laboratoire pour déterminer les paramètres physico-mécaniques
des matériaux rencontrés.
Ces essais concernent :
 La granulométrie et la sédimentométrie ;
 Les limites d’Atterberg ;
 L’essai Proctor normal ;
 Les essais de cisaillement UU.
Lithologie des affleurements et des puits de reconnaissance
les affleurement visible au niveau du site de l’oued Oummilha ont des grés tendres
jaunâtres à blanc grisâtres parfois nettement individualisées et parfois alternant avec des lits
argiles grises à gris verdâtres.
L’appui rive gauche est recouvert par des éboulis de pente constitués de sable fins silteux, de
cailloutis et de blocs rocheux de 40 cm de diamètre. L’épaisseur de ces éboulis peut atteindre
70 cm. Sables fins silteux et grès dolomitique reposant sur des argiles rouges
Au niveau du lit mineur, les alluvions sont constituées de sables fins à moyen, de galets
et de blocs de dolomies grises et de grès roux très durs parfois métrique. Le lit mineur est
constitué de sables fins à moyen, peu graveleux sur au moins 1 m de profondeur. Ces
alluvions surmontent le substratum constitué par des alternances de grès tendres et d’argiles.
On n’a pas pu excaver le puit au-delà de 1 m de profondeur pour présence de nappe
rencontrée à 1 m de profondeur.
Sur la berge rive droite, au niveau de l’axe du site, nous pouvons observer un
conglomérat à éléments dolomitiques arrondis dans un ciment de silt argileux à calcaire

98
induré reposant en discordance sur des grès fins jaunâtres. L’épaisseur peut atteindre 1.5m.
Eboulis à blocaux et sables fins peu graveleux reposant sur des argiles brunâtres.
D’une façon général, les matériaux constituant la fondation de la future digue sont
hétérogènes et varient d’une rive à une autre.
Le substratum, constitué par des alternances de grés tendres et d’argiles fait suite
directement à ces alluvions et ces matériaux d’altération superficielle, se situe entre 0.8 et 1.6
m de profondeur.
Le substratum, relativement imperméable de part cette alternance et la présence de ces
lits d’argile intercalés entre les bancs de grés, ne posera aucun problème d’infiltration sous la
digue. Cependant, il est surmonté de passées graveleuses ou sableuses et particulièrement au
niveau du lit mineur et des matériaux superficielles d’altération sur les rives, il est conseillé de
mettre en œuvre une tranchée d’ancrage jusqu’à ce niveau.
Les essais de laboratoire
Les essais effectués au laboratoire sont les suivants :
- La granulométrie et la sédimentométrie ;
- Les limites d’Atterberg ;
- L’essai Proctor normal ;
- Les essais de cisaillement UU.
 L’analyse granulométrique
L’analyse granulométrique permet de caractériser la texture d’un sol. Elle est définie par
la courbe représentant les pourcentages cumulés des partie des matériaux formés de grains
ayant un diamètre inférieur à une dimension donnée. La courbe est constituée par points à
partir d’essais en laboratoire. Ces essais consistent en un tamisage pour les éléments de
diamètre supérieur à 80 microns et en mesure du temps de sédimentation d’une suspension
pour les éléments fins (Rolly et al, 1977)
 Les limites d’Atterberg
Les argiles ont la propriété d’absorber des quantités d’eau très importante ou, au
contraire, de se dessécher, ceci en fonction des conditions d’humidité aux quelle elles sont
soumises. L’argile a un comportement liquide lorsqu’elle est malaxée avec des quantités
d’eau de plus en plus importantes. Par contre, si elle est desséchée, les grains sont très serrés
et les liaisons deviennent intenses définissant ainsi un comportement solide. Entre ces deux
états, l’argile est malléable et ait un comportement plastique. La transition d’un état à un autre
est très progressive, c’est pourquoi toute tentative pour fixer la limite entre deux états
comporte une d’arbitraire :

99
 - la limite de liquidité WL est la teneur en eau correspondant au passage entre
l’état liquide et l’état plastique ;
- la limite de plasticité WP est la teneur en eau correspondant au passage entre
l’état plastique et l’état solide.
- l’indice de plasticité qui s’exprime en % et qui est d’autant plus important que
l’échantillon contient plus d’argile : IP= WL- WP (Rolly et al, 1977)
 L’essai Proctor normal – essai de compactage
Le compactage est l’ensemble des opérations mécaniques permettant d’augmenter la
résistance d’un sol en réduisant les vides. Cette action augmente la compacité du sol, réduit
les possibilités de déformation et améliore sa capacité portante. Cet essai résulte de la
 
variation de la densité sèche  d  en fonction de la teneur en eau (W). Ainsi on obtient une
w 
courbe en cloche qui présente un point haut que l’on appelle optimum Proctor. Cet optimum
s’explique lorsque :
- la teneur en eau est élevée, l’eau absorbe une importante partie de l’énergie de
compactage et la densité diminue ;
- par contre lorsque la teneur en eau est faible, l’eau a un rôle lubrifiant. Ainsi la
densité sèche augmente avec la teneur en eau. (M’hiri, 2004)
 L’essai de résistance au cisaillement
La résistance au cisaillement d’un sol est fonction de l’angle de frottement interne des
matériaux qui le constituent et de leur cohésion (φ, C). L’étude mécanique élémentaire de la
résistance d’un sol s’effectue sous deux contraintes principales pour une déformation plane :
une contrainte normale au plan de cisaillement (n) et une contrainte de cisaillement (  )
provoquant le mouvement de glissement. (Rolly et al, 1977)
Les parmètres sur la résistance au cisaillement d’un sol sont données par la notion de
courbe intrinsèque :
- Dans le cas d’un sable la contrainte (  ) nécessaire pour provoquer le
glissement est proportionnel à la contrainte (n). On a :  =   tg (58)
- Dans le cas d’un sol argileux, on a :  = C+   tg (59)
La mesure de C et φ s’effectue avec l’appareil triaxial. Pour notre cas on a réalisé l’essai
non consolidé non drainé qui nous donne une idée sur le comportement de notre ouvrage à
court terme.

100
 Résultats et interprétations des essais de laboratoire
Au total, sept échantillons remaniés représentatifs des différents horizons ont été
prélevés et soumis aux essais de laboratoire.
 Appui rive droite (PA3)
Les caractéristiques géotechniques de ces argiles brunâtres qui constitueront les
fondations de l’appui rive droite sont les suivantes :

Tableau 1- Les caractéristiques géotechniques de l’appui rive droite

Caractéristiques géotechniques Résultats
Granulométrie Diamètres (mm) Pourcentage (%)
<5mm 100%
<1mm 98.6%
<80μ 53.2%
<2µ 20.3%
Limites d’Atterberg Limite de liquidité (LL Limite de plasticité
%) (LP %)
25.4 % 13 %
Poids spécifique « γs » 2.68 kN/m3
Coefficient de perméabilité « K » 10-6 m/s

Les fondations rive gauche sont étanches étant donné qu’elles sont constituées par des
silts sableux peu argileux moyennement plastique (IP=13%). Ces argiles sont de type CL
d’après la classification USCS (El Ouni, 2003). Le pourcentage des éléments ayant un
diamètre <2µ est de 20.3% ce qui confère à la fondation une étanchéité satisfaisante,
moyennant une ancrage de la digue dans ce substratum sur une profondeur de 1m.
 Lit mineur (PA2)
Les caractéristiques géotechniques du lit mineur de l’oued Oummilha sont les suivants :

Tableau 2- Les caractéristiques géotechniques du lit mineur

Caractéristiques géotechniques Résultats
Granulométrie Diamètres (mm) Pourcentage (%)
<20mm 100%
<10mm 73.7%
<5mm 66.6%
<1mm 48.3%

<80µ 1.6%
Poids spécifique « γs » 2.7 kN/m3
Coefficient de permiabilité « K » 10-3m/s

Les matériaux de la fondation au niveau du lit mineur sont constitués des sables fins
moyens peu graveleux, propre de type SP d’après la classification USCS.

101
 Le pourcentage des éléments fins inférieurs à 80µ est de 1.6 %. La présence de ces sables
et sables graveleux situés au niveau du lit mineur nous recommandent de prévoir une tranchée
d’étanchéité pour éviter toute circulation d’eau soit au niveau des interlits gravelo-sableux,
soit au niveau du contact (sablo-graveleux)-substratum pouvant entraîner des phénomènes de
renard.
 Appui rive droite
Les matériaux de la fondation au niveau de l’appui rive droite sont constitués par des
sables fins silteux d’une épaisseur de 0.4m. Ces sables reposent sur des grès dolomitiques
jaunâtres jusqu’à 1.6m de profondeur.

Tableau 3- L es caractéristiques géotechniques de l’appui rive gauche

Caractéristiques géotechniques Résultats

Granulométrie Diamètres (mm) Pourcentage (%)

<20mm 100%
<10mm 98.2%
<5mm 96.5%
<1mm 87.4%
<80µ 11.6%
Poids spécifique « γs » 2.73 kN/m3
Coefficient de perméabilité « K » 10-3m/s
 Zones d’emprunts
Les travaux de reconnaissances ont intéressé cinq zones d’emprunts puisque les trois
premiers situés en amont et en aval du site se sont avérés sableuses graveleuses sur une
épaisseur de 2.4 m, a cause de ça on a poussé les recherches vers d’autres zones à l’aval du
site et les travaux ont consisté en l’excavation de deux puits PE4 et PE5 de 2m de profondeur.
Les investigations par puits ont montré que moyennant un décapage de 0.3m de
profondeur, la zone d’emprunts située à l’amont de l’axe du site aussi bien en rive droite
qu’on rive gauche, la zone d’emprunt peut être exploités jusqu’à 2m de profondeur.
Les caractéristiques géotechniques des matériaux sont les suivantes :

Tableau 7- Les caractéristiques géotechniques des matériaux de construction

Caractéristiques géotechniques Résultats
Granulométrie Diamètres (mm) Pourcentage (%)
<5mm 100%
<1mm 99.1-99.2%

<80μ 60.4-80.7%
<2.8µ 23.2-44.6%
Limites d’Atterberg Limite de liquidité Limite de
(LL %) plasticité (LP %)
36.2-45.5 % 19.8-24.9 %

102
 Poids spécifique « γs » 2.68 kN/m3
Coefficient de permiabilité « K » 10-7 m/s
Cohésion « C » 1.26-1.59 kg/m3
Angle de frottement « φ » 13-20.8°

Il ressort de ces analyses que les matériaux d’emprunts rencontrés sont des argiles
sableuses plastiques, type CL selon la classification USCS, et dont l’indice de plasticité varie
entre 19.8 et 24.9. Ces matériaux, une fois compactés à l’optimum Proctor donneront une
densité sèche de 1.75 et une teneur en eau entre 14% et 15% et peuvent avoir un coefficient de
perméabilité de 10-7 m/s ce qui assurera une étanchéité acceptable de la digue.

Chapitre VI : Etude de la digue

Généralités
Suivant la nature et la quantité des matériaux disponibles sur les zones d’emprunts du
site du lac, la topographie du site et la nature des fondations on peut concevoir le type de
digue convenable. Généralement, les lacs collinaires sont réalisés avec des digues en terre
présentant trois schémas principaux : Une digue homogène, une digue à noyau central et une
digue à écran d’étanchéité

1.1 Digue homogène

Elle est formée par du remblai en matériaux suffisamment imperméables (sable

limoneux), par un drain horizontal d’épaisseur constante et de largeur 1/3 à 1/4 de la largeur
de la digue et par une butée drainante qui joue le rôle de mur de pied contrebutant l’extrémité
aval du massif, elle est ancrée dans l’assise de hauteur entre 1/4 à 1/6 de la hauteur de la
digue.(Gassem, 2002)

1.2 Digue zonée ou à noyau

Dans le cas où les matériaux disponibles ne permettent pas la réalisation d’une digue
homogène, on peut concevoir une digue zonée qui se compose :
 d’une partie centrale étanche (noyau) réalisée avec des matériaux imperméables du type
limono-argileux et de forme trapézoïdale ayant une largeur en crête minimale de 3 m,
 de recharges (massifs) constitués de matériaux plus perméables et stables permettant
d’abaisser la ligne de saturation,
 de zone filtrante de transition (entre le noyau et les massifs) dont l’épaisseur peut varier
entre 1 et 2 m. (Gassem, 2002)

1.3 Digue avec masque étanche

103
 S’il y a un manque de matériaux convenables pour la réalisation d’un noyau étanche, on
peut réaliser une digue homogène à masque amont étanche. Le masque amont est une paroi
étanche plaquée sur le talus amont du barrage. Il peut être en béton de ciment ou béton
bitumineux, chapes préfabriquées, membranes souples, etc. Le masque repose en général sur
une couche d’éléments fins drainant et peut nécessiter, suivant sa nature, une couche de
protection. (Gassem, 2002)

Conception générale de la future digue

Généralement, les lacs collinaires sont réalisés avec des digues homogènes du fait que
c’est le type de digue le plus économique et le plus facile à mettre en œuvre lorsqu’on dispose
en quantité suffisante de matériaux ayant une granulométrie étendue et permettant d’obtenir
après compactage des conditions d’étanchéité et de stabilité.
En effet, compte tenu des résultats des reconnaissances géotechniques, des conditions
topographiques du site, et des caractéristiques des matériaux disponibles, ainsi que de la taille
modeste de l’ouvrage, une digue homogène en terre compactée s’avère la mieux adaptée pour
le lac collinaire d’oued Oummilha .
La conception de l’ouvrage a été faite dans l’esprit de simplifier l’exécution. L’ouvrage
doit remplir les deux principales fonctions suivantes :
 Une fonction poids qui exige que les caractéristiques mécaniques des matériaux compactés
doivent assurer la stabilité de l’ouvrage.
 Une fonction écran qui exige que le corps de la digue soit assez imperméable pour éviter
une circulation d’eau pouvant entraîner des sous-pressions importantes et des renards.
Ainsi le corps de l’ouvrage sera constitué d’une digue homogène en terre compactée qui
sera munie d’un dispositif de drainage et d’une tranchée d’ancrage, ainsi que d’une protection
du talus amont avec de l’enrochement contre les effets des vagues.

Dimensionnement de la digue
Le dimensionnement de la digue consiste à déterminer ses caractéristiques géométriques
permettant l’évaluation de son volume à savoir la hauteur, la largeur en crête et les pentes des
talus amont et aval.
Le calage de la crête de la digue, correspondant à sa hauteur Hd se fait après la
détermination de :
 La côte de la retenue normale RN, correspondant à la hauteur utile de la digue Hu ;

104
 Les hauteurs des plus hautes eaux, correspondant à la charge maximale « h » au dessus du
seuil déversant de l’évacuateur de crues ;
 La revanche R.
La détermination de ces paramètres permettra de déterminer la côte crête de la digue
C.C par la relation suivante : C.C=RN+h+R
Ainsi la hauteur de la digue sera donnée par la relation suivante : Hd = Hu+h+R
 La hauteur utile Hu est déterminée par l’étude d’optimisation de la retenue normale ;
 La charge maximale h est déterminée lors du dimensionnement de l’évacuateur de crues ;
 La revanche R est déterminée en fonction de la longueur d’eau et de la hauteur et la vitesse
des vagues.
Optimisation de la hauteur de la retenue
L’optimisation de la hauteur de la digue dépend de plusieurs paramètres notamment :
 La topographie du site et de la cuvette ;
 L’hydrologie du bassin versant ;
 Le coût d’investissement du mètre cube d’eau mobilisée.
Ainsi le volume utile Vu de la retenue est lié à l’apport annuel A du bassin versant par la
relation suivante : Vu=K*A
Où K est un coefficient compris entre 1 et 1.5 pour les lacs collinaires.
Par ailleurs, la hauteur optimale de la digue, doit correspondre au meilleur rendement
topographique du site et au coût d’eau mobilisée le plus faible. Ainsi pour déterminer la
hauteur optimale de la retenue, on a procédé à l’analyse de la variation de ces différents
paramètres en fonction de la hauteur de la digue et celle de la retenue.
En effet, pour chaque hauteur de la digue on calcule :
 Le volume de la retenue ;
 Le volume de la digue ;
 Le rendement topographique ;
 Le coefficient K ;
 Le prix unitaire du m3 de remblai ;
 Le coût global de l’ouvrage ;
 Le coût de m3 d’eau mobilisé.
Les résultats de ce calcul sont donnés dans le tableau suivant :

105
 Tableau 1- Les résultats de ce calcul d’optimisation de la hauteur de la retenue
Cote RN Hu Vu Hd Vd Rt A (m3) K PU Coût Coût/m3
(mNGT) (m) (m3) (m) (m3) (DT) global d’eau (DT)
(DT)
246 10 46500 12 13500 3.44 182000 0.26 11 148000
247 11 58000 13 17000 3.41 182000 0.32 10.5 178500
248 12 71000 14 21000 3.38 182000 0.39 10 210000
249 13 86000 15 26000 3.31 182000 0.47 9.5 247000
250 14 104000 16 31500 3.30 182000 0.57 9 283500
251 15 122000 17 39000 3.13 182000 0.67 8.5 331500
252 16 143000 18 48500 2.95 182000 0.79 8 388000
Sachant que PU est le prix unitaire du m3 de remblai donné par la formule suivante :
PU m3 de remblai= 121. Vd -0.2511

Et le Coût/m3d’eau est donné par la formule suivante :

Coût/m3d’eau= Coût global de la digue/ Vu
En analysant les résultats de tableau ci-dessus :
 La capacité de stockage de la cuvette est assez faible par rapport à l’apport moyen annuel
du bassin versant. En effet, le coefficient K est toujours inférieur à 1 même pour une hauteur
très élevée de la digue. Il est de 0.79 pour une hauteur de l’ouvrage de 18m, d’autant plus que
les matériaux disponibles en dehors de la cuvette ne permettent pas la réalisation d’une digue
dont le volume est au-delà de 50 000 m3.
 La variation du rendement topographique du site en fonction de la hauteur de la digue est
peu significative. En effet, pour une variation de la hauteur de la digue de 50% passant de 12
à 18 m, le rendement topographique diminue de 14.25% passant de 3.44 à 2.95, cette
réduction du rendement topographique est négligeable pour les hauteurs inférieures ou égales
à 16 m, ce qui fait que sur le plan rendement topographique les hauteurs de la digue entre 12
et 16 m sont presque équivalentes.
 De même la variation du coût du m3 d’eau mobilisé en fonction de la hauteur de la digue se
réduit de 15%. Le passage de la hauteur de la digue de 16 à 18 m, n’améliore pas trop la
situation sauf au niveau du volume de la retenue qui augmente de 38% et le coefficient K qui
passe de 0.57 à 0.79. En effet, il engendre une réduction du rendement topographique de 16%
avec une légère variation du coût du m3 d’eau mobilisée. Cependant le volume de la digue
augmente de 56% en passant de 32 000 à 49 000 m3, ce qui peut poser des problèmes au
niveau de la disponibilité des matériaux d’emprunts et du coût de l’ouvrage qui serait de
l’ordre de 400 000 DT. Ainsi, la hauteur de 16m constitue une limite à ne pas dépasser pour
rester dans le cadre des matériaux disponibles et le coût du lac collinaire.
 Le passage de la hauteur de 15 à 16m, permet d’une part l’augmentation du volume de la
retenue de 21% qui passe de 86 000 à 104 000 m3 et le coefficient K qui passe de 0.47 à 0.57,

106
et d’autre part une légère diminution de 5% du coût du m3 d’eau mobilisée avec le même
rendement topographique. Ainsi la hauteur de 16 m présente certains avantages par rapport à
celle de 15 m, ce qui nous conduit à adopter pour le lac collinaire de l’oued Oummilha une
hauteur de la digue de 16 m correspondant à :
 Une côte de la retenue normale de l’ordre de 250 m NGT ;
 Une hauteur utile de 14m ;
 Un volume de la retenue Vu de 104 000 m3 ;
 Un volume de la digue Vd de 32 000 m3 ;
 Un rendement topographique de 3.30 ;
 Un coût de l’ouvrage de 283 500 DT ;
 Un coût du m3 d’eau mobilisée de 2.5 DT ;
 Un coefficient K de 0.57.
La hauteur des plus hautes eaux
La hauteur des plus hautes eaux est la charge maximale au-dessus du déversoir de
l’évacuateur des crues, correspondante à la crue du projet, fixée par le maître de l’ouvrage,
après son laminage. En effet, pour les lacs collinaires, la crue du projet est la crue
cinquantennale tout en vérifiant la capacité de transit de l’évacuateur de crues pour la crue
centennale pour garantir la sécurité de la digue.
La hauteur des plus hautes eaux est calculée lors de l’étude de laminage qui a donné les
valeurs suivantes :
PHE= Hu +1.29 m pour T=50 ans ;
PHE= Hu +1.50 m pour T= 100 ans.
La revanche
La revanche correspond à la tranche de la digue au-dessus des PHE, elle est destinée à la
protection contre le batillage des vagues en vue d’éviter la submersion de la digue. Elle
constitue une sécurité supplémentaire en cas de sous-estimation des crues et compense le
tassement de la digue. La revanche dépend de la longueur du plan d’eau et de la hauteur des
vagues influencée par la vitesse du vent, ainsi elle constituera une tranche supplémentaire
d’amortissement des crues centennale. La revanche est donnée par la formule de Gaillard
(Zannou.B, 2004):
V2
R= 0.75  H 
2g

L
H= 0.5 
3

107
 V= 1.5  2  H
R : revanche de la digue en m ;
H : hauteur des vagues en m, H=0.65 m ;
L : longueur de la retenue (Flech) en Km, L=0.2 Km ;
V : vitesse de propagation des vagues en m/s, V=2.8 m/s.
Tout calcul fait on aura R=0.7 m.

Le calage de la crête de la digue

Le calage de la crête de la digue consiste à la sommation de :
 La côte de la retenue normale CRN = 250 mNGT ;
 La hauteur de plus hautes eaux PHE = Hu +1.3 m ;
 La revanche R = 0.7 m.
La crête de la digue sera calée à 252 m NGT, ce qui correspond à une hauteur de 16m et
une longueur en crête L=140m.
La largeur en crête
La largeur en crête doit être suffisante pour permettre la circulation des engins pour les
travaux de finition et d’entretien. Elle est fonction de la hauteur de la digue, pour des hauteurs
H d supérieures à 10 m, on adopte une largeur en crête de l’ordre de 1/3 de la hauteur de la
digue, ou on la calcule par la formule suivante : LC  1.65 H d

Pour notre ouvrage H d=16m, la formule donne une largeur de 6.6m, ce qui assez élevée
pour un lac collinaire et pour réduire le volume de la digue on adopte une largeur en crête de
5m correspondant au 1/3 de la hauteur de la digue, donc Lc=5m.
La tranchée d’ancrage
La tranchée d’ancrage et d’étanchéité sert à éviter les fuites à travers la fondation, se
prémunir contre tout risque d’érosion régressive et assurer la stabilité de l’ouvrage. La
détermination des caractéristiques de la tranchée d’ancrage s’est basée sur les critères
géologiques et lithologiques des puits de reconnaissances réalisés sur l’axe de la digue, soit un
ancrage dans les argilites gris verdâtre. La fouille centrale d’étanchéité sera aménagée afin
d’enlever la couche superficielle du terrain et s’étendra tout le long de l’axe de la digue sur
134 m. La tranchée d’ancrage aura une section trapézoïdale, une largeur à la base de 3.5 m,
une profondeur de 3 m et des talus latéraux de fruit 1/3.
Etude d’infiltration
Dans les ouvrages en terre, constitués de matériaux plus ou moins perméables, l’eau
stockée dans la retenue a tendance à s’infiltrer dans la digue et la fondation pour venir resurgir
à l’aval. Si certaines précautions ne sont pas prises, ces infiltrations peuvent mettre en jeu la

108
sécurité de l’ouvrage et compromettre sa stabilité. Ainsi il est nécessaire d’étudier ces
phénomènes dans l’objectif de réduire à des valeurs tolérables les débits d’infiltration et de les
contrôler par la mise en place de dispositif de drainage permettant d’éviter l’érosion
régressive provoqué par le phénomène de renard et d’empêcher la saturation du talus aval de
la digue et par conséquent les résurgences à l’aval des eaux infiltrées.
Positionnement de la ligne de saturation
La ligne d’infiltration, improprement appelée ligne de saturation, est la ligne située dans
le plan vertical coupant le massif, et au dessus de laquelle la pression hydrostatique est nulle.
Elle est aussi appelée ligne phréatique. La ligne d’infiltration sépare, dans le massif, la zone
sèche (ou humide) de la zone saturée Pour tracer la ligne de saturation on se base sur le cas
théorique de Kozeny, qui a montrer que pour un massif homogène démuni de drain, la ligne
de saturation peut être assimilée à une courbe parabolique d’axe horizontal et de foyer situé au
pied aval de la digue (Zahaf, 2004).
Pour une digue munie d’un drain horizontal, qui va rabattre la ligne de saturation, la
courbe de Kozeny a pour foyer le point O situé à l’extrémité amont du drain et la ligne de
saturation qui a pour point de départ le point B, sera rabattue au point C. la figure suivante
présentera la ligne de saturation.

d Y

A B

H
C
y0
x
0.7b 0.3b Oa

L’équation de cette parabole est :

Y= 2Y0 X  Y02 , avec Y0= h 2  d 2  d (68) et d= E-0.7p-Ld

h : la hauteur de la retenue normale en m ;

d : distance horizontale entre les points A et A1 en m ;
E : empattement de la digue en m ; E=Lc+H  (Fam+Fav)
Avec : Lc: largeur en crête en m;

109
 H: hauteur de la digue ;
Fam , Fav: fruit amont et aval des talus de la digue.
p : largeur mouillée rapportée à l’horizontale ;
Ld : largeur maximale du drain
Les points remarquables de cette parabole, facilitant le traçage de la ligne de saturation sont :
A (d, h) : point d’intersection de la parabole avec le plan d’eau ;
A1 (0, Y0) : point d’intersection de la parabole avec l’axe des ordonnées y ;
Dans le cas d’une digue homogène munie d’un drain horizontal, le point C se trouve sur l’axe
des abscisses et a pour coordonnées (-b, 0), où :
Y0
OC=b= C  Pour 60°<  <180°
1  cos
 : Angle entre l’horizontale et l’axe du drain ;  =180° (drain horizontal)
C : Coefficient qui dépend de  ; C =1

Y0
D’où b=
2
Y0
C (- , o) : point d’intersection de la parabole avec l’axe des abscisses x.
2
Le traçage de la ligne de saturation par la méthode de la parabole théorique de Kozeny
suppose que le massif est homogène et isotrope, ce qui n’est pas notre cas. Afin de valider
cette méthode dans le cas d’un massif anisotrope on multiplie les dimensions horizontales

KV
(dans l’axe des x) par le facteur correctif suivant :  0.32
Kh

(Généralement on prend KV=[Link])

KV et Kh : perméabilités verticale et horizontale.
Le tableau ci-dessous présente le calcul des différents paramètres de la ligne de
saturation :

Tableau 2- Les paramètres de la ligne de saturation

Paramètres Valeurs calculées Valeurs corrigées
h 14 m 14 m
E 290.62 m 93 m
P 100 m 32 m
Ld 96.87 m 31 m
d 9m 32.6 m
Y0 4.5 m 2.88 m
Ainsi les points remarquables de la parabole auront comme coordonnées :
A (32.6, 14) ;

110
 A1 (0, 2.88) ; C (-1.44, 0).
Dispositif de drainage
1.3.1 Généralités
La prévision d’un dispositif drainant à l’intérieur de la digue sert à protéger l’ouvrage
contre les effets des eaux d’infiltration, rabattre la ligne de saturation et intercepter les eaux
infiltrées. Ce dispositif est généralement constitué par un drain et un filtre.
Les drains sont constitués par du gravier, alors que les filtres sont constitués de couches
successives de sables de granulométrie de plus en plus fine assurant la transition entre le drain
et les matériaux constituant le remblai de la digue.
Pour les lacs collinaires, le drain est généralement constitué par un tapis filtrant, ou drain
horizontal, qui s’étend sur environ le 1/3 aval de l’emprise de la digue. destiné à évacuer les
débits de fuite, le tapis filtrant est muni par un réseau de conduites en PVC perforé qui
débouche dans un petit regard permettant ainsi le contrôle du débit d’infiltration.
1.3.2 Dimensionnement
Afin d’éviter la dégradation et le colmatage du filtre, le sable constitutif doit avoir un
D60
coefficient d’uniformité supérieur à 2.
D10
D’après Terzaghi, la condition à respecter pour une granulométrie étendue est :
F15
5
S15
F : la dimension du grain du filtre ;
S : la dimension du matériau de base.
D’autre part, l’épaisseur de chaque couche du tapis filtrant, doit être au moins de 20 à 30
cm et dans tous les cas elle doit être supérieure à 50 fois le diamètre F15.
Ainsi le tapis filtrant, sera constitué par deux couches de sable de 20 cm chacune
intercalées d’une couche de gravier de 30 cm d’épaisseur.
Débit de fuite
Le débit spécifique de fuite q est donné par la formule suivante (Zannou, 2004) :


q= K m  H 2  d 2   d
0.5

Avec :
Km : la perméabilité de remblai compacté en m/s ; Km =10-7 m/s
H : la hauteur de l’eau dans la retenue en m ; H =14 m
d : distance horizontale entre les points A et A1 en m ; d=32.6 m
q=2.88 10-7 l/s/ml

111
 Protection des talus
Le talus amont de la digue est soumis à l’attaque des vagues sous l’action des vents,
donc il y a lieu de prévoir une protection par de l’enrochement pour éviter l’érosion de la
digue. Le talus aval de la digue est soumis aux effets de l’érosion due au ruissellement des
eaux de pluie, sa protection se fait par de la plantation qui parait suffisante.
L’épaisseur minimale de la couche d’enrochement est donnée par la relation suivante :
e= c*v2
e : L’épaisseur minimale de la couche d’enrochement en m ;
c : coefficient dépendant de la pente du talus et du poids spécifique du matériau
d’enrochement.
Pour un talus amont de 1/3 et un poids spécifique moyen de l’enrochement de 2.5, on a
c= 0.028. V=2.8m/s déterminée lors du calcul de la revanche, ce qui donne une épaisseur de
l’ordre de 22 cm. Pour les lacs collinaires, on adopte pour la couche d’enrochement une
épaisseur constructive de 40 cm.
Le tassement
La construction d’un lac collinaire engendre le tassement des terrains de fondation dont
l’importance dépend essentiellement de la nature de ces terrains, à ce tassement s’ajoute
l’effet du remblai en cours et après constructions. Ainsi le tassement global d’un ouvrage en
terre vaut le tassement du remblai, augmenté de celui des terrains de fondation.
Pour les lacs collinaires, on adopte généralement un tassement résiduel global de 3%de
la hauteur de la digue, ce qui donne pour notre site un tassement de l’ordre de 0.5 m. pour
compenser ce tassement résiduel, on prévoit la surélévation de la crête de la digue de 0.5 m,
contre flèche, au niveau de sa partie centrale.
La stabilité de la digue
La stabilité d’un remblai est déterminée par la capacité à résister, le long des surfaces de
glissement potentielles, aux efforts de cisaillement( poussée de l’eau) et des charges
permanentes dépendant du poids mort du remblai et des pentes des talus (Projet RAB/80/011,
1987).
Le calcul de la stabilité de la digue se fait selon plusieurs méthodes dont la plus utilisée est
celle de Fellenius qui s’effectue par un logiciel STAB.
Principe et Hypothèses de la méthode de Fellenius
La méthode de Fellenius consiste au découpage du terrain en tranches verticales des
faibles épaisseurs juxtaposées et à l’étude d’équilibre de l’ensemble à la limite du glissement
le long du cercle. Il faut calculer les moments des forces stabilisatrices (force de frottement et
de cohésion) et des forces déstabilisatrices (composante tangentielle du poids propre).

112
 La méthode de Fellenius suppose que :
 La surface de rupture pour les talus amont et aval au contact de la quelle il peut y avoir
une surface de glissement est une surface cylindrique à axe horizontal qui apparaît comme un
cercle de glissement dans un profil en travers de la digue.

Talus avant
rupture

Talus après
rupture

Surface de glissement

Figure 2- Principe de la méthode de Fellenius (Gassem , 2002)

 L’interaction entre les tranches est nulle.
Le coefficient de sécurité est donnée par la relation suivante (Zahaf, 2004) :
m
1 b
FS    Ci   Pn  cos  tgi
m
cos
 P  sin 
n 1
n
n 1

Ci : la cohésion des matériaux de construction de la couche « i » correspondant à la tranche n ;

φi : l’angle de frottement des matériaux de construction de la digue ;
Pn : poids de la tranche n ;
 : angle orienté que fait le rayon du cercle passant par le milieu de la base de la tranche avec
la verticale ;
m : nombre total de tranches.
Le calcul se fait pour plusieurs surfaces de rupture sur le profil en travers type de la
digue dans la zone de la plus importante hauteur. En effet, c’est ce profil qui présente le plus
de risque d’instabilité.
Le calcul de stabilité a été effectué en utilisant un logiciel basé sur la méthode de
Fellenius et ce pour différents états de fonctionnement de la digue soient :
 Court terme : talus aval ;
 Long terme : talus aval ;

113
 Vidange rapide : talus amont
Les résultats de calcul de la stabilité
Suite à l’introduction des données suivantes (concernant la géométrie de la section
étudiée, les propriétés géotechnique des matériaux de construction et le niveau d’eau de la
retenue normale) :
 un talus amont de 1/3 et un talus aval de 1/2.5 ;
 hauteur de la digue de 16 m ;
 cohésion de 1.59 kg/m3 ;
 angle de frottement de 20.8° ;
 hauteur de la retenue normale de 14 m.
Le logiciel permet l’élaboration des coefficients de sécurité pour différents rayons pour
différents cas :
 Court terme (fin de construction) FS>1.5 : la digue est soumise à l’effet du poids des
matériaux de construction qui sont des charges permanentes;
 Long terme (après mise en eau) FS>1.5 : la digue est soumise à l’effet du poids des
matériaux de construction et à la pression hydrostatique
Dans ces deux cas le talus aval est le plus sollicité ;
 Vidange rapide FS>1.2 : c’est le cas le plus défavorable pour talus amont, lorsque le niveau
de l’eau libre contenue dans la digue se baisse rapidement, créant un régime transitoire et le
sens des lignes équipotentielles s’inverse. (Philipponnat et Hubert, 2000)
Les valeurs des paramètres rentrant dans le calcul de la stabilité ainsi que les résultats de
calcul sont présentés en détail en annexe.
Le tableau suivant présente les résultats de calcul de la stabilité par la méthode de
Fellenius :

Tableau 3- Les résultats de calcul de la stabilité par la méthode de Fellenius

Cas étudié Coefficient de sécurité FS
Court terme : talus aval 4.30
Long terme : talus aval 1.98
Vidange rapide : talus amont 1.43
D’après les résultats issus du calcul de la stabilité on peut conclure que l’ouvrage est
stable pour les différents cas étudiés.

114
Chapitre VII : Etude des ouvrages annexes

L’évacuateur de crues

Définition et conception
On nomme évacuateurs de crues ou déversoirs les ouvrages complémentaires aux
ouvrages de retenues, qui permettent la restitution des débits de crues excédentaires à l’aval
du lac collinaire. D’une importance primordiale pour la sécurité du lac, l’évacuateur doit être
en mesure d’empêcher le débordement de l’eau par-dessus de la digue et l’apparition de
phénomène d’érosion à l’aval de la digue dans la zone de rejet dans l’oued. Il y a deux
groupes d’évacuateurs de crues ; en charge et de surface, ce dernier est le plus communément
utilisé grâce à son coût de réalisation moins élevé et à sa marge de sécurité plus importante.
L’évacuateur de surface débute par un seuil, ce seuil débite dans un chenal à faible pente qui
amène l’eau à l’aval de la digue, l’eau emprunte ensuite le coursier dont la pente permet de
rattraper la différence de cote entre le niveau de la retenue et celui du lit de l’oued à l’aval
(Alouini, 2005).
L’emplacement de l’évacuateur de crues
L’emplacement de l’évacuateur de crues dépend de plusieurs paramètres :
 La topographie des deux rives en particulier en aval de l’axe de la digue ;
 La pente et la longueur du canal en terre ;
 Le volume des terrassements ;
 L’accès pour l’approvisionnement en eau.
Pour le cas de notre site, la prise en compte de ces paramètres, a conduit à placer
l’évacuateur des crues sur la rive droite qui présente un petit talweg qui pourra servir pour la
mise en place du canal en terre pour réduire les travaux de terrassement.
Dimensionnement
1.3.1 Le seuil déversant- Laminage de la crue
Le dimensionnement du seuil déversoir de l’évacuateur de crues se fait sur la base de la
crue de projet, pour les lacs collinaire c’est la crue cinquantennale, et de manière à pouvoir
transiter en pleine section et sans risque de débordement la crue centennale après son
laminage au niveau de la retenue, ce qui nous amène à l’étude de laminage des crues. Le
laminage d’une crue est la diminution du débit de pointe de son hydrogramme par des moyens
naturels ou artificiels.

115
 Par ailleurs, si l’onde de crue trouve au cours de son déplacement un réservoir muni
d’un système de vidange quelconque, un certain volume de la crue servira à remplir le
réservoir jusqu’à la côte de déversement. Ensuite, le débit sortant suivra la courbe
caractéristique du réservoir ou des vannes de fond, et l’hydrogramme de sortie présentera une
pointe plus faible et décalée dans le temps par rapport à celle de l’hydrogarmme d’entrée.
Dans les petits bassins versants, les débits de crue instantanés sont souvent
relativement élevés, mais le volume des crues est en général assez faible et les crues peuvent
être laminées de façon appréciable par stockage dans la tranche disponible entre la côte du
déversoir et celle des plus hautes eaux.
L’étude du laminage consiste soit à déterminer la côte maximale atteinte par le plan
d’eau pendant la crue, la largeur du déversoir étant fixée à priori, soit la longueur, la hauteur
au-dessus du seuil étant fixée à priori(Alouini, 2005).
Le mécanisme du laminage est basé sur l’équation de continuité:

Qe(t). dt = Qs(z)dt + A(z). dz

Qe(t). dt : élément de volume entrant

Qs(z)dt : élément de volume sortant

A(z). dz : variation de volume de stock dans la retenue

t : temps en s ;
z : côte du plan d’eau en mNGT ;
Qe : débit de crue entant dans la retenue en m3/s ;
Qs : débit sortant de la retenue par l’évacuateur de crues en m3/s ;
A(z) : aire du plan d’eau à la côte z en m2.
L’étendue du plan d’eau supérieur de la retenue et l’hydrogramme de crue, déterminés
lors de l’étude hydrologique, permet donc de procéder à l’étude du laminage.

L’hypothèse de départ est qu’à l’instant zéro, le réservoir est plein, c’est à dire que le
niveau d’eau est à la côte de retenue normale. On fixe alors, une largeur de déversoir et un pas
de temps ∆t assez court permettant de supposer que le changement de l’hydrogramme
d’entrée et linéaire pendant la période de laminage ∆t.
Nous devons pour chaque pas de temps et pour la largeur fixée déterminer le débit
sortant par le déversoir et le volume stocker au niveau de la retenue.

116
 Pour ce faire, nous avons procédé à une programmation très simple sous EXCEL. Cette
programmation fonctionne d’une manière itérative. En faisant entrer le débit de crue
et l’équation de la courbe de volume on calcul successivement :
 Le volume à l’entrée qui n’est autre que le débit d’entrée par le temps.
 La côte maximale atteinte par le plan d’eau.
 La hauteur déversée.
 Le débit sortant par le déversoir qui est donnée par la formule suivante:

Q = μ(2g)1/2Lseuil.h3/2
L : largeur déversante en m ;
h : charge sur le déversoir en m ;
µ : coefficient du débit du seuil, pour un déversoir à seuil large µ=0.38
Toutefois le calcul est limité par la fixation d’une charge maximale au-dessus du seuil
déversant : 0.6 <h<1.5 m environ (projet RAB/80/011, 1987)
L’étude du laminage de la crue cinquentenale et centennale a conduit aux résultats
suivants :
Tableau - L étude du laminage de la crue cinquentenale et centennale
L seuil (m) T=50 ans T=100ans
Qe (m3/s) Qs(m3/s) h (m) Qe (m3/s) Qs(m3/s) h (m)
24 69.44 66.14 1.39 87 83.19 1.61
25 69.44 66.28 1.35 87 83.43 1.57
26 69.44 66.41 1.32 87 83.66 1.54
27 69.44 66.59 1.29 87 83.87 1.50
28 69.44 66.76 1.26 87 84.06 1.47

Les figures suivantes présentent le laminage des crues cinquantennale et centennale :

117
 80
70
Débit(m3/s)
60
50
40
30
20
10
0
0 0,25 0,5 0,75 1 1,25 1,5 1,75 2 2,25 2,5 2,75 3 3,25
Temps(h)

Q50 Qe

Figure 20- Laminage de la crue cinquantennale

100
80
Débit (m3/s)

60
40
20
0
0 0,25 0,5 0,75 1 1,25 1,5 1,75 2 2,25 2,5 2,75 3 3,25
Temps (h)

Q100 Qe

Figure 21- Laminage de la crue centennale

Interprétation des résultats de l’étude de laminage
D’après l’étude de laminage on remarque que le volume laminé est assez faible ceci
s’explique par la topographie de la cuvette qui est assez profonde et peu large présentant de
talus assez fort, en plus le lit de l’oued présente une forte pente, ce qui lui donne une retenue
peu étendue et assez courte et ayant une capacité de stockage assez limitée, et par conséquent
l’élargissement de l’évacuateur du crue au-delà de 27 m ne diminue pas trop la hauteur d’eau
au dessus du seuil déversoir.
Description de l’évacuateur de crue

L’évacuateur de crue est dimensionné pour évacuer la crue cinquentennale laminée

(Q=66.14 m3/s).La solution proposée est un évacuateur de crue à dégâts maîtrisable. Il est

118
constitué de trois éléments à savoir, le seuil déversoir, le canal coursier et le bassin
d’amortissement à ressaut hydraulique.
La configuration topographique des pentes de la rive droite nous parait plus propice à
l’implantation de l’évacuateur de crue. Dans l’étude de laminage de crue nous avons pris une
largeur de seuil déversoir est de 27 m, le niveau de l’eau atteint dans la retenue est de 1.29 m
au dessus de seuil déversoir.
Canal coursier

Le canal en terre se situe à l’aval immédiat du seuil déversant. Il sera constitué en terre
et aura une forme trapézoïdale avec de talus de 2/3 et une largeur à la base de 27 m jusqu’au
bassin de dissipation. Compte tenu des conditions topographique du terrain, le canal aura des
pentes successives de 15.7 et 13% et une longueur total de 125 m. Le passage en écoulement
torrentiel au niveau du coursier est assuré par le rétrécissement formé au niveau du
changement de pente (section de contrôle).

Bassin d’amortissement

Le canal en terre débouche dans un bassin de dissipation qui sert à amortir l’énergie des
crues arrivant sur le coursier avec une grande vitesse. Cette énergie dissipée peut être
exprimée en fonction de l’état critique du courant à la section initiale du ressaut ou, plus
précisément, en fonction de la valeur du nombre de Froude F en cette même section.
V
F=
g y

V : vitesse de l’écoulement en m/s ;

y : le tirant d’eau en m ;

g : l’accélération du pesanteur ; g=9.81m2/s.

Soient L la longueur du bassin et l sa largeur et soient y1 et y2 respectivement les tirants

d’eau amont et aval dans le bassin de dissipation.
La vitesse d’eau dans le bassin de dissipation est donnée par la relation suivante :
V=0.8  2 g  ( H  y1 )

Avec :
H : hauteur d’eau par rapport au fond du bassin de dissipation en m ;
y1 : tirant d’eau en amont du bassin de dissipation en m.

119
 Le débit Q à restituer par le bassin de dissipation est celui de l’évacuateur de crue pour
la crue de projet, soit Q50= 69.44 m3/s.
 Calcul de y1 ; tirant d’eau en amont du bassin de dissipation :
La détermination de y1 se fait par un calcul itératif qui consiste à la modification de la
valeur de y1, ensuite on détermine à chaque fois la vitesse V et le débit Q jusqu’à obtenir
Q=Q50 .
Pour le lac collinaire de l’oued Oummilha, la largeur du bassin de dissipation est de
l’ordre de 27 m, le débit du projet Q50= 69.44 m3/s.
 Calcul de y2 ; tirant d’eau en aval du bassin de dissipation :
Le tirant d’eau aval y2 est déterminé en fonction de y1 par la relation suivante :
y1
Y2=  1 8 F 2 1
2
 Calcul de la longueur L du bassin de dissipation :
La longueur L du bassin de dissipation correspond à celle du ressaut. Par définition c’est
la distance entre sa face amont et la zone atteinte lorsque toute énergie est pratiquement
dissipée et ne provoque plus d’érosion que l’écoulement fluvial. L est donnée par la relation
suivante :
L=6  (y2- y1)
 Calcul de y3 ; tirant d’eau au niveau du lit de l’oued :
Soit y3, le tirant d’eau au niveau du lit de l’oued, donné par la formule de Manning-
Strickler (Lebdi, 2004) :
Q50= KS  S  R2/3  I
Avec :
Q50 : le débit du projet ;
KS : coefficient de Strickler ; KS = 42 m1/3/s ;
S : la section mouillée (m2) ;
P : le périmètre mouillé (m) ;
R : rayon hydraulique (m) ;
I : pente du tronçon du cours d’eau ; I= 1.5 %
En réalisant des itérations sur EXCEL, qui consiste à la modification de la valeur de y3,
ensuite on détermine le débit Q jusqu’à obtenir Q=Q50.
 Calcul de p la profondeur du bassin de dissipation :
La profondeur du bassin de dissipation est donnée par l’expression suivante :

120
 p= y2  y3  Z
p : profondeur du bassin (m) ;
y2 : tirant d’eau en aval du bassin de dissipation (m) ;
y3 : tirant d’eau au niveau du lit de l’oued (m) ;
 : coefficient de majoration qui vari entre 1.1 et 1.3 ;
Z : perte de charge à la sortie du bassin, généralement elle est de l’ordre de 0.2 m.
Le calcul hydraulique relatif au bassin de dissipation du lac est récapitulé dans le tableau
suivant :

Le calcul hydraulique du bassin de dissipation

Désignation Valeur
F : nombre de Froud 11.7
V : vitesse d’écoulement 14.2 m/s
y1 : tirant d’eau en amont du bassin 0.15 m
y2 : tirant d’eau en aval du bassin 2.4 m
L : longueur du bassin de dissipation 13.5 m (Pour plus de sécurité on prend L=15 m)
y3 : tirant d’eau au niveau du lit de l’oued 0.88 m
p : profondeur du bassin 1.8 m

Ouvrages de prise et de vidange

Pour les lacs collinaires, les ouvrages de prise et de vidange sont généralement
combinés. Ils sont constitués par une conduite placée sous le corps de la digue, et munie d’une
tour de prise à l’amont et d’une chambre de vannes à l’aval, permettant le contrôle des
différents usages d’eau.
Ces ouvrages ont pour objectifs essentiels :
 Le prélèvement des quantités d’eau nécessaires à l’irrigation en fonction de la demande ;
 Le dévasement de la retenue ;
 La vidange rapide de la retenue en cas d’incident ;
 L’alimentation en eau d’abreuvoir et éventuellement de la population.
Conduite de vidange
Conception
Le long de la conduite de vidange il y a des risques de formations de renard, liés aux
écoulements préférentiels susceptibles de développer à la surface de celle-ci. Ainsi du fait que
la conduite fonctionne sous pression, il faut s’assure de son étanchéité absolue, car toute fuite
peut introduire des pressions interstitielles et peut provoquer la formation du renard. D’autre
part pour rallonger les lignes de fuites et éviter les phénomènes de renard, des écrans anti-

121
renards seront placés perpendiculairement à la conduite de vidange et réparties le long de
celle-ci avec des écarts de 10 à 15 m. par ailleurs la conduite de vidange sera soumise à
certaines sollicitations mécaniques qui peuvent provoquer sa rupture. Ces sollicitations sont
dues essentiellement aux passages des engins de compactage et au tassement. L’enrobage de
la conduite avec du béton armé coulé en pleine fouille afin d’assurer un bon contact entre la
conduite et le remblai de la digue pourra éviter les risques de rupture, en lui permettant de
résister au différentes sollicitations.
Pour éviter les risques de décollement du sol sous la conduite enrobée et les risques de
sa rupture du fait qu’elle travaille en poutre, dans le cas des fondations compressibles sous
l’effet des tassements différentiels, l’enrobage sera réalisé en tronçon réunis par des joints
water stop et du polystère. La conduite sera constituée de 7 m en béton fretté reliés par de
joints souples, ce qui lui permettra d’avoir plus de souplesse.
Dimensionnement
Le débit de la conduite de vidange varie en fonction du temps du fait qu’il dépend de la
charge d’eau au dessus de la conduite, en effet il est maximum pour une retenue pleine et
minimum pour une retenue presque vide. Ainsi, pour déterminer le diamètre optimal de la
conduite de vidange on calcule les vitesses d’abaissement du plan d’eau dans la retenue par
tranche de 0.5 m, et ce pour différents diamètres.
Le diamètre de la conduite de vidange sera déterminé en fonction des conditions de
vidange, compte tenu des pertes de charge linéaires et singulières dont les expressions sont les
suivantes :
 Les pertes de charge linéaires :
Vc2
JL= Lc 
K S2  Rh4 / 3
JL : perte de charge linéaire (m) ;
Lc : longueur de la conduite (m) ;
Vc : vitesse de la conduite (m/s) ;
Ks : coefficient de rugosité de Strickler (m1/3/s) ; Ks =70 m1/3/s pour les conduite en béton ;
D
Rh : le rayon hydraulique (m) ; Rh = pour les conduites ayant une section circulaire.
4
 Les pertes de charge singulières :
Vc2
Js= (0.55+C) 
2g
Vc : vitesse de la conduite (m/s) ;

122
Js : pertes de charge singulières (m) ;
Vc : vitesse de la conduite (m/s) ;
C : coefficient de débit, dépendant du degré d’ouverture de la vanne.

Les coefficients des débits en fonction du degré d'ouverture de la vanne (Touihri, 2004)
% ouverture de la vanne 30 50 70 80 100
C 120 19 3.5 1.5 0.2

Et des contraintes suivantes :

 La vitesse d’abaissement du plan d’eau lors de la vidange doit être assez faible pour éviter
de créer des dégâts dans le corps de la digue. Généralement on se limite à une vitesse de
1.5m/jour.
 Une marge de vitesse maximale et minimale dans la conduite de vidange de l’ordre de 1 et
6 m/s.
La vitesse dans la conduite de vidange ne doit pas être ni très élevée ni trop faible. En
effet, on admet généralement une vitesse maximale de 6 m/s et une vitesse minimale de 1m /s
permettant d’assurer l’autocurage de la conduite en faible charge.
D’après le théorème de Bernouilli, assimilant la retenue à un déversoir alimentant une
conduite qui débouche à l’air libre on a :
2
VRN V2
Z RN  PRN   Z C  PC  C  J L  J S
2g 2g
PRN : la pression à la surface de la retenue normale ; PRN = Patm ;
PC: la pression en aval de la retenue normale ; PC = Patm (débouchant à l’air libre);
ZRN : cote de la retenue normale (mNGT) ;
ZC : cote avale de la conduite de vidange ;
JL : perte de charge linéaire (m) ;
Lc : longueur de la conduite (m) ;
Vc : vitesse de la conduite (m/s) ;
Js : pertes de charge singulières (m) ;
2 g  Z RN  Z C 
VC 
2 g  LC
1  0.55  C   4
K C2  RH3

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

ZRN : cote de la retenue normale (mNGT) ;

123
ZC : cote avale de la conduite de vidange ;
Lc : longueur de la conduite (m) ;
Ks : coefficient de rugosité de Strickler (m1/3/s) ; Ks =70 m1/3/s pour les conduite en béton ;
Rh : le rayon hydraulique (m) ;
C : coefficient de débit, dépendant du degré d’ouverture de la vanne.

Ainsi pour déterminer le diamètre optimal de la conduite de vidange on calcule les

vitesses d’abaissement du plan d’eau et celles dans la conduite par tranche de 0.5 m, et ce
pour différents diamètres et compte tenu des caractéristiques de la retenue et de la longueur de
la conduite.
Pour le lac collinaire de l’oued Oummilha, la conduite de vidange aura une longueur
d’environ de 75 m, calée aux cotes amont 242 m NGT et avale 240.95 m NGT, envisageant
ainsi une pente de l’ordre de 1.5%. Les résultats de calcul hydraulique effectuer pour les
diamètres de la conduite de 250 à 450 mm sont donnés dans le tableau suivant (le détail de
calcul est dans l’annexe) :

Le dimensionnement de la conduite de vidange

Diamètre Temps de vidange Vitesse dans la conduite (m/s) Vitesse
(mm) (jour) 50% 100% d’abaissement
(m/j)
300 7.88 1.16 3.97 0.62
350 5.31 1.19 4.33 0.88
400 3.79 1.22 4.65 1.18
450 1.28 1.24 4.94 1.54
D’après le tableau précédent qui synthétise les calculs hydrauliques effectués pour
différents diamètres et degré d’ouverture de la vanne de la conduite de vidange, détaillés dans
l’annexe, on remarque que la marge des vitesse dans la conduite de vidange est respectée et la
vitesse moyenne d’abaissement du plan d’eau varie entre 0.62 et 1.54 m/s, or la limite pour
laquelle la vitesse moyenne d’abaissement ne présente aucun risque pour la digue est de
l’ordre de 1.5 m/s. Ainsi, les diamètre Ф 350 et Ф 400 sont les plus convenables, mais du fait
que la conduite Ф 350 présente des risques de colmatage et que le vidange des lacs se fait
généralement lors des crues et pour les premiers mètres de la retenue, la conduite Ф400
fonctionnera mieux dans les conditions hydrauliques exigées ; faible vitesse d’abaissement du
plan d’eau de 1.18 m/s et fortes vitesses de vidange de 1.22 (50%) et 4.65 (100%) m/s.
Ecrans anti-renard

124
 La mise en place des écrans anti-renard le long de la conduite de vidange permettent de
rallonger les li les chemin d’écoulement et donc de diminuer le risque de percolation et
d’éviter la formation des renards.
Leur nombre dépend de la charge hydrostatique à laquelle est soumise la conduite, est
déterminé en respectant la règle de LANE (Zannou, 2004) :
Lv+Lh/3  C.H
Lv  C.H-Lh/3
Lv : longueur totale des cheminements verticaux en m ;
Lh : longueur totale des cheminements horizontaux en m ; Lh=75m
H : hauteur maximale d’eau à l’amont de la conduite ; H=8m
C : coefficient de LANE fonction de la nature des matériaux ; C=6.
Pour le lac collinaire de l’oued Oummilha on a Lv  23m, on adoptant un écran
dépassant le diamètre de la conduite de 2 m, on aura un cheminement vertical par écran de
4m. Ainsi la règle de LANE pourra être respectée par la mise en place de six écrans anti-
renard avec un écartement de 10 à 15m.
Tour de prise
Lors de calage de la tour de prise il faut tenir compte de la tranche morte afin de ne pas
se trouver dans l’obligation de la surélever après quelques années de fonctionnement.
Généralement on adopte une tour de prise de hauteur comprise ente 1/3 et 1/2 de la hauteur de
la retenue normale, ce qui permettra d’être en sécurité d’une part et favoriser l’exploitation de
la tranche supérieure de la retenue.
Pour permettre l’exploitation de la tranche inférieure de la retenue, la tour de prise sera
munie latéralement d’un certain nombre d’orifices pourvus de grilles métalliques disposées à
différentes hauteurs. Pour notre cas on adopte une tour de prise circulaire en maçonnerie de
2m de diamètre intérieur et aura une hauteur totale de 6m. Elle sera munie de 6 orifices carrés
de 0.5 m de coté. En effet, la crête de la tour de prise est calée à la cote 248 mNGT, dépassant
d’au moins de 1 m la tranche morte, qui atteindra la cote 247 mNGT au bout de 20 ans de
fonctionnement. La tête de la tour sera protégée par une grille métallique sous forme conique.
La tour de prise reposera sur un radier en béton armé qui sera calé à 242 mNGT et au niveau
duquel sera raccordée la conduite de vidange. Il est à signale qu’il a été prévu la mise en place
d’un raccordement souple entre la tour de prise et la conduite de vidange pour amortir les
effets de tassement différentiels.

125
 Chambre de vannes
La chambre de vanne abritera un dispositif de vannage permettant de contrôler les
opérations de vidange et d’alimentation des différents usages. Ce dispositif sera constitué de
deux vannes papillon Ф 400 en série dont une d’arrêt, et d’un robinet-vanne Ф200 en
dérivation pour servir à l’alimentation d’un abreuvoir et d’une potence.
La chambre de vannes sera dimensionnée de manière à faciliter la manipulation des
vannes avec des ouvertures à panneaux amovibles aménagés au plafond. Un bassin
d’amortissement est contigu à cette chambre de vanne. Il est à signaler qu’un canal en
enrochement sera réalisé à l’aval de la chambre de vannes, permettant de restituer les
écoulements de vidange dans l’oued sans risque d’érosion.

EXAMEN de PROTECTION CONTRE LES INONDATIONS

2 HAR , 2 IBR , 2 GS , Documents non autorisés , durée 2 heures
Mardi 22 Juin 2004 de 8h à 10h

Questions de Cours (10 points ):

1- Donner les circonstances et les hypothèses d’utilisation de la méthode du

gradex ainsi que les différentes étapes de son application (3 points )

2- Donner le principe , les hypothèses , les paramètres et les procédures

pratiques d’utilisation de la méthode Socose ( 2 points )

3- Décrire brièvement les différentes études préliminaires nécessaires pour

la création et la construction d’un lac collinaire ( 2 points )

4- Donner et définir brièvement les différents éléments constitutifs d’un

évacuateur de crue ( 2 points )

Questions à réponse courtes ( maximum 5 lignes ) : 1 point par question

1- citer quelques pays souvent touchés par les inondations

2- citer quelques dates repères des inondations en Tunisie et les villes

touchées

3- citer quelques ouvrages de protection contre les inondations en Tunisie

126
 4- citer quelques formules ( sans donner la formule ) utilisées en Tunisie
pour l’estimation de l’apport moyen

5- citer quelques formules ( sans donner la formule ) utilisées en Tunisie

pour l’estimation de l’apport fréquentiel
Surface Indice Hmax Hmin H50% Pmm E(mm)
Nom IG
(km2) kc m m
255 10 65 470 1100
Gueriana 84.5 1.27 8.72

Sidi 185 10 58 475 1050

Hcine 52.5 1.35 11.40
Sijoumi
Fouchana 150 10 45 475 1100
31 1.36 12.20
Mhamdia
80 10 64 470 1050
Mourouj 19.5 1.54 4.27

6- citer quelques formules ( sans donner la formule ) utilisées en Tunisie

pour l’estimation du débit fréquentiel

7- citer quelques grands bassins versants en Tunisie et leurs principaux cours

d’eau
8- Comment utiliser la méthode Socose pour estimer des crues centennales

9- Quels sont les paramètres qui aggravent les conséquences des inondations

10- Quels sont les moyens amont et aval pour se protéger contre les
inondations
EXAMEN DE Protection Contre les Inondations : Session de Rattrapage : 2004/2005
Documents non autorisés
Bilan hydrologique de sebkhat Essijoumi

127
En utilisant les formules empiriques appropriées calculer les apports moyens annuels et
fréquentiels ( T=10,20,50 et 100 ) ainsi que les débits fréquentiels des différents bassins
versants dont les cours d’eau se jettent à Sebkhat Essijoumi. Déterminer le bilan annuel et
fréquentiel du [Link] et discuter les résultats
p3
Formule de Tixeront : Le= , Le en m, P en m , E en m
3.E 2
Formule de Fersi : Le = 163,9. P. I g . 10-4 le en mm, P en mm
5/3 1/ 3
P  Hmed 
Formule de Turc : Le =     , P en mm, Hmed ( H50%) en m
 34   100 
Formule de Kallel : Q(T) = 5,5.S0,5.T0, 41
QT max
Formule de Ghorbel : RT,Q= avec Qmax (moy) = 2,86 .S 0.8 , Q(T) en m3/s, S en
Q moy max

km2 et pour la zone d’étude on RT,Q = 1,47 . T 0.4 -1,35

A(T )
Formule de Ghorbel :  RTV , avec le Nord : RTV = 0,92 Ln(T) + 0,66
Amoyen
Formule de Fersi A(T) = 0,21 Amoyen T0,89.

EXAMEN de PROTECTION CONTRE LES INONDATIONS

2 HAR , 2 GS , Documents non autorisés , durée 2 heures
Session de rattrapage , juillet 2006

Partie I ( 10 points )

5- Donner le principe , les hypothèses , les paramètres et les procédures

pratiques d’utilisation de la méthode Socose ( 3 points )
6- Décrire brièvement les différentes études préliminaires nécessaires pour
la création et la construction d’un lac collinaire ( 4 points )
7- Donner et définir brièvement les différents éléments constitutifs d’un lac
collinaire ( 3 points )

128
Partie II ( 10 points ) : 2 points par question
11- citer quelques ouvrages utilisés dans la protection contre les
inondations en Tunisie
12- citer quelques formules ( sans donner l’expression ) utilisées en
Tunisie pour l’estimation de l’apport moyen
13- citer quelques formules (sans donner l’expression) utilisées en
Tunisie pour l’estimation de l’apport fréquentiel
14- citer quelques formules (sans donner l’expression) utilisées en
Tunisie pour l’estimation du débit fréquentiel
15- Comment utiliser la méthode Socose pour estimer des crues
centennales
Examen de Protection Contre les Inondations, 2 HAR+GS,
Durée 2 Heures , Documents non autorisés , Juin 2006

Donner les différentes étapes d’étude pour l’ implantation d’un lac collinaire sur un
cours d’eau drainant un bassin versant.
Adopter le schéma suivant : 1- Introduction Générale ( 1 point ) , 2- Etude
Géomorphologique ( 2 points) , 3- Etude Géologique et Géotechnique( 2 points ) ,
4- Etude Hydrologique ( 6 points) , 5- Dimensionnement de la digue ( 4 points ),
6- Dimensionnement de l’évacuateur des crues ( 4 points ) , 7- Conclusions et
recommandations ( 1 point )

Examen de protection contre les inondations , 2 HAR , 2006/2007

Le 19 mai 2007

On dispose de 32 ans de mesure de pluie maximales annuelles de 3 heures et on a

16 ans de mesure de débit maximum annuel au niveau de l’exutoire d’un bassin

versant de superficie 20 km2 et ayant un temps de concentration de 3 heures

Année Pluie (mm) Année Pluie (mm) Année Débits en m3/s

1 10 17 29 1 8,30
2 19 18 35 2 14,8
3 33 19 15 3 3,70
4 8 20 48 4 14,8

129
5 46 21 54 5 7,4
6 52 22 25 6 6,40
7 28 23 43 7 4,6
8 22 24 43 8 12
9 5 25 31 9 9,2
10 25 26 35 10 12,8
11 19 27 29 11 11,10
12 42 28 52 12 13,80
13 40 29 33 13 6,40
14 12 30 12 14 10,7
15 37 31 27 15 5,8
16 16 32 46 16 10,5

On cherche à déterminer les débits de période de retour 50 ans et 100 ans en utilisant
- L’ajustement direct de la loi de Gumbel sur la série des débits
- L’ajustement par la méthode de Gradex
1- Donner rapidement le principe , les hypothèses et les étapes de la méthode du
Gradex ( 3 points)
2- Ajuster la loi de Gumbel sur la série des débits et déterminer Qd10 et Qd100 ( 6
points )
3- Déterminer par la méthode du Gradex Qg10 et Qg100 ( 10 points)
4- Comparer les deux estimations et commenter ( 2 points)
Rappel : Les paramètres de la loi de Gumbel par la méthode des moments :
a 0.78 Sx X0  X 0.577 a

130