La surverse d’une digue fluviale (levée) peut conduire au développement d’une brèche par érosion externe, suivie d’une rupture brutale de la digue. Cela provoquerait une vague de submersion dans la plaine protégée, laquelle peut présenter des enjeux humains, économiques et financiers majeurs. La gestion et la prévention du risque d’inondation passe par une détermination précise de l’aléa. Pour ce faire, il est nécessaire d’avoir une estimation fiable du débit au travers de la brèche, donc du mécanisme de formation de la brèche et de sa dynamique d’expansion. Les approches existantes sont souvent adaptées pour les digues frontales (barrage et remblai en terre) soumises aux surverses. La transposition de ces approches pour les digues fluviales demeure peu fiable. Les processus qui régissent la formation des brèches dans les digues fluviales restent donc encore du domaine de la recherche. Un programme expérimental visant à améliorer notre compréhension des processus physiques qui régissent la rupture graduelle des digues fluviales par surverse a été mené conjointement par le Laboratoire National d’Hydraulique et Environnement (LNHE) de la division R&D d’EDF et le groupe de recherche Hydraulics in Environmental and Civil Engineering (HECE) de l’Université de Liège. Les travaux ont été conduits sur deux dispositifs expérimentaux distincts, chacun constitué d’un canal principal et d’une plaine d’inondation, séparés par une digue fluviale. Nous nous sommes focalisés sur les surverses localisées de digues homogènes non-cohésives. Une métrologie adaptée, incluant la mesure détaillée de l’évolution de la géométrie de la brèche en continu, par une technique non intrusive (profilométrie laser), a été développée et exploitée dans les travaux de cette thèse. Les tests, réalisés sous conditions contrôlées, ont permis d’investiguer l’évolution de la bèche et des débits sortants pour différentes conditions hydrauliques (débits d’entrée dansle canal principal, régulation du débit sortant en aval du canal principal, confinement de la plaine inondable). Les effets des dimensions du canal principal, de la taille des sédiments et de la cohésion apparente ou encore de la mobilité des fonds au pied de la digue ont fait également l’objet d’étude. En exploitant les mesures, l’évolution des écoulements au voisinage de la brèche a été simulée avec le code de calcul bidimensionnel TELEMAC-2D, permettant d’évaluer les performances de ce code pour des cas d’écoulements, en rupture de digue fluviale, hautement transitoires et complexes. Le couplage avec le code morpho dynamique 2-D SISYPHE a permis d’apprécier l’apport d’une modélisation hydro-morpho dynamique détaillée à l’étude des brèches dans les digues fluviales / Overtopping of fluvial dikes (dykes or embankment levees) can promote external erosion, leading to the initiation of breaching and potentially brutal dike failure and inundation of the protected area. This can generate major human, economic, and financial losses. Flood risk management and prevention require precise hazard quantification. Accurate estimate of the flow through the breach is paramount, for which a precise understanding of the breach formation and expansion is required. Existing methods are often the result of investigations conducted on overtopping of frontal dikes (embankment dams). The application of such approaches to fluvial dikes is not reliable and processes underpinning breach expansion are still under research. An innovative experimental program was conducted to fill this gap by investigating the physical processes involved in overtopping induced fluvial dike gradual breaching. Experiments were conducted in the framework of collaboration between the National Laboratory for Hydraulics and Environment (LNHE) of the R&D division of EDF and the research team Hydraulics in Environmental and Civil Engineering (HECE) of University of Liège. Experiments were conducted on two distinct experimental setups, each consisting of a main channel and floodplain area separated by an erodible fluvial dike. The focus was made on overtopping induced spatial erosion of homogenous, non-cohesive dikes. Measurements included continuous scanning of the dike geometry using a non-intrusive method (Laser Profilometry Technique), which was designed and developed specifically for the present works. Tests conducted under controlled flow and dike configurations allowed assessing the effects of channel inflow discharge, downstream channel regulation system, and floodplain confinement on the breach development and outflow. Effects of main channel size, dike material size, apparent cohesion, and bottom erodibility were studied as well. Using the experimental data, the flow features near the breach area was simulated using the two-dimensional depth-averaged hydrodynamic code TELEMAC-2D, which allowed assessing the performance of the code for highly transient and complex flows such as involved in dike breaching. Coupling TELEMAC-2D with the morphodynamic model SISYPHE enabled investigating the interest of a detailed hydro-morphodynamic modeling for fluvial dike breaching studies

