Délimitation Et Caractérisation D

Mini projet : Hydrogéologie
Délimitation et caractérisation d’un bassin versant sous ARCGIS

Oued Loukkos
Réalisé par : AIT BAALI OUIAME

ZINE EL ABIDINE Abderrahman
ALAOUI HICHAMI Hiba
Année universitaire: 2023/2024
Page 1 sur 38
14/12/2023
Table de Matière:
1-Introduction
2- Généralité sur le bassin.
3- Présentation des outils de travail.
4- Etapes préliminaires.
5-Préparation du MNT en ArcMap.
6-Délimitation du bassin.
7- Les caractéristiques du bassin.
8-Calcul des caractéristiques du bassin.
9-Caractéristiques du sol et de végétation.
10- Conclusion.
Page 2 sur 38
14/12/2023
Introduction :
L'analyse topographique du terrain, en particulier celle des bassins versants, constitue

une composante cruciale de l'hydrologie et de l'évaluation d'un site de projet pour un
ingénieur spécialisé dans l'hydraulique. La digitalisation des données a
considérablement simplifié cette tâche, avec des logiciels tels qu'ArcGIS et d'autres,
qui permettent d'extraire des paramètres essentiels pour la conception d'infrastructures
hydrologiques. Ces paramètres incluent le calcul du débit et de l'intensité des
écoulements superficiels, ainsi que la sélection optimale de la retenue du barrage et de
son type.
Le mini-projet d'hydrologie se présente comme une opportunité de se pencher les

étudiants en génie hydraulique aux outils des Systèmes d'Information Géographique
(SIG) à travers une approche d'auto-apprentissage. L'objectif est de mener une étude
exhaustive sur un cas concret, permettant aux futurs ingénieurs de développer leurs
compétences dans l'utilisation de ces logiciels pour l'analyse et la modélisation des
systèmes hydrauliques. Ce projet offre ainsi une expérience pratique visant à renforcer
la compréhension des méthodes hydrologiques appliquées dans le domaine de
l'ingénierie hydraulique.
Page 3 sur 38
14/12/2023
 Objectifs du mini-projet :
1. Acquisition de Compétences SIG
2. Analyse Hydrologique
3. Utilisation d'ArcGIS et Autres Logiciels
4. Interprétation des Résultats
5. Auto-Apprentissage
6. Dimensionnement des Ouvrages
Page 4 sur 38
14/12/2023
Généralité sur le bassin :
1- Le bassin versant de Moulouya :

Le bassin hydraulique du Loukkos est situé dans la partie nord-ouest du Maroc et
s'étend sur une superficie de 22 000 km²,
représentant environ 3% du territoire national.
Cette zone englobe principalement la région
administrative de Tanger-Tétouan-Al Hoceima,
comprenant les provinces de Tanger-Asilah,
Larache, et Chefchaouen.
Le bassin hydrologique du Loukkos couvre la
majeure partie de la zone d'étude, soit environ
18 000 km². Il englobe également des sous-
bassins tels que ceux de l'Oued Laou, Oued
Martil, et Oued Moulay Bousselham,
contribuant à la diversité hydrologique de la
région.
Cette région hydraulique est d'une grande
importance pour le Maroc, tant sur le plan environnemental qu'économique, en
abritant des écosystèmes variés et en soutenant des activités agricoles. Les
provinces incluses dans le bassin du Loukkos sont essentielles pour la gestion
durable des ressources en eau et la préservation de la biodiversité dans le contexte
de la planification hydraulique nationale.
2- Ressources en eau :
La zone du Loukkos s'étend sur environ 4 771 km2 et se divise en deux sous-bassins:
Dradère–Soueire-Mda : Ce sous bassin d'une superficie de 1 041 km2 regroupe les
bassins versants des oueds Dradère, Soueire et la partie amont de l'oued Mda.
Page 5 sur 38
14/12/2023
Loukkos : Ce sous bassin d une superficie totale de 3 730 km2 recouvre les bassins
hydrologiques des oueds Loukkos, Ouarour et Makhazine qui prennent naissance
dans la chaîne rifaine.
La situation sur la façade atlantique et au Nord-ouest du pays et la variabilité de
son relief confèrent à la zone du Loukkos un climat méditerranéen à influence
océanique. La pluviométrie moyenne interannuelle varie de 600 mm au niveau de
la bande côtière et la basse vallée du Loukkos à 1800 mm sur les hauts reliefs du
bassin du Loukkos.
Les principaux cours d'eau traversant la zone du Loukkos sont :
Oued Loukkos avec un apport moyen de 1200 Mm3/an et dont les principaux
affluents sont Ouarour (160 Km2) et El Makhazine (630 Km2)
Oued Dradère-Souiere avec un apport moyen de 260 Mm3/an
3- Oued Loukkos :
a- Caractéristiques :
- Longueur : 176 km
- Bassin : 3730 km²
- Bassin collecteur : Loukkos
- Débit moyen : 2221 µs/cm
b- Cours
- Source : Rif
- Localisation : au nord du Tanaqoub
- Altitude : 0 m
- Coordonnées : 35° 12′ 33″ N, 5° 28′ 06″ O
- Embouchure : Océan Atlantique
- Localisation : Larache
- Altitude : 0 m
- Coordonnées : 35° 12′ 15″ N, 6° 08′ 34″ O
Page 6 sur 38
14/12/2023
c- Définition :
Le Loukkos est un cours d'eau au Maroc qui prend sa source à la confluence des
massifs du Moyen et du Haut Atlas, dans la région d’Almssid, province de Midelt.
S'étirant sur une distance de 600 kilomètres, il se déverse dans la mer
Méditerranée, irriguant les régions du Rif et des Béni-Snassen, notamment les
plaines de Kebdana et Triffa, à l’extrême nord-ouest du Maroc. L'embouchure du
Loukkos se trouve à une distance de 14 kilomètres de la frontière algéro-marocaine.
 Mise en situation :
