Numerical modeling of seismic field and soil interaction : application to the sedimentary basin of Argostoli (Greece)
Institution:
université Paris-SaclayDisciplines:
Directors:
Abstract EN:
A prediction of the seismic response of civil engineering structures that requires a high level of safety (i.e. nuclear power plants or dams) faces several major difficulties given the complexity of the problem being treated. To this end, the source simulation, seismic wave propagation and site effects have been studied by different approaches over the last two decades. Recently, numerical methods, such as the spectral element, combined with massively parallel computers have proved a good efficiency in modelling the seismic wave propagation from source to site in complex three-dimensional geological environments. However, the accuracy of these predictions remains limited due to the multi-scale nature of the problem and the large uncertainties in the data to be introduced into the model (i.e. the geometric and kinematic characterization of the seismic source, the detailed geological and numerical model of the source-to-site propagation path). Therefore, the use of a regional numerical model able to simulate the seismic phenomenon from source to site would allow a better analysis and classification of the origin of the associated uncertainties.This work aims to study numerically the effect of local and regional geology on the seismic response of a basin and precisely the Argostoli site located on the island of Kefalonia (Greece). Firstly, the numerical code used in this work (SEM3D) is verified using three canonical cases. The simulated results showed a good agreement with the reference solutions. In the cases of a numerical model with different scales or with complex geologies, one of the important difficulties is the conformity of the numerical meshes with the geological interfaces, that will result an increase in the numerical cost. Because of the characteristics of the used numerical method, one possible approach was to use the “non-conforming” or “not-honoring” meshes to overcome this difficulty. A parametric study on the applicability of this approach was then carried out in order to highlight the influence of some numerical parameters on the obtained results.Thereafter, parametric studies on several seismic scenarios in the Argostoli site were conducted. Concerning seismic loading, two types of source have been studied: point sources and extended faults. The study with point sources clearly revealed a phenomenon of amplification and trapping of waves in the basin, leading to a complex and elongated signal, with significant energy compared to a study with a simplified geology. For the second type of source, the modeled fault is close to the surface. It allows to study, in addition to the influence of the basin, the effect of the near-field on the seismic response of the site. Indeed, near-fault ground motion can be significantly different from ground motion observed far from the seismic source. Based on the results obtained, the basin effect is even more pronounced but with higher amplification and different resonance frequencies. In addition, the near-field effect has been highlighted, marked by strong velocity pulses at some locations in the basin. The order of magnitude of the obtained response spectra is comparable to the ones obtained during the earthquake sequence of similar magnitude that took place in 2014.In the last part, a parametric study (allowed by the development of computing power) on the numerical aspects related to the computational accuracy was carried out. With this study, it is possible to increase the frequency resolution from 7 Hz to 10 Hz with soft soil mechanical characteristics while keeping the same domain size. These simulations open even more questions on the interdependence of the fineness of resolution of physical data and meshes for numerical simulations. In conclusion, this thesis corresponds to a first step in the numerical characterization of the seismic response of the Argostoli basin and the effects due to the basin, the type of source and their interactions.
Abstract FR:
Une prédiction de la réponse sismique des structures de génie civil (telles que les centrales nucléaires ou les barrages), doit faire face à plusieurs difficultés majeures compte tenu de la complexité du problème traité. Pour cela, la simulation de la source, la propagation des ondes sismique et les effets de site, ont été étudiés par différentes approches dernièrement. Récemment, des méthodes combinées avec des ordinateurs massivement parallèles se sont avérées très efficaces pour modéliser la propagation des ondes sismiques depuis la source jusqu’au site, dans des environnements géologiques tridimensionnels complexes. Cependant, la précision reste limitée en raison du caractère multi-échelle du problème et des grandes incertitudes sur les données à introduire dans le modèle (i.e. la caractérisation géométrique et cinématique de la source sismique, le modèle géologique et numérique du chemin de propagation). Pour cela, l’utilisation d'un modèle numérique régional capable de simuler le phénomène sismique de la source au site permettrait de mieux comprendre et de classer l'origine de ces incertitudes.Ce travail vise à étudier numériquement l'effet de la géologie locale et régionale sur la réponse sismique d’un bassin et plus particulièrement le site d’Argostoli situé sur l'île de Céphalonie (Grèce). Dans un premier temps, le code numérique utilisé dans ce travail (SEM3D) est vérifié à l'aide de trois cas canoniques. Les résultats obtenus montrent un bon accord avec les solutions de référence. Dans les cas d’un modèle numérique avec diverses échelles ou avec des géologies complexes, l’un des points clefs est la conformité du maillage numérique avec les interfaces géologiques, ce qui se traduit par l’augmentation du coût numérique. Grâce aux caractéristiques de la méthode numérique utilisée, pour contourner cette difficulté, une approche possible est l’utilisation des maillages “non-conformes” ou “not-honouring”. Une étude paramétrique sur la faisabilité de cette approche est ainsi réalisée afin de mettre en évidence l'influence de certains paramètres numériques sur les résultats obtenus.Par la suite, des études paramétriques de plusieurs scénarii sismiques sur le site d'Argostoli ont été réalisées. Concernant le chargement sismique, deux types de sources ont été étudiés : des sources ponctuelles et des failles étendues. L’étude avec les sources ponctuelles révèle un phénomène d'amplification et de piégeage des ondes dans le bassin, conduisant à un signal complexe et allongé, avec une énergie importante par rapport à une étude avec une géologie simplifiée. Pour le second type de source, la faille modélisée est proche de la surface. Cela permet d'étudier l’effets du bassin et du champ-proche sur la réponse sismique du site. En effet, le mouvement du sol à proximité d'une faille peut être différent du mouvement du sol observé loin de la source sismique. D’après les résultats obtenus, l'effet du bassin est plus marqué mais avec une amplification plus forte et des fréquences de résonance différentes. De plus, l'effet du champ-proche a été mis en évidence, marqué par de fortes impulsions de vitesse à certains endroits du bassin. L’ordre de grandeur des spectres de réponses obtenues est comparable à ceux obtenus lors des séismes de magnitude semblable qui ont eu lieu en 2014 au même endroit.Dans la dernière partie, une étude paramétrique sur des aspects numériques liées à la précision du calcul a été réalisée. Cette étude a permis d’augmenter la résolution fréquentielle de 7 Hz à 10 Hz avec des caractéristiques mécaniques de sols mous avec la même taille de domaine. Ces simulations ouvrent plus de questions sur l’interdépendance de la finesse de la résolution des données physiques et des maillages pour les calculs numériques. En conclusion, cette thèse correspond à une première étape dans la caractérisation numérique de la réponse sismique du bassin d’Argostoli et les effets dus au bassin, au type de source et à leurs interactions.