Analyse des Dammbruchs der Tailings Storage Facility (TSF): Schlammfluss mit hohem Feststoffgehalt und Routenführung von Schlammüberschwemmungen
Computergestützte hydraulische Modellierung und numerische Hochwassersimulation haben revolutionäre Paradigmenwechsel erfahren. Die Kombination aus hochauflösenden, von LiDAR abgeleiteten digitalen Höhenmodellen (DEMs), GPU-beschleunigten parallelen Lösungsarchitekturen und gekoppelten hydrodynamisch-hydrologischen Formulierungen ermöglicht es Wasserressourceningenieuren, Überschwemmungsdynamiken mit beispielloser physikalischer Genauigkeit zu simulieren. Beim AnalysierenAnalyse des Dammbruchs der Tailings Storage Facility (TSF): Schlammfluss mit hohem Feststoffgehalt und Routenführung von SchlammüberschwemmungenModellierungsteams müssen eine Brücke zwischen theoretischer Strömungsmechanik, numerischen Stabilitätskriterien und praktischen Regulierungsrahmen schlagen.
Formulierung nicht-Newtonscher viskoplastischer Flüssigkeitsparameter, der Bingham-Viskosität und der volumetrischen Quellung in Modellen für Dammbrüche bei katastrophalem Bergbauabfall.
In diesem umfassenden technischen Leitfaden für Ingenieure werden die maßgeblichen Gleichungen, numerischen Diskretisierungsschemata, Netzoptimierungsstrategien, Parametrisierungen hydraulischer Strukturen und reale Kalibrierungsabläufe aufgeschlüsselt, die für die Durchführung robuster hydrodynamischer Hochwassersimulationen erforderlich sind.
1. Mathematische Formulierungen und hydrodynamische Grundgleichungen
Im Mittelpunkt der modernen Modellierung von Fluss-, Stadt- und Überlandfluten stehen die zweidimensionalen Flachwassergleichungen (2D SWE) – auch als 2D-Saint-Venant-Gleichungen bekannt –, die durch Tiefenintegration der Navier-Stokes-Gleichungen unter der Annahme einer hydrostatischen Druckverteilung abgeleitet werden:
Die konservative Form der 2D-Flachwassergleichungen wird ausgedrückt als:
Zweidimensionale Flachwasser-Konservierungsmatrixformulierung:
∂U/∂t + ∂E/∂x + ∂G/∂y = S
Wo der konservierte Vektor
Uund gerichtete FlussvektorenEundGsind definiert als:
U = [ h, h·u, h·v ]ᵀ
E = [ h·u, h·u² + ½·g·h², h·u·v ]ᵀ
G = [ h·v, h·u·v, h·v² + ½·g·h² ]ᵀ
Und der kombinierte Quellbegriff
Sberücksichtigt die Bodenneigung und die Bodenreibungsscherspannung:
S = [ 0, -g·h·(∂z_b/∂x + S_fx), -g·h·(∂z_b/∂y + S_fy) ]ᵀ
Wo:
-
hstellt die Wassertiefe dar undz_bist Betthöhe.-
uundvsind tiefengemittelte Strömungsgeschwindigkeiten in der x- und y-kartesischen Richtung.-
gist die Erdbeschleunigung (9,81 m/s²).-
S_fxundS_fysind Mannings Reibungssteigungen (S_fx = (n²·u·√(u²+v²)) / h^(4/3)).
In Überschwemmungsgebieten mit geringem Gefälle, in denen konvektive Beschleunigungsterme vernachlässigbar sind, vereinfacht die Diffusionswellen-Approximation die Berechnung durch Weglassen von Trägheitstermen, was die Laufzeiten drastisch verkürzt und gleichzeitig die Tiefengenauigkeit beibehält. In städtischen Straßenkorridoren mit schnellen Wellenfronten, Dammbruchwellen oder überkritischen Kontraktionen sind jedoch vollständige Impulsgleichungen unbedingt erforderlich, um eine Unterschätzung der Spitzenwasserstände zu vermeiden. Detaillierte Informationen zur Mechanik der geotechnischen Brucherosion finden Sie in unserer Analyse Wahrscheinlichkeitssimulation von Dammbrüchen und Deichbrüchen .
2. Integrierter numerischer Workflow: Gefahr von Gelände bis Überschwemmung
Die Entwicklung vertretbarer Hochwassermodelle erfordert eine strenge Datenvorverarbeitung, Netzgenerierung, Randbedingungserzwingung und Qualitätsüberprüfung nach der Simulation.
Die End-to-End-Pipeline für die hydrodynamische Modellierung ist wie folgt aufgebaut:
<!--ec:block {"_type": "Flussdiagramm", "Quelle": "Flussdiagramm TD
Daten[LiDAR- und Bathymetrie-Datenerfassung] -> HydroEnforce[Hydro-Enforcement: Bridge & Culvert Notching]
HydroEnforce --> MeshGen[Computergestützte Netzgenerierung: Unstrukturiert / Untergitter]
MeshGen --> ParamRough[Räumlich verteiltes Manning-Rauheitskartierung]
ParamRough --> BCSetup[Randbedingungen: Ganglinien, Gezeiten und Niederschlag]
BCSetup --> Solvers{Numerische Engine: 2D-SWE vs. Diffusionswelle}
Solver --> Simulation[GPU-beschleunigtes Zeitschrittverfahren und CFL-Courant-Stabilität]
Simulation -> PostProcess[Wasseroberflächenhöhe und Scherspannungsinterpolation]
PostProcess --> HazardClass[Flood Hazard Index & Velocity-Depth Slicing]
HazardClass --> Verification{Messgerät- und Hochwassermarkenkalibrierung?}
Verifizierung – Kalibriert –> Exportieren[FEMA / Regulatory Inundation Boundary Delivery]
Verifizierung – Residual High –> TuneParam[Manning n & Infiltration Rekalibration] –> Simulation
Durch die Durchsetzung strikter Hydroflattening-Bruchlinien und die Kalibrierung der hydraulischen Reibung anhand beobachteter Pegelganglinien eliminieren Modellierer unerwünschte Trockenzellen-Poolings und liefern zertifizierte technische Karten. Zur Konvergenz von Küsten- und Flussgewässern mit mehreren Gefahren lesen Sie unsere umfassende Studie zum Thema gekoppelte zusammengesetzte Hochwassermodellierung (Fluss-, Pluvial- und ADCIRC-Flutflut) .
3. Computergestützte Netzschemata und numerische Diskretisierung
Die Auswahl der Computernetztopologie regelt direkt die numerische Diffusion, die Erhaltung der Wellenfront und die Ausführungseffizienz:
| Mesh Topology Scheme | Spatial Resolution | Computational Demand | Optimal Application Domain |
| Structured Cartesian Grid | Uniform 1m - 10m cells | Low to Moderate | Broad regional flat terrain and rural floodplains |
| Unstructured Triangular / Poly | Variable 0.5m - 50m cells | High (requires GPU parallel) | Complex urban geometry, street corridors, bridge piers |
| Sub-Grid High-Res Bathymetry | Detailed hydraulic tables | Optimized High Performance | Natural channel conveyance with large floodplain cells |
| 1D-2D Coupled Hybrid Mesh | 1D channels + 2D overbanks | Minimal | Long river networks where 2D channel meshes are costly |
Kernnetz- und Solver-Richtlinien
- Courant-Friedrichs-Lewy (CFL)-Zustand: Behalten Sie in expliziten Zeitschrittschemata die Courant-Zahl bei
C = (√(g·h) + |u|) · (Δt / Δx) ≤ 1,0um eine katastrophale numerische Divergenz zu verhindern. - Bathymetrische Sub-Grid-Probenahme: Moderne Solver berechnen Bühnenflächen- und Bühnentransportkurven aus zugrunde liegenden 0,5-Meter-DEMs vor und erfassen die Mikrotopographie sogar innerhalb von 10-Meter-Berechnungszellen.
- Bruchkantenausrichtung: Richten Sie die Flächen der Berechnungszellen immer entlang der Deichkämme, Straßenkronen und Kanalufer aus, um eine künstliche numerische Überschreitung zu verhindern. Benchmarks zur Hochgeschwindigkeits-GPU-Ausführung finden Sie in unserer Forschung zu GPU-beschleunigte Flachwasser-Gleichungs-Engines .
4. Hydraulic Roughness & Manning's n Parametrisierungs-Benchmarks
Die genaue Zuordnung der Oberflächenrauheit ist der empfindlichste Parameter bei der hydrodynamischen Kalibrierung:
| Land Cover Classification | Standard Manning n | Calibration Range | Hydraulic Impact on Wave Velocity |
| Concrete Flume / Paved Road | 0.013 – 0.015 | 0.011 – 0.018 | Minimal resistance; high supercritical velocity |
| Main River Channel (Clean) | 0.030 – 0.035 | 0.025 – 0.045 | Governs base channel conveyance and stage |
| Cultivated Farmland / Turf | 0.040 – 0.050 | 0.035 – 0.060 | Moderate floodplain attenuation; hydrograph lag |
| Dense Riparian Forest / Brush | 0.080 – 0.120 | 0.060 – 0.160 | Significant stage surcharge; energy dissipation |
| Urban Building Block Footprint | 0.150 – 0.300 | 0.100 – 0.500 | Approximates building blockage and macro-roughness |
5. Fünf technische Regeln für eine vertretbare hydrodynamische Modellierung
Um die Reproduzierbarkeit des Modells, die Einhaltung gesetzlicher Vorschriften und die physikalische Gültigkeit sicherzustellen, sollten Wasserressourcenexperten die folgenden Grundprinzipien befolgen:
- Überprüfen Sie das Massenerhaltungsgleichgewicht: Überprüfen Sie immer den kumulativen Volumenfehler (
Volumenfehler % = (Zufluss – Abfluss – ΔSpeicher) / Gesamtzufluss · 100). Stellen Sie sicher, dass der globale Kontinuitätsfehler streng unter 1,0 % bleibt. - Führen Sie Qualitätsaudits zur Wasserdurchsetzung durch: Überprüfen Sie manuell alle Brücken, Durchlässe und Straßenböschungen im DEM. Durch ungekerbte Straßenkreuzungen entstehen künstliche Dämme, die flussaufwärts gelegene Gemeinden fälschlicherweise überschwemmen.
- Simulieren Sie die Empfindlichkeit variabler Rauheit: Führen Sie automatisierte Sensitivitätsdurchläufe durch, bei denen Mannings n-Werte um ±20 % variiert werden, um das Band der Wasseroberflächenunsicherheit für die Risikokommunikation zu quantifizieren.
- Vermeiden Sie verkürzte Grenznähe: Positionieren Sie die hydraulischen Grenzen flussaufwärts und flussabwärts ausreichend weit vom Untersuchungsgebiet entfernt (mindestens 5 bis 10 Kanalbreiten), um künstliche Grenzreflexionsartefakte zu verhindern.
- Kalibrieren Sie anhand unabhängiger historischer Benchmarks: Kalibrieren Sie Modelle sowohl für kontinuierliche Ganglinienaufzeichnungen (Nash-Sutcliffe-Effizienz > 0,80) als auch für vermessene physikalische Hochwassermarken (mittlerer Fehler < 0,15 Meter).
6. Häufig gestellte Fragen (FAQ)
Wann sollte 2D-SWE mit vollem Impuls anstelle der Diffusionswellen-Näherung verwendet werden?
2D-SWE mit vollem Impuls muss immer dann genutzt werden, wenn Trägheitskräfte, dynamische Wellenstoßfronten, schnelle Strömungsbeschleunigungen oder hydraulische Übergänge dominieren. Dazu gehören Dammbrüche, steile städtische Sturzfluten, Kontraktionen von Brückenpfeilern, Strömungen um scharfe Straßenecken und Überläufe mit hoher Geschwindigkeit. Diffusionswellengleichungen sind nur für Flussauen mit geringem Gefälle und langsamem Anstieg geeignet.
Wie werden unterirdische Regenwasserkanalrohre mit 2D-Überlandströmungsgittern gekoppelt?
Gekoppelte 1D-2D-Motoren (z. B. EPA SWMM, verbunden mit 2D-Hydrodynamiklösern) tauschen Masse und Impuls an Straßenauffangbecken und Mannlöchern aus. Wenn die Kapazität der unterirdischen Leitungen überschritten wird, strömt Wasser durch Schachtabdeckungen auf das 2D-Oberflächengitter als Überlandwasser; Wenn die Oberflächenkapazität es zulässt, läuft das Wasser basierend auf den Formeln zum Auffangen von Öffnungen und Wehren wieder nach unten ab.
Was ist die Hauptursache für numerische Instabilität in 2D-Überschwemmungsmodellen mit endlichem Volumen?
Die häufigste Ursache für numerische Instabilität sind Nass-Trocken-Zellgrenzenübergänge in Verbindung mit übermäßig großen Zeitschritten, die die Courant-Bedingung überschreiten. Moderne Schockerfassungslöser verwenden positivitätserhaltende Riemann-Löser (z. B. HLLC- oder Roe-Schemata) und adaptive Zeitschrittverfahren, um nicht negative Wassertiefen an steilen Benetzungsfronten aufrechtzuerhalten.
Fazit: Aufbau widerstandsfähiger Gemeinschaften durch Computerhydraulik
Da Urbanisierung und Klimainstationarität extreme Niederschlagsereignisse verstärken, dient eine hochauflösende hydrodynamische Simulation als Grundlage für die Katastrophenvorsorge, die Gestaltung der zivilen Infrastruktur und die lebenssichere Evakuierungsplanung. Durch die Anwendung strenger mathematischer Formeln, die Validierung anhand empirischer Überschwemmungsaufzeichnungen und die Nutzung moderner Parallelrechnungen können Wasserbauingenieure Überschwemmungsgefahren genau kartieren und gefährdete Gemeinden weltweit schützen.