L'exutoire est localisé près de la région du Gharb-Chrarda-Beni Hssen, avec des
coordonnées Lambert :
- Coordonnée X : 459.8 km
- Coordonnée Y : 482.5 km
-Superficie du bassin versant (bv) : 3 730 km²
La contribution annuelle moyenne à l'embouchure du Loukkos est d'environ

920 millions de mètres cubes par an (pour la période de 1970 à 1998). Le
climat de la zone du bassin est de type méditerranéen, caractérisé par une
pluviométrie annuelle relativement basse et irrégulière, avec des cumuls
annuels variant de 200 à 400 mm.
Page 7 sur 38
14/12/2023
Présentation des outils de travail :

1- USGS :
L'USGS (United States Geological Survey) est une
qui joue un rôle essentiel dans la collecte de
données hydrologiques pour comprendre
et surveiller les ressources en eau. Dans le domaine de l'hydrologie et des

barrages, l'USGS fournit des informations cruciales sur les régimes
hydrologiques, les débits des cours d'eau et les caractéristiques des bassins
versants, contribuant ainsi à la planification, à la conception et à la gestion
efficace des projets de barrages pour assurer une utilisation durable des
ressources hydriques.
2- Arc Gis :
Permet l’extraction des données, le
Tracé de la courbe hypsométrique et
De la délimitation du bassin et la
Détermination des caractéristiques,
Avec l’aide de son extension
Arc Hydro.
3- Excel
Permet de tracer les courbes et faire les calculs
Page 8 sur 38
14/12/2023
Etapes préliminaires :
1-Transformation des cordonnées :
Il faut convertir les coordonnées Lamberts (X, Y) fournies pour cette étude, en
coordonnées du système WGS84 adopté par Google Earth, pour pouvoir
déterminer l’emplacement du bassin.
On utilise le site : GEOFREE.
X = 459.800 km
Y = 482.500 km
On entre les coordonnées (X, Y) via le site

https://tool-online.com/conversion-coordonnees.php
on obtient la latitude et la longitude en système de coordonnées Merchich/Long-Lat.
Page 9 sur 38
14/12/2023
2- Localisation de la Zone :
Après avoir convertir dans un premier temps les coordonnés Lambert en
coordonnées Merchich : (longitude et latitude du point)
On fait entrer ces coordonnés dans le site : https://earthexplorer.usgs.gov/
afin de générer le mnt :
o Le MNT :
Page 10 sur 38
14/12/2023
o Résultat :
4-préparation du MNT en ArcMap :

