京都大学 / 新制・課程博士 / 博士(工学) / 甲第23863号 / 工博第4950号 / 新制||工||1773(附属図書館) / 京都大学大学院工学研究科都市社会工学専攻 / (主査)教授 小池 克明, 教授 須崎 純一, 准教授 柏谷 公希 / 学位規則第4条第1項該当 / Doctor of Philosophy (Engineering) / Kyoto University / DGAM

The distributed point source method, or DPSM, developed in the last decade has been used for solving various engineering problems-such as elastic and electromagnetic wave propagation, electrostatic, and fluid flow problems. Based on a semi-analytical formulation, the DPSM solution is generally built by superimposing the point source solutions or Green's functions. However, the DPSM solution can be also obtained by superimposing elemental solutions of volume sources having some source density called the equivalent source density (ESD). In earlier works mostly point sources were used. In this paper the DPSM formulation is modified to introduce a new kind of ESD, replacing the classical single point source by a family of point sources that are referred to as quantum sources. The proposed formulation with these quantum sources do not change the dimension of the global matrix to be inverted to solve the problem when compared with the classical point source-based DPSM formulation. To assess the performance of this new formulation, the ultrasonic field generated by a circular planer transducer was compared with the classical DPSM formulation and analytical solution. The results show a significant improvement in the near field computation.

abstract: The dynamic Earth involves feedbacks between the solid crust and both natural and anthropogenic fluid flows. Fluid-rock interactions drive many Earth phenomena, including volcanic unrest, seismic activities, and hydrological responses. Mitigating the hazards associated with these activities requires fundamental understanding of the underlying physical processes. Therefore, geophysical monitoring in combination with modeling provides valuable tools, suitable for hazard mitigation and risk management efforts. Magmatic activities and induced seismicity linked to fluid injection are two natural and anthropogenic processes discussed in this dissertation. Successful forecasting of the timing, style, and intensity of a volcanic eruption is made possible by improved understanding of the volcano life cycle as well as building quantitative models incorporating the processes that govern rock melting, melt ascending, magma storage, eruption initiation, and interaction between magma and surrounding host rocks at different spatial extent and time scale. One key part of such models is the shallow magma chamber, which is generally directly linked to volcano’s eruptive behaviors. However, its actual shape, size, and temporal evolution are often not entirely known. To address this issue, I use space-based geodetic data with high spatiotemporal resolution to measure surface deformation at Kilauea volcano. The obtained maps of InSAR (Interferometric Synthetic Aperture Radar) deformation time series are exploited with two novel modeling schemes to investigate Kilauea’s shallow magmatic system. Both models can explain the same observation, leading to a new compartment model of magma chamber. Such models significantly advance the understanding of the physical processes associated with Kilauea’s summit plumbing system with potential applications for volcanoes around the world. The unprecedented increase in the number of earthquakes in the Central and Eastern United States since 2008 is attributed to massive deep subsurface injection of saltwater. The elevated chance of moderate-large damaging earthquakes stemming from increased seismicity rate causes broad societal concerns among industry, regulators, and the public. Thus, quantifying the time-dependent seismic hazard associated with the fluid injection is of great importance. To this end, I investigate the large-scale seismic, hydrogeologic, and injection data in northern Texas for period of 2007-2015 and in northern-central Oklahoma for period of 1995-2017. An effective induced earthquake forecasting model is developed, considering a complex relationship between injection operations and consequent seismicity. I find that the timing and magnitude of regional induced earthquakes are fully controlled by the process of fluid diffusion in a poroelastic medium and thus can be successfully forecasted. The obtained time-dependent seismic hazard model is spatiotemporally heterogeneous and decreasing injection rates does not immediately reduce the probability of an earthquake. The presented framework can be used for operational induced earthquake forecasting. Information about the associated fundamental processes, inducing conditions, and probabilistic seismic hazards has broad benefits to the society. / Dissertation/Thesis / Doctoral Dissertation Geological Sciences 2018

