Le Liban n'a pas souffert de grands tremblements de terre destructeurs depuis près de deux siècles. Il est toutefois traversé par la faille transformante majeure du Levant, séparant sur1000 km de longueur la plaque Arabique à l’est de la plaque Africaine à l’ouest. Ses principales branches au Liban sont la faille de Yammouneh qui traverse le pays du sud au nord, les failles de Serghaya et Rachaya dans sa partie Est, la faille de Roum et les failles inverses du Mont Liban dans la partie Ouest. Ces failles ont généré des séismes destructeurs dans la longue histoire connue de la région, parmi lesquels les plus importants sont ceux deJuillet 551 sur la faille du Mont-Liban, de mai 1202 sur la faille de Yammouneh, d’Octobre1759 sur la faille de Rachaya et de Novembre 1759 sur la faille de Serghaya. L'évaluation del'aléa et du risque sismique local est donc de première importance pour l'ensemble du pays.L'objectif du travail effectué dans cette thèse est d'appliquer au Liban les avancées réalisées ces dernières années dans le développement de nouveaux outils à la fois fiables et économiquement abordables pour l’évaluation de l’aléa sismique, en commençant par les grandes villes et en particulier la capitale Beyrouth. L'objectif est de mieux appréhender et comprendre le risque sismique sur le territoire libanais, pour pouvoir ensuite commencer à élaborer des politiques de prévention et des codes parasismiques qui puissent le réduire à terme.Un réseau sismologique temporaire composé de 10 stations a été installé dans Beyrouth etune partie de sa banlieue sur des sites représentatifs des principales unités géologiques présentes. Plusieurs dizaines de séismes locaux et régionaux ont pu y être enregistrés, et leur réponse sismique a été évaluée par la méthode du rapport spectral site sur référence (SSR),comparé au rapport spectral de la composante horizontale sur la composante verticale (H/V)calculé sur les tremblements de terre et sur le bruit ambiant. Les mêmes enregistrements ont également été utilisés pour prédire, par la technique des Fonctions de Green empiriques(FGE), le mouvement sismique correspondant à un événement majeur (Mw7.5) sur la faille de Yammouneh. Cet exercice de prédiction a toutefois été réalisé en deux étapes en raison des limitations dans l'application de la technique FGE en champ proche, avec deux techniques complémentaires: l'enregistrement d’un petit événement a été d'abord utilisé pour simuler un séisme de Mw6.5 sur la faille de Yammouneh, en parallèle à l'utilisation de plusieurs équations de prédiction de mouvement du sol (GMPE), soigneusement sélectionnées, pour effectuer une prédiction similaire. La comparaison FGE/GMPE a alors permis de calibrer la prédiction du mouvement du sol par GMPE à différents sites de Beyrouth pour l'événement cible de Mw7.5. Ces résultats ponctuels ont ensuite été étendus à l'ensemble de la municipalité de Beyrouth et de sa banlieue proche, en vue de mieux cerner les contours d'une future carte de microzonage, au travers d'une vaste campagne de mesures de bruit ambiant sur615 sites. Leur traitement H/V a permis d’obtenir une carte de la fréquence de résonance pour l'ensemble de la zone, carte dont la robustesse a été testée et prouvée. Des mesures sismiques actives et passives ont en outre été menées sur les principales unités géologiques à proximité des 10 sites préalablement sélectionnés et instrumentés, permettant ainsi d'obtenir les premières estimations directes de la vitesse des ondes de cisaillement (via les courbes de dispersion des ondes de Rayleigh). La comparaison de ces mesures avec les estimations -très dispersées - issues de la compilation des paramètres géologiques/géotechniques disponibles et des équations de corrélation existantes avec les valeurs N des SPT, montre tout l'intérêt de ces mesures simples et fiables. / Lebanon is one of the countries that have not suffered from large destructive earthquakes foralmost two centuries. It is however lying on the 1000 km long, left lateral Levant fault thatseparates the Arabic plate in the east from the African plate in the west. Its main branches inLebanon are the Yammouneh fault that crosses the country from south to north, the Serghayaand Rachaya faults in its Eastern part, the Roum and Mount Lebanon Thrust faults in itsWestern part. These faults have generated destructive earthquakes in the long known historyof the area. The largest events are: The July 551 earthquake on the Mount Lebanon Thrustfault, the May 1202 earthquake on the Yammouneh fault, the October 1759 on the Rachayafault and the November 1759 on the Serghaya fault. From all above, one can conclude thatLebanon is exposed to a significant seismic hazard. Assessing the local seismic hazard andrisk is therefore of primary importance for the whole country.The objective of the work undergone in this PhD is to take advantage of the latest advancesachieved worldwide to promote rather inexpensive, though reliable, seismic hazardassessment tools, to try to apply them in Lebanon starting with the big cities and specificallythe capital Beirut. These studies will help to understand the Lebanese seismic risk andsubsequently to start to elaborate seismic policies and codes that may help reducing this risk.A temporary seismological network consisting of 10 stations has been installed in Beirut anda part of its suburbs. Several tens of local and regional earthquakes could be recorded, andallowed to estimate the site response at selected sites in Beirut through the standard site toreference spectral ratio method ("SSR") on earthquakes, compared to the horizontal tovertical ratio ("H/V") calculated on earthquakes and on ambient noise. The same recordingscould also be used via the empirical Green Function’s technique ("EGF") to predict theseismic ground motion corresponding to a Mw7.5 on the Yammouneh fault. However, due tolimitations in near-field applications of the EGF technique, this prediction exercise wasperformed in two steps and with two complementary techniques: a weak event recording wasfirst used to simulate a Mw6.5 earthquake on the Yammouneh fault, while several, carefullyselected ground motion prediction equations (GMPE) were used to perform a comparativeprediction for the same earthquake. This EGF/GMPE comparison then allowed tuning theGMPE prediction of ground motion at various sites within Beirut for the target Mw7.5 event.The results were then extended in view of proposing a framework for a future microzonationviimap. A comprehensive campaign of ambient noise measurements was achieved for 615 sitesof Beirut municipality and close suburbs, the H/V processing of which allowed to derive arobust map of resonance frequency for the whole area. In addition, active and passive seismicmeasurements were conducted on different geological units near the 10 formerlyinstrumented sites, which provided quantitative estimates of the shallow S-wave velocitythrough the Rayleigh wave dispersion curves. These geophysical measurements permitted toprovide direct estimates of the shear waves velocity, which prove much more reliable thanthe highly scattered estimates derived from the compilation of the availablegeological/geotechnical parameters and the use of existing correlations equations betweenSPT N-value and S-wave velocity Vs.

