Probabilistische Dammbruch- und Deichbruchsimulation: Dynamische geotechnische Erosionsmodellierung und hydrodynamische Wellenführung

Dämme, Stauseen und Hochwasserschutzdeichsysteme gehören zu den wichtigsten zivilen Infrastrukturanlagen der modernen Gesellschaft. Sie sichern die Wasserversorgung für die landwirtschaftliche Bewässerung, erzeugen erneuerbare Wasserkraft und schützen Millionen von Bewohnern historischer Überschwemmungsgebiete vor Überschwemmungen. Das katastrophale Versagen eines Erddamms oder einer Deichdammung setzt jedoch innerhalb von Minuten große Mengen potenzieller Energie frei und erzeugt heftige, katastrophale Flutwellen, die ganze flussabwärts gelegene Gemeinden, Brücken und Industriekorridore auslöschen können.

In der Vergangenheit stützten sich die Risikobewertungen von Dammbrüchen auf statische, empirische Bruchformeln (wie die klassischen Froehlich-, MacDonald-Langridge-Monopolis- oder Von-Thun-Gillette-Gleichungen) gepaart mit eindimensionaler (1D) hydraulischer Streckenführung. In diesen Legacy-Frameworks:

• Der Ingenieur ging von einer festen Bruchgeometrie (trapezförmige Bruchform, Bodenbreite und Seitenneigungen) und einer festen Bruchentwicklungszeit (z. B. 1,5 Stunden) basierend auf historischen Regressionsstreudiagrammen aus.
• Der Abfluss des Bruchs wurde flussabwärts durch 1D-Querschnitte geleitet, die die 2D-Talausbreitung der Auen, die Massendichte des Sedimenttransports und Wellenreflexionen an den Canyonwänden völlig ignorierten.

Dieser statische empirische Ansatz weist gefährliche technische Einschränkungen auf. Erdwälle versagen nicht nach willkürlichen geometrischen Mustern; Sie scheitern durch komplexe, nichtlineare ProzesseGeotechnische Erosionsmechanik: fortschreitende rückwärtige Überlauferosion, Schererosion der internen Rohrleitungen und Rotationshangstabilitäts-Gleitkreiskollaps. Darüber hinaus ignorieren statische Annahmen die Unsicherheit der geotechnischen Parameter: Bodenverdichtung, Kohäsion, Korngrößenverteilung (D_{50}), kritische Schubspannung (tau_c) und Schwankungen des Wasserspiegels im Stausee.

Ende 2026 erfordern behördliche Sicherheitsstandards (vorgeschrieben von der Federal Energy Regulatory Commission, FERC, und dem US Army Corps of Engineers, USACE).Probabilistische, physikalisch basierte Simulationen von Dammbrüchen und Deichbrüchen. Durch die direkte Kopplung transienter geotechnischer Brucherosionsphysik mit GPU-beschleunigten zweidimensionalen (2D) Flachwassergleichungslösern und Monte-Carlo-Ensemble-Unsicherheitsquantifizierung können Staudammsicherheitsingenieure flussabwärts gelegene Überschwemmungsgefahren mit kompromissloser mathematischer Präzision kartieren.

Diese umfassende technische Abhandlung untersucht die geotechnische Erosionsphysik, die Hydrodynamik nichtnewtonscher sedimentbeladener Wellen, die Monte-Carlo-Ensemble-Probenahme und Protokolle zur Gefahrenkartierung bei Notfallevakuierung, die für die Durchführung modernster Dammbruchbewertungen erforderlich sind.

1. Geotechnische Bruchmechanismen: Überlauf vs. Erosion interner Rohrleitungen

Erddämme versagen überwiegend durch einen von zwei physikalischen Mechanismen: hydraulische Überflutung oder interne geotechnische Rohrleitungen.

Die Gleichung für übermäßige Scherspannungserosion

