1.

CARACTERISTIQUES PHYSIQUES DES BASSINS VERSANTS

1.1
Dlimitation et superficie des bassins versants tudis
La dlimitation du bassin versant est effectue sur la base des cartes topographiques au 1 /25000me et
1/50000me
Les lignes de partage des eaux sont dtermines compte tenu aussi bien de la topographie que des
limites imposes par lurbanisation (route, construction, talus), le systme de collecte des eaux pluviales
et le sens de lcoulement.
Ltude des diffrentes caractristiques physiques des bassins versants devra nous fournir les lments
ncessaires une bonne apprciation du ruissellement tant sous son aspect spatial (coefficient de
ruissellement) que sous son aspect temporel (temps de concentration).
1.2
Pente moyenne
La pente moyenne dun bassin versant est assimile celle de son cours deau principal. Elle est
dtermine par la formule suivante :

LK

I L/

IK

Avec,
L : longueur totale du cours deau en km, L= LK
LK : longueur des diffrents tronons du cours deau, en km
lk : Pente de chaque tronon du cours deau.
1.3
Coefficient de compacit
Pour caractriser la forme d'un bassin versant, on utilise le coefficient de compacit Kc, dtermin par la
formule suivante : Kc = 0,28 P /

Avec
P : Primtre du bassin en km
A : Superficie en km, du bassin versant

2.

EVALUATION DES DEBITS DE CRUES

Les dbits de crue des bassins versants seront valus selon la taille des bassins versants et sur la
base dune analyse comparative des rsultats obtenus par diffrentes mthodes applicables pour la zone
du projet.
2.1 Choix de la priode de retour
Les dgts conomiques causs par une interruption de trafic routier ont toujours prsent des valeurs
considrables. Ils constituent probablement le poste le plus coteux dans les dommages dus aux
inondations aprs les pertes de vies humaines.
Aussi, un drainage efficace, cohrent et rationnel de la route implique un dimensionnement judicieux des
ouvrages de drainage. Ce qui appelle une plus grande vigilance quant au choix de la priode de retour
de l'vnement pluvieux. Ce choix dpend de l'importance de l'coulement intercept par l'ouvrage, du
type d'ouvrage projet et du niveau d'amnagement adopt pour la route tudie.
Contrairement ce que l'on peut penser, en effet, une averse de priode de retour donne donne trs
rarement un dbit de la mme priode de retour.
La plupart des mthodes de prdtermination des dbits sont d'ailleurs bases sur le concept d'une
priode de retour du dbit gale celle de l'averse, seul le bon sens de l'hydrologue ou de l'ingnieur et
son analyse fine des divers aspects contribuant la formation du ruissellement permettent de mesurer
jusqu' quel degr on peut admettre cette hypothse dans les calculs.
2.2

Mthodes de calculs des dbits

2.2.1 Formules empiriques (S<100 km2)

A- la mthode rationnelle :
Q(T) = Kr Ka I(tc) S/3,6 (m3/s)
Avec
S : surface du bassin versant en km2
I(T) : intensit de la pluie en mm/heure
Kr : coefficient de ruissellement.
Ka : Coefficient d'abattement de la pluie.
Temps de concentration tc :
Le temps de concentration est calcul par la formule Empirique de KIRPICH :
tc =

1
L1,15
0,38
52 H

Avec :
L : Distance en mtre entre l'exutoire et le point le plus loign du bassin (m)
H : dnivele en mtre entre l'exutoire et le point le plus loign du bassin (m)

S en Km

Tableau donnant les valeurs de Ka

S en km2
25 >
25-50
25 Ka
> 50 25 1 100 50 0,95
150 100

Ka

1
S en Km 25 >
Ka

0.95
50 25
0.95

0.9
100 50
0.9

0.85
150 100
0.85

50-100
150
0,9

100-150
0,85
200

200

250

0.8
150

0.75
200

200

250

0.8

0.75

150-200
0,8

200-250
0,75

Coefficient
de
S en Km
ruissellement
Le coefficient de ruissellement diminue
lorsque la superficie du bassin versant
augmente. D'autre part, il varie dans le
Ka
mme sens que la pente du bassin. Un
S en Km
terrain permable, une couverture
vgtale
dense
abaissent
considrablement le coefficient de
Ka
ruissellement.
On pourra dterminer les coefficients
de ruissellement partir du tableau suivant :