Lors de séismes forts, la propagation des ondes sismiques dans les sols met en jeu des non linéarités de comportement qui se manifestent différemment selon le niveau de sollicitation. En effet, pour de faibles déformations (généralement <10^{-6}), une loi de comportement linéaire (i.e. module et amortissement indépendants du niveau de sollicitation) permet de reproduire les observations expérimentales sur site. En revanche, pour des déformations plus élevées, une loi de comportement non linéaire hystérétique est nécessaire pour décrire l'évolution de la rigidité et des dissipations énergétiques au cours de la sollicitation sismique. De plus, comme les séismes forts sont caractérisés par des amplitudes et des durées plus importantes, le rôle de la pression interstitielle ne peut pas être négligé pour les sols saturés sous fortes sollicitations (mobilité cyclique et liquéfaction). Ces phénomènes peuvent conduire à l'annulation des contraintes effectives et devenir cause de dommages sévères pour les structures et les ouvrages. L'analyse proposée élargit l'applicabilité des modèles de calcul actuels pour une analyse plus fine du risque sismique. En partant d'une formulation aux éléments finis décrivant la propagation des ondes sismiques suivant la direction verticale en prenant en compte le chargement en 3D (l'approche « 1D-3C» : une direction-trois composantes) dans les sols nonlinéaires secs, de nouvelles stratégies pour la prise en compte du rôle de l'eau sont développées. Le modèle de comportement est basé sur la relation entre la pression interstitielle et le travail de la contrainte de cisaillement. Ce modèle décrivant l'évolution de la pression interstitielle considère l'état de contrainte tridimensionnel du matériau. Le modèle est validé par comparaison avec des résultats expérimentaux. L'approche « 1D-3C » a été utilisée pour modéliser la réponse des sols pour 4 séismes réels : le séisme de Superstition Hills en 1987 aux États-Unis (M_{w}=6.7); le séisme de Tohoku en 2011 au Japon (M_{w}=9.1 ); le séisme de Kushiro en 1993 au Japon (M_{w}=7.8) et le séisme d'Emilie Romagne en 2012 en Italie (M_{w}=5.9). Pour les trois premiers séismes, des enregistrements en profondeur et en surface sont disponibles. L'étude de ces trois premiers cas rend possible la validation du modèle par comparaison des données enregistrées et calculées. Le modèle peut donc être considéré comme un outil fiable pour la prédiction de la réponse sismique des sols saturés / During strong earthquakes, the seismic wave propagation in soils involves nonlinear behaviors strongly depending on the strain level. Indeed, for small strain (typically <10^{-6}), a linear constitutive law (modulus and damping independent on the load level) can reproduce the experimental observations on site. However, for larger strains, a nonlinear hysteretic constitutive law is needed to describe the evolution of stiffness and energy dissipation during seismic loading. In addition, as strong earthquakes are characterized by larger amplitudes and durations, the role of pore pressure cannot be neglected for saturated soils. Indeed pore water pressure controls phenomena such as cyclic mobility and liquefaction due to the loss of soil strength. This can lead to a fast decrease of effective stresses and permanent deformations in the soil causing severe damage to structures. This work extends the applicability of existing calculation models for a more detailed analysis of seismic risk. Starting from a FEM approach describing the propagation of seismic waves in the vertical direction, taking into account 3D loading (so-called "1D-3C" approach: 1 direction - 3 components) in nonlinear dry soils, new strategies to consider the role of water are developed. The model is based on the relationship between the pore pressure and the shear work. The three-dimensional stress state of the material is considered. The model is validated by comparison with experimental results. The "1D-3C" approach was used to model the response of soils for four real earthquakes: the Superstition Hills earthquake in 1987 in the United States (M_{w}=6.7), the Tohoku earthquake in 2011 in Japan (M_{w}=9.1), the Kushiro earthquake in Japan in 1993 (M_{w}=7.8) and the Emilia Romagna earthquake in Italy in 2012 (M_{w}=5.9). For the first three earthquakes, records at depth and on the surface are available. The study of the first three cases makes possible the validation of the model by comparing the calculated accelerations on the surface with the available records. The model can then be considered as an advanced tool for the prediction of the seismic soil response

