Spatio-Temporal Performance of 2D Local Inertial Hydrodynamic Models for Urban Drainage and Dam-Break Applications

Cette étude démontre que le modèle hydrodynamique 2D HydroPol2D, basé sur l'approximation inertielle locale, offre un compromis efficace entre précision et rapidité de calcul (23 fois plus rapide que HEC-RAS 2D) pour la modélisation des inondations urbaines et des ruptures de barrage, bien que l'omission des infrastructures urbaines et l'absence d'inertie convective puissent entraîner des écarts significatifs dans les débits de pointe et la vitesse de propagation des ondes de crue.

L'essentiel

🌊 Le Grand Défi : Prévoir les Inondations sans se noyer dans les calculs

Imaginez que vous devez prédire comment l'eau va inonder une ville après une forte pluie ou la rupture d'un barrage. C'est comme essayer de suivre le mouvement d'une foule immense dans un métro bondé.

Les scientifiques ont deux façons de faire :

1. La méthode "Super-Héros" (Modèles complets) : Ils calculent chaque détail, chaque tourbillon, chaque collision de l'eau. C'est très précis, mais c'est comme si vous deviez compter chaque goutte d'eau individuellement. Le problème ? Cela prend tellement de temps de calcul que, quand le résultat arrive, l'inondation est déjà terminée. C'est trop lent pour une alerte en temps réel.
2. La méthode "Estimation Rapide" (Modèles locaux-inertiels) : C'est une version simplifiée. On ignore certains détails complexes (comme les petits tourbillons) pour aller plus vite. C'est comme regarder la foule de loin et dire "ça va aller vers la gauche". C'est rapide, mais parfois, si la foule se comporte bizarrement (eau très rapide), l'estimation peut être fausse.

L'objectif de cette étude était de tester si cette méthode rapide (le modèle HydroPol2D) est assez fiable pour des situations dangereuses, comme la rupture d'un barrage, et si elle peut gérer les "obstacles" urbains (égouts, ponts, pompes) sans avoir besoin de connaître chaque tuyau de la ville.

---