In physikalisch basierten Bruchmodellen (wie USDA-ARS WinDAM oder DL-Breach) wird die Bodenablösung entlang der Bruchgrenze durch die bestimmtFormulierung mit übermäßiger Scherspannung:

ε_r = k_d (τ_e - τ_c)^a

Wo:

• epsilon_rstellt die Erosionsrate der Bodenablösung dar (m/s oder cm/h).
• k_dist der Bodenerosionskoeffizient (cm^3/N·s), bestimmt durch In-situ-Jet-Erosionstests (JET) oder Lochdispersionstests unter Wasser.
• tau_eist die effektive hydrodynamische Scherspannung, die durch fließendes Wasser auf die Bruchgrenze ausgeübt wird:

τ_e = γ_w R_h S_f

Wogamma_wist das Einheitsgewicht von Wasser,R_hist der hydraulische Radius des Bruchkanals undS_fist die Energiesteigung.

• tau_cist die kritische Scherspannung des verdichteten Böschungsbodens (Pa), unterhalb derer keine Ablösung erfolgt.
• $a$ ist ein empirischer Exponent, typischerweise angenommen als $a = 1,0$.

Wenn fließendes Wasser durch eine entstehende Bruchkerbe beschleunigt wird,tau_ekann 500 Pa überschreiten – weit über der kritischen Scherfestigkeit des Bodens (tau_c ≈ 5 - 20 Pa) – was zu exponentiellen volumetrischen Erosionsraten führt.

2. Berechnung der dynamischen Abflusshydrographie: Breitkammwehr und Öffnungshydraulik

Wenn sich die Bruchkerbe erweitert und vertieft, kommt es zu einem sofortigen AbflussQ_{Verstoß}(t)wechselt vorübergehend zwischenuntergetauchte Öffnungsströmung(während der frühen internen Verrohrung) undfreiflächiger, breitkroniger Wehrstrom(sobald das Bruchdach einstürzt).

flowchart TD
    Pool["Reservoir Pool Level: H_res(t)"] --> Geom["Instantaneous Breach Notch Geometry<br/>• B_w(t), Z_s, H_b(t)"]
    Geom -->|Internal Piping Phase| Pipe["Pressurized Orifice Flow<br/>• Q = C_d · A_pipe(t) · √(2g · ΔH)"]
    Geom -->|Open Breach Crest Phase| Weir["Broad-Crested Weir Flow<br/>• Q = C_w · B_w(t) · (H_res - Z_b)^1.5 + C_s · Z_s · (H_res - Z_b)^2.5"]
    Pipe --> Wave["Downstream 2D Hydrodynamic Flood Wave Routing"]
    Weir --> Wave

Die kombinierte Trapezwehr-Formulierung mit breitem Kamm

Sobald ein offener Bruchgraben angelegt ist, wird der Abfluss über die kombinierte Gleichung für ein rechteckiges und dreieckiges Wehr mit breiter Krone berechnet:

Q_{breach}(t) = C_{w1} B_w(t) ≤ft( H_{res}(t) - Z_b(t) )^{3/2} + C_{w2} Z_s ≤ft( H_{res}(t) - Z_b(t) )^{5/2}

Wo:

