Calcul D Un Bassin Versant Qgis

Cette page vous permet d’estimer rapidement les paramètres hydrologiques clés d’un bassin versant extrait dans QGIS à partir d’un MNT, d’un exutoire et d’outils de traitement raster. Le calculateur ci-dessous synthétise l’aire, la pente moyenne simplifiée, l’indice de compacité de Gravelius, le temps de concentration et un débit de pointe théorique selon la méthode rationnelle.

Calculateur hydrologique

Guide expert : comment réaliser le calcul d’un bassin versant dans QGIS

Le calcul d’un bassin versant dans QGIS consiste à délimiter l’aire qui contribue à un exutoire donné, puis à dériver des indicateurs morphométriques et hydrologiques utiles pour l’analyse des écoulements. Dans un projet réel, cette opération sert aussi bien à la gestion des inondations qu’au dimensionnement d’ouvrages hydrauliques, à l’étude d’un sous-bassin agricole, à l’aménagement urbain ou à l’expertise environnementale. QGIS, combiné aux algorithmes GDAL, SAGA ou GRASS, fournit un environnement robuste pour travailler à partir d’un modèle numérique de terrain, repérer les lignes de flux et calculer des paramètres quantifiables.

Le principe général est simple : l’eau s’écoule de l’amont vers l’aval en suivant les pentes du relief. Un bassin versant est donc la surface topographique qui alimente un point de sortie. Une fois ce polygone extrait, l’analyste peut calculer son aire, son périmètre, la longueur du cours principal, sa dénivelée, sa pente moyenne et divers indices de forme. Ces données alimentent ensuite des formules d’estimation du temps de concentration, du débit de pointe ou de la sensibilité au ruissellement. Le calculateur de cette page reprend cette logique de travail et propose un résultat immédiat à partir de valeurs déjà dérivées dans QGIS.

Pourquoi utiliser QGIS pour la délimitation d’un bassin versant

QGIS est particulièrement adapté à cette tâche pour plusieurs raisons. D’abord, il permet de gérer facilement des rasters de terrain, des couches vectorielles d’hydrologie, des points d’exutoire et des données d’occupation du sol. Ensuite, les boîtes à outils de traitement intègrent des workflows complets pour corriger les dépressions du MNT, calculer la direction d’écoulement, l’accumulation de flux, puis délimiter automatiquement le bassin versant. Enfin, QGIS est transparent sur les unités, les projections et les attributs, ce qui réduit les risques d’erreur lorsque l’on passe de la géomorphologie à l’hydrologie appliquée.

Point essentiel : la qualité du résultat dépend d’abord de la qualité du MNT, du bon choix de projection métrique et du positionnement précis de l’exutoire. Une erreur de quelques pixels au point d’exutoire peut modifier notablement l’aire contributive calculée.

Données nécessaires avant le calcul

• Un modèle numérique de terrain cohérent, sans trous majeurs ni artefacts significatifs.
• Une projection en mètres, typiquement un système UTM ou un Lambert adapté à la zone d’étude.
• Le point d’exutoire, idéalement calé sur le réseau hydrographique ou sur la cellule de plus forte accumulation.
• Des données d’occupation du sol si vous souhaitez affecter un coefficient de ruissellement.
• Des données pluviométriques IDF pour estimer une intensité de pluie pertinente.

Étapes de traitement dans QGIS

1. Importer le MNT et vérifier sa projection. Une projection géographique en degrés n’est pas adaptée aux mesures de longueur et de surface.
2. Corriger les dépressions à l’aide d’un outil de remplissage des cuvettes. Cette étape évite les blocages artificiels de l’écoulement.
3. Calculer la direction d’écoulement pour connaître, pixel par pixel, la direction de transfert vers l’aval.
4. Calculer l’accumulation de flux afin d’identifier les chenaux potentiels et aider au positionnement de l’exutoire.
5. Snapper l’exutoire sur la cellule d’accumulation appropriée si nécessaire, afin d’éviter les décalages.
6. Délimiter le bassin versant à partir du point d’exutoire et de la grille de direction d’écoulement.
7. Vectoriser le raster du bassin puis calculer les champs géométriques : aire, périmètre, altitudes, longueur du cours principal.
8. Extraire les statistiques raster sur le bassin pour obtenir les altitudes min, max, moyenne et éventuellement les pentes.