Wireless has proved to be one of the most important and fastest growing fields of communications especially during last few decades. To achieve reliable communication, we model a wireless system to analyze its performance and to find ways to improve the reliability of a particular system. Extensive research is being done to accurately model wireless systems, and to achieve better performance. Simulation techniques have been in use for many years to support the design and evaluation of electronic communication systems. Over the past few decades, Computer Aided Design (CAD) techniques (including both computerized analytical techniques and simulation) have matured, and are now usually applied at some point in the system design/development process. The aim of this thesis is to find efficient algorithms that can model third generation wireless channels in a discrete sense. For modeling these channels, mathematical tools known as hidden Markov models are used. These models have proved themselves to be very efficient in many areas of electrical engineering including speech recognition, pattern recognition, artificial intelligence, wavelets and queuing theory. Wideband Code Division Multiple Access (W-CDMA) wireless communication parameters including channels fading statistics, Bit Error Rate (BER) performance and interval distribution of errors are modeled using different Markov models, and their results are tested and validated. Four algorithms for modeling error sources are implemented, and their results are discussed. Both hidden Markov models and semi-hidden Markov models are used in this thesis, and their results are validated for the W-CDMA environment. The state duration distributions for these channels are also approximated using Phase-Type (PH) distribution. / Master of Science

The Finite Difference Time Domain (FDTD) Method has been a powerful tool in numerical simulation of electromagnetic (EM) problems for decades. In recent years, it has also been applied to biomedical research to investigate the interaction between EM waves and biological tissues. In Wireless Body Area Networks (WBANs) studies, to better understand the localization problem within the body, an accurate source/receiver model must be investigated. However, the traditional source models in FDTD involve effective volume and may cause error in near field arbitrary direction. This thesis reviews the basic mathematical and numerical foundation of the Finite Difference Time Domain method and the material properties needed when modeling a human body in FDTD. Then Coincident Phase Centers (CPCs) point sources models have been introduced which provide nearly the same accuracy at the distances as small as 3 unit cells from the phase center. Simultaneously, this model outperforms the usual sources in the near field when an arbitrary direction of the electric or magnetic dipole moment is required.