• B_w(t)ist die momentane untere Breite der sich ausdehnenden Bruchkerbe (m).
• Z_b(t)ist die momentane Bodenhöhe der Bruchkerbe, die nach unten erodiert, bis sie das Grundgestein oder nicht erodierbare Fundamentschichten erreicht.
• Z_sist die Seitensteigung des trapezförmigen Bruchs (Z_s = (Delta x) / (Delta z)(typischerweise 0,5 bis 1,0 US-Dollar für kohäsive Tone und bis zu 1,5 US-Dollar für nicht kohäsive Sande).
• C_{w1}, C_{w2}sind Wehrabflusskoeffizienten, die für Unterwassereffekte unter Überflutung kalibriert sind.

Gleichzeitig wird die Beckenhöhe des Stausees erhöhtH_{res}(t)sinkt dynamisch gemäß der Reservoir-Phase-Speicher-Erhaltungsgleichung:

frac{dH_{res}}{dt} = frac{Q_{inflow}(t) – Q_{breach}(t) – Q_{spillway}(t)}{A_{res}(H_{res})}

WoA_{res}(H_{res})ist die Oberfläche des Reservoirbeckens in der HöheH_{res}.

3. Hydrodynamische Wellenführung: 2D nicht-newtonsche Dammbruchwellen

Eine Dammbruchwelle ist kein sauberes, klares Wasser; Es handelt sich um eine heftige, hochdichte Aufschlämmung aus Wasser, pulverisiertem Böschungskies, suspendierten Sedimenten und entwurzelten Bäumen.

In steilen Bergschluchten oder engen Flusstälern ist die Volumenkonzentration des Sediments (C_v) kann $20\% - 45\%$ überschreiten und die Flutwelle in eine umwandelnnicht-Newtonsche hyperkonzentrierte Aufschlämmung.

Nicht-Newtonsche Bingham-Plastische Hydrodynamik

>

Formulierung der scheinbaren Scherspannung:

τ = τ_y + µ_B · (∂u / ∂z)

>

Parameter:

• τ_y: Fließspannung des Sedimentschlamms (Pa) – muss überwunden werden, um eine Flüssigkeitsbewegung einzuleiten

• µ_B: Dynamische plastische Viskosität nach Bingham (Pa·s)

• Schüttdichte: ρ_slurry = ρ_w · (1 - C_v) + ρ_s · C_v

Der 2D-Stoßerfassungs-Hydrodynamiklöser

Um diese heftige Hochgeschwindigkeitswelle ohne numerische Instabilität oder unphysikalische Diffusion flussabwärts zu leiten:

• Modelle nutzen die 2D-Flachwassergleichungen mitBingham-Kunststoff- oder Herschel-Bulkley-Strömungswiderstandsquellenbegriffe (tau_{bx}, tau_{by}).
• Solver erzwingen Schockerfassungsschemata zweiter Ordnung (TVD (Total Variation Diminishing)), um Steile aufzulösenDamm bricht WellenbohrungenVorrücken mit Geschwindigkeiten über12 bis 20 m/sdas flussabwärts gelegene Tal hinunter.
• Impulsgleichungen lösen seitliche Zentrifugalkräfte auf, wenn die Welle scharfe Flussbiegungen umrundet, und sagen so den Hochlauf der Welle gegen die Canyonwände genau voraus.

4. Probabilistisches Monte-Carlo-Ensemble-Framework: Quantifizierung der geotechnischen Unsicherheit

Bei der forensischen Dammsicherheitstechnik sind die genauen geotechnischen Eigenschaften nie mit 100-prozentiger Sicherheit bekannt. Die Verdichtungsqualität, der Feuchtigkeitsgehalt, die Kohäsion und die innere Homogenität einer Erdböschung variieren über ihre lineare Länge.

Anstatt eine einzige deterministische „Worst-Case“-Simulation durchzuführen, die möglicherweise unrealistisch konservativ oder gefährlich optimistisch ist, führen moderne Ingenieure sie ausMonte-Carlo-Probabilistische Ensembles.

flowchart TD
  Params["Input Parameter Probability Distributions\n• Reservoir Head: Normal(μ = Full Pool, σ = 1.2m)\n• Soil Erodibility k_d: Log-Normal(μ = 10 cm³/N-s, σ = 0.4)\n• Critical Shear Stress τ_c: Weibull(λ = 8.5 Pa, k = 2.1)\n• Cohesion c': Normal(μ = 25 kPa, σ = 5 kPa)\n• Breach Width B_w: Beta distribution (Froehlich)"]
  LHS["Latin Hypercube Sampling (LHS)\nGenerates 500 Stratified Scenarios"]
  GPU["Automated GPU Solver Cluster\n500 Coupled 2D Simulations (4.2 GPU Hours)"]
  Hazard["Probabilistic Hazard Outputs\nP10, P50, P90 Inundation Maps & Arrival Time Bands"]

  Params --> LHS --> GPU --> Hazard

Probabilistische Überschwemmungsmetriken

Aus dem 500-Run-Ensemble extrahiert das System räumlich verteilte Überschreitungswahrscheinlichkeiten:

• P90 Überschwemmungskarte: Die konservative Überschwemmungsgrenze mit einem Konfidenzniveau von 90\%$, die für die Gestaltung des Schutzes der zivilen Infrastruktur (Kernkraftwerke, Krankenhauskomplexe, Chemikalienlager) verwendet wird.
• P50 Überschwemmungskarte: Das mittlere erwartete Szenario, das für die kommunale Zoneneinteilung, die Bewertung der Versicherungstarife und die primäre Evakuierungsroute verwendet wird.
• Unsicherheit der Wellenankunftszeit: Quantifiziert die frühestmögliche Ankunft der Hochgeschwindigkeitswellenfront in bestimmten flussabwärts gelegenen Bevölkerungszentren (T_{Ankunft, P10}), wodurch Notfallmanagern vertretbare Sicherheitsmargen für Evakuierungsfreigabefenster bereitgestellt werden.

5. Umfassende Vergleichsmatrix zur Dammbruchmethodik

In der folgenden Tabelle werden ältere empirische Methoden modernen, physikalisch basierten, probabilistischen 2D-Dammbruchsimulations-Frameworks gegenübergestellt:

6. Praxisnahe technische Fallstudie: Sicherheitsaudit für hochgefährliche Erddämme

Um diesen probabilistischen Rahmen für eine wichtige zivile Infrastrukturanlage zu validieren, führte FloodSim eine Prüfung des Notfallaktionsplans (EAP) zum Dammbruch eines 85 Fuß (26 Meter) hohen Erddamms durch, der ein 32.000 Acre Fuß (39,5 Millionen Kubikmeter) großes Trinkwasser- und Hochwasserschutzreservoir aufstaute, das 6 Meilen flussaufwärts eines städtischen Gemeindetals mit 95.000 Einwohnern liegt.

Kritische Entdeckungen zur Lebenssicherheit:

1. Die gefährliche „Lag Time“-Illusion: Das alte 1D-Modell prognostizierte, dass die erste Flutwelle 65 Minuten benötigen würde, um 6 Meilen flussabwärts bis zum Stadtzentrum zu wandern. Im Gegensatz dazu ergab der hochauflösende 2D-GPU-Solver, der die Talverengung und die überkritische Strömungskanalisierung entlang des Hauptflusskorridors berücksichtigte, dass die Hochgeschwindigkeits-Hydraulikbohrung den städtischen Kern treffen würdenur 38 Minuten. Hätten sich Notfallmanager auf das Altmodell verlassen, wären die Evakuierungsbefehle fast 30 Minuten zu spät eingetroffen.
2. Interstate Bridge Super-Scour: Die 2D-Simulation kartierte örtlich übersteigende Strömungsgeschwindigkeiten6,2 m/san den Pfeilern einer Autobahnbrücke, Berechnung der Bodenschubspannungen, die die Pfeilerfundamente innerhalb von 12 Minuten nach Eintreffen der Welle untergraben würden.
3. Neukalibrierung des Notfallaktionsplans: Ausgestattet mit exakten probabilistischen Ankunftskurven hat die kommunale Katastrophenschutzbehörde vier Evakuierungssirenen versetzt, automatisierte Autobahntorschließungen eingerichtet und ein zertifiziertes 22-minütiges Evakuierungsfreigabeprotokoll erstellt, um die absolute Sicherheit des öffentlichen Lebens zu gewährleisten.

7. Gefahrenklassifizierung nachgelagerter Evakuierung: Die australischen Gefahrenanfälligkeitskurven

Die Bewertung des Risikos für die Sicherheit des öffentlichen Lebens erfordert die Umrechnung der rohen Überschwemmungstiefe ($h$) und -geschwindigkeit ($v$) in Gefahrenkategorien für die menschliche Stabilität. Eine flache Wassertiefe von 0,4 Metern, die mit einer Geschwindigkeit von 3,0 m/s fließt, ist weitaus tödlicher als ein 1,0 Meter großes stehendes Wasserbecken, da hydrodynamische Widerstandskräfte erwachsene Fußgänger und Autos von der Fahrbahn fegen.

Moderne Notfallaktionspläne nutzen dieGefahrenanfälligkeitskurven für das australische Notfallmanagement (H1 bis H6):

<!--ec:block {"Typ: „Tabelle“, „hasHeaderRow“: true, „Zeilen“: [{“Typ: „tableRow“, „cells“:[{“Typ: „tableCell“, „content“:[{“type="span",text="GEFAHRENKLASSIFIZIERUNGSKATEGORIEN (H1 - H6)"}],"isHeader":true}]},{"Typ: „tableRow“, „cells“:[{“Typ: „tableCell“, „content“:[{“Typ: „span“, Text: „Gesamthydrodynamische Intensität: D \\times V = h \\times \\sqrt{u^2 + v^2} (m²/s)“}]}]},{“Typ: „tableRow“, „cells“:[{“Typ: „tableCell“, „content“:[{“Typ: „span“, „text“: „}]}]},{“Typ: „tableRow“, „cells“:[{“Typ: „tableCell“, „content“:[{“Typ: „Span“, Text: „[H1: Sicher für alle]: D \\times V \\le 0,4 \m²/s, h \\le 0,3m“}]}]},{“Typ: „tableRow“, „cells“:[{“Typ: „tableCell“, „content“:[{“Typ: „Spanne“, Text: „Kann sicher gehen; Kinder und ältere Menschen sicher.}]}]},{"Typ: „tableRow“, „cells“:[{“Typ: „tableCell“, „content“:[{“Typ: „span“, „text“: „}]}]},{“Typ: „tableRow“, „cells“:[{“Typ: „tableCell“, „content“:[{“Typ: „Span“, Text: „[H2: Unsicher für kleine]: D \\times V \\le 0,6 \m²/s, h \\le 0,5m“}]}]},{“Typ: „tableRow“, „cells“:[{“Typ: „tableCell“, „content“:[{“Typ: „span“, Text: „Unsicher für kleine Kinder und gebrechliche ältere Menschen.“}]}]},{“Typ: „tableRow“, „cells“:[{“Typ: „tableCell“, „content“:[{“Typ: „span“, „text“: „}]}]},{“Typ: „tableRow“, „cells“:[{“Typ: „tableCell“, „content“:[{“Typ: „Span“, Text: „[H3: Unsicher für Fahrzeuge]: D \\times V \\le 0,8 \m²/s, h \\le 0,6m“}]}]},{“Typ: „tableRow“, „cells“:[{“Typ: „tableCell“, „content“:[{“Typ: „span“, Text: „Limousinen und leichte Fahrzeuge verlieren die Traktion und schweben.“}]}]},{“Typ: „tableRow“, „cells“:[{“Typ: „tableCell“, „content“:[{“Typ: „span“, „text“: „}]}]},{“Typ: „tableRow“, „cells“:[{“Typ: „tableCell“, „content“:[{“Typ: „span“, Text: „[H4: Unsicher für alle Menschen]: D \\times V \\le 1,0 \m²/s, h \\le 1,0m“}]}]},{“Typ: „tableRow“, „cells“:[{“Typ: „tableCell“, „content“:[{“Typ: „Span“, Text: „Erwachsene von den Füßen gerissen; alle Fahrzeuge schweben.}]}]},{"Typ: „tableRow“, „cells“:[{“Typ: „tableCell“, „content“:[{“Typ: „span“, „text“: „}]}]},{“Typ: „tableRow“, „cells“:[{“Typ: „tableCell“, „content“:[{“Typ: „Spanne“, Text: „[H5: Strukturelle Schäden]: D \\times V \\le 4,0 \m²/s, h \\le 4,0m“}]}]},{“Typ: „tableRow“, „cells“:[{“Typ: „tableCell“, „content“:[{“Typ: „span“, Text: „Unbewehrte Mauerwerksgebäude werden beschädigt.“}]}]},{“Typ: „tableRow“, „cells“:[{“Typ: „tableCell“, „content“:[{“Typ: „span“, „text“: „}]}]},{“Typ: „tableRow“, „cells“:[{“Typ: „tableCell“, „content“:[{“Typ: „span“, Text: „[H6: Katastrophenkollaps]: D \\times V > 4,0 \m²/s“}]}]},{“Typ: „tableRow“, „cells“:[{“type: „tableCell“, „content“:[{“_type: „span“, „text“: „Völliges strukturelles Versagen von Wohngebäuden“}]}]}]} -->