a- Définir le repère du travail :
Dans ce travail, le système choisi est celui des coordonnées Lamberts puisque
les coordonnées de l’exutoire sont données dans ce dernier et le système
géodésique adopté au Maroc est Merchich (projection Lambert Conforme).
Pour le bassin étudié, il est situé dans la zone 1 =Maroc-Nord.
Pour ce faire il faut définir un nouveau système de projection : changer le
système de coordonnées du MNT car la transformation WGS84 en
coordonnées Merchich n’est pas prédéfinie. Il faut la créer grâce à l’outil :
Data management tools , projections and transformations , create custom

geographic transformation.
Page 11 sur 38
14/12/2023
Dans le Logiciel Arc Gis on ouvre le MNT téléchargé sur USGS :
On a choisi la résolution du MNT de 90m :
Page 12 sur 38
14/12/2023
Page 13 sur 38
14/12/2023
b- Application de la transformation sur le MNT :
Il faut ensuite appliquer la transformation au projet à travers :

Projections and transformations > Raster >Project Raster.
Ensuite, il faut projeter Le raster sur le système de coordonnées Nord Maroc

(degré).
Page 14 sur 38
14/12/2023
c- Changement de l’environnement du travail :
Cette étape est très primordiale puisqu’on travaillait avec les coordonnées
WGS 84 et on compte par la suite de retourner à Lambert, ce qui nous
pousse à relancer un nouveau projet sur ArcMap et copier le nouveau
calque transformé.
Page 15 sur 38
14/12/2023
Page 16 sur 38
14/12/2023
Délimitation du bassin :
Par la suite, il faut travailler avec l’outil : ArcHydro.
Arc Hydro est un système ArcGIS adapté pour prendre en charge les
applications de ressources hydrologiques. Afin de délimiter le bassin
versant, on exécute les commandes suivantes dans l’ordre.
1- Correction des erreurs :
Pour bien exécuter les autres commandes, cette étape est très importante
pour corriger les erreurs d’altitude aberrantes ou extrêmes qui risquant
d’influencer la détermination du sens d’écoulement : cuvettes ou crêtes, et
ce par une méthode d’interpolation.
Terrain Processing > DEM Manipulation > Fill Sinks.
Page 17 sur 38
14/12/2023
Résultat :
2- Orientation du sens de l’écoulement :
Afin de déterminer le sens de l’écoulement dans chaque cellule élémentaire de la carte

on procède par l’outil flow direction :
Terrain Processing > Flow Direction.

Page 18 sur 38
14/12/2023
- Résultat :
3- Détermination du réseau Hydraulique :
Page 19 sur 38
14/12/2023
Cette commande utilise le sens de l’écoulement déterminé précédemment pour définir le réseau
hydraulique de la carte :
Terrain Processing > Flow Accumulation.
- Résultat :
Page 20 sur 38
14/12/2023
4- Affinement du réseau Hydrique :
Cette étape permet de distinguer entre les cours d’eau principaux et les cours d’eau
secondaires. Ceci est réalisé en utilisant deux méthodes :
Terrain Processing > Stream Definition.
Terrain Processing > Stream Segmentation.

Ces deux méthodes servent respectivement à affiner le réseau hydraulique précèdent
et puis segmenter le réseau d’écoulement en tronçons à identifiants uniques.
4- Définition et Transformation des surfaces de drainage :
Cette étape permet d’associer à chaque segment d’écoulement définit

précédemment sa surface de drainage en lui donnant un identifiant. Ensuite on
transforme ces surfaces de drainages en polygones. Ceci est réalisé en exécutant
ces deux commandes :
Terrain Processing > Catchment Grid Delineation.
Terrain Processing > Catchment Polygon Processing.
Page 21 sur 38
14/12/2023
Page 22 sur 38
14/12/2023
- Résultat :
Page 23 sur 38
14/12/2023
6- Détermination des lignes de drainage :
Dans cette étape, on convertit les cours d’eau de la représentation en grille en
représentation vectorielle. Ensuite, on définit pour chaque surface, définit un
polygone de la surface drainée en amont.
Ceci est réalisé en exécutant ces deux commandes :
Terrain Processing > Drainage line Processing.
- Résultat :
Page 24 sur 38
14/12/2023
Terrain Processing > Adjoint Catchment Processing.
7- Positionnement de l’exutoire :
Dans cette étape, on cherche le positionnement de nos exutoires à l’aide du point
donné. Tout d’abord, on précise le point donné en coordonnées Lambert, puis on
cherche le cours d’eau le plus proche à ce point, la projection de ce dernier sur ce
cours d’eau constitue l’exutoire.
On utilise l’outil : Batch Point et l’option Go to XY.
Page 25 sur 38
14/12/2023
8- Délimitation du bassin :
On délimite le bassin par rapport au exutoire définit.
On utilise la commande :
Watershed Processing > Batch Watershed Processing.
Page 26 sur 38
14/12/2023
Résultat:
Page 27 sur 38
14/12/2023
7- Les caractéristiques du Bassin :