Surface displacement at volcanic edifices is related to subsurface processes associated with magma movements, fluid transfers within the volcano edifice and gravity-driven deformation processes. Understanding of associated ground displacements is of importance for assessment of volcanic hazards. For example, volcanic unrest is often preceded by surface uplift, caused by magma intrusion and followed by subsidence, after the withdrawal of magma. Continuous monitoring of the surface displacement at volcanoes therefore might allow the forecasting of upcoming eruptions to some extent. In geophysics, the measured surface displacements allow the parameters of possible deformation sources to be estimated through analytical or numerical modeling. This is one way to improve the understanding of subsurface processes acting at volcanoes. Although the monitoring of volcanoes has significantly improved in the last decades (in terms of technical advancements and number of monitored volcanoes), the forecasting of volcanic eruptions remains puzzling. In this work I contribute towards the understanding of the subsurface processes at volcanoes and thus to the improvement of volcano eruption forecasting. I have investigated the displacement field of Llaima volcano in Chile and of Tendürek volcano in East Turkey by using synthetic aperture radar interferometry (InSAR). Through modeling of the deformation sources with the extracted displacement data, it was possible to gain insights into potential subsurface processes occurring at these two volcanoes that had been barely studied before. The two volcanoes, although of very different origin, composition and geometry, both show a complexity of interacting deformation sources. At Llaima volcano, the InSAR technique was difficult to apply, due to the large decorrelation of the radar signal between the acquisition of images. I developed a model-based unwrapping scheme, which allows the production of reliable displacement maps at the volcano that I used for deformation source modeling. The modeling results show significant differences in pre- and post-eruptive magmatic deformation source parameters. Therefore, I conjecture that two magma chambers exist below Llaima volcano: a post-eruptive deep one and a shallow one possibly due to the pre-eruptive ascent of magma. Similar reservoir depths at Llaima have been confirmed by independent petrologic studies. These reservoirs are interpreted to be temporally coupled. At Tendürek volcano I have found long-term subsidence of the volcanic edifice, which can be described by a large, magmatic, sill-like source that is subject to cooling contraction. The displacement data in conjunction with high-resolution optical images, however, reveal arcuate fractures at the eastern and western flank of the volcano. These are most likely the surface expressions of concentric ring-faults around the volcanic edifice that show low magnitudes of slip over a long time. This might be an alternative mechanism for the development of large caldera structures, which are so far assumed to be generated during large catastrophic collapse events. To investigate the potential subsurface geometry and relation of the two proposed interacting sources at Tendürek, a sill-like magmatic source and ring-faults, I have performed a more sophisticated numerical modeling approach. The optimum source geometries show, that the size of the sill-like source was overestimated in the simple models and that it is difficult to determine the dip angle of the ring-faults with surface displacement data only. However, considering physical and geological criteria a combination of outward-dipping reverse faults in the west and inward-dipping normal faults in the east seem to be the most likely. Consequently, the underground structure at the Tendürek volcano consists of a small, sill-like, contracting, magmatic source below the western summit crater that causes a trapdoor-like faulting along the ring-faults around the volcanic edifice. Therefore, the magmatic source and the ring-faults are also interpreted to be temporally coupled. In addition, a method for data reduction has been improved. The modeling of subsurface deformation sources requires only a relatively small number of well distributed InSAR observations at the earth’s surface. Satellite radar images, however, consist of several millions of these observations. Therefore, the large amount of data needs to be reduced by several orders of magnitude for source modeling, to save computation time and increase model flexibility. I have introduced a model-based subsampling approach in particular for heterogeneously-distributed observations. It allows a fast calculation of the data error variance-covariance matrix, also supports the modeling of time dependent displacement data and is, therefore, an alternative to existing methods. / Oberflächenverschiebungen an Vulkanen können einerseits durch unterirdische Magmen- oder Fluidbewegungen oder andererseits durch Gravitation verursacht werden. So sind insbesondere vor Eruptionen oft Aufwölbungen an Vulkanen zu beobachten, verursacht durch Magmenintrusion in die Erdkruste. Nach Eruptionen hingegen sinkt das Vulkangebäude aufgrund von Magmenextrusion wieder. Kontinuierliche Messungen an Vulkanen ermöglichen es, Eruptionen teilweise bis auf wenige Tage vorherzusagen. Die gemessenen Oberflächenverschiebungen können in analytischen oder numerischen Modellierungen genutzt werden, um Parameter eines möglichen Quellprozesses abzuschätzen. Auf diese Art und Weise kann das Verständnis über die unterirdischen Prozesse, die an Vulkanen stattfinden, verbessert werden. Obwohl es in den letzten Jahrzehnten eine enorme Entwicklung und Verbesserung der Überwachung von Vulkanen gab, sind viele Vorhersagen sehr vage und ungenau. Mit dieser Arbeit möchte ich einen Beitrag zum Verständnis von unterirdischen Prozessen an Vulkanen und auf lange Sicht gesehen, zur Vorhersage von Eruptionen leisten. Ich habe die Vulkane, Llaima in Chile und Tendürek im Osten der Türkei, mit Hilfe der Interferometrie von Radardaten (InSAR) untersucht. Die somit gemessenen Verschiebungen an der Erdoberfläche ermöglichen es, durch Modellierung der möglichen Deformationsquellen, Informationen über die Untergrundstrukturen dieser beiden bisher kaum erforschten Vulkane zu bekommen. Obwohl unterschiedlich in Aufbau, Gesteinszusammensetzung und Entstehung, zeigen beide Vulkane Anzeichen dafür, dass jeweils mehrere interagierende Deformationsquellen im Untergrund existieren. Am Vulkan Llaima war es schwierig, aufgrund der starken Dekorrelation des Radarsignals zwischen den Satellitenaufnahmen, die InSAR Methode anzuwenden. Ich entwickelte eine Methode um die doppeldeutigen relativen Phasenwerte der Interferogramme modellbasiert in eindeutige relative Phasenwerte umzurechnen. Die damit erzeugten Oberflächenverschiebungskarten am Vulkan eigneten sich nun für eine anschließende Modellierung der Deformationsquelle. Die Modellierungsergebnisse zeigen signifikante Unterschiede zwischen den Parametern der präeruptiven- und posteruptiven Deformationsquellen. Demzufolge könnten zwei unterschiedliche, interagierende Magmenkammern unter Llaima existieren, eine tiefe, posteruptiv aktive Kammer und eine flache, durch den Aufstieg von Magma präeruptiv aktive Kammer. Am Vulkan Tendürek ist eine langfristige, kontinuierliche Senkung des Vulkangebäudes zu beobachten, die mit einem großen, aufgrund von Kühlung sich kontrahierenden, magmatischen Sill, erklärbar ist. Unter Hinzunahme von hochauflösenden, optischen Daten jedoch, sind bei genauerer Untersuchung bogenförmige Strukturen an der Erdoberfläche sichtbar. Diese sind Anzeichen dafür, dass Verwerfungen existieren, die das gesamte Vulkangebäude in einem elliptischen Ring umgeben. Dabei ist zu beobachten, dass die Ringstörungen über Jahrtausende, möglicherweise sogar kontinuierlich, geringe Magnituden von Versatz aufweisen. Bei langer, kontinuierlicher Aktivität über mehrere zehntausende von Jahren, könnte dies ein weiterer Mechanismus zur Entstehung von Calderastrukturen an Vulkanen darstellen, der jedoch sehr langsam verläuft. Im Gegensatz dazu ist die heutige weit verbreitete Auffassung, dass Calderen als Folge katastrophaler Einstürze von Vulkangebäuden entstehen. Um zu untersuchen welche Geometrie die vorgeschlagenen Strukturen Sill und Ringstörungen an Tendürek im Untergund haben könnten, vollführte ich eine weitaus komplexere numerische Modellierung. Diese zeigt, dass die Größe des Sills ohne Berücksichtigung der Ringstörung um ein Vielfaches überschätzt ist. Die Orientierung und Geometrie der Ringstörungen ist jedoch nicht eindeutig nur mit Oberflächenverschiebungsdaten auflösbar. Unter der Berücksichtigung von geologischen und physikalischen Gesichtspunkten sind nach Außen einfallende Aufschiebungen im Westen und nach Innen einfallende Abschiebungen im Osten die plausibelste Erklärung. Außerdem habe ich eine Methode zur Datenreduzierung entwickelt. Abhängig vom zu untersuchenden Prozess sind für die Modellierung von unterirdischen Deformationsquellen verhältnismäßig wenige gut verteilte Messpunkte an der Erdoberfläche ausreichend. Satelliten gestützte Radaraufnahmen haben jedoch oft mehrere Millionen dieser Punkte. Deshalb müssen diese riesigen Datensätze auf eine Art und Weise reduziert werden, dass keine oder nur möglichst wenige Informationen verloren gehen. Für diesen Zweck habe ich, ausgehend von einem existierenden Algorithmus, eine modellbasierte Methode zur Reduzierung von besonders heterogen verteilten Oberflächendaten entwickelt. Diese Methode ist besonders gut auf Zeitreihendatensätze anwendbar und stellt somit eine Alternative zu existierenden Algorithmen dar.