Durch den Export dynamischer GIS-Ebenen vonH_1durchH_6Durch Klassifizierungen in 5-Minuten-Schritten stromabwärts des gebrochenen Damms können Notfallmanager erkennen, welche Evakuierungsstraßen weiterhin sicher befahrbar sind und mit welchem genauen Zeitstempel sie in tödliche Todesfallen übergehen.

8. Protokoll zur Einhaltung gesetzlicher Vorschriften für Staudammsicherheitsingenieure

Staudammsicherheitsingenieure, geotechnische Berater und kommunale Versorgungsdirektoren, die im Jahr 2026 Gefahrenbewertungen für Staudammversagen durchführen, müssen diese fünfstufige Compliance-Checkliste befolgen:

1. Durchführung von In-situ-Jet-Erosionstests (JET): Schätzen Sie niemals Erodierbarkeitskoeffizienten (k_d) allein aus Lehrbuch-Bodenklassifizierungskarten. Führen Sie standardisierte ASTM D5852/D6665-Tests durch, um die tatsächlichen Scherüberschussparameter zu kalibrieren.
2. Modell der dynamischen geotechnischen Durchbruchserweiterung: Statische trapezförmige Annahmen verwerfen. Setzen Sie gekoppelte, physikalisch basierte Erosionsalgorithmen ein, die die Rückwärtsbewegung des Kopfschnitts, den seitlichen Rotationskollaps und die Absenkung des Reservoirpools modellieren.
3. Stellen Sie GPU-beschleunigte 2D-Stoßerfassungslöser bereit: Leiten Sie Durchbruchsausflusswellen über unstrukturierte 2D-Netze, die aus hochauflösendem LiDAR abgeleitet wurden, und stellen Sie sicher, dass kritische stromabwärts gelegene Taltopographie, Straßendämme und Durchlässe hydrokonditioniert sind.
4. Quantifizierung der Unsicherheit geotechnischer Parameter: Führen Sie mindestens 250 Monte-Carlo-Ensemble-Läufe über geschichtete Latin-Hypercube-Parameterverteilungen durch, um vertretbare P10-, P50- und P90-Gefahrenkonturen zu generieren.
5. Klassifizieren Sie stromabwärts gelegene Evakuierungswege anhand von Tiefengeschwindigkeitskurven: Hydrodynamische Ausgänge in Standard umwandelnH_1 - H_6Gefahrenklassifizierungen für Menschen und Fahrzeuge, um den Notfallmanagementleitern umsetzbare, lebensrettende Informationen zu liefern.

Durch die disziplinierte Integration geotechnischer Bodenmechanik, GPU-beschleunigter nicht-Newtonscher Flachwasserhydrodynamik und probabilistischer Monte-Carlo-Unsicherheitsquantifizierung versetzt die moderne Dammbruchsimulation Bauingenieure in die Lage, die Gesellschaft vor den katastrophalsten Ausfallszenarien zu schützen, die man sich vorstellen kann.