L'évaluation de l’aléa sismique doit tenir compte des différents aspects qui interviennent dans le processus sismique et qui affectent le mouvement du sol en surface. Ces aspects peuvent être classés en trois grandes catégories : 1) les effets de source liés au processus de rupture et à la libération d'énergie sur la faille. 2) les effets liés à la propagation de l'énergie sismique à l'intérieur de la Terre. 3) l'influence des caractéristiques géotechniques des couches peu profondes ; appelé effet de site. Les effets de site sont pris en compte dans la mitigation des risques par l'évaluation de la réponse sismique du sol. Lors de sollicitations cycliques, le sol présente un comportement non-linéaire, ce qui signifie que la réponse dépendra non seulement des paramètres du sol mais aussi des caractéristiques du mouvement sismique (amplitude, contenu en fréquence, durée, etc.). Pour estimer la réponse non-linéaire du site, la pratique habituelle consiste à utiliser des simulations numériques avec une analyse linéaire équivalente ou une approche non-linéaire complète. Dans ce document, nous étudions l'influence du comportement non-linéaire du sol sur la réponse du site sismique en analysant les enregistrements sismiques des configurations des réseaux de forages. Nous utilisons les données du réseau Kiban Kyoshin (KiK-Net). Les 688 sites sont tous équipés de deux accéléromètres à trois composantes, l'un situé à la surface et l'autre en profondeur. À partir de ces données, nous calculons les amplifications du mouvement du sol depuis la surface jusqu'aux enregistrements en fond de puit à l'aide des rapports spectraux de Fourier. Une comparaison entre le rapport spectral pour le faible et le fort mouvement du sol est alors réalisée. Le principal effet du comportement non-linéaire du sol sur la fonction de transfert du site est un déplacement de l'amplification vers les basses fréquences. Nous proposons une nouvelle méthodologie et un nouveau paramètre appelé fsp pour quantifier ces changements et étudier les effets non-linéaires. Ces travaux permettent d'établir une relation site-dépendante entre le paramètre fsp et le paramètre d'intensité du mouvement du sol. La méthode est testée sur les données accélérométriques du séisme de Kumamoto (Mw 7.1, 2016). Nous proposons ensuite d’utiliser des corrélations entre moment seismic et la duration de la faille (Courboulex et al., 2016), obtenues à partir d’une base de données globale de fonctions source et une méthode basée sur l’approche des fonctions de Green empiriques (EGF) stochastiques pour simuler les mouvements forts du sol dus à un futur séisme. Cette méthodologie est appliquée à la simulation d’un séisme de subduction en Équateur et comparée aux données réelles du séisme de Pedernales (Mw 7.8, 16 avril 2016) dans la ville de Quito. Nous proposons enfin de combiner la méthode de simulation de mouvements forts par EGF et la prise en compte des effets non-linéaires proposée dans les premiers chapitres. La méthode est testée sur les données accélérométriques d’une réplique du séisme de Tohoku (Mw 7.9). / Seismic hazard assessments must consider different aspects that are involved in an earthquake process and affect the surface ground motion. Those aspects can be classified into three main kinds. 1) the source effects are related to the rupture process and the release of energy. 2) the path effects related to the propagation of energy inside Earth. 3) the influence of the shallow layers geotechnical characteristics; the so-called site-effects. The site effects are considered in risk mitigation through the evaluation of the seismic soil response. Under cyclic solicitations the soil shows a non-linear behavior, meaning that the response will not only depend on soil parameters but also on seismic motion input characteristics (amplitude, frequency content, duration, …). To estimate the non-linear site response, the usual practice is to use numerical simulations with equivalent linear analysis or truly non-linear time domain approach. In this document, we study the influence of the nonlinear soil behavior on the seismic site response by analyzing the earthquake recordings from borehole array configurations. We use the Kiban Kyoshin network (KiK-Net) data. All 688 sites are instrumented with two 3-components accelerometers, one located at the surface and the another at depth. From these data, we compute the ground motion amplifications from the surface to downhole recordings by the computing Fourier spectral ratios for the aim to compare between the spectral ratio for weak and strong ground motion. The main effect of the non-linear behavior of the soil on the site transfer function is a shift of the amplification towards lower frequencies. We propose a new methodology to quantify those changes and study the nonlinear effects. This work results in a site-dependent relationship between the changes in the site response and the intensity parameter of the ground motion. The method is tested analyzing the records of the earthquake of Kumamoto (Mw 7.1, 2016). Posteriorly, we propose to integrate a correlation between seismic moment and the duration of the fault (Courboulex et al., 2016) in the empirical Green’s function method. This methodology was applied to simulate one seduction event in Equator, and we compare the results with the records of the Pedernales earthquake (Mw 7.8, 2016) in the city of Quito. We attempt to take in account the nonlinear effects in the empirical Green’s function method. We use the methodologies of the first part of this document based on the frequency shift parameter. The procedure could be implemented in other methodologies that can predict an earthquake at a rock reference site, such as the stochastic methods. We test the procedure using the accelerometric records for one of the aftershocks o the Tôhoku earthquake (Mw 7.9).

