Results from two-station surface-wave inversions across the Archean Kaapvaal craton of southern Africa are compared with seismic velocities estimated from approximately 100 mantle xenoliths brought to the surface in kimberlite pipes. As the xenoliths represent a snapshot of the mantle at the time of their eruption, comparison with recently recorded seismic data provides an opportunity to compare and contrast the independently gained results. These cratonic xenoliths from the southern Kaapvaal, all less than 100Ma in age, have been analyzed geothermobarometrically to obtain the equilibrium P-T conditions of the cratonic mantle to about 180km depth [James et al 2004]. Seismic velocity-depth and density-depth profiles calculated on the basis of these P-T data and the mineral modes of the xenoliths are used to produce theoretical surface-wave dispersion curves and to generate roughly the upper 200km of a starting/reference model. A regionally-developed crustal structure [Niu and James 2002] was used for the crust and 300km of mantle values taken from PREM filled in down to 500km depth. This composite model was used as the starting/reference model for a Neighbourhood Algorithm surface-wave inversion using fundamental-mode Rayleigh-wave phase velocities for 16 paths within the Kaapvaal Craton from five events. The velocity structures found by that inversion are consistent with those derived from the xenolith data. Hence the velocity structure (i.e. thermal structure) of the mantle to a depth of 180km beneath the Kaapvaal craton is basically the same today as it was 80-90Ma. Further, synthetics runs show that for this surface-wave dataset, there is no strong low-velocity zone at depths shallower than at least 200km. / Master of Science

The seismic structure (Vp, Vp/Vs, and Qp anomalies) contributes to the physical properties and the lithology of rocks and possible fluid distribution in the region. The Vp model images the geometry of the Tarutung and the Sarulla basins. Both basins have a depth of around 2.0 km. High Vp/Vs and high attenuation (low Qp) anomalies are observed along the Sarulla graben associated with a weak zone caused by volcanic activities along the graben. Low Vp/Vs and low conductivity anomalies are found in the west of the Tarutung basin. This anomaly is interpreted as dry, compact, and rigid granitic rock in the region as also found by geological observations. Low Vp, high Vp/Vs and low Qp anomalies are found at the east of the Tarutung basin which appear to be associated with the three big geothermal manifestations in Sipoholon, Hutabarat, and Panabungan area. These anomalies are connected with high Vp/Vs and low Qp anomalies below the Tarutung basin at depth of around 3 - 10 km. This suggests that these geothermal manifestations are fed by the same source of the hot fluid below the Tarutung basin. The hot fluids from below the Tarutung basin propagate to the more dilatational and more permeable zone in the northeast. Granite found in the west of the Tarutung basin could also be abundant underneath the basin at a certain depth so that it prevents the hot fluid to be transported directly to the Tarutung basin. High seismic attenuation and low Vp/Vs anomalies are found in the southwest of the Tarutung basin below the Martimbang volcano. These anomalies are associated with hot rock below the volcano without or with less amount of partial melting. There is no indication that the volcano controls the geothermal system around the Tarutung basin. The geothermal resources around the Tarutung basin is a fault-controlled system as a result of deep circulation of fluids. Outside of the basin, the seismicity delineation and the focal mechanism correlate with the shape and the characteristics of the strike-slip Sumatran fault. Within the Tarutung basin, the seismicity is distributed more broadly which coincides with the margin of the basin. An extensional duplex system in the Tarutung basin is derived from the seismicity and focal mechanism analysis which is also consistent with the geological observations. The vertical distribution of the seismicity suggests the presence of a negative flower structure within the Tarutung basin. / Indonesien zählt zu den weltweit führenden Ländern bei der Nutzung von geothermischer Energie. Die geothermischen Energiequellen sind im Wesentlichen an den aktiven Vulkanismus gebunden, der durch die Prozesse an der indonesischen Subduktionszone verursacht wird. Darüber hinaus sind geotektonische Strukturen wie beispielsweise die Sumatra-Störung als verstärkende Faktoren für das geothermische Potenzial von Bedeutung. Bei der geophysikalischen Erkundung der indonesischen Geothermie-Ressourcen konzentrierte man sich bisher vor allem auf die Magnetotellurik. Passive Seismologie wurde dahingegen ausschließlich für die Überwachung von im Betrieb befindlichen Geothermie-Anlagen verwendet. Jüngste Untersuchungungen z.B. in Island und in den USA haben jedoch gezeigt, dass seismologische Verfahren bereits in der Erkundungsphase wichtige Informationen zu den physikalischen Eigenschaften, zum Spannungsfeld und zu möglichen Fluid- und Wärmetransportwegen liefern können. In der vorgelegten Doktorarbeit werden verschiedene moderne Methoden der passiven Seismologie verwendet, um beispielhaft ein neues, von der indonesischen Regierung für zukünftige geothermische Energiegewinnung ausgewiesenes Gebiet im nördlichen Teil Sumatras (Indonesien) zu erkunden. Die konkreten Ziele der Untersuchungen umfassten (1) die Ableitung von 3D Strukturmodellen der P- und S-Wellen Geschwindigkeiten (Parameter Vp und Vs), (2) die Bestimmung der Absorptionseigenschaften (Parameter Qp), und (3) die Kartierung und Charakterisierung von Störungssystemen auf der Grundlage der Seismizitätsverteilung und der Herdflächenlösungen. Für diese Zwecke habe ich zusammen mit Kollegen ein seismologisches Netzwerk in Tarutung (Sumatra) aufgebaut und über einen Zeitraum von 10 Monaten (Mai 2011 – Februar 2012) betrieben. Insgesamt wurden hierbei 42 Stationen (jeweils ausgestattet mit EDL-Datenlogger, 3-Komponenten, 1 Hz Seismometer) über eine Fläche von etwa 35 x 35 km verteilt. Mit dem Netzwerk wurden im gesamten Zeitraum 2568 lokale Erdbeben registriert. Die integrierte Betrachtung der Ergebnisse aus den verschiedenen Teilstudien (Tomographie, Erdbebenverteilung) erlaubt neue Einblicke in die generelle geologische Stukturierung sowie eine Eingrenzung von Bereichen mit einem erhöhten geothermischen Potenzial. Das tomographische Vp-Modell ermöglicht eine Bestimmung der Geometrie von Sedimentbecken entlang der Sumatra-Störung. Für die Geothermie besonders interessant ist der Bereich nordwestlich des Tarutung-Beckens. Die dort abgebildeten Anomalien (erhöhtes Vp/Vs, geringes Qp) habe ich als mögliche Aufstiegswege von warmen Fluiden interpretiert. Die scheinbar asymetrische Verteilung der Anomalien wird hierbei im Zusammenhang mit der Seismizitätsverteilung, der Geometrie der Beben-Bruchflächen, sowie struktur-geologischen Modellvorstellungen diskutiert. Damit werden wesentliche Informationen für die Planung einer zukünftigen geothermischen Anlage bereitgestellt.

