TY - THES A1 - Ziegler, Moritz O. T1 - The 3D in-situ stress field and its changes in geothermal reservoirs T1 - Das 3D in-situ Spannungsfeld und seine Änderungen in Geothermiereservoiren N2 - Information on the contemporary in-situ stress state of the earth’s crust is essential for geotechnical applications and physics-based seismic hazard assessment. Yet, stress data records for a data point are incomplete and their availability is usually not dense enough to allow conclusive statements. This demands a thorough examination of the in-situ stress field which is achieved by 3D geomechanicalnumerical models. However, the models spatial resolution is limited and the resulting local stress state is subject to large uncertainties that confine the significance of the findings. In addition, temporal variations of the in-situ stress field are naturally or anthropogenically induced. In my thesis I address these challenges in three manuscripts that investigate (1) the current crustal stress field orientation, (2) the 3D geomechanical-numerical modelling of the in-situ stress state, and (3) the phenomenon of injection induced temporal stress tensor rotations. In the first manuscript I present the first comprehensive stress data compilation of Iceland with 495 data records. Therefore, I analysed image logs from 57 boreholes in Iceland for indicators of the orientation of the maximum horizontal stress component. The study is the first stress survey from different kinds of stress indicators in a geologically very young and tectonically active area of an onshore spreading ridge. It reveals a distinct stress field with a depth independent stress orientation even very close to the spreading centre. In the second manuscript I present a calibrated 3D geomechanical-numerical modelling approach of the in-situ stress state of the Bavarian Molasse Basin that investigates the regional (70x70x10km³) and local (10x10x10km³) stress state. To link these two models I develop a multi-stage modelling approach that provides a reliable and efficient method to derive from the larger scale model initial and boundary conditions for the smaller scale model. Furthermore, I quantify the uncertainties in the models results which are inherent to geomechanical-numerical modelling in general and the multi-stage approach in particular. I show that the significance of the models results is mainly reduced due to the uncertainties in the material properties and the low number of available stress magnitude data records for calibration. In the third manuscript I investigate the phenomenon of injection induced temporal stress tensor rotation and its controlling factors. I conduct a sensitivity study with a 3D generic thermo-hydro-mechanical model. I show that the key control factors for the stress tensor rotation are the permeability as the decisive factor, the injection rate, and the initial differential stress. In particular for enhanced geothermal systems with a low permeability large rotations of the stress tensor are indicated. According to these findings the estimation of the initial differential stress in a reservoir is possible provided the permeability is known and the angle of stress rotation is observed. I propose that the stress tensor rotations can be a key factor in terms of the potential for induced seismicity on pre-existing faults due to the reorientation of the stress field that changes the optimal orientation of faults. N2 - Kenntnis des derzeitigen in-situ Spannungszustandes der Erdkruste ist essenziell für geotechnische Anwendungen und seismische Gefährdungsabschätzungen, welche auf physikalischen Beobachtungen basieren. Jedoch sind die Spannungsinformationen jedes Datenpunktes unvollständig und die Menge an vorhandenen Datenpunkten ist normalerweise nicht groß genug, um schlüssige Ergebnisse zu erzielen. Daher ist eine eingehende Untersuchung des in-situ Spannungsfeldes, welche durch 3D geomechanisch-numerische Modellierung geleistet wird, erforderlich. Jedoch ist die räumliche Auflösung der Modelle begrenzt und der resultierende Spannungszustand ist großen Unsicherheiten unterworfen, welche die Aussagekraft der Ergebnisse beschränken. Zusätzlich gibt es zeitliche Änderungen des Spannungsfeldes, welche durch natürliche Prozesse bedingt oder menschengemacht sind. In meiner Dissertation behandle ich diese Herausforderungen in drei Manuskripten, welche (1) die Orientierung des derzeitigen Spannungszustandes, (2) die 3D geomechanisch-numerische Modellierung des in-situ Spannungszustandes und (3) das Phänomen injektionsinduzierter zeitlicher Rotationen des Spannungstensors zum Thema haben. In dem ersten Manuskript präsentiere ich die erste umfassende Spannungsdatensammlung von Island mit insgesamt 495 Einträgen. Dafür habe ich Bilddatenlogs aus 57 Bohrlöchern in Island auf Indikatoren der maximalen horizontalen Spannungsorientierung hin untersucht. Diese Studie ist die erste ganzheitliche Spannungsuntersuchung, welche sich auf verschiedene Indikatoren der Spannungsorientierung stützt und in einem geologisch sehr jungen und tektonisch aktiven Gebiet auf einem Mittelozeanischen Rücken an Land liegt. Es zeigt sich, dass selbst sehr nahe an der Plattengrenze eine tiefenunabhängige, eindeutige Spannungsorientierung existiert. In dem zweiten Manuskript präsentiere ich einen kalibrierten 3D geomechanisch-numerischen Modellierungsansatz des in-situ Spannungszustandes des bayrischen Molassebeckens welches den regionalen (70x70x10km³) und den lokalen (10x10x10km³) Spannungszustand untersucht. Um diese zwei Modelle zu verbinden, habe ich ein mehrstufigen Modellansatz entworfen, welcher eine zuverlässige und effiziente Methode darstellt um Randbedingungen und Initialbedingungen für das kleinere Modell aus dem größeren Modell abzuleiten. Des Weiteren quantifiziere ich die Unsicherheiten in den Modellergebnissen, welche im Allgemeinen durch geomechanisch-numerische Modellierung und im Speziellen durch den Mehrstufenansatz entstehen. Ich zeige, dass die Signifikanz der Modellergebnisse hauptsächlich durch die Unsicherheiten in den Materialeigenschaften sowie der geringen Anzahl vorhandener Spannungsmagnitudendaten zur Kalibrierung reduziert wird. In dem dritten Manuskript untersuche ich das Phänomen injektionsinduzierter zeitlicher Rotationen des Spannungstensors und deren kontrollierende Parameter. Ich führe eine Sensitivitätsanalyse mit einem generischen 3D thermo-hydro-mechanischen Modell durch. Darin zeige ich, dass die Schlüsselparameter, welche die Rotationen des Spannungstensors kontrollieren, die Permeabilität des Reservoirgesteins als der entscheidende Faktor, die Injektionsrate und die initiale Differenzspannung sind. Insbesondere für geothermische Systeme nach dem Hot-Dry-Rock-Verfahren mit einer geringen Permeabilität weisen die Ergebnisse auf große Rotationen des Spannungstensors hin. Gemäß diesen Ergebnissen kann die initiale Differenzspannung in einem Reservoir abgeschätzt werden, sollte die Permeabilität bekannt sein und der Winkel der Spannungsrotation beobachtet werden. Ich schlage vor, dass Spannungsrotationen ein Schlüsselfaktor in Bezug auf das Potenzial für induzierte Seismizität sind, welche auf prä-existierenden Störungen entsteht, die durch die Reorientierung des Spannungsfeldes optimal orientiert werden. KW - stress KW - stress changes KW - induced seismicity KW - geothermal KW - geomechanical modelling KW - Spannung KW - Spannungsänderungen KW - induzierte Seismizität KW - Geothermie KW - geomechanische Modellierung Y1 - 2017 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus4-403838 ER - TY - THES A1 - Sen, Ali Tolga T1 - Inversion of seismic source parameters for weak mining-induced and natural earthquakes T1 - Inversion von seismischen Quellparametern für schwache Bergbau-induzierte und natürliche Erdbeben N2 - The purpose of this thesis is to develop an automated inversion scheme to derive point and finite source parameters for weak earthquakes, here intended with the unusual meaning of earthquakes with magnitudes at the limit or below the bottom magnitude threshold of standard source inversion routines. The adopted inversion approaches entirely rely on existing inversion software, the methodological work mostly targeting the development and tuning of optimized inversion flows. The resulting inversion scheme is tested for very different datasets, and thus allows the discussion on the source inversion problem at different scales. In the first application, dealing with mining induced seismicity, the source parameters determination is addressed at a local scale, with source-sensor distance of less than 3 km. In this context, weak seismicity corresponds to event below magnitude MW 2.0, which are rarely target of automated source inversion routines. The second application considers a regional dataset, namely the aftershock sequence of the 2010 Maule earthquake (Chile), using broadband stations at regional distances, below 300 km. In this case, the magnitude range of the target aftershocks range down to MW 4.0. This dataset is here considered as a weak seismicity case, since the analysis of such moderate seismicity is generally investigated only by moment tensor inversion routines, with no attempt to resolve source duration or finite source parameters. In this work, automated multi-step inversion schemes are applied to both datasets with the aim of resolving point source parameters, both using double couple (DC) and full moment tensor (MT) models, source duration and finite source parameters. A major result of the analysis of weaker events is the increased size of resulting moment tensor catalogues, which interpretation may become not trivial. For this reason, a novel focal mechanism clustering approach is used to automatically classify focal mechanisms, allowing the investigation of the most relevant and repetitive rupture features. The inversion of the mining induced seismicity dataset reveals the repetitive occurrence of similar rupture processes, where the source geometry is controlled by the shape of the mined panel. Moreover, moment tensor solutions indicate a significant contribution of tensile processes. Also the second application highlights some characteristic geometrical features of the fault planes, which show a general consistency with the orientation of the slab. The additional inversion for source duration allowed to verify the empirical correlation for moment normalized earthquakes in subduction zones among a decreasing rupture duration with increasing source depth, which was so far only observed for larger events. N2 - Ziel dieser Dissertation ist es ein automatisches Inversionsschema zur Ableitung von Parametern für Punkt- und finite Quellen von "schwachen" Erdbeben im Sinne von aussergewöhnlich kleinen Erdbeben an oder unterhalb der Magnitudenschwelle für Standard-Inversionsroutinen zu entwickeln. Der verwendete Inversionsansatz baut auf bestehender Inversionssoftware auf, die methodische Arbeit zielt hauptsächlich auf die Entwicklung und das Tuning des optimierten Inversionsablaufs. Das resultierende Inversionsverfahren ist für verschiedene Datensätze getestet und erlaubt somit eine Diskussion über das Herdinversionsproblem auf unterschiedlichen Skalierungen. In der ersten Anwendung, die sich mit induzierter Seismizität in Minen beschäftigt, werden die Herdparameter im lokalen Bereich mit einer Herd-Empfänger Entfernung von weniger als 3 km bestimmt. In diesem Zusammenhang bezieht sich die "schwache" Seismizität auf Ereignisse mit Magnituden unter Mw 2.0, welche sehr selten das Ziel von automatisierten Herdinversionsroutinen sind. Die zweite Anwendung berücksichtigt den regionalen Datensatz der Nachbebesequenz des Maule (Chile) Erdbebens von 2010, aufgezeichnet durch Breitband-Stationen in einer maximalen Entfernung von 300 km. In diesem Fall umfasst der Magnitudenbereich Ereignisse mit Magnituden kleiner als Mw 4.0. Auch dieser Datensatz wird hier als "schwache" Seismizität interpretiert, da die Analyse solch moderater Seismizität normalerweise nur durch Routinen für Momententensorinversion erfolgt, welche keinerlei Auflösung für die Herddauer oder finite Herdparameter beinhaltet. In dieser Arbeit wird ein automatisiertes, mehrstufiges Inversionsverfahren auf beide Datensätze angewendet mit dem Ziel sowohl Punktquellen-Herdparameter für double-couple (DC) und volle Momententensor (MT) Modelle als auch für Herddauer und finite Herdparameter zu bestimmen. Ein Hauptergebnis dieser Analyse von schwachen Ereignissen ist die Vergrößerung von Momententensor-Katalogen, deren Interpretation wichtig werden könnte. Aus diesem Grund wurde ein neuartiger Cluster-Ansatz für Herdmechanismen angewendet um diese automatisch zu klassifizieren, der somit die Untersuchung der wichtigsten und sich wiederholenden Bruchprozessen erlaubt. Die Inversion des Datensatzes von induzierter Seismizität in Minen enthüllt das sich wiederholende Auftreten von ähnlichen Bruchprozessen, wobei die Herdgeometrie von der Form der Minentäfelung kontrolliert wird. Darüber hinaus deuten die Momententensorlösungen einen signifikanten Beitrag von tensilen Prozessen an. Auch die zweite Anwendung zeigt einige charackteristische, geometrische Besonderheiten der Bruchflächen auf, welche prinzipiell konsistent zur Orientierung der subduzierten Platte sind. Die zusätzliche Inversion nach der Herdzeit erlaubt die Verifizierung von empirischen Korrelationen von momentennormalisierten Erdbeben in Subduktionszonen zwischen einer abnehmenden Bruchdauer mit zunehmender Herdtiefe, welche bisher nur für größere Ereignisse beobachtet wurde. KW - Seismologie KW - induzierten Seismizität KW - Quellenumkehr KW - Momententensor KW - Aufenthaltsdauer KW - seismology KW - induced seismicity KW - source inversion KW - moment tensor KW - source duration Y1 - 2014 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus-71914 ER - TY - THES A1 - Niemz, Peter T1 - Imaging and modeling of hydraulic fractures in crystalline rock via induced seismic activity T1 - Charakterisierung und Modellierung hydraulischer Brüche in Kristallingestein mit Hilfe induzierter Seismizität N2 - Enhanced geothermal systems (EGS) are considered a cornerstone of future sustainable energy production. In such systems, high-pressure fluid injections break the rock to provide pathways for water to circulate in and heat up. This approach inherently induces small seismic events that, in rare cases, are felt or can even cause damage. Controlling and reducing the seismic impact of EGS is crucial for a broader public acceptance. To evaluate the applicability of hydraulic fracturing (HF) in EGS and to improve the understanding of fracturing processes and the hydromechanical relation to induced seismicity, six in-situ, meter-scale HF experiments with different injection schemes were performed under controlled conditions in crystalline rock in a depth of 410 m at the Äspö Hard Rock Laboratory (Sweden). I developed a semi-automated, full-waveform-based detection, classification, and location workflow to extract and characterize the acoustic emission (AE) activity from the continuous recordings of 11 piezoelectric AE sensors. Based on the resulting catalog of 20,000 AEs, with rupture sizes of cm to dm, I mapped and characterized the fracture growth in great detail. The injection using a novel cyclic injection scheme (HF3) had a lower seismic impact than the conventional injections. HF3 induced fewer AEs with a reduced maximum magnitude and significantly larger b-values, implying a decreased number of large events relative to the number of small ones. Furthermore, HF3 showed an increased fracture complexity with multiple fractures or a fracture network. In contrast, the conventional injections developed single, planar fracture zones (Publication 1). An independent, complementary approach based on a comparison of modeled and observed tilt exploits transient long-period signals recorded at the horizontal components of two broad-band seismometers a few tens of meters apart from the injections. It validated the efficient creation of hydraulic fractures and verified the AE-based fracture geometries. The innovative joint analysis of AEs and tilt signals revealed different phases of the fracturing process, including the (re-)opening, growth, and aftergrowth of fractures, and provided evidence for the reactivation of a preexisting fault in one of the experiments (Publication 2). A newly developed network-based waveform-similarity analysis applied to the massive AE activity supports the latter finding. To validate whether the reduction of the seismic impact as observed for the cyclic injection schemes during the Äspö mine-scale experiments is transferable to other scales, I additionally calculated energy budgets for injection experiments from previously conducted laboratory tests and from a field application. Across all three scales, the cyclic injections reduce the seismic impact, as depicted by smaller maximum magnitudes, larger b-values, and decreased injection efficiencies (Publication 3). N2 - Hydraulisch-stimulierte tiefengeothermale Systeme (Enhanced Geothermal systems, EGS) gelten als einer der Eckpfeiler für die nachhaltige Energieerzeugung der Zukunft. In diesen geothermalen Systemen wird heißes Tiefengestein durch Fluidinjektionen unter hohem Druck aufgebrochen, um Wegsamkeiten zur Erwärmung von Wasser oder anderen Fluiden zu schaffen. Beim Aufbrechen werden zwangsläufig kleine seismische Ereignisse ausgelöst (induzierte Seismizität), die in sehr seltenen Fällen an der Oberfläche spürbar sind, jedoch in extremen Fällen auch Schäden verursachen können. Die Kontrolle bzw. die Reduzierung der seismischen Aktivität in EGS ist daher ein entscheidender Punkt, damit diese Art der Energiegewinnung eine breite gesellschaftliche Akzeptanz findet. Grundlage dieser Dissertation ist eine Serie von kontrollierten, hydraulischen Bruchexperimenten mit Bruchdimensionen von einigen Metern. Die Experimente wurden in einer Tiefe von 410 m in kristallinem Gestein eines Versuchsbergwerks (Äspö Hard Rock Laboratory, Schweden) mit unterschiedlichen Injektionsstrategien durchgeführt. Die detaillierte Auswertung der Bruch-Experimente in dieser Dissertation zielt darauf ab, die Nutzbarkeit von hydraulischen Stimulationen (hydraulic fracturing, HF) in EGS zu untersuchen und das Verständnis von Bruchprozessen sowie der hydromechanischen Beziehung zur induzierten Seismizität zu verbessern. Um die Schallemissionsaktivität (acoustic emissions, AE), die durch 11 piezoelektrische AE-Sensoren kontinuierlich aufgezeichnet wurde, zu extrahieren und zu charakterisieren, wurde ein halbautomatischer, wellenformbasierter Detektions-, Klassifizierungs- und Lokalisierungsworkflow entwickelt. Mit Hilfe des resultierenden Katalogs von 20000 AEs wurde das Bruchwachstum detailliert kartiert und charakterisiert. Das Experiment mit der neuartigen, zyklischen Injektionsstrategie (HF3) weist einen geringeren seismischen Fußabdruck auf als die Standard-Injektionsstrategie. HF3 induzierte weniger AEs und eine kleinere Maximalmagnitude. Außerdem hatte das Experiment einen signifikant höheren b-Wert, was einer verringerten Anzahl von großen AEs relativ zur Anzahl der kleineren AEs entspricht. Darüber hinaus zeigte HF3 eine erhöhte Komplexität im Bruchmuster mit mehreren Brüchen bzw. einem Netzwerk von Brüchen. Im Gegensatz dazu entwickelten die Standard-Injektionen einzelne, ebene Bruchzonen (Publikation 1). Zusätzlich zu den induzierten AEs wurden transiente, langperiodische Signale auf den horizontalen Komponenten von zwei Breitband-Seismometern, die wenige Meter von den Brüchen installiert waren, ausgewertet. Diese Signale wurden als Neigungssignale interpretiert und mit modellierten Neigungssignalen verglichen. Der Vergleich zeigt unabhängig, dass hydraulische Brüche geöffnet wurden und bestätigt, dass die AE-basierte Analyse die Bruchgeometrie verlässlich kartieren kann. Die gemeinsame Betrachtung von AEs und Neigungssignalen offenbart verschiedene Phasen des Bruchprozesses: das (wiederholte) Öffnen des Bruches, das Bruchwachstum und das weitere Wachsen des Bruches nach dem Ende der Injektion. Außerdem liefert die Analyse Hinweise auf die Reaktivierung einer natürlichen Bruchzone in einem der Experimente (Publikation 2). Eine neuentwickelte und hier präsentierte Wellenform-Ähnlichkeitsanalyse, die Informationen des gesamten Sensornetzwerkes nutzt und zum ersten Mal auf einen umfangreichen AE-Katalog angewendet wurde, unterstützt diese Interpretation. Um zu validieren, ob die verringerte Seismizität während der zyklischen Injektion auf der Meter-Skala (Bergwerk) auf andere Maßstäbe übertragbar ist, wurden Energie-Budgets für Injektionsexperimente aus zuvor durchgeführten Laborversuchen und aus einem Tiefengeothermie-Projekt berechnet. Über alle drei Skalen hinweg zeigen die zyklischen Injektionen einen verringerten seismischen Fußabdruck mit kleineren Maximalmagnituden, größeren b-Werte und einem kleineren Verhältnis von seismisch-abgestrahlter zu injizierter Energie (Publikation 3). KW - induced seismicity KW - hydraulic fracturing KW - enhanced geothermal systems (EGS) KW - injection KW - deformation KW - acoustic emissions KW - fracture growth KW - injection scheme KW - basement rock KW - Schallemissionen KW - Grundgestein KW - Deformation KW - verbesserte geothermische Systeme KW - Bruchausbreitung KW - hydraulisches Aufbrechen KW - Induzierte Seismizität KW - Injektion KW - Injektionsschema Y1 - 2022 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus4-556593 ER - TY - THES A1 - Grigoli, Francesco T1 - Automated seismic event location by waveform coherence analysis T1 - Die automatische Lokalisierung von Mikroerdbeben durch die Wellenform Kohärenzanalyse N2 - Automated location of seismic events is a very important task in microseismic monitoring operations as well for local and regional seismic monitoring. Since microseismic records are generally characterised by low signal-to-noise ratio, such methods are requested to be noise robust and sufficiently accurate. Most of the standard automated location routines are based on the automated picking, identification and association of the first arrivals of P and S waves and on the minimization of the residuals between theoretical and observed arrival times of the considered seismic phases. Although current methods can accurately pick P onsets, the automatic picking of the S onset is still problematic, especially when the P coda overlaps the S wave onset. In this thesis I developed a picking free automated method based on the Short-Term-Average/Long-Term-Average (STA/LTA) traces at different stations as observed data. I used the STA/LTA of several characteristic functions in order to increase the sensitiveness to the P wave and the S waves. For the P phases we use the STA/LTA traces of the vertical energy function, while for the S phases, we use the STA/LTA traces of the horizontal energy trace and then a more optimized characteristic function which is obtained using the principal component analysis technique. The orientation of the horizontal components can be retrieved by robust and linear approach of waveform comparison between stations within a network using seismic sources outside the network (chapter 2). To locate the seismic event, we scan the space of possible hypocentral locations and origin times, and stack the STA/LTA traces along the theoretical arrival time surface for both P and S phases. Iterating this procedure on a three-dimensional grid we retrieve a multidimensional matrix whose absolute maximum corresponds to the spatial and temporal coordinates of the seismic event. Location uncertainties are then estimated by perturbing the STA/LTA parameters (i.e the length of both long and short time windows) and relocating each event several times. In order to test the location method I firstly applied it to a set of 200 synthetic events. Then we applied it to two different real datasets. A first one related to mining induced microseismicity in a coal mine in the northern Germany (chapter 3). In this case we successfully located 391 microseismic event with magnitude range between 0.5 and 2.0 Ml. To further validate the location method I compared the retrieved locations with those obtained by manual picking procedure. The second dataset consist in a pilot application performed in the Campania-Lucania region (southern Italy) using a 33 stations seismic network (Irpinia Seismic Network) with an aperture of about 150 km (chapter 4). We located 196 crustal earthquakes (depth < 20 km) with magnitude range 1.1 < Ml < 2.7. A subset of these locations were compared with accurate locations retrieved by a manual location procedure based on the use of a double difference technique. In both cases results indicate good agreement with manual locations. Moreover, the waveform stacking location method results noise robust and performs better than classical location methods based on the automatic picking of the P and S waves first arrivals. N2 - Die automatische Lokalisierung seismischer Ereignisse ist eine wichtige Aufgabe, sowohl im Bereich des Mikroseismischen Monitorings im Bergbau und von Untegrund Aktivitäten,  wie auch für die lokale und regionale Überwachung von natürlichen Erdbeben. Da mikroseismische Datensätze häufig ein schlechtes Signal-Rausch-Verhältnis haben müssen die Lokalisierungsmethoden robust gegen Rauschsignale und trotzdem hinreichend genau sein. Aufgrund der in der Regel sehr hochfrequent aufgezeichneten Messreihen und der dadurch sehr umfangreichen Datensätze sind automatische Auswertungen erstrebenswert. Solche Methoden benutzen in der Regel automatisch gepickte und  den P und S Phasen zugeordnete Ersteinsätze und Minimieren die Summe der quadratischen Zeitdifferenz zwischen den beobachteten und theoretischen Einsatzzeiten. Obgleich das automatische Picken der P Phase in der Regel sehr genau möglich ist, hat man beim Picken der S Phasen häufig Probleme, z.B. wenn die Coda der P Phase sehr lang ist und in den Bereich der S Phase hineinreicht. In dieser Doktorarbeit wird eine Methode vollautomatische, Wellenform-basierte Lokalisierungsmethode entwickelt, die Funktionen des Verhältnisses "Short Term Average / Long Term Average"  (STA/LTA) verwendet und keine Pickzeiten invertiert. Die STA/LTA charakteristische Funktion wurde für unterschiedliche Wellenform Attribute getestet, um die Empfindlichkeit für P und S Phasen zu erhöhen. Für die P Phase wird die STA/LTA Funktion für die Energie der Vertikalkomponente der Bodenbewegung benutzt, wohingegen für die S Phase entweder die Energie der horizontalen Partikelbewegung oder eine optimierte Funktion auf Basis der Eigenwertzerlegung benutzt wird. Um die Ereignisse zu lokalisieren wird eine Gittersuche über alle möglichen Untergrundlokalisierungen durchgeführt. Für jeden räumlichen und zeitlichen Gitterpunkt werden die charakteristischen Funktionen entlang der theoretischen Einsatzkurve aufsummiert. Als Ergebnis erhält man eine 4-dimensionale Matrix über Ort und Zeit des Ereignisses, deren Maxima die wahrscheinlichsten Lokalisierungen darstellen. Um die Unsicherheiten der Lokalisierung abzuschätzen wurden die Parameter der STA/LTA Funktionen willkürlich verändert und das Ereignis relokalisiert. Die Punktwolke aller möglichen Lokalisierungen gibt ein Maß für die Unsicherheit des Ergebnisses. Die neu entwickelte Methode wurde an einem synthetischen Datensatz von 200 Ereignissen getestet und für zwei beobachtete Datensätze demonstriert. Der erste davon betrifft induzierte Seismizität in einem Kohlebergbau in Norddeutschland. Es wurden 391 Mikrobeben mit Magnituden zwischen Ml 0.5 und 2.0 erfolgreich lokalisiert und durch Vergleich mit manuell ausgewerteten Lokalisierungen verifziert.Der zweite Datensatz stammt von einem Anwednung auf des Regionale Überwachungsnetz in der Region Campania-Lucania (Süditalien)  mit 33 seismischen Stationen und einer Apertur von etwa 150 km. Wir konnten 196 Erdbeben mit Tiefen < 20 km und Magnituden zwischen Ml 1.1 und 2.7 lokalisieren. Eine Untergruppe der eigenen Lokalisierungen wurde mit den Lokalisierungen einer Standard Lokalisierung sowie einer  hochgenauen Relativlokalisierung verglichen. In beiden Fällen ist die Übereinstimmung mit den manuellen Lokalisierungen groß. Außerdem finden wir, dass die Wellenform Summations Lokalisierung ronbust gegen Rauschen ist und bessere Ergebnisse liefert als die Standard Lokalisierung, die auf dem automatischen Picken von Ersteinsatzzeiten alleine basiert. KW - earthquake location KW - induced seismicity KW - microseismicity Y1 - 2014 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus-70329 ER -