Ce travail de recherche se situe dans le cadre du développement de la méthode des éléments de frontière (BEM) pour les milieux poreux multiphasiques. À l'heure actuelle, l'application de la BEM aux pr oblèmes des milieux poreux non-saturés est encore limitée, car l'expression analytique exacte de la solution fondamentale n'a pas été obtenue, ni dans le domaine transformé ni dans le domaine réel. Ceci provient de la complexité du système d'équations régissant le comportement des milieux poreux non-saturés. Les développements de la BEM pour les sols non-saturés effectués au cours de cette thèse sont basés sur les modèles thermo-hydro-mécanique (THHM) et hydro-mécanique (HHM) présentés dans la première partie de ce mémoire. Ces modèles phénoménologiques basés sur la théorie de la poromécanique et les acquis expérimentaux sont obtenus dans le cadre du modèle mathématique présenté par Gatmiri (1997) et Gatmiri et al. (1998). Après avoir présenté les modèles THHM et HHM, on établit pour la première fois les équations intégrales de frontière et les solutions fondamentales associées pour un milieu poreux non-saturé sous chargement quasi-statique pour les deux cas isotherme (2D dans le domaine de Laplace) et non-isotherme (2D et 3D dans les domaines de Laplace et temporel). Aussi, les équations intégrales de frontière ainsi que les solutions fondamentales 2D et 3D (dans le domaine de Laplace) pour le modèle dynamique couplé des sols non-saturés sont obtenues. Ensuite, les formulations d'éléments de frontière (BEM) basées sur la méthode quadrature de convolution (MQC) concernant les milieux poreux saturé et non-saturé sous chargement quasi-statique isotherme et dynamique sont implémentées dans le code de calcul « HYBRID ». Ayant intégrées les formulations de BEM pour les problèmes de propagation d'ondes ainsi que pour les problèmes de consolidation dans les milieux poreux saturés et non-saturés, il semble que nous ayons fourni à l'heure actuelle le premier code de calcul aux éléments de frontière (BEM) qui modélise les différents problèmes dans les sols secs, saturés et non-saturés. Une fois le code vérifié et validé, des études paramétriques portant sur des effets de site sismiques sont effectuées. Le but recherché est d'aboutir à un critère simple, directement exploitable par les ingénieurs, combinant les caractéristiques géométriques et les caractéristiques du sol, permettant de prédire l'amplification du spectre de réponse en accélération dans des vallées sédimentaires aussi bien que vides / The purpose of this dissertation is to develop a boundary element method (BEM) for multiphase porous media. Nowadays, the application of the BEM for solving problems of unsaturated porous media is still limited, because no fundamental solution exists in the published literature, neither in the frequency nor time domain. This fact rises from the complexity of the coupled partial differential equations governing the behaviour of such media. The developments of the BEM for the unsaturated soils carried out during this thesis are based on the thermo-hydro-mechanical (THHM) and hydro-mechanical (HHM) models presented in the first part of this dissertation. These phenomenological models are presented based on the experimental observations and with respect to the poromechanics theory within the framework of the suction-based mathematical model presented by Gatmiri (1997) and Gatmiri et al. (1998). After having presented the THHM and HHM models, for the first time, one establishes the boundary integral equations (BIE) and the associated fundamental solutions for the unsaturated porous media subjected to quasi-static loading for both isothermal (2D in the Laplace transform domain) and non-isothermal (2D and 3D in Laplace transform and time domains) cases. Also, the boundary integral equations as well as the fundamental solutions (2D and 3D in the Laplace transform domain) are obtained for the fully coupled dynamic model of unsaturated soils.In the next step, the boundary element formulations (BEM) based on the convolution quadrature method (CQM) regarding the saturated and unsaturated porous media subjected to isothermal quasi-static and dynamic loadings are implemented via the computer code HYBRID. Having integrated the BEM formulations for the wave propagation, as well as the consolidation problems in the saturated and unsaturated porous media, it seems that now the first boundary element code is obtained that can model the various problems in dry, saturated and unsaturated soils. Once the code is verified and validated, parametric studies on seismic site effects are carried out. The aim is to achieve a simple criterion directly usable by engineers, combining the topographical and geological characteristics of the soil, to predict the amplification of acceleration response spectra in sedimentary as well as hollow valleys

