Die Förderung von Erdgas aus diversen Rotliegend-Gaslagerstätten in den Niederlanden und Norddeutschland hat Erdbeben ausgelöst, die mit einer signifikanten Porendruckabsenkung im Untergrund korrelieren. Für die betroffene Region in Norddeutschland wird die Erdbebengefährdung als vernachlässigbar eingeschätzt. Entsprechend haben die induzierten Erdbeben, ungeachtet ihrer sehr geringfügigen Oberflächenauswirkungen, große Besorgnis in der Bevölkerung hervorgerufen. Ein vereinfachter Erklärungsansatz für die induzierten Erdbeben sieht Porendruck-induzierte Spannungsänderungen vor, die mit dem vorherrschenden tektonischen Spannungszustand interferieren und dabei die Integrität von vorhandenen Störungen infrage stellen. Analytische Ansätze beschreiben erstens, wie die Porendruckabnahme zu einer differentiellen Zunahme der effektiven vertikalen und horizontalen Spannungskomponenten innerhalb des Reservoirs führen und zweitens, wie die kompressiven Verformungstendenzen des Reservoirs Kompensationsspannungen in umliegenden Gesteinsformationen hervorrufen. Die analytischen Modelle können die Reaktivierung von Störungen durch Porendruckabsenkungen grundsätzlich erklären, realistische mechanische Stratigraphien und die strukturelle Komplexität der Lagerstätten aber nur unzureichend abbilden. Die vorliegende Arbeit abstrahiert grundlegende geologischen Gegebenheiten der norddeutschen Rotliegend-Gasfelder in zweidimensionale Finite Elemente (FE) Modelle. Mithilfe der numerischen Modellierung soll das quantitative Verständnis der induzierten Erdbeben verbessert und der Einfluss verschiedener Parameter, die das Auftreten von induzierten Erdbeben in Norddeutschland potenziell begünstigen, evaluiert werden. Ein Parameterraum, der die Tiefe und Mächtigkeit des Reservoirs, die mechanischen und hydraulischen Modelleigenschaften sowie die Geometrie der Reservoirhorizonte entlang kompartmentbildender Verwerfungen variiert, wurde im Hinblick auf das jeweilige Reaktivierungspotenzial der Störung als Abschiebung analysiert. Charakteristische geologische Merkmale wie die lateral variierende Mächtigkeit von Evaporiten und kompartmentbildende Verwerfungen innerhalb der Gasfelder wurden in drei unterschiedlichen Modell-settings untersucht und der Grabenrand als Referenzszenario herangezogen. Die bodenmechanische Konsolidationsanalyse des FE-Programms ABAQUS wurde für die Lösung der voll gekoppelten poroelastischen Gleichungen verwendet und dafür zunächst ein Modellierungsansatz konzipiert. Hierbei ist der Biot-Koeffizient zentraler Inputparameter und nimmt für Festgesteine normalerweise Werte kleiner eins an. In der numerischen Prozedur war er in der Berechnung von Versagensspannungen und der Ausgabe von Terzaghi-Spannungen zunächst scheinbar nicht berücksichtigt. Es wird aber gezeigt, dass der Biot-Koeffizient den totalen Spannungszustand determiniert und damit die Magnituden beider abgeleiteter Effektivspannungskonzepte bestimmt. Biot-Willis Effektivspannungen beschreiben das Spannungs-Verformungs-Verhalten des poroelastischen Materials in direkter Weise, Terzaghi Effektivspannungen determinieren die Festigkeits- bzw. Fließgrenze von porösen Gesteinen unter hohem Umgebungsdruck. Eine weitere zentrale Bedeutung kommt der deterministischen Darstellung von Störungen als Kontaktflächen im numerischen Modell zu. Kontaktflächen bilden den diskontinuierlichen Charakter von realen Störungen ab und erlauben kleinräumige Relativverschiebungen innerhalb des FE-Netzes. In der poroelastischen Modellierung führt die Repräsentation einer abdichtenden Störung durch Kontakt zusätzliche Kontinuitätsbedingungen ein, die zusammen mit der kontroversen numerischen Versagensspannung die Einsetzung der Störungsbelastungsgröße SSR erforderten. Der Parameter SSR quantifiziert die Belastung der Störung im Sinne des Mohr-Coulomb Kriteriums und erlaubt eine vergleichende Bewertung verschiedener Parametervariationen. Ein zweiter Simulationslauf in ABAQUS explicit analysiert die Dissipation der vorhergehend identifizierten, kritischen Belastungsfälle in Form von dynamischem Versatz auf der Modellstörung. Das implementierte Reibungsgesetz verschreibt der Störung einen abnehmenden Reibungswiderstand mit zunehmender Gleitgeschwindigkeit und bildet in Kombination mit den dynamischen Spannungen Grundzüge der dynamischen Bruchausbreitung bei der seismischen Reaktivierung von Störungen ab. In der Modellierung wurden verschieden Faktoren identifiziert, die das Auftreten von produktions-induzierten Erdbeben begünstigen. Einzelne Parameter, wie eine geringe Tiefenlage des Reservoirs, eine große Reservoir-Mächtigkeit, große Biot-Koeffizienten sowie ein großer Steifigkeitskontrast zwischen Deckgestein und Reservoir begünstigen die Reaktivierung von Randstörungen des Reservoirhorizontes. Viskoelastische Evaporite im Hangenden spielen eine besondere Rolle für die Kritikalität der Störung, da sie eine lateral inhomogene Auflast auf das Reservoir ausüben und die räumliche Dissipation der produktions-induzierten Spannungen einschränken. Im Vergleich zum Grabenrand wurden Verwerfungen, die die Lagerstätte in interne Kompartimente unterteilen, als wahrscheinlichere Lokation für Erdbeben identifiziert. Hier erfahren Störungen, vor allem durch Produktionsszenarien, die den Porendruck entweder in beiden gegeneinander versetzten Kompartimenten oder nur im Liegenden Kompartiment absenken, eine besondere Belastung. Dabei kontrolliert der Versatz der Kompartimente entlang der Störung eine, durch die Geometrie bedingte, differenzielle Kompaktion der Kompartimente. Die bevorzugte Reaktivierung von steil einfallenden Störungen (> 60°) spiegelt in allen Modellen einen essentiellen Beitrag der elastischen Reservoirkompaktion zur kritischen Belastung der Störung wider, die im Kontrast zu tektonischen Modellen steht. Das größte Reaktivierungspotenzial der Grabenstörung wurde in den Modellen im oberen Bereich des Reservoirs und einem schmalen Abschnitt darüber beobachtet, an den kompartmentbildenden Verwerfungen tritt die Bruchinitiierung im unteren Bereich des Liegenden Reservoirhorizontes auf. Die Bruchausbreitung erfolgt entweder abwärts oder gleichzeitig auf- und abwärts entlang der Störung, wobei der übergeordnete heterogene Spannungszustand und die initiale Beschleunigung die Bruchausbreitung entweder begünstigen oder dämpfen können. Der isotrope Spannungszustand im Salz und die Aufschiebungstendenz unterhalb des Reservoirs stellen ultimative geologische Barrieren für die Bruchausbreitung dar. Die Bruchlänge und der maximale Versatz hängen dabei vorwiegend vom tektonischen Spannungszustand und seiner Modifikation durch die Porendruckabsenkung ab und spiegeln nur sekundär die simulierte Porendruckabsenkung wider. Die mögliche Abbildung der Reservoirverformung durch zwei unterschiedliche Effektivspannungskonzepte sowie die starke Sensitivität der Ergebnisse gegenüber verschiedenen Möglichkeiten der Störungsdarstellung heben fachspezifische Herausforderungen der numerischen Modellierung hervor. Aus den Modellierungsergebnissen lässt sich zusammenfassend kein einzelner, isolierter Parameter als Indikator für eine wahrscheinliche produktions-induzierten Reaktivierung von Störungen ableiten, weil das Reaktivierungspotenzial vor allem durch eine komplexe Wechselwirkung aus Anfangsspannungszustand, der Porendruckänderung und der Festigkeit von Störungen determiniert wird. Um die beobachtete Wechselwirkung verschiedener Parameter im Sinne einer Verstärkung oder Abschwächung des Reaktivierungspotenzials zu limitieren, muss der Parameterraum für detaillierte Vorhersagemodelle durch Struktur- und fallspezifische Daten weiter eingeschränkt werden. Nichtsdestotrotz liefert die vorgestellte FE- Modellierung ein grundlegendes Prozess- und Parameterverständnis und ermöglicht die Betrachtung von Szenarien und Grenzfällen in der Planung und dem Betrieb von Gaslagerstätten. Der präsentierte zweistufige Ansatz erlaubt zunächst die Identifizierung von kritischen Störungen und Belastungsfällen und evaluiert in einem zweiten Schritt die potenziellen Konsequenzen des kritischen Belastungsfalles in Form von dynamischer Bruchausbreitung. Hierbei stellt der vorgestellte Modellierungsansatz auch für andere geomechanische Problemstellungen, wie der Integrität von Deckgesteinen, ein vielversprechendes Werkzeug dar.