Une fois le bassin et bien délimité, on commence par les calculs.
1- Surface et périmètre du bassin :
L’outil Drainage Area Characterization qui permet de calculer les

caractéristiques physiographiques du bassin délimité à savoir : la surface, le
périmètre, ainsi que les altitudes minimales et maximales.
On utilise la commande :
Terrain Morphology > Drainage Area Characterization.
Ces résultats nous permet ainsi le traçage de la courbe hypsométrique en utilisant le

tableau hypsométrique de la couche DrainEAV, et exécutant ensuite la commande
Create Graph pour les deux colonnes Top Elev en vertical et Cum Area en horizontal.
List By source > DrainEAV > Open Table
Page 28 sur 38
14/12/2023
Page 29 sur 38
14/12/2023
Page 30 sur 38
14/12/2023
3- Cours d’eau principal :
Pour calculer la longueur de plus long chemin on procède par la méthode Longest
Flow Path :
Watershed Processing > Longest Flow Path
On obtient la longueur du cours d’eau à partir du tableau :
Layer Longest Flow Path > Open Attribute
Page 31 sur 38
14/12/2023
8- Calcul des caractéristiques du Bassin :
D’après les données produites par l’outil ArcGIS, nous pouvons

maintenant calculer les paramètres qui caractérise notre
bassin et influencent son temps et la forme de sa réponse qui
sont à la fois :
●  Les caractéristiques physiques (forme, relief, pente).

●  Les caractéristiques du réseau de drainage.
●  Les caractéristiques du sol et de végétation.
Les caractéristiques physiques :
1- Les caractéristiques de forme :
On détermine maintenant les deux indices :
Compacité de Gravelius et indice de forme de Horton.

P P
KG= =0.28
2 √ πA √A
=1.8
Indice de forme de Horton :
A
K H= 2
L
=0.29314
La forme du bassin :
Page 32 sur 38
14/12/2023
On a 1.5< K G < 1.8 et K H < 1 , donc le bassin est allongé. (Ce qui est
conforme avec la délimitation).
1- Les caractéristiques du relief :
On commence par le tracé de : la courbe hypsométrique :

-Nous avons essayé de la tracer tous d’abord en utilisant les
données produites par ArcGis.
- Ensuite, en utilisant les données fournis en Mini-Projets.
Courbe Hypsométrique (Données ArcGIS) :
On utilise ces deux colonnes pour le tracé :
Page 33 sur 38
14/12/2023
On obtient la courbe Hypsométrique :
Page 34 sur 38
14/12/2023
Après qu’on corrige 2 points on obtient:
L’altitude médiane : point d’abscisse 50% sur la courbe Hypsométrique :
H 50 %=1486.4 m (Courbe Hypsométrique).
L’altitude minimale : point de contrôle se situe à l’exutoire du bassin :
H min =1449 m
L’altitude maximale : point cumulant du bassin, l’altitude la plus élevée :
H max =1544.94 m
1
❑
h i+ hi+1
L’altitude moyenne : h moy= =∑ ❑ Si =¿1494,30m
S ❑ 2
Page 35 sur 38
14/12/2023
1469.2+ 1479
Le mode : M= 2
=1474.1m
2-Rectangle équivalent :
Il s’agit d’une transformation purement géométrique en vertu de

laquelle on assimile le bassin à un rectangle ayant le même périmètre et la
même superficie du bassin.
3-Les indices de pente :
Page 36 sur 38
14/12/2023
Leur connaissance est d’une grande importance car il est
évident que les eaux ruissellent d’autant plus que la pente des
versants est grande.
La pente moyenne du bassin :
Indice de pente de roche :
Indice du pente globale :
Du dénivelée utile : Du=H5%-H95%=1519.7-1463.4=56.3