Les études de risque sismique en Iran sont dans leur phase initiale. Ainsi, il n'existait jusqu'à présent aucun catalogue d'enregistrements accélérométriques en Iran, et aucune étude sur l'ensemble des données iraniennes es n'avait été réalisée. Ce mémoire présente donc la première étude d'ensemble sur les mouvements forts en Iran. Un premier objectif a été d'établir un catalogue des mouvements forts d'Iran en faisant correspondre à chaque enregistrement de qualité suffisante, un événement sismique de localisation et magnitude connues ues, pour 279 enregistrements accélérométriques, analogues et digitaux. Cette correspondance a pu se faire sur la base des catalogues sismologiques internationaux (ISC, NEIC, ... etc.) ou nationaux. Pour 189 autres enregistrements numériques, pour lesquelles aucune information télésismique ou locale n'était disponible, les distance hypocentrales et les moments sismiques ont été estimés sur la base des enregistrements eux-mêmes. Une attention particuliere a été accordée aux effets de sites, en choisissant notamment 50 stations accélérométriques du réseau national oilles enregistrements ont été nombreux. Une nouvelle classification pour les effets de sites est proposée -basée sur les rapports H/V - qui s'avère relativement simple dans son principe, et plus fiable que les classifications déjà utilisées. Finalement, les lois d'atténuation pour différents paramètres de mouvement fort (dont Arms, Intensité Arias, PGA, PGV, et PGD, énergie, durée et valeurs spectrales) sont obtenues et discutées sur la base d'un total de 468 enregistrements accélérométriques en 3 composantes.