1
25 >
1

50

25

0.95
50 25

50

100

0.9
100 50

0.95

150

100

0.85
150 100

0.9

150

Plates-formes et chausss
Terrains dnuds ou
vgtation non couvrante
Terrains attaqus par l'rosion
Labours frais
Cultures couvrantes, crale
haute
Terrains de parcours
Prairies
Brousse dense, savane
sous bois
Fort ordinaire
Sous bois touffus
Grande fort primaire

<5%
0.95

0<S<10 ha
5<i<10% 10<i<30%
-

>30%
-

<5%
-

200

200

250

0.8
150

0.75
200

200

250

0.8

0.75

0.85

Valeurs de Kr

Nature de couverture
vgtale

25 >

10<S<400 ha
5<i<10% 10<i<30%
-

>30%
-

0.80

0.85

0.90

0.95

0.70

0.75

0.80

0.85

0.75

0.80

0.85

0.90

0.52

0.60

0.72

0.80

0.70

0.75

0.80

0.85

0.30

0.36

0.12

0.50

0.30

0.50

0.60

0.70

0.13

0.20

0.25

0.30

0.20

0.25

0.30

0.40

0.15

0.18

0.22

0.25

Intensit de la pluie :

Lintensit de la pluie est dtermine laide des courbes intensit-dure-frquence (IDF). Elle est exprime par la
relation suivante :

I = a tc b
tc : temps de concentration en heures
I : intensit moyenne est en mm/h

B- Mthode SOGREAH (SPEED)

La mthode SPEED a t dfinie et applique dans le cadre du projet national de protection contre les inondations.
Elle se base sur la formule suivante :
0 , 75
Q (T) = S P(T ) P0 / 12 , (m3/s)
Avec :
QT: Dbit de pointe de crue de priode de retour T, en m 3/s
PT : Pluie journalire de priode de retour T (en mm)
P0 : Seuil de ruissellement (mm),
S : Superficie du bassin versant tudi (km2)
La mthode SPEED ne tient pas compte des caractristiques morphologiques : pente et compacit des bassins
versants ; sont utilisation pour des B.V de trs petite taille ou reprsentant une morphologie particulire (pente trs
faible) devrait tre effectue avec prudence.

Cette formule ne peut tre utilise que dans la zone ayant fait lobjet de ltude du projet national de protection
contre les inondations et par utilisation des cartes de pluies tablies dans le cadre de cette mthode .
i)
Les pluies
L'tude probabiliste des pluies journalires conduit la formulation rgionale des pluies de temps de retour T qui
suit une loi de Gumbel et qui est donc reprsente par une droite dans le plan Y, P
Variable de Gumbel

YT = - Log Log (1 )
T

Pluie journalire

PT = a [YT - Y0 ]

Cette droite est bien dfinie si l'on fournit les deux valeurs :
P10 = pluie dcennale journalire, correspondant Y10 = 2,25
P100 = pluie centennale journalire, correspondant Y100 = 4,60
Ce qui permet de calculer :
PT = P10 +

P100 P10
(YT Y10) pour T < 100 ans
Y100 Y10

ii)
Les dbits de pointe de crue
Les dbits de pointe de crue sont calculs par la formule :
QT =

S 0, 75
(PT P0)
12

Avec :
S

Superficie du bassin versant (km2)

PT

Pluie journalire de priode de retour T, en (mm)

Po

Seuil de ruissellement (mm)

Ces dbits suivent donc, comme les pluies journalires, deux lois de Gumbel, l'une pour les temps de retour jusqu'
100 ans, l'autre pour les temps de retour au-del de 100 ans.
On dfinit trois dbits caractristiques :
Q10 =A S0,75

avec A =

Q100 =B S0,75

P10 P 0
12
avec B =

P100 P 0
12

Pour les autres temps de retour, on a :

QT = Q10 +

Q100 Q10
Y100 Y10

(YT Y10)