Donc
Indice de pente classique :
Dénivelée spécifique :
Pente moyenne de cours d’eau principal : 7.6391
Page 37 sur 38
14/12/2023
Temps de concentrations avec les différentes méthodes:
Conclusion:
Au terme de ce travail réalisé on a pu mettre à l’épreuve notre bagage théorique acquis
Durant les cours de l’hydrologie et de voir une application concrète des lois et des
Expressions du cours. Ainsi, l’outil d’information géographique ArcGIS et son extension
ArcHydro nous ont permis, loin des méthodes numériques, de délimiter le bassin et de
Calculer leur caractéristiques physiques et géométriques.
Page 38 sur 38
14/12/2023

Délimitation Et Caractérisation D

Transféré par

Droits d'auteur :

Formats disponibles

Délimitation Et Caractérisation D

Transféré par

Informations du document

Titre original

Copyright

Formats disponibles

Partager ce document

Partager ou intégrer le document

Options de partage

Avez-vous trouvé ce document utile ?

Ce contenu est-il inapproprié ?

Droits d'auteur :

Formats disponibles

Délimitation Et Caractérisation D

Transféré par

Droits d'auteur :

Formats disponibles

Mini projet : Hydrogéologie

Délimitation et caractérisation d’un bassin versant sous ARCGIS

Réalisé par : AIT BAALI OUIAME

Année universitaire: 2023/2024

2- Généralité sur le bassin.

3- Présentation des outils de travail.

5-Préparation du MNT en ArcMap.

7- Les caractéristiques du bassin.

8-Calcul des caractéristiques du bassin.

9-Caractéristiques du sol et de végétation.

L'analyse topographique du terrain, en particulier celle des bassins versants, constitue

Le mini-projet d'hydrologie se présente comme une opportunité de se pencher les

1. Acquisition de Compétences SIG

3. Utilisation d'ArcGIS et Autres Logiciels

4. Interprétation des Résultats

6. Dimensionnement des Ouvrages

Généralité sur le bassin :

1- Le bassin versant de Moulouya :

La contribution annuelle moyenne à l'embouchure du Loukkos est d'environ

Présentation des outils de travail :

qui joue un rôle essentiel dans la collecte de

données hydrologiques pour comprendre

et surveiller les ressources en eau. Dans le domaine de l'hydrologie et des

On entre les coordonnées (X, Y) via le site

4-préparation du MNT en ArcMap :

Data management tools , projections and transformations , create custom

Dans le Logiciel Arc Gis on ouvre le MNT téléchargé sur USGS :

On a choisi la résolution du MNT de 90m :

Il faut ensuite appliquer la transformation au projet à travers :

Ensuite, il faut projeter Le raster sur le système de coordonnées Nord Maroc

c- Changement de l’environnement du travail :

1- Correction des erreurs :

Terrain Processing > DEM Manipulation > Fill Sinks.

2- Orientation du sens de l’écoulement :

Afin de déterminer le sens de l’écoulement dans chaque cellule élémentaire de la carte

Terrain Processing > Flow Direction.

3- Détermination du réseau Hydraulique :

Terrain Processing > Flow Accumulation.

Terrain Processing > Stream Definition.

Terrain Processing > Stream Segmentation.

Cette étape permet d’associer à chaque segment d’écoulement définit

Terrain Processing > Catchment Grid Delineation.

Terrain Processing > Catchment Polygon Processing.

Terrain Processing > Drainage line Processing.

7- Les caractéristiques du Bassin :

1- Surface et périmètre du bassin :

L’outil Drainage Area Characterization qui permet de calculer les

Terrain Morphology > Drainage Area Characterization.

Ces résultats nous permet ainsi le traçage de la courbe hypsométrique en utilisant le

Watershed Processing > Longest Flow Path

On obtient la longueur du cours d’eau à partir du tableau :

Layer Longest Flow Path > Open Attribute

D’après les données produites par l’outil ArcGIS, nous pouvons

●  Les caractéristiques physiques (forme, relief, pente).

Les caractéristiques physiques :

1- Les caractéristiques de forme :

On détermine maintenant les deux indices :

Compacité de Gravelius et indice de forme de Horton.

1- Les caractéristiques du relief :