Les conditions locales de sol sont un des facteurs influençant le mouvement du sol lors d'un séisme. Cette thèse se propose de déterminer l'effet de ces conditions locales dans le cas de topographies et de vallées sédimentaires. Nous avons ainsi utilisé trois jeux de données acquis dans des contextes différents, que nous avons analysé par des méthodes expérimentales : les rapports spectraux par rapport à une station de référence, les rapports spectraux entre les composantes horizontales et verticales calculés sur des enregistrements sismiques et sur le bruit de fond. Nous avons aussi utilisé la méthode de fonctions de Green empiriques pour simuler le mouvement du sol lors d'un séisme majeur ainsi que des méthodes numériques pour calculer les fonctions de transfert des sites et les comparer aux résultats expérimentaux. Nous avons traité d 'une part les effets d'une topographie de grandes dimensions sur les mouvements sismiques, à partir d'une expérience qui s'est déroulée à Corinthe en 1993. Cette étude a montré que l'effet de topogra phies lisses et cylindriques sur le mouvement du sol est faible et que les fortes amplifications qui peuvent être constatées sont dues à des hétérogénéités de la géologie. D'autre part nous avons étudié les effets d 'une couche sédimentaire à partir de deux expériences à Grenoble en i995 et à Annecy en 1906. L'expérience de Grenoble a montré que le mouvement du sol dans la vallée est amplifié sur une bande de fréquence de 0.5Hz à 5Hz avec des amplifications de 10. L'étude fine de microzonage a pu mettre en évidence une variation de la fréquence de résonance avec l'emplacement de la station. La simulation de mouvements forts a montré que des séismes de magnitude 6.4 à des distances supérieures à 80km pouvaient amener des accélérations de l'ordre de 5% de g. L'étude à Annecy, suite au séisme d 'Epagny du 15 juillet 1996, a montré une corrélation entre les résultats des méthodes expérimentales et les dégats observés lors du choc principal.