iii)
Les hydrogrammes de crue
Le volume de la crue de temps de retour T est fonction du dbit de pointe QT(m3/s) de mme temps de retour et du
temps de base TB (heures) suppos indpendant du temps de retour (thorie de l'hydrogramme unitaire) :

VT =

1,8 QT TB
1000

Cette formule est celle de la crue "triangulaire quivalente", c'est--dire d'un triangle ayant mme dbit de pointe et
mme volume que la crue relle. Elle sert la fois dfinir TB sur les mesures relles et calculer VT en fonction de
QT.
Le temps de base est fonction de la superficie S (km2) du bassin versant : TB = S0,385
Par ailleurs, la valeur de est la mme pour toute la zone d'tude ( = 1,5).
Le temps de monte de la crue est pris gal Tm = 0,40 TB.
iv)
Mise en uvre
Pour pouvoir appliquer les formules ci-dessus un bassin versant particulier, situ dans une zone donne, il faut
fournir les valeurs utilisables dans cette zone :
soit directement les valeurs A, B et C par ajustement statistique des crues observes, si l'on dispose d'une station
de jaugeage avec une srie d'observations suffisamment longue,
soit indirectement partir des donnes pluviomtriques de P10 et P100, utilisables dans cette zone.
En labsence de donnes hydromtriques au niveau des ouvrages existants, les calculs hydrologiques sont effectus
partir des donnes pluviomtriques. La synthse des rsultats de l'tude pluviomtrique mene dans le cadre du
Projet National de Protection contre les Inondations est matrialise sous la forme de trois cartes.
la carte rgionalise des pluies journalires dcennales, P10,
la carte rgionalise des pluies journalires centennales, P100,
la carte des seuils de ruissellement, P0.
Il faut souligner que les rsultats obtenus sont valables grande chelle et ne dispensent pas d'tudes particulires
locales destines affiner les valeurs utiliser pour tel ou tel projet. En effet, la Tunisie centrale est particulirement
complexe sur le plan de l'hydrologie :
les pluies journalires varient considrablement suivant la position du pluviomtre par rapport aux reliefs :
l'utilisation des sries pluviomtriques locales, mmes courtes, est indispensable pour l'tude fine d'un petit
bassin,
le comportement des oueds peut galement varier beaucoup suivant la nature du bassin versant (paisseur de la
couche permable) et surtout suivant les pertes dans le lit plus ou moins permable. Cela signifie que P0
peut prendre localement des valeurs trs loignes de la valeur rgionale cartographie.

2.2.2 Formules utilisant uniquement la superficie du bassin versant (S>100 km2)

A- Mthode de Franco-Rodier:
Cest un moyen assez commode pour lvaluation des crues exceptionnelles pour les grands basins versants.
Franco et Rodier ont tudi les crues exceptionnelles dans le monde et se sont aperu que, pour une mme
rgion (caractrise par un coefficient k). Les points reprsentant le dbit, Q, en fonction de la surface dun
basin versant, S, reports sur un graphique bi logarithmique saligne et les diffrentes droites convergent vers
le point (So, Qo).
k

Q Q0 (

S 110
)
S0

Avec :
So: superficie de rfrence gale 10 8 km2
Qo: dbit de rfrence gale 10 6 m3 / s
K : coefficient rgional, fonction de la priode de retour T.
Cette mthode a t extrapole pour le calcul des crues de diffrentes priodes de retour, en donnant une
signification frquentielle au coefficient K. ce coefficient est dtermin partir des rsultats des observations
hydromtriques.

T(ans)

10

20

T(ans)

K
T(ans)

3.80
10

3.85
20

k 3.98
50

10
100
3.80
4.15
100

3.80

3.85

3.98

4.15

50

20
3.85

B- Mthode Ghorbel
A partir d'une analyse statistique des dbits maximas

T(ans)
50
k
3.98
T(ans)

10 100 20
3.80
3.85
10 4.15 20

50

100

3.98
50

4.15
100

k
T(ans)

3.80
10

3.85
20

3.98
50

4.15
100

K
T(ans)

3.80
10

3.85
20

3.98
50

4.15
100

3.80

3.85

3.98

4.15

observs dans la rgion, l'auteur a tabli la relation rgionale donnant le dbit de crue de retour T et
qui s'crit comme suit
Q (T) = R (T) Qmax(moy)

(m3 /s)

Avec :
Qmax(moy) = dbit maximum moyen en m3/s, exprim par la relation suivante :
Qmax(moy) = S 0,8 [1.075 (

P H
L

I c ) - 0.232]