Near surface measurements of seismic velocity and Q are useful in a number of situations, for example for when carrying out re-datuming and migration for depth images, or when analysing ground conditions for building. This thesis concentrates on the estimation of surface wave group and phase velocity as well as Q structure through the use of cross correlations of ambient noise recordings. Linearised tomography estimates are made for the British Isles, the Permain Basin of Texas and New Mexico, the Ekofisk Life of Field Seismic (LoFS) array and the Aquistore CO2 storage site. The results correspond well with the known geological structure and/or structure observed in velocity maps by other researchers. For the Ekofisk array a non-linearised non-linear method was also applied and the results estimated by these two methods for the Ekofisk LoFS array are compared. By non-linearised non-linear it is meant that the inversion method is fully non-linear and no linearisation has taken place in the method, this term will be used throughout this thesis for all methods which fall into this category. The tomography results from the two methods had similarities in their general structure but differences in the finer details, and so suggest that the substantial increase in time required for the non-linearised non-linear method is not justified. Linearised tomographic inversion of the Aquistore array was used to determine the potential of using ambient noise tomography for monitoring of CO2 storage sites. It was found that the repeatability of the tomography at the Aquistore site was not good enough to allow ambient noise tomography to be used for monitoring; however, it may be possible to apply the method at other sites. A Q and phase velocity inversion of the Ekofisk array is also presented, with results mostly showing excellent correlation with known geological features. It is shown that the higher frequencies are more sensitive to the effects of sea floor subsidence at the site, while lower frequencies are more sensitive to the effects of faulting. A final near surface method called ambient noise gradiometry was applied, this concentrates on estimating locations of sources of seismic energy within receiver arrays. Ambient noise gradiometry is applied to synthetic and real data for this purpose. It was found that using ambient noise gradiometry allows internal sources of energy to be identified but they produce a bias in the phase velocity tomography result. Two methods of reducing this bias are presented, both of which also provide an estimate of the source term for different sections of time of the recording.

京都大学 / 新制・課程博士 / 博士(工学) / 甲第23861号 / 工博第4948号 / 新制||工||1773(附属図書館) / 京都大学大学院工学研究科都市社会工学専攻 / (主査)教授 清野 純史, 教授 三村 衛, 准教授 古川 愛子 / 学位規則第4条第1項該当 / Doctor of Philosophy (Engineering) / Kyoto University / DGAM