The aim of this dissertation is to develop a new algorithm to embed a watermark in JPEG compressed images, using encoding methods. This encompasses the embedding of proprietary information, such as identity and authentication bitstrings, into the compressed material. This watermark encoding scheme involves combining entropy coding with homophonic coding, in order to embed a watermark in a JPEG image. Arithmetic coding was used as the entropy encoder for this scheme. It is often desired to obtain a robust digital watermarking method that does not distort the digital image, even if this implies that the image is slightly expanded in size before final compression. In this dissertation an algorithm that combines homophonic and arithmetic coding for JPEG images was developed and implemented in software. A detailed analysis of this algorithm is given and the compression (in number of bits) obtained when using the newly developed algorithm (homophonic and arithmetic coding). This research shows that homophonic coding can be used to embed a watermark in a JPEG image by using the watermark information for the selection of the homophones. The proposed algorithm can thus be viewed as a ‘key-less’ encryption technique, where an external bitstring is used as a ‘key’ and is embedded intrinsically into the message stream. The algorithm has achieved to create JPEG images with minimal distortion, with Peak Signal to Noise Ratios (PSNR) of above 35dB. The resulting increase in the entropy of the file is within the expected 2 bits per symbol. This research endeavor consequently provides a unique watermarking technique for images compressed using the JPEG standard. / Dissertation (MEng)--University of Pretoria, 2008. / Electrical, Electronic and Computer Engineering / unrestricted