🛠️ Les Trois Outils de Test (Les "Recettes" de l'eau)

Les chercheurs ont testé trois versions différentes de leur modèle rapide, comme si c'étaient trois chefs cuisiniers essayant de faire le même plat :

1. La recette originale (lim) : La version de base.
2. La recette "Centrée" (s-centered) : Une version qui lisse les erreurs.
3. La recette "Vent arrière" (s-upwind) : Une version qui anticipe mieux la direction de l'eau.

Ils ont comparé ces trois recettes à la méthode "Super-Héros" (le logiciel HEC-RAS) qui sert de référence absolue.

---

🧪 Les Quatre Scénarios de Test

Pour voir si leurs outils fonctionnaient, ils ont joué à quatre jeux différents :

1. La vague sur un tapis roulant (Test théorique)

Imaginez une goutte d'eau glissant sur un plan parfaitement lisse. C'est le test de base.

• Résultat : La recette "Vent arrière" (s-upwind) a été la plus précise et la plus stable, même quand l'eau glissait très vite. La recette originale a parfois "trébuché" (instabilité).

2. Le Bassin de rétention (Le robinet et la vanne)

Ils ont simulé un bassin qui reçoit de l'eau et la relâche par un tuyau (culvert) et une décharge (spillway).

• Le problème : Souvent, on ne connaît pas les plans exacts des tuyaux dans les vieilles villes.
• L'astuce : Au lieu de modéliser chaque tuyau en détail, ils ont utilisé une "formule magique" (une courbe de débit) basée sur la taille du trou.
• Résultat : Cette astuce a fonctionné à merveille ! L'erreur était inférieure à 5% par rapport au modèle complexe. C'est comme deviner le débit d'un robinet en regardant juste la taille de la poignée, et avoir raison presque à tous les coups.

3. La Ville sans égouts (Le piège de la carte)

Ils ont pris une ville réelle (São Paulo) avec une carte très précise (LiDAR), mais sans information sur les égouts.

• Sans les égouts (modèle simple) : L'eau s'accumule partout, comme si la ville était un champ de boue. Le modèle pense qu'il y a une inondation énorme, alors que l'eau s'écoule normalement dans les tuyaux. L'erreur de débit à la sortie est de 17,5% ! C'est énorme.
• Avec les égouts (modèle intelligent) : En ajoutant les "trous" (égouts, tunnels) comme des vannes dans le modèle, l'eau s'écoule correctement.
• Le gain : Le modèle avec les vannes a été deux fois plus rapide à calculer car il n'avait pas besoin de simuler des inondations inutiles qui n'existent pas en réalité.

4. La Rupture de Barrage (Le scénario catastrophe)

C'est le test ultime : un barrage de 42 mètres de haut qui cède soudainement, envoyant une vague de 49 000 m³ d'eau par seconde vers une ville de 200 000 habitants.

• La vitesse : Le modèle rapide a été 23 fois plus rapide que le modèle complet. Au lieu de 43 heures de calcul, il a fallu moins de 2 heures. C'est crucial pour une alerte en temps réel !
• La précision :
  • Pour savoir où l'eau va aller (la zone inondée), le modèle rapide est excellent (95% de réussite).
  • Pour savoir quand l'eau arrive exactement, il y a un petit décalage. Comme le modèle ignore certains effets de vitesse, la vague arrive un tout petit peu plus tôt dans la simulation que dans la réalité.
  • Près du barrage (où l'eau va très vite), les erreurs sont plus grandes, mais plus loin dans la ville, c'est très précis.

---

💡 Les Grandes Leçons à retenir

1. La simplicité a du sens : On n'a pas besoin d'un super-ordinateur pour tout savoir. Une version simplifiée, bien calibrée, peut prédire les inondations urbaines avec une grande précision.
2. Les "trous" dans la ville comptent : Si vous modélisez une ville sans tenir compte des égouts et des ponts, vous allez paniquer pour rien (fausses inondations) ou sous-estimer le danger. Il faut ajouter ces "vannes" dans le modèle, même sans avoir les plans exacts.
3. Le temps, c'est de l'argent (et des vies) : Parce que ce modèle est 23 fois plus rapide, on peut l'utiliser pour faire des centaines de simulations en même temps. Par exemple, on peut simuler 100 scénarios de rupture de barrage différents pour voir quelle est la pire situation possible et mieux préparer les secours.

En résumé : Les chercheurs ont créé un outil de prédiction d'inondation qui est comme un moteur de voiture hybride. Il est assez puissant pour gérer les situations extrêmes (rupture de barrage) et assez économe en énergie (calcul rapide) pour être utilisé en temps réel, à condition de bien régler les "vannes" de la ville. C'est une avancée majeure pour protéger les populations.

Résumé technique

1. Problématique

La modélisation précise et rapide des inondations est cruciale pour l'analyse des risques et la prévision, notamment dans les zones urbaines denses et lors d'événements extrêmes comme les ruptures de barrages.

• Dilemme Complexité/Performance : Les modèles hydrodynamiques complets (résolvant toutes les équations de Saint-Venant, y compris l'inertie convective) offrent une grande précision mais nécessitent des temps de calcul prohibitifs pour les grandes zones à haute résolution, les rendant inadaptés à la prévision en temps réel ou aux analyses probabilistes.
• Limites des modèles simplifiés : Les modèles à faible complexité, tels que les approximations « inertielle locale » (qui négligent l'inertie convective), sont plus rapides mais peuvent manquer de précision dans les régimes d'écoulement supercritiques (nombre de Froude > 1) ou dans des terrains plats.
• Défi Urbain : En milieu urbain, l'absence de données détaillées sur les réseaux d'assainissement (égouts, déversoirs, pompes) conduit souvent à une mauvaise représentation de la dynamique de l'eau. Les modèles standards traitent souvent les bâtiments comme des obstacles topographiques, ce qui crée des artefacts numériques et surestime les zones inondées, sans tenir compte des mécanismes de drainage réels.

L'objectif de l'étude est d'évaluer la performance spatio-temporelle d'un modèle inertielle local 2D (HydroPol2D) par rapport à un solveur complet (HEC-RAS 2D) et d'analyser l'impact de l'intégration de conditions aux limites internes représentant les infrastructures urbaines.

2. Méthodologie

Les auteurs ont développé et testé le modèle HydroPol2D en utilisant trois schémas numériques différents basés sur l'approximation inertielle locale :

1. Formulation originale (lim) : Celle proposée par Bates et al. (2010).
2. Schéma centré (s-centered) : Introduit pour ajouter une diffusion numérique artificielle.
3. Schéma amont (s-upwind) : Une autre variante pour stabiliser les écoulements rapides.

Approche des infrastructures urbaines :
Au lieu de coupler un modèle 1D complexe de réseau d'égouts (qui nécessite des données topographiques détaillées souvent indisponibles), les auteurs ont implémenté des conditions aux limites internes basées sur des lois hydrauliques simples (courbes de tarage). Ces conditions modélisent :

• Les conduits (culverts) et déversoirs.
• Les orifices et les pompes.
• Ces structures sont définies par des paramètres géométriques (ex: surface de la section) dérivés de plans ou de données GIS (Google Earth), réduisant ainsi le nombre de paramètres à calibrer.

Études de cas :
Quatre cas d'étude ont été simulés et comparés au solveur complet HEC-RAS 2D :

1. Validation théorique : Propagation d'une onde non cassante sur un plan horizontal (comparaison avec une solution analytique).
2. Bassin de rétention (São Paulo) : Simulation d'un bassin recevant un hydrogramme de crue centennale, avec contrôle par un conduit et un déversoir.
3. Bassin urbain (São Paulo) : Comparaison d'une simulation urbaine avec et sans représentation des infrastructures de drainage (conduits, tunnels) dans un bassin de 2,36 km² soumis à une pluie de 50 ans.
4. Rupture de barrage (Pirapama, Brésil) : Scénario de rupture de barrage impactant la ville de Cabo de Santo Agostinho (200 000 habitants), avec un débit de pointe de ~49 000 m³/s et des zones de terrain très plat.

3. Contributions Clés

• Intégration d'infrastructures simplifiée : Démonstration qu'il est possible de représenter fidèlement le drainage urbain (conduits, déversoirs) via des conditions aux limites internes basées sur des lois de tarage, sans nécessiter de données de réseau d'égouts exhaustives.
• Comparaison de schémas numériques : Évaluation approfondie de trois schémas inertiels locaux (lim, s-centered, s-upwind) dans des conditions variées, y compris des écoulements supercritiques (rupture de barrage).
• Gain de temps de calcul : Validation que le modèle inertielle local peut être utilisé pour des scénarios complexes (rupture de barrage) avec un temps de calcul considérablement réduit par rapport aux modèles complets.
• Analyse spatio-temporelle : Au-delà de la simple comparaison des zones inondées maximales, l'étude analyse la performance du modèle au cours du temps (arrivée de l'onde, hydrogrammes).

4. Résultats Principaux

A. Performance des schémas numériques :

• Onde non cassante : Le schéma s-upwind présente la meilleure concordance avec la solution analytique, tandis que la formulation originale (lim) montre des instabilités près des fronts mouillés/séchés.
• Rupture de barrage :
  • Le modèle HydroPol2D prédit correctement les zones inondées maximales avec un indice de succès critique (CSI) élevé : 0,95 (lim), 0,92 (s-centered) et 0,89 (s-upwind).
  • Cependant, l'absence d'inertie convective entraîne une avance plus rapide de l'onde de crue par rapport au modèle complet HEC-RAS, surtout dans les zones très plates et près de la rupture.
  • Le schéma s-upwind offre le meilleur compromis global (meilleurs scores NSE, KGE, RMSE) pour prédire à la fois les étendues et les profondeurs, bien que le schéma centré ait un CSI légèrement supérieur pour la profondeur maximale.

B. Impact des infrastructures urbaines (Cas 2) :

• Sans infrastructures : Le modèle surestime les profondeurs d'eau de près de 1,8 m et les zones inondées, car l'eau s'accumule artificiellement dans les dépressions topographiques du LiDAR. Les hydrogrammes de sortie sont totalement décalés (différence de pic de 17,5 % et décalage temporel important).
• Avec infrastructures : L'intégration des conditions aux limites (conduits, tunnels) permet de drainer correctement l'eau, réduisant les erreurs de profondeur et corrigeant les hydrogrammes.
• Efficacité : Paradoxalement, la simulation avec les infrastructures est deux fois plus rapide que celle sans, car elle évite les itérations numériques complexes nécessaires pour résoudre les écoulements bloqués dans les dépressions topographiques.

C. Performance temporelle :

• Pour le scénario de rupture de barrage, HydroPol2D est 23 fois plus rapide que HEC-RAS 2D (114 min contre 2 610 min sur un CPU standard), rendant possible les analyses probabilistes et les prévisions en temps réel.

5. Signification et Conclusion

Cette étude démontre que les modèles hydrodynamiques 2D à inertie locale, lorsqu'ils sont couplés à une représentation adéquate des infrastructures urbaines via des conditions aux limites internes, constituent une alternative robuste et efficace aux modèles complets.

• Pour la prévision : Le gain de temps de calcul (facteur >20) permet d'envisager des ensembles probabilistes (Monte Carlo) pour les ruptures de barrages, intégrant les incertitudes sur la géométrie de la brèche et les paramètres hydrauliques.
• Pour l'urbanisme : La méthode proposée offre une solution pratique pour les pays en développement ou les zones où les données de réseaux d'assainissement sont absentes, permettant de simuler des écoulements réalistes avec des données topographiques brutes.
• Limites : Bien que performants pour les étendues d'inondation, ces modèles simplifiés doivent être utilisés avec prudence pour la prédiction précise des temps d'arrivée des crues dans des régimes supercritiques ou sur des terrains très plats, où l'absence d'inertie convective introduit un biais d'avance.

En conclusion, l'approche proposée équilibre efficacement précision et rapidité, ouvrant la voie à une utilisation plus large de la modélisation hydrodynamique dans les systèmes d'alerte précoce et la gestion des risques urbains.