Comprendre les indicateurs calculés

L’aire du bassin est la mesure la plus fondamentale. Elle détermine le volume potentiel d’eau mobilisable lors d’un épisode pluvieux. Plus le bassin est grand, plus le débit total susceptible d’atteindre l’exutoire peut être élevé, toutes choses égales par ailleurs.

Le périmètre intervient notamment dans le calcul de l’indice de compacité de Gravelius. Cet indice compare la forme du bassin à celle d’un cercle de même aire. Une valeur proche de 1 traduit un bassin compact, généralement plus réactif, tandis qu’une valeur plus élevée indique un bassin plus allongé, où les écoulements tendent à arriver de manière plus étalée dans le temps.

La pente moyenne simplifiée peut être approchée par la dénivelée entre l’amont et l’exutoire divisée par la longueur du talweg principal. Cette pente influence fortement la vitesse de ruissellement et donc le temps de concentration.

Le temps de concentration représente le temps nécessaire pour que l’eau provenant du point hydraulique le plus éloigné atteigne l’exutoire. Dans les approches simplifiées, il peut être estimé avec des formules empiriques comme celle de Kirpich pour les petits bassins. C’est une grandeur essentielle, car l’intensité de pluie à retenir dans la méthode rationnelle doit être associée à une durée proche de ce temps.

Le débit de pointe est souvent approximé par la méthode rationnelle, surtout pour les petits bassins urbains ou ruraux. La formule Q = 0,278 × C × i × A exprime le débit en m³/s lorsque A est en km², i en mm/h et C sans unité. Ce calcul reste une estimation initiale et non un modèle hydrologique complet.

Formules utilisées dans ce calculateur

• Pente moyenne du talweg (%) = ((Altitude max – Altitude exutoire) / Longueur du talweg en mètres) × 100
• Indice de Gravelius = P / (2 × √(π × A)) avec P en km et A en km²
• Temps de concentration de Kirpich (minutes) = 0,01947 × L^0,77 × S^-0,385 avec L en mètres et S pente décimale
• Débit de pointe rationnel (m³/s) = 0,278 × C × i × A

Valeurs de référence pour le coefficient de ruissellement

Le coefficient C varie selon l’usage du sol, le taux d’imperméabilisation, l’état de saturation et la pente. Les valeurs ci-dessous sont indicatives pour un calcul préliminaire. Dans un projet réglementaire, il faut s’appuyer sur des références locales, des guides nationaux ou des études antérieures comparables.

Occupation du sol dominante	Coefficient C typique	Réactivité hydrologique attendue
Forêt dense	0,10 à 0,30	Ruissellement limité, infiltration élevée, réponse plus amortie.
Agricole mixte	0,30 à 0,50	Réponse modérée, variable selon les pratiques culturales et la saison.
Prairie, sols nus, terrains dégradés	0,40 à 0,60	Réponse rapide si sols compactés ou érodés.
Périurbain	0,50 à 0,75	Ruissellement plus marqué, réseau de drainage souvent plus efficace.
Urbain dense	0,70 à 0,95	Très forte réactivité, pics de débit rapides et élevés.

Exemple de statistiques morphométriques selon la taille du bassin

Les valeurs suivantes sont des ordres de grandeur utiles pour situer un bassin versant obtenu dans QGIS. Elles ne remplacent pas l’analyse de terrain mais aident à interpréter la réactivité hydrologique probable.