Les travaux de cette thèse concernent la modélisation et la simulation numérique de la propagation d’ondes acoustiques en présence d’un écoulement. Le modèle retenu pour ces études est l’équation de Galbrun. Les travaux faits sur l’équation de Galbrun ont essentiellement porté sur le régime harmonique. En revanche, la plupart des études mathématiques et numériques du problème de l’aéroacoustique est en régime transitoire. C’est pourquoi, il est intéressant pour nous d’étudier l’équation de Galbrun en régime transitoire. Pour résoudre cette équation en régime transitoire, notre approche a reposé sur la transformée de Laplace, qui nous permet de faire l’échange entre le domaine harmonique et le domaine réel. Un autre sujet abordé dans cette thèse est celui du traitement des conditions aux limites non réfléchissantes en écoulement uniforme et non-uniforme. Nous proposons la méthode PML pour l’équation de Galbrun. Inspirée par la méthode de Hu, nous proposons un nouveau modèle PML associé à l’équation de Galbrun, qui a toujours conduit à une solution exponentiellement décroissante dans la couche, même en présence d’ondes inverses. Les simulations acoustiques montrent étonnamment d’erreur de convergence pour les deux modèles classiques et nouveaux. Nous validons notre modèle PML à travers plusieurs exemples numériques dans l’écoulement uniforme et non-uniforme. Le dernier objectif est de proposer des modèles de sources aéroacoustiques associées à l’équation de Galbrun. Après une présentation en détail des modèles existants, on adapte une méthode hybride (EIF) à l’équation de Galbrun. Pour assurer la validité de l’approche globale, certains tests classiques sont choisis parmi la littérature et les résultats sont comparés avec les approches existantes et les solutions analytiques. / The work of this thesis is about the numerical modeling and simulation of the propagation of acoustic waves in the presence of a flow. The model used for these studies is the equation of Galbrun. The work done on the Galbrun equation focused on the harmonic regime. In contrast, most of the mathematical and numerical studies of the aeroacoustics problems are in the transient regime. That is why it is interesting for us to study the Galbrun equation in the transient regime. To solve this equation in the transient regime, our approach is based on the Laplace transform, which allows us to exchange between the frequency domain and the real domain. Another topic discussed in this thesis is the treatment of non-reflecting boundary conditions in uniform and non-uniform flow. We propose the Perfectly Matched Layer method for the Galbrun equation. Inspired by the Hu’s method, we propose a new PML model associated with the Galbrun equation, which always leads to an exponentially decreasing solution in the layer, even in the presence of reverse waves. Acoustic simulations show surprisingly error convergence for both classic and new models. We validate our PML model through several numerical examples in uniform and non-uniform flow. The final objective is to propose models for aeroacoustics sources associated with the Galbrun equation. After presenting in detail the existing models, we adapt a hybrid method (Expansion about Incompressible Flow) in Galbrun equation. To ensure the validity of the overall approach, some classical tests are selected from the literature and the results are compared with existing approaches and analytical solutions.

Bolu-Ilgaz region was damaged by several large earthquakes in the last century and the structural damage was substantial especially after the 1944 and 1999 earthquakes. The objective of this study is to build the seismic source characterization model for the rupture zone of 1944 Bolu-Gerede earthquake and perform probabilistic seismic hazard assessment (PSHA) in the region. One of the major improvements over the previous PSHA practices accomplished in this study is the development of advanced seismic source models in terms of source geometry and reoccurrence relations. Geometry of the linear fault segments are determined and incorporated with the help of available fault maps. Composite magnitude distribution model is used to properly represent the characteristic behavior of NAF without an additional background zone. Fault segments, rupture sources, rupture scenarios and fault rupture models are determined using the WG-2003 terminology. The Turkey-Adjusted NGAW1 (G&uuml / lerce et al., 2013) prediction models are employed for the first time on NAF system. The results of the study is presented in terms of hazard curves, deaggregation of the hazard and uniform hazard spectrum for four main locations in the region to provide basis for evaluation of the seismic design of special structures in the area. Hazard maps of the region for rock site conditions and for the proposed site characterization model are provided to allow the user perform site-specific hazard assessment for local site conditions and develop site-specific design spectrum. The results of the study will be useful to manage the future seismic hazard in the region.

The North Anatolian Fault System is one of the most important active strike-slip fault systems in the world. The August 17, 1999 and November 12, 1999 earthquakes at Kocaeli and D&uuml / zce are the most recent devastating earthquakes. The study area lies in the Eastern Marmara Region and is bounded by the 28.55-33.75 E and 40.00-41.20 N, latitude and longitude coordinates, respectively. There are numerous studies conducted in the study area in terms of active tectonics and seismicity, however studies are scale dependent. Therefore, a comprehensive literature survey regarding active tectonics of the region was conducted and these previous studies were combined with the lineaments extracted from 10 ASTER images via principle component analysis manual extraction method. Therefore, a line seismic source model for the Eastern Marmara region was compiled mainly based on major seismic events of instrumental period. The seismicity of these line segments were compared with the instrumental period earthquake catalogue compiled by Kandilli Observatory and Earthquake Research Institute with a homogeneous magnitude scale between 1900 and 2005. Secondary event and completeness of this catalogue was checked. The final catalogue was matched with the compiled seismic source for historical seismicity and source-scenario-segment-weight relationships were developed. This developed seismic source model was tested by a probabilistic seismic hazard assessment for D&uuml / zce city center by utilizing four different ground motion prediction equations. It was observed that Gutenberg-Richter seismicity parameter &lsquo / b&rsquo / does not have significant effect over the model, however change in the segmentation model have a low but certain influence.