L'objectif de cette thèse est, l'étude fine des effets d'amplification suite à la géologie locale, notamment dans des bassins sédimentaires, ainsi que l'évaluation des avantages et Iimites des méthodes d'estimations. La partie expérimentale est basée sur des données acquises au site test EURO-SEISTEST, le premier site-test en Europe occidentale proche de Thessalonique (Grèce). Après une évaluation des effets d'amplification au site même et une étude de leur stabilité en fonction de plusieurs paramètres, les estimations obtenues par cinq méthodes différentes (la méthode des rapports spectraux classique, la méthode des fonctions récepteurs, une inversion généralisée, une méthode de la Coda et la méthode de Nakamura) sont comparées entre elles. Les avantages et limites de chaque méthode sont évaluées. Les résultats sont dans la suite confrontés et comparés avec des modélisations uni - et bidimensionnelles. La variabilité spatiale à travers la vallée examinée est étudiée par la cohérence spatiale. La grande couverture spatiale par des stations installées au site a permis de séparer deux zones d'un comportement différent de la cohérence entre deux signaux dépendant de la distance qui les sépare. Ceci n'a jamais été remarqué dans des études précédents à cause d 'une couverture spatiale trop faible. Finalement, l'étude paramétrique des effets d'amplifications dans un milieu hétérogène en trois dimensions montre que l'effet géométrique est important et peut provoquer des amplifications significatives. La réponse d'un milieu plus complexe n'est difficilement à prévoir, dépend de l 'azimut de l'onde incidente et montre des amplifications importantes en présence d'un forte contraste d'impedence présent.

Ce travail est consacré à une étude numérique des effest produits sur les signaux sismiques par la structure géométrique des formations superficielles, dans une optique essentiellement tournée vers le risque sismique. Pour cela nous avons appliqué la méthode d'Aki et Larner ( milieux viscoélastiques linéaires homogènes isotropes) au calcul de la réponse de reliefs topographiques et de remplissages sédimentaires bidimentionnels ( 2D) à des ondes de volume planes incidentes. Npous nous sommes attachés à faire le lien entre les phénomènes physiques de base, les caractéristiques des spectres de réponse calculés en surface et les observations expérimentales. Le calcul montre clairement l'importance des phénomènes de diffraction et de formation d'ondes de surface sur toute variation latérale de la structure superficielle; de plus, pour les vallées alluviales, les ondes ainsi diffractées sont au moins partiellement piégées dans la couverture, et il en résulte sous certaines conditions géométriques des phénomènes de résonance bidimensionnelle extrêmement prononcés. Quant aux spectres de réponse, une étude paramétrique portant sur les facteurs sismiques, géométriques et mécaniques fait apparaître les résultats généraux suivants: - pour les reliefs topographiques, les effets, relativement large bande, sont maximaux pour des longueurs d'onde comparables à leur dimension horizontale. Les sommets des reliefs subissent alors une amplification, leur base une déamplification, et leurs versants d'importants mouvements différentiels. Cela est conforme aux résultats expérimentaux, bien que les amplifications observées soient souvent plus importantes que celles obtenues ici. - les vallées alluviales sont toujours sujettes à de fortes amplifications et prolongations des signaux sismiques, et leurs bords à d'importants mouvements différentiels. On peut cependant distinguer deux types de comportements: lorsqu'elles sont peu profondes, Ieur fréquence fondamentale est identique à celle prévue par l'approche 1 D, mais les amplificatIons correspondantes peuvent être 2 fois plus importantes, et les ondes de surface donnent lieu à des pics spectraux secondaires nombreux et rapprochés. Au contraire, la résonance 2D dans les vallées alluviales profondes survient à des fréquences reliées tant à leur largeur qu'à leur épaisseur, et l'amplification peut atteindre 4 fois celle prédite par l'approche 1D. La comparaison tant quantitative que qualitative avec les résultats expérimentaux s'avère tres bonne. Enfin, on discute l'applicabilité pratique de l'approche 2D et des resultats présentés ici au domaine du risque sismique. Comme le souligne en particulier un exemple d'étude de micro* zonage du site de Nice, les modèles 1D couramment utilisés en génie sismique peuvent non seulement sous-estimer l'ampleur des effets de site, mais en oublier "une des composantes essentielles, les mouvements différentiels .