The structure of fault zones (FZs) plays an important role in understanding fault mechanics, earthquake rupture and seismic hazards. Fault zone seismic guided waves (GW) carry important information about internal structure of the low-velocity fault damage zone. Numerical modeling of observed FZGWs has been used to construct models of FZ structure. However, the depth extent of the waveguide and the uniqueness of deep structure in the models have been debated. Elastic finite-difference synthetic seismograms were generated for FZ models that include an increase in seismic velocity with depth both inside and outside the FZ. Strong GWs were created from sources both in and out of the waveguide, in contrast with previous homogenous-FZ studies that required an in-fault source to create GW. This is because the frequency-dependent trapping efficiency of the waveguide changes with depth. The near-surface fault structure efficiently guides waves at lower frequencies than the deeper fault. Fault structure at seismogenic depth requires the analysis of data at higher frequencies than the GWs that dominate at the surface. Adapting a two-station technique from surface wave studies, dispersive differential group arrival times between two earthquakes can be used to solve for FZ structures between the earthquakes. This method was tested with synthetic data and shallow events recorded in the SAFOD borehole in the San Andreas Fault. A pair of deep earthquakes recorded in the SAFOD borehole indicate a ~150 m wide San Andreas Fault waveguide with >20% velocity contrast at 10-12 km depth. With additional earthquakes, the full FZ structure at seismogenic depth could be imaged. Subsurface FZ structure can also be derived from a surface source and receiver array analogous to a body-wave refraction survey. Synthetic seismograms for such source-receiver geometry were generated and verified that FZGWs are refracted by the increase in velocity with depth. Synthetic data from a surface array were successfully inverted to derive FZ structure in the subsurface. The new methods presented in this dissertation extend the potential of FZGWs to image deeper FZ structure than has been uniquely constrained in the past. / Ph. D.

Expanding the standard, single-well recording geometry used to monitor seismicity during hydro-fracture treatments could provide more accurate hypocenter locations and seismic velocities, improving general reservoir characterization. However, for the real, two-well data set obtained for this project, only S-wave picks were available, and testing resulted in anomalous hypocenter location behavior. This study uses a hypocenter location algorithm and both real and synthetic data sets to investigate how the accuracy of the velocity model, starting hypocenter location, recording geometry, and arrival-time picking error affect final hypocenter locations. Hypocenter locations improved using a velocity model that closely matched the observed sonic log rather than a smoothed version of this model. The starting hypocenter location did not affect the final location solution if both starting and final locations were between the wells. Two solutions were possible when the true solution was not directly between the wells. Adding realistic random picking errors to synthetic data closely modeled the dispersed hypocenter error pattern observed in the real data results. Adding data from a third well to synthetic tests dramatically reduced location error and removed horizontal geometric bias observed in the two-well case. Seismic event data recorded during hydro-fracture treatments could potentially be used for three-dimensional joint hypocenter-velocity tomography. This would require observation wells close enough to earthquakes to record P- and S-wave arrivals or wells at orientations sufficient to properly triangulate hypocenter locations. Simulating results with synthetic tests before drilling could optimize survey design to collect data more effectively and make analysis more useful. / Master of Science

Les séismes de type longue période (LP) sont aujourd’hui enregistrés sur la plupart des volcans dans le monde entier. Malgré cela, le mécanisme à leur source n’est encore que très peu compris. A l’heure actuelle les modèles proposés pour expliquer leur origine sont : 1) la résonance d’une fracture remplie de gaz ou de fluides excités par des instabilités dans l’écoulement des fluides ou par la rupture fragile du magma ; 2) la fracturation lente des sédiments peu consolidés à la surface des volcans, dans des conditions de transition entre le ductile et le fragile. L’outil le plus utilisé pour comprendre leur nature est aujourd’hui l’inversion du tenseur des moments. Au cours des dernières années, les inversions du tenseur des moments se concentraient principalement sur la compréhension du mécanisme physique à l’origine des séismes LP qui souvent supposaient des milieux géologiques très simples, voire homogènes. Des études récentes ont montré l’influence des sédiments peu consolidés à la surface des volcans sur la propagation des ondes à basse fréquence et en conséquence, sur l’inversion du tenseur des moments quand ils ne sont pas pris en compte dans le processus d’inversion. Le but de cette thèse est de mieux comprendre les processus physiques qui génèrent les séismes LP et de quantifier les incertitudes liées à leur interprétation. / Long-period (LP) seismic events are abundantly recorded during rest and unrest periods at many volcanoes worldwide. However, their source mechanism is still poorly understood. Models which have been proposed so far to explain their origin are: 1) the resonance of a fluid-filled cavity triggered by fluid instabilities or the brittle failure of magma; 2) slow-rupture earthquakes occurring in the low consolidated materials composing the shallow portion of the volcanic edifice. Nowadays the main tool used to get insights into their nature is moment tensor (MT) inversion. MT inversions carried out in the past years focused mainly on the understanding of the physical origin of LP events and often supposed a relative simple geological structure of the medium. Recent studies highlighted the strong influence of shallow unconsolidated materials on the retrieved MT solutions and the importance of considering geological inhomogeneity in the inversion process. The principal aim of this thesis is to gain a better understanding of the source processes that generate LP events and to quantify the uncertainties related to the MT inversion process.