Avec
P : Pluviomtrie moyenne annuelle sur le bassin en m ;
Ic : Coefficient de compacit ;
H : Diffrence entre l'altitude Mdiane et l'altitude l'exutoire en m ;
L : Longueur du cours d'eau principal, en km ;
R (T) : Paramtre rgional.

Cette quation est valable pour les bassins des oueds qui ont une pente infrieure 5,5%.
Pour les bassins versants de la dorsale fortes pentes (suprieures 5,5%) :
Q max(moy) = 2,86 S0,8

T (ans)
Zone 1
Zone 2
Zone 3
Zone 4
Zone 5

2
0.86
0.70
0.59
0.5
0.3

5
1.39
1.33
1.45
1.6
1.0

10
1.79
1.98
2.34
2.5
2.2

20
2.19
2.84
3.52
3.5
3.7

50
2.72
4.40
5.68
5.1
6.7

100
3.12
6.04
7.93
6.2
9.2

Zone 1 : LIchkeul, lextrme Nord et les affluents rive gauche de le Medjerda

Zone 2 : La Medjerda avec ses affluents rive droites, le Cap Bon et le Zroude Khanget Zazia
Zone 3 : les oueds Mliane, Merguelil et la branche Nord du Zruod
Zone 4 : le sahel et la rgion de sfax
Zone 5 : le Sud

C- Mthode FRIGUI :
Les formules de Frigui sont des formules adaptes aux diffrentes rgions tunisiennes :
Pour la Tunisie du Nord

,2
QT T ( A26
A
1) 0.47
,5
QT T ( A53
A
1) 0 , 53

Pour Medjerda
Pour le Cap Bon et Mliane
Pour le centre et le Sud

,4
QT T ( A38
A
1) 0 , 44
,7
QT T ( A76
A
1) 0.44

Avec

QT Dbit de crue de la priode de retour T en m3/s

T Coefficient dajustement dpendant de la rgion et de la priode de retour

A superficie du basin versant en km

Le coefficient de correction a t tabli comme suit suivant la priode de retour T

Rgion

100

Nord
Medjerda
Cap Bon et Mliane
Centre et Sud

1
1
1
1

D- Mthode de Kallel :

50

0,80
0,78
0,77
0,74

20

10

0,58
0,54
0,50
0,48

0,45
0,38
0,35
0,33

A partir dun inventaire des dbits spcifiques maximums, lauteur a tabli des quations rgionales donnant le
dbit maximal (Qt) de priode de retour (T) pour un bassin versant de superficie (S).
Pour le noyau de la Tunis septentrionale :

q( T ) q0 .S .T .S
O :
q(T) : Dbit spcifique, en m3/s/km2 ;
S : Superficie du bassin versant en km ;
T : Priode de retour en ans.
q0 , et : Constantes rgionales.

Le Dbit maximum de priode de retour T :

Q = q(T) . S

Pour la Tunisie du Nord

Q( T ) 5,5 S .T 0.41

Pour le noyau de la Dorsale tunisienne

Q( T ) 2,6.S 0, 3 . S .T 0.41

Pour la Tunisie centrale

T=10, 20 ans

Q( T ) 12,3. S .T 0.41

T=50, 100 ans

Q( T ) 24,7. S .T 0.41

Pour la Tunisie mridionale

Q( T ) 12.35. S .T 0.41

E- Mthode Fersi
A partir des sries tendues et homognises, des dbits maximas annuels, des bassins versants controls
et de lame coule moyenne, l'auteur a pu dterminer pour le Sud Tunisien la relation suivante :
Qmax (moy) = 0,266 Hr S0.5 log S, en m3/s,
Avec :
S : surface du bassin versant en km
Qmax (moy): dbit maximum annuel moyen, en m3/s
Hr : lame coule moyenne, en mm
Hr = 16,4 P Ig , en mm, avec
P : pluviomtrie moyenne interannuelle sur le bassin versant, en m
Ig : indice de pente global, en m/km.
Q(T) = Y(T) Qmax (moy), en m3/s,
Avec :
Y(T) : coefficient rgional, pour le sud tunisien et une priode de retour centennale, Y (100) = 9,2 et Y (50) = 6,7.
Q(T) : le dbit maximum de crue de retour T en ans, en m3 / s.