Le comportement non-linéaire des sols observé lors des mouvements sismiques forts est maintenant bien établi et le déploiement des puits accélérométriques a permis des analyses détaillés de la propagation des ondes ainsi qu'une évaluation quantitative des paramètres physiques tels que la vitesse de cisaillement et de compression des ondes et les facteurs d'amortissements en fonction de la déformation. En dépit du nombre grandissant d'´etudes sur ce phénomène, sa connaissance est encore récente et les recherches sur les données de puits accélérométriques restent une étape importante vers la compréhension du comportement complexe in-situ des sédiments soumis à des mouvements sismiques forts. L'objectif de ces travaux est triple. Premièrement, un code d'inversion par algorithme génétique est développé afin d'inverser des données de puits accélérométriques via la théorie des matrices de propagation de Thomson- Haskell. Cette technique nous permet dans un premier temps de valider la structure en une dimension (1D) (e.g., vitesse des ondes de cisaillement, facteurs d' amortissements) d'un puits accélérométrique dans le domaine linéaire et dans un second temps de mettre en évidence de manière quantitative le comportement non-linéaire des sédiments lors du séisme de Fukuoka, 2005, Japon. Deuxièmement, les résultats de l'inversion sont utilisés pour tester des lois de comportement simples et avancées en utilisant la Méthode des Eléments Finis. Les résultats montrent clairement que l'hypothèse bilinéaire de la loi de comportement simple produit des séries temporelles non réalistes en vitesse et en accélération. L'utilisation d'une loi de comportement avancée mène à de meilleurs résultats, cependant, le nombre de paramètres ajustables pour obtenir des résultats consistants avec l'observation est un obstacle inévitable. Troisièmement, afin d'étendre l'étude des effets de site à des dimensions supérieures, des codes 2D et 3D de la Méthode en Eléments Spectraux sont développés et validés en comparant leurs résultats dans le domaine linéaire avec ceux obtenus théoriquement ou via d'autres méthodes numériques.

Ce travail est une étude des effets de site observés à Mexico, lors du séisme du 19 septembre 1985 (séisme de Michoacan), qui ont causés la mort de milliers de personnes et de effondrement de centaines de structures. Si l'amplification du mouvement est raisonnablement expliquée par le contraste d'impédance entre la couche d'argile superficielle et son substratum, la différence de durée entre les enregistrements obtenus dans les différentes zones de la ville n'a pas encore été expliquée de façon satisfaisante. Ce mémoire comprend deux parties: la première fait un bref rappel de l'histoire géologique du bassin de Mexico ainsi qu'une revue rapide des données géotechniques et géophysiques relatives à cette vallée et des données de mouvement fort. La deuxième partie est consacrée à la simulation numérique (en utilisant les méthodes d'Aki-Larner et des Différences Finies) des effets de diffraction et d'amplification dans ce bassin sous trois aspects différents: les effets à grande et petite échelle, et l'influence de la gravité. Nous obtenons une première estimation de l'amplification due à la structure à grande échelle et montrons que l'accroissement de la durée ne peut être expliquée en termes d'ondes de surface pour la valeur du facteur de qualité généralement admise. Concernant la petite éehelle, nous montrons que les varlatlons d'épaisseur de la couche d'argile, la présence d'une couche de remblai irrégulière en surface ainsi que des variations continues de vitesse au sein de la couche d'argile sont à l'origine d'amplifications locales significatives qui peuvent expliquer les différences observées entre stations proches, mais pas l'accroissement de la durée observée. . Finalement, nous montrons qu'il n'est pas possible de générer des ondes gravitaires dans un solide élastique, même dans le cas cl 'un coefficient de Poisson très proche de 0,5. Par contre, des ondes gravitaires sont générées de façon efficace dans une couche fluide non visqueuse mais sont rapidement atténuées dans le cas visqueux.