Catégorie de bassin	Surface typique	Temps de concentration souvent observé	Usage principal de l’analyse
Micro-bassin urbain	< 1 km²	5 à 30 min	Gestion pluviale, ouvrages de rétention, dimensionnement local.
Petit bassin rural	1 à 20 km²	20 min à 3 h	Hydrologie d’événement, érosion, franchissements routiers.
Bassin intermédiaire	20 à 200 km²	2 à 12 h	Études d’inondation locales, hiérarchisation des affluents.
Grand bassin	> 200 km²	Plusieurs heures à plusieurs jours	Gestion intégrée, modélisation hydrologique distribuée.

Bonnes pratiques pour fiabiliser le calcul d’un bassin versant dans QGIS

• Travaillez toujours dans un système projeté métrique avant de calculer surfaces et longueurs.
• Contrôlez visuellement la cohérence du réseau extrait avec l’hydrographie connue.
• Évitez de placer l’exutoire à côté du chenal principal : utilisez le snapping sur accumulation maximale.
• Vérifiez la résolution du MNT. Un pixel trop grossier peut lisser les vallons et fausser la délimitation.
• Pour les zones plates, croisez le résultat avec des orthophotos, des levés locaux ou une couche hydrographique de référence.
• Ne confondez pas pente moyenne du bassin et pente du talweg principal, qui n’ont pas le même usage.
• Si le bassin est karstique, fortement anthropisé ou régulé, l’interprétation purement topographique peut être insuffisante.

Limites de la méthode rationnelle et du temps de concentration simplifié

La méthode rationnelle est utile pour un premier ordre de grandeur mais elle repose sur plusieurs hypothèses fortes : pluie uniforme sur l’ensemble du bassin, intensité constante sur la durée critique, coefficient C global et réponse hydrologique simplifiée. Ces hypothèses sont souvent acceptables pour les petits bassins, mais elles deviennent plus fragiles à mesure que la surface augmente ou que l’hétérogénéité spatiale s’accroît. De même, la formule de Kirpich n’est qu’une approximation empirique. Pour des études détaillées, des modèles événementiels ou continus plus complets peuvent être nécessaires.

Malgré ces limites, la combinaison QGIS + calculs morphométriques reste extrêmement puissante dans les phases d’avant-projet, de diagnostic territorial et de pré-dimensionnement. Elle permet de produire rapidement une base quantitative cohérente et cartographiée, à condition de documenter soigneusement les hypothèses retenues.

Ressources institutionnelles et universitaires recommandées

Pour approfondir la délimitation des bassins versants, les MNT et l’analyse hydrologique, consultez des sources faisant autorité comme le U.S. Geological Survey, les jeux d’élévation de la NOAA pour le contexte topographique et climatique, ainsi que les ressources pédagogiques d’universités comme Carleton College sur l’analyse SIG et les bassins versants.

Comment interpréter les résultats du calculateur

Lorsque vous lancez le calcul, la page affiche plusieurs indicateurs. Si l’indice de Gravelius est proche de 1,2 ou inférieur, le bassin est relativement compact et peut produire des pointes plus synchrones. Si la pente moyenne dépasse quelques pourcents, la vitesse de transfert augmente généralement. Un temps de concentration court signale une réactivité potentiellement forte et justifie l’utilisation d’intensités de pluie élevées à courte durée. Le débit de pointe obtenu doit être comparé aux sections hydrauliques existantes, aux capacités de buses ou de fossés, ou à des événements de référence dans des bassins similaires.

En pratique, l’expert confronte toujours ces résultats à des observations de terrain : présence de zones humides, de drains agricoles, d’ouvrages de ralentissement, d’endiguements, d’imperméabilisation locale, de sols argileux ou d’affleurements rocheux. Le SIG ne remplace pas l’hydrologie de terrain, mais il fournit un cadre spatial rigoureux pour structurer l’analyse et documenter les choix techniques.

Avertissement : les résultats affichés sont des estimations simplifiées destinées à l’analyse préliminaire. Pour un dossier réglementaire, un dimensionnement sensible ou une expertise d’inondation, faites valider les hypothèses, les données pluviométriques et les formules par un hydrologue ou un bureau d’études qualifié.