River reaches protected by dikes exhibit high damage potential due to strong value accumulation in the hinterland areas. While providing an efficient protection against low magnitude flood events, dikes may fail under the load of extreme water levels and long flood durations. Hazard and risk assessments for river reaches protected by dikes have not adequately considered the fluvial inundation processes up to now. Particularly, the processes of dike failures and their influence on the hinterland inundation and flood wave propagation lack comprehensive consideration. This study focuses on the development and application of a new modelling system which allows a comprehensive flood hazard assessment along diked river reaches under consideration of dike failures. The proposed Inundation Hazard Assessment Model (IHAM) represents a hybrid probabilistic-deterministic model. It comprises three models interactively coupled at runtime. These are: (1) 1D unsteady hydrodynamic model of river channel and floodplain flow between dikes, (2) probabilistic dike breach model which determines possible dike breach locations, breach widths and breach outflow discharges, and (3) 2D raster-based diffusion wave storage cell model of the hinterland areas behind the dikes. Due to the unsteady nature of the 1D and 2D coupled models, the dependence between hydraulic load at various locations along the reach is explicitly considered. The probabilistic dike breach model describes dike failures due to three failure mechanisms: overtopping, piping and slope instability caused by the seepage flow through the dike core (micro-instability). The 2D storage cell model driven by the breach outflow boundary conditions computes an extended spectrum of flood intensity indicators such as water depth, flow velocity, impulse, inundation duration and rate of water rise. IHAM is embedded in a Monte Carlo simulation in order to account for the natural variability of the flood generation processes reflected in the form of input hydrographs and for the randomness of dike failures given by breach locations, times and widths. The model was developed and tested on a ca. 91 km heavily diked river reach on the German part of the Elbe River between gauges Torgau and Vockerode. The reach is characterised by low slope and fairly flat extended hinterland areas. The scenario calculations for the developed synthetic input hydrographs for the main river and tributary were carried out for floods with return periods of T = 100, 200, 500, 1000 a. Based on the modelling results, probabilistic dike hazard maps could be generated that indicate the failure probability of each discretised dike section for every scenario magnitude. In the disaggregated display mode, the dike hazard maps indicate the failure probabilities for each considered breach mechanism. Besides the binary inundation patterns that indicate the probability of raster cells being inundated, IHAM generates probabilistic flood hazard maps. These maps display spatial patterns of the considered flood intensity indicators and their associated return periods. Finally, scenarios of polder deployment for the extreme floods with T = 200, 500, 1000 were simulated with IHAM. The developed IHAM simulation system represents a new scientific tool for studying fluvial inundation dynamics under extreme conditions incorporating effects of technical flood protection measures. With its major outputs in form of novel probabilistic inundation and dike hazard maps, the IHAM system has a high practical value for decision support in flood management. / Entlang eingedeichter Flussabschnitte kann das Hinterland ein hohes Schadenspotential, aufgrund der starken Akkumulation der Werte, aufweisen. Obwohl Deiche einen effizienten Schutz gegen kleinere häufiger auftretende Hochwässer bieten, können sie unter der Last hoher Wasserstände sowie langer Anstaudauer versagen. Gefährdungs- und Risikoabschätzungsmethoden für die eingedeichten Flussstrecken haben bisher die fluvialen Überflutungsprozesse nicht hinreichend berücksichtigt. Besonders, die Prozesse der Deichbrüche und deren Einfluss auf Überflutung im Hinterland und Fortschreiten der Hochwasserwelle verlangen eine umfassende Betrachtung. Die vorliegende Studie setzt ihren Fokus auf die Entwicklung und Anwendung eines neuen Modellierungssystems, das eine umfassende Hochwassergefährdungsanalyse entlang eingedeichter Flussstrecken unter Berücksichtigung von Deichbrüchen ermöglicht. Das vorgeschlagene Inundation Hazard Assessment Model (IHAM) stellt ein hybrides probabilistisch-deterministisches Modell dar. Es besteht aus drei laufzeitgekoppelten Modellen: (1) einem 1D instationären hydrodynamisch-numerischen Modell für den Flussschlauch und die Vorländer zwischen den Deichen, (2) einem probabilistischen Deichbruchmodell, welches die möglichen Bruchstellen, Breschenbreiten und Breschenausflüsse berechnet, und (3) einem 2D raster-basierten Überflutungsmodell für das Hinterland, das auf dem Speiherzellenansatz und der Diffusionswellengleichung basiert ist. Das probabilistische Deichbruchmodell beschreibt Deichbrüche, die infolge von drei Bruchmechanismen auftreten: dem Überströmen, dem Piping im Deichuntergrund und dem Versagen der landseitigen Böschung als Folge des Sickerflusses und der Erosion im Deichkörper (Mikro-Instabilität). Das 2D Speicherzellenmodell, angetrieben durch den Breschenausfluss als Randbedingung, berechnet ein erweitertes Spektrum der Hochwasserintensitätsindikatoren wie: Überflutungstiefe, Fliessgeschwindigkeit, Impuls, Überflutungsdauer und Wasseranstiegsrate. IHAM wird im Rahmen einer Monte Carlo Simulation ausgeführt und berücksichtigt die natürliche Variabilität der Hochwasserentstehungsprozesse, die in der Form der Hydrographen und deren Häufigkeit abgebildet wird, und die Zufälligkeit des Deichversagens, gegeben durch die Lokationen der Bruchstellen, der Zeitpunkte der Brüche und der Breschenbreiten. Das Modell wurde entwickelt und getestet an einem ca. 91 km langen Flussabschnitt. Dieser Flussabschnitt ist durchgängig eingedeicht und befindet sich an der deutschen Elbe zwischen den Pegeln Torgau und Vockerode. Die Szenarioberechnungen wurden von synthetischen Hydrographen für den Hauptstrom und Nebenfluss angetrieben, die für Hochwässer mit Wiederkehrintervallen von 100, 200, 500, und 1000 Jahren entwickelt wurden. Basierend auf den Modellierungsergebnissen wurden probabilistische Deichgefährdungskarten generiert. Sie zeigen die Versagenswahrscheinlichkeiten der diskretisierten Deichabschnitte für jede modellierte Hochwassermagnitude. Die Deichgefährdungskarten im disaggregierten Darstellungsmodus zeigen die Bruchwahrscheinlichkeiten für jeden betrachteten Bruchmechanismus. Neben den binären Überflutungsmustern, die die Wahrscheinlichkeit der Überflutung jeder Rasterzelle im Hinterland zeigen, generiert IHAM probabilistische Hochwassergefährdungskarten. Diese Karten stellen räumliche Muster der in Betracht gezogenen Hochwasserintensitätsindikatoren und entsprechende Jährlichkeiten dar. Schließlich, wurden mit IHAM Szenarien mit Aktivierung vom Polder bei extremen Hochwässern mit Jährlichkeiten von 200, 500, 1000 Jahren simuliert. Das entwickelte IHAM Modellierungssystem stellt ein neues wissenschaftliches Werkzeug für die Untersuchung fluvialer Überflutungsdynamik in extremen Hochwassersituationen unter Berücksichtigung des Einflusses technischer Hochwasserschutzmaßnahmen dar. Das IHAM System hat eine hohe praktische Bedeutung für die Entscheidungsunterstützung im Hochwassermanagement aufgrund der neuartigen Deichbruch- und Hochwassergefährdungskarten, die das Hauptprodukt der Simulationen darstellen.