Dissertations / Theses on the topic 'Mécanismes aux foyers'

La première partie de ce mémoire est une mise au point sur la notion de contraintes en tectonique cassante. Nous présentons une classification schématique des différents types d'états de contraintes que l'on trouve dans la croûte supérieure, et des déformations qui leur sont associées. Pour chaque état de contraintes, on analyse la cinématique des plans néoformés et des plans réactivés. Nous discutons ensuite à l'aide d'exemples, les limites d'application des méthodes inverses pour la détermination des états de contraintes (exemple de la méthode Etchecopar et coll.). On montre également que la prise en compte de la typologie des stries (stries pénétratives avec dissolution ou stries avec recristallisation), permet de mieux contraindre la détermination du tenseur. Dans la deuxième partie, la tectonique récente des Alpes du Sud est abordée par l'étude des paléocontraintes. 150 stations microtectoniques réparties dans les nappes de Digne, de Castellane et de Nice sont analysées. Les résultats permettent de reconstituer l'évolution du champ de contraintes depuis l'Oligocène. Celle-ci dépend des conditions aux limites propres à chaque domaine. Ainsi, la variation rapide des directions de contraintes dans les Arcs de Digne et de Nice est liée à la présence de bordures libres situées au front de ces deux nappes (respectivement les bassins de Valensole et du Var). Ces bassins ont permis l'expression de phénomènes gravitaires dans les parties internes des nappes, soumises à un uplift. La nappe de Castellane bloquée à son front par le massif des Maures montre une plus grande stabilité des directions de contraintes. La reconstruction du champ de contraintes pour la période quaternaire nous permet de déterminer la cinématique et les mécanismes aux foyers possibles des accidents qui affectent les différents domaines. Ces résultats sont comparés aux données de la sismicité historique et instrumentale.

Etablissement d'une methode de determination des mecanismes au foyer avec une station unique appliquee a des seismes euro-mediterraneens. Un seul enregistrement 3 composantes longue periode suffit a inverser le tenseur des moments sismiques ekstrom et al (1986). La principale difficulte est l'estimation des termes de propagation. Une analayse de dispersion des ondes de love et rayleigh pour des periodes (30-100 sec) subit un filtrage multiple par analyse temps-frequence. L'integration de la vitesse de groupe apparente permet d'obtenir la vitesse de phase. Les 5 composantes du tenseur des moments sismiques sont inversees; les parametres de faille: plongement, glissement et azimuth sont calcules. Un sismographe 3-composantes large-bande de type geotech a ete utilise. On utilise les ondes de love et de rayleigh par comparaison aux resultats du national earthquake information service. Cette methode s'applique a des seismes de faible intensite (mo10#2#4 dynes?cm) et a des donnees anciennes

L'exploitation pétrolière ou toute injection de fluide dans un milieu géologique peut imposer des variations de contraintes pouvant générer une micro-sismicité. Nous examinons si l'utilisation des modèles de frottement dynamique permettrait une meilleure détermination de la mécanique sur le plan de faille à partir de sismomètres placés en profondeur. Nous avons cherché à identifier le modèle de loi de frottement à partir des sismogrammes. Avant de traiter ce problème inverse, nous avons examiné le problème direct en construisant des sismogrammes synthétiques sur des stations virtuelles, sur la base d'une loi frottement postulée, en utilisant le code de simulation de séismes FD3D et le code de rayonnement AXITRA. L'utilisation de deux lois différentes, dont une variation originale issue de la loi Dieterich-Ruina, conduit à des sismogrammes semblables. Ceci pourrait indiquer qu'une inversion systématique à partir de l'unique observation de sismogrammes pourrait s'avérer difficile.

Le développement de nouvelles technologies pour le transport aérien moins polluant est de plus en plus basé sur la simulation numérique, qui nécessite alors une description fiable de la chimie. Pour la plupart des carburants, la description de la combustion nécessite des mécanismes détaillés mais leur utilisation dans une simulation numérique de combustion turbulente est limitée par le coût calcul. Des mécanismes cinétiques réduits et des méthodes de tabulation ont été proposés pour surmonter ce problème. Ces descriptions chimiques simplifiées ayant été développées dans le cadre de configurations laminaires, cette thèse propose de les évaluer dans des configurations turbulentes: une DNS de flamme prémélangée méthane/air de type Bunsen et une LES d'un brûleur expérimental. Les mécanismes sont analysés en termes de structure de flamme, paramètres de flamme globaux, longuer de flamme, prediction des concentrations en espèces majoritaires et des émissions polluantes. Une méthodologie pour évaluer a priori la capacité d'un mécanisme à prédire correctement des phénomènes chimiques tridimensionnels est proposée en se basant sur les résultats de flammes laminaires monodimensionnelles non étirées et étirées. Il ressort que, d'une part, pour construire un mécanisme réduit, il est nécessaire de faire un compromis entre coût calcul, robustesse et qualité des résultats. D'autre part, la qualité des résultats de DNS et LES de configurations tridimensionnelles turbulentes peut être anticipée par une analyse du comportement des schémas réduits dans des configurations simplifiées de flammes monodimensionnelles laminaires non étirées et étirées.

Nous proposons une méthode d'inversion non linéaire des amplitudes des ondes directes P, SV et SH basée sur l'algorithme du recuit simulé, afin de déterminer à partir d'un nombre limité de stations, le mécanisme au foyer de séismes induits en contexte de réservoir. Cette méthode permet de déterminer aussi bien les paramètres du plan de faille (azimut, pendage et angle de glissement) décrivant une source double-couple, que les six composantes du tenseur des moments décrivant une source plus générale.
L'inversion double-couple et l'inversion du tenseur des moments sont testées sur quatre séismes induits dans le réservoir géothermique de Soultz-sous-Forêts. Les mécanismes obtenus sont en accord avec ceux déterminés par Charléty et al. 2007.
La méthode d'inversion est appliquée à la sismicité de la saline de Vauvert enregistrée par un réseau permanent de deux stations 3-composante. Dans un premier temps, la méthode est testée sur 15 séismes enregistrés durant le déploiement temporaire d'une antenne de quatre capteurs 3-composante. La comparaison entre l'inversion utilisant les deux stations permanentes et l'antenne temporaires et l'inversion utilisant seulement les deux stations permanentes montre des mécanismes au foyer double-couple identiques pour les séismes localisés entre les deux stations permanentes. Dans un deuxième temps, la méthode est appliquée à un essaim de sismicité. Les mécanismes au foyer obtenus pour 532 évènements, indiquent pour la majorité une rupture le long de fractures sub-verticales NE-SW, interprétée comme de probables ruptures sur les plans stratigraphiques des bancs d'insolubles intercalés dans la formation de sel.

ANALYSE SISMOTECTONIQUE DE LA ZONE SISMIQUE SUD-ISLANDAISE. La Zone Sismique Sud-Islandaise (ZSSI) est une zone transformante sénestre E-W, d’environ 70-80 km de long et 20-30 km de large, localisée dans le sud de l’Islande à la latitude 64°N. La ZSSI est le siège d’une importante sismicité, avec notamment l’occurrence de deux séismes majeurs le 17 juin (Mw= 6. 5), et le 21 juin 2000 (Mw= 6. 4) plus à l’ouest. Une inversion des mécanismes au foyer a été réalisée sur 33878 séismes, correspondant à la période comprise entre le 01 juillet 1991 et le 31 décembre 2004, afin d’identifier les régimes de contraintes dans la ZSSI depuis l’installation du réseau sismique islandais en juillet 1991 (réseau SIL). Cette inversion révèle la présence de deux régimes de contraintes opposés. Le régime principal, qui se définit par une compression NE-SW et une extension NW-SE. Le régime secondaire est perpendiculaire et correspond à une compression NW-SE et une extension NE-SW. L’étude spatio-temporelle de la sismicité montre que ces deux régimes sont liés de façon complexe, ce qui suggère que la réponse élastique de la croûte à la déformation peu profonde, incluant le rebond élastique et la chute des contraintes, joue un rôle prépondérant dans la permutation des axes de contraintes σ1 et σ3. La reconstitution du déplacement co-sismique moyen, faite à partir d’une étude minutieuse des «push-ups», avec mesures de terrain, micro-gravimétrie et essai de modélisation indiquent des valeurs de déplacements comparables comprises entre 1. 1 et 1. 6 m. L’évaluation des magnitudes des séismes responsables de ces différentes ruptures a été faite par comparaison avec les traces des séismes instrumentés récents. Elle démontre que par le passé des séismes de magnitude égale à 7 ou supérieure se sont produits dans la ZSSI. Le cisaillement simple profond, lié au fonctionnement de la transformante et affectant les couches les plus ductiles, prévaut dans la ZSSI le long d’un corridor d’environ 20 km de large. En surface, il semblerait que ce cisaillement s’exprime par une rotation de blocs. La distribution des failles et les dimensions reconstituées des blocs sont cohérentes avec un modèle d’évolution structurale simple de la Zone Sismique Sud-Islandaise, permettant d’apporter une interprétation tectonique précise dans un type de zone transformante océanique généralement inaccessible à l’observation terrestre
In southern Iceland, the South Iceland Seismic Zone (SISZ) is an E-W trending, left-lateral active transform zone at latitude 64°N. The SISZ is approximately 20-30 km wide and 70-80 km long. The last major seismic crisis recorded in the SISZ occurred in June 2000. A first major event occurred on June 17th, 2000 (Mw= 6. 5) and was followed by a second event on June 21st, 2000. The mass inversion aimed at to identifying the seismotectonic stress regimes that best account for the whole set of data. The inversion reveals consistent orientations, for both a primary strike-slip regime with NE-SW compression and NW-SE extension and a secondary strike-slip regime with NW-SE compression and NE-SW extension. A simple permutation, or switch, occurs between the extreme stress axes of these two regimes (maximum and minimum principal stresses), so that they should be considered as mechanically ‘opposite’. The study of the space and time evolution of the two stress regimes, at the regional scale of the SISZ, indicates that the corresponding earthquakes are highly intricate in space and time. This suggests that the secondary regime does not correspond to distinct regional seismotectonic mechanisms but mainly results from the effects of elastic rebound, stress drop and fluid migration. However, despite homogeneity revealed in the first approximation by the stress analysis of the whole SISZ, the study of stress tensors at a local scale around the faults reveals some interesting spatial variations. The results indicate a larger average deviation of the stress between 5 to 10 km, up to about 40 degrees instead of about 20 degrees at shallower depths. This difference in angular dispersion may be attributed to an effect of the higher fluid pore pressure beneath the depth of about 5 km. The three-dimensional structural analysis of the push-up structures reveals the amount of displacement along earthquake strike-slip faults and indicates similar values for each fault in a range of 1. 1-1. 6 m. The determinations of the magnitudes from co-seismic displacement evaluation demonstrate that in the past, the earthquakes may have reached magnitude equal to 7 or more in the SISZ. The proximity of the Icelandic Mantle Plume, marked by a positive thermal anomaly, certainly affected the rheology of the upper lithosphere and hence the structural development of the transform deformation zone. Simple shear probably prevails within a 20 km wide corridor, with a typical pattern of brittle deformation in the uppermost crust and a more pervasive deformation in the underlying, partly viscous layer. The block faulting at the surface is quite consistent with a Riedel type distribution. Both the fault distribution and the effective block dimensions are consistent with a simple structural evolution model

La caractérisation de la source d’un séisme se fait à partir de l’analyse des mesures des déplacements transitoires et statiques du sol, et dépend de la quantité et de la qualité de ces mesures. Nous avons travaillé sur la détermination des mécanismes au foyer des répliques du séisme de Saintes (MW 6.4, 2004), et sur la détermination de la distribution spatio-temporelle du glissement cosismique des séismes de L’Aquila (Mw 6.3, 2009), et de Miyagi-Oki (Mw 7.2, 2005) et de Sanriku-Oki (Mw 7.3, 2011). Ces travaux se sont basés sur des méthodes d’inversions, et différents jeux de données (accélérométriques, large-bandes, GPS et InSAR) accessibles ou non selon le séisme considéré. La seule diversité des mesures n’est pas suffisante pour décrire la rupture. La modélisation des données se confronte à des difficultés, comme par exemple la pertinence des modèles de vitesses sismiques pour la modélisation des données accélérométriques. Une autre problématique récurrente est la non-unicité de la meilleure solution déterminée par les méthodes d’inversions pour décrire les données. Pour répondre à ces deux problématiques, nous avons d‘une part développé une procédure d’exploration de modèles de vitesse pour déterminer les valeurs optimales capables de décrire au mieux les données accélérométriques du séisme de L’Aquila. D’autre part, nous avons développé une procédure de construction d’un modèle de source moyen que nous avons appliqué pour la détermination du glissement cosismique des séismes de L’Aquila, de Miyagi-Oki, et de Sanriku-Oki. L’ensemble de ces travaux et les réponses aux problèmes soulevés sont présentés dans ce travail de thèse
Studies of the earthquake source are based on observations of seismic ground motions. They also depend on the quality and the density of measurements. In this present work we will present studies of the determination of focal mechanism of main aftershocks of the Les Saintes (MW 6.4, 2004) earthquake, and the determination of the coseismic slip of the L’Aquila (MW 6.3, 2009), the Miyagi-Oki (MW 7.2, 2005), ant the Sanriku-Oki (MW 7.3, 2011) earthquakes. These studies were based on two inversion methods. Different kinds of data were available (strong motion, broadband teleseismic, GPS and InSAR) depending on the earthquake studied. But the multiplicity of data is not sufficient to well describe rupture process. There are others difficulties as the data modeling of strong motion. Seismic velocity models are used to describe the characteristics of layers crossed by seismic waves. The quality of the modeling is depending on the pertinence of these seismic velocity models. The description of the rupture process is also depending on the non-uniqueness of the best solution given by global inversion methods. We propose two procedures in order to take into account these two classic issues. First, we developed a velocity model exploration procedure to obtain optimized 1D velocity models in order to improve the strong motion modeling of the L’Aquila earthquake. Then we developed a procedure to build an average rupture model from the combined results of several joint inversions, which was applied to the L’Aquila, the Miyagi-Oki, and the Sanriku-Oki earthquake. This thesis presents all these works and answers to the raised issues

Le détroit de la sonde constitue une zone de transition entre deux modes différents de subduction: une subduction frontale au niveau de Java et une subduction oblique au niveau de Sumatra. Par conséquent, le détroit de la Sonde est c0nsidéré comme une, zone clé pour la compréhension des processus géodynamiques de la région. Dans le but de mieux contraindre la distribution des séismes et l'état actuel des contraintes dans cette région, une campagne de microsismicité a été effectuée dans le détroit de la Sonde pendant l'été 1984. Les résultats de cette étude montrent clairement que la sismicité crustale est concentrée dans trois essaims: 1) sous le volcan Krakatoa où les séismes enregistrés ont une origine tectonique, 2) dans le graben situé à l'ouest du détroit et 3) au sud de Sumatra où la concentration est plus ou moins diffuse. Les mécanismes focaux individuels des séismes enregistrés dans le détroit ainsi que les mécanismes composites révèlent un régime en extension. La détermination d'un tenseur de contraintes à partir des séismes du groupe du Krakatoa montre que la direction d'extension est N130°E. Ce résultat confirme que le détroit de la Sonde est actuellement en régime tectonique extensif du fait du mouvement vers le Nord-Ouest de la plaque d'avant arc de Sumatra le long de la faille de Semangko. Deux réservoirs magmatiques sous le complexe du Krakatoa ont été détectés grâce à l'étude des anomalies de propagation des ondes S: le premier est situé à une profondeur de l'ordre de 9 km et est probablement composé de plusieurs petits réservoirs; le deuxième est plus profond (environ 22 km de profondeur) et doit avoir une grande extension latérale. L'existence de ces deux réservoirs magmatiques ainsi que le contexte extensif du détroit peuvent expliquer les variations dans la composition des produits volcanique émis lors des différentes éruptions du volcan.

La caractérisation de la source d’un séisme se fait à partir de l’analyse des mesures des déplacements transitoires et statiques du sol, et dépend de la quantité et de la qualité de ces mesures. Nous avons travaillé sur la détermination des mécanismes au foyer des séismes modérés de Tadjena (Mw 5.0, 2006), Béni-Ilmane (Mw 5.5, 2010), Hammam Melouane (Mw 5.0, 2013), Bordj-Ménaïel (Mw 4.1, 2014), Alger (Mw 5.7, 2014) et M’ziraa (Mw 5.1, 2016). Pour cela, nous avons inversés les formes d’ondes en champ proche et régional des stations large bandes et accélérométriques du réseau sismologique algérien (ADSN). Nous avons déterminé la distribution spatio-temporelle du glissement cosismique du séisme de Boumerdes-Zemmouri (Mw 6.8, 2003) en inversant conjointement l’ensemble des données disponibles, sismologiques (télésismiques, accélérometriques) et géodésiques (GPS, InSAR, soulèvement côtier). Nous avons travaillé sur la relation qui peut exister entre le séisme de Boumerdes-Zemmouri et les séismes de Hammam Melouane, Bordj Ménaïel et Alger (Mw 5.7, 2014) en termes de transfert de contrainte (CFF). Nous avons calculé le champ de contrainte dans différentes régions de l’Algérie par inversion des mécanismes aux foyers disponibles dans chaque région
Studies of the earthquake source are based on observations of seismic and static ground motions, which depend on the quality and the density of measurements. In this thesis, we present studies of the determination of the focal mechanism of the Tadjena (Mw 5.0, 2006), Béni-Ilmane (Mw 5.5, 2010), Hammam Melouane (Mw 5.0, 2013), Bordj-Ménaïel (Mw 4.1, 2014), Alger (Mw 5.7, 2014) and M’ziraa (Mw 5.1, 2016) earthquakes, by modeling waveforms in the near and regional field with broadband and strong motion stations of the Algerian Digital Seismic Netwok (ADSN). In addition, we determined the coseismic slip distribution of the Boumerdes-Zemmouri earthquake (Mw 6.8, 2003) by inverting a most comprehensive set of data (teleseismic, strong motion, GPS, InSAR, coastal uplift). We calculated the Coulomb Failure Function between the Boumerdes-Zemmouri earthquake (source fault) and the Hammam Melouane, Bordj Ménaïel and Algiers (Mw 5.7, 2014) events (receiver faults). We computed the stress tensor in different regions of Algeria obtained from the inversion of the available focal mechanisms

La localisation de séismes induits ainsi que les mécanismes au foyer associés sont des outils fréquemment utilisés afin, entre autres, d’imager la structure d’un réservoir. Cette thèse présente une technique permettant de quantifier les erreurs associées à ces deux paramètres. Par cette méthode, incertitudes et imprécisions sont distinguées. La méthode a été appliquée aux sites de Soultz et de Rittershoffen pour étudier l’impact de plusieurs critères sur la localisation de la sismicité induite. Ainsi, il a été montré que l’utilisation de capteurs installés profondément dans des puits et qu’une bonne couverture sismique azimutale réduit sérieusement les incertitudes de localisation. Les incertitudes du modèle de vitesse, représentées par une distribution gaussienne des modèles avec un écart de 5% autour du modèle de référence, multiplient les incertitudes de localisation par un facteur 2 à 3. Des simplifications utilisées pour le calcul ou une mauvaise connaissance du milieu peuvent mener à des imprécisions de l’ordre de 10% spatialement non isotropes. Ainsi, les structures du sous-sol peuvent être déformées dans les interprétations. L’application d’un tir de calibration peut néanmoins corriger ce fait en grande partie. L’étude d’erreurs associées aux mécanismes au foyer ne semble cependant pas conduire aux mêmes conclusions. Le biais angulaire peut certes être augmenté par l’omission de la faille dans le modèle de vitesse, mais dans plusieurs cas il est le même que dans le cas idéal voire diminué. En outre, une meilleure couverture sismique améliorerait toujours le mécanisme au foyer obtenu. Ainsi, il n’est pas conseillé d’imager un réservoir en n’utilisant que la localisation de séismes, mais une combinaison de plusieurs paramètres sismiques pourrait s’avérer efficace. La méthode appliquée dans le cadre de cette thèse pourra servir pour d’autres sites à condition d’en avoir une bonne connaissance a priori
Induced seismicity location and focal mechanisms are commonly used to image the sub-surface designin reservoirs among other tasks. In this Ph.D. the inaccuracies and uncertainties on earthquake location and focal mechanisms are quantified using a three-step method. The technique is applied to the geothermal sites of Soultz and Rittershoffen to investigate the effect of several criteria on thee arthquake location. A good azimuthal seismic coverage and the use of seismic down-hole sensors seriously decrease the location uncertainty. On the contrary, velocity model uncertainties, represented by a 5% Gaussian distribution of the velocity model around the reference model, will multiply location uncertainties by a factor of 2 to 3. An incorrect knowledge of the sub-surface or the simplifications performed before the earthquake location can lead to biases of 10% of the vertical distance separating the source and the stations with a non-isotropic spatial distribution. Hence the sub-surface design maybe distorted in the interpretations. To prevent from that fact, the calibration shot method was proved to be efficient. The study on focal mechanism errors seems to lead to different conclusions. Obviously, the angular bias may be increased by neglecting the fault in the velocity. But, it may also be the same as or even smaller than the bias calculated for the case simulating a perfect knowledge of the medium of propagation. Furthermore a better seismic coverage always leads to smaller angular biases. Hence,it is worth advising to use more than only earthquake location in order to image a reservoir. Other geothermal sites and reservoirs may benefit from the method developed here
Die korrekte Lokalisierung von induzierter Seismizität und den dazugehörigen Herdflächenlösungensind sehr wichtige Parameter. So werden zum Beispiel die Verteilung der Erdbeben und die Orientierung ihrer Herdflächenlösungen dazu benutzt um in der Tiefe liegende Reservoirs zulokalisieren und abzubilden. In dieser Doktorarbeit wird eine Technik vorgeschlagen um diemethodisch bedingten Fehler zu quantifizieren. Mit dieser Methode werden die verschiedenen Fehlerquellen, die Unsicherheiten und die Fehler im Modell getrennt. Die Technik wird für die geothermischen Felder in Soultz und in Rittershoffen benutzt um den Einfluss verschiedener Parameter (Annahmen) auf die Lokalisierung der induzierten Seismizität zu bestimmen. Es wurde festgestellt, dass Bohrlochseismometer und eine gute azimutale Verteilung der seismischen Stationen die Unbestimmtheiten verkleinern. Die Geschwindigkeitsunbestimmheiten, die durch eine Gauss-Verteilung mit 5% Fehler dargestellt werden, vervielfachen die Lokalisierungsungenauigkeiten um einen Faktor 2 bis 3. Eine ungenaue Kenntnis des Untergrunds oder die verwendete vereinfachte Darstellung der Geschwindigkeitsverhältnisse im Untergrund (notwendig um die synthetischen Rechnungen durchführen zu können) führen zu anisotropen Abweichungen und Fehlern in der Herdtiefe von bis zu 10%. Diese können die Interpretationen des Untergrunds deutlich verfälschen. Ein “calibration shot” kann diese Fehler korrigieren. Leider können die Fehler für die Herdflächenlösungen nicht in derselben Weise korrigiert werden. Es erscheint daher als keine gute Idee, ein Reservoir nur über die Lokalisierung von Erdbeben zu bestimmen. Eine Kombination mehrerer seismischer Methoden scheint angezeigt. Die hier besprochene Methode kann als Grundlage dienen für die Erkundung anderer (geothermischer)

La sismicité induite par le gisement de gaz de Lacq a débuté en 1969, dans une région tectoniquement considérée comme asismique auparavant (M<3). L'exploitation est située au Sud du bassin Aquitain, à 30 km ad Nord de la chaîne pyrénéenne. Un réseau local a été installé en 1974 afin de suivre l'évolution spatio-temporelle des séismes *Plus de 1000 événements (1.0 <=MS <=4.2) ont été enregistrés entre 1974 et 1992. La très bonne connaissance de la structure géomécanique qui piège le réservoir (forages, profils sismiques, carottages ... ) a permis de construire un modèle de vitesse 3D permettant de relocaliser de manière fine la sismicité sur la période 1976- 1989. La position des hypocentres est étroitement corrélée avec la structure en dôme du réservoir ainsi qu'avec les discontinuités préexistantes. La construction des mécanismes au foyer (133 événements sélectionnés sur la période 1976-1989) montre une grande diversité dans les mécanismes de rupture. Cependant, la méthode des dièdres droits a permis de voir une évolution temporelle globale à l'échelle du gisement. Dans un premier temps (1976-1986) nous avons eu un régime en compression maximale verticale, suivi dans la seconde période (1987-1989) par un régime en compression maximale horizontale. L'inversion du tenseur des contraintes a confirmé ces résultats. La sismicité induite de Lacq étant très localisée (du fait de l'extension limitée du gisement), cela nous a permis de considérer le gisement comme un objet géologique isolé, à une échelle intermédiaire entre les processus tectoniques et les expériences de laboratoire. Nous n'avons pas retrouvé les résultats expérimentaux obtenus sur le coefficient b (relation fréquence-magnitude des séismes) car la présence d'événements asismiques (non décelables) sur Lacq semble jouer un rôle majeur dans la libération des contraintes. La géométrie fractale nous a permis de distinguer deux classes de séismes: (1) des séismes diffus répartis sur tout le gisement et (2) des nids de sismicité corrélés avec la structure en dôme du réservoir. Bien qu'il n'y ait pas de corrélation temporelle entre D et b, nous proposons cependant qu'il existe un lien géomécanique entre ces deux paramètres.

L'interprétation d'un secteur des -Alpes occidentales externes en terme de paléomarge continentale téthysienne nous a amené à considérer le secteur du haut-Var sous deux angles différents: -d'une part son évolution stratigraphique et structurale locale. -d'autre part son incorporation dans le schéma général du bassin subalpin et du domaine téthysien. L'analyse des variations latérales de faciès permet de reconstituer un schéma du secteur du haut-Var cohérent à condition de placer un accident de direction N 20°-30° à vergence Ouest ayant joué en faille normale au cours du Lias et du Dogger et limitant un domaine de bassin à l'Ouest et une plate-forme à l'Est. Les manifestations tectonosédimentaires observées suggèrent, outre les déformations dues à la tectonique distensive, l'existence d'une halocinèse précoce le long du faisceau de failles N 20°-30°. Ainsi au Jurassique, la région du Dôme de Barrot serait limitée à l'Ouest par cette faille "Ouest-Barrot", et correspondrait au symétrique de la marge cévenole par rapport au bassin subalpin. C'est donc plus à l'E et au SE qu'il faut rechercher les blocs basculés de la marge téthysienne proprement dite. La géophysique: La néotectonique de la région Sud des Alpes Occidentales et de la mer ligure, évaluée grâce au calcul de 20 mécanismes au foyer de séismes récents nous permettent de localiser les zones de concentration de la déformation. Ces résultats, associés à d'autres données géophysiques, nous permettent de constater la présence d'une bande étroite à comportement crustal anormal, orientée N 120°- 130° et qui s'étend depuis le massif du Pelvoux jusqu'à Savone (Ligurie occidentale), en passant au NE du massif de l'Argentera. Cette zone est considérée comme un grand linéament structural issu d'orogénèses antérieures et qui fonctionnerait actuellement en décrochement dextre. La géodynamique: La confrontation des données paléogéographiques mésozoïques d'une part, et des résultats géophysiques et structuraux d'autre part, nous incite à envisager l'existence d'une zone transformante à jeu dextre au Jurassique. Dans ce schéma, le secteur étudié appartiendrait à la bordure orientale du bassin subalpin intramarge tout en ayant un regard sur la marge transformante téthysienne située plus au Nord.

Un réseau local de 44 stations sismologiques portables installé dans le sud de la Grèce nous a permit d'enregistrer plus de 800 séismes. Le but de cette campagne était d'étudier la tectonique de l'arc hellénique et la géométrie de la subduction. Les séismes crustaux sont localisés en majorité entre les fosses helléniques et l'arc égéen méridional, sans montrer de lacune de sismicité, au niveau du gap d'Anticythère. Les fosses de Pline et de Strabon ont très peu de sismicité. Nous avons recensé une centaine de séismes intermédiaires, les séismes subcrustaux se répartissent sous l'arc égéen méridional, tandis que les séismes profonds sont situés à l'est du Péloponnèse et sous Astipaléa. La Mer de Crète est une région asismique alors que la plupart des modèles indiquent qu'elle est la zone où l'extension est maximale. On ne peut pas préciser la géométrie de la subduction puisqu'il n'y a pas continuité d'est en ouest. 98 mécanismes au foyer ont pu être déterminés dans le sud du domaine égéen. Les mécanismes des séismes crustaux confirment les données de microtectonique et de mécanismes au foyer de séismes majeurs; on distingue en plus de la compression E-W au niveau de l'île de Samos. Les séismes subcrustaux montrent un champ de contraintes différent et très complexe. Au contraire, les séismes les plus profonds ont tous des mécanismes avec des axes T plongeant vers le nord parallèles au plongement supposé de la plaque subduite. Les vecteurs glissement sont N-S, parallèles à la direction de convergence entre l'Afrique et l'Eurasie.

Ce mémoire d'habilitation décrit ma recherche effectuée en vue de mieux comprendre le mécanisme des tremblements de terre. Les ondes sismiques sont les principales données qui donnent accès à cette information, et mes travaux sont donc en lien direct avec les observations des capteurs sismiques (ou géodésiques) capables de détecter ces ondes. Les qualité, quantité et type de données disponibles m'ont conduit à aborder cette thématique sous trois angles différents, qui recoupent les trois chapitres de ce mémoire. A l'échelle mondiale, les capteurs large-bande du réseau global permettent une étude systématique et homogène de tous les séismes dont la magnitude est suffisante pour générer des ondes clairement analysables. Cette magnitude seuil est de l'ordre de 5.5-6 lorsque l'on analyse les ondes de volume se propageant dans la Terre. Par ailleurs, ces données sont accessibles en temps réel, ce qui permet d'obtenir des informations rapides sur les séismes, à condition que des techniques d'analyse automatisées soient mises en place. Ce double intérêt de l'approche globale -systématisme et rapidité- est développé dans le Chapitre I de ce mémoire. L'échelle mondiale trouve ses limites lorsqu'on s'intéresse aux détails du processus de rupture. Alors que les données mondiales sont suffisantes pour imager les caractéristiques moyennes de source (mécanisme au foyer, profondeur, magnitude de moment, fonction source), elles ne permettent généralement pas d'extraire avec précision les informations internes à la rupture sismique (distribution de glissement, vitesse de rupture locale). Pour ce faire, il est nécessaire d'utiliser des données plus proches de l'événement sismique, et de développer des techniques d'analyse adaptées à chaque configuration (analyses en réseau, analyses comparées entre le séisme principal et l'un de ses précurseurs...). Par ailleurs, et contrairement à l'échelle globale - où les données large-bande sont les données " reines " de l'analyse de la source -, l'analyse à distance plus proche permet et/ou requiert de diversifier les données utilisées. Cela peut passer par l'utilisation des capteurs accélérometriques, ou par la nouvelle utilisation des capteurs GPS en tant que sismomètres. Le Chapitre II s'intéresse à cette thématique de résolution fine des propriétés de source, avec un intérêt particulier pour la détermination de la vitesse de rupture. La mise en évidence récente des vitesses de rupture supershear (plus rapides que les ondes de cisaillement) et de leurs conséquences constituent un point important de ce chapitre. Enfin, il est intéressant de replacer la rupture sismique dans le cadre plus large des conditions qui stimulent ou inhibent son déclenchement. Dans cet esprit, je présente dans le Chapitre III des éléments de réponse venant de l'observation continue de la zone de subduction Equatorienne. Cette observation - incluant capteurs large-bande, accéléromètres, et GPS - se déroule depuis 2008 dans le cadre du projet ANR ADN, soutenu par l'IRD et en collaboration avec l'Institut de Géophysique à Quito. Nous montrerons dans ce chapitre comment l'occurrence de séismes lents sur l'interface de subduction est un facteur déclenchant de la sismicité.

La sismicité de la terminaison Nord-Ouest de l'Arc égéen est analysée à partir de données microsismiques. Celles-ci ont été recueillies par 52 stations qui ont fonctionné, en été 1989, sur les îles Ioniennes, en Acarnanie et en Epire. 656 microséismes ont été localisés et 128 solutions focales construites. La région étudiée constitue la zone de transition entre subduction océanique égéenne (au Sud de Céphalonie) et subduction continentale adriatique (au Nord de Zante). L'activité microsismique y est superficielle (0-40 km). Les microséismes observés à l'Est de Zante (z = 30-40 km) sont associés à la subduction océanique égéenne. Au Nord de Zante, la subduction continentale adriatique est sismiquement peu marquée ; les microséismes profonds (15- 40 km) jalonnent. en particulier, les failles N160° sénestres du couloir étolien et l'accident N20° dextre de Lixourion-Leucade. Le jeu de cette structure se manifeste sur Céphalonie par la présence de microséismes le long de la faille d'Ano-Angon ; ce jeu est, par ailleurs, associé à l'activité microsismique de la faille de l'Ainos. L'analyse des mécanismes focaux microsismiques par la méthode des dièdres droits révèle l'existence d'un raccourcissement horizontal profond N65° - N80° affectant la zone étudiée. Dans les premiers kilomètres de la couverture sédimentaire, les solutions focales indiquent la présence d'étirements de nature et d'orientations variées. En profondeur, au raccourcissement ENE-WSW, s'ajoute un étirement horizontal perpendiculaire. Tout se passe comme si un bloc apparemment asismique limité par le décrochement de Lixourion-Leucade, le couloir étolien et les failles normales du golfe Ambracique, subissait un mouvement vers le Sud. Le raccourcissement se traduit par les jeux dextre inverse du décrochement de Lixourion-Leucade et sénestre des failles du couloir étolien. D'autre part, les solutions focales microsismiques indiquent l'ouverture N-S du fossé d'Agrinion (lac Trikhonis) qui est connectée, par le jeu sénestre des failles du couloir étolien, à celle (apparemment asismique) du golfe Ambracique.

L'importante sismicité qui caractérise la partie Est du bassin Méditerranéen est l'expression de la déformation intense de ce domaine qui se trouve pris entre la convergence des plaques d'Afrique, d'Europe et de Turquie . Le rapprochement de l'Afrique par rapport à l'Europe est une rotation anti-horaire autour d'un pôle situé dans l'Atlantique . Ce déplacement qui est continu depuis à peu près 70 Ma est de l'ordre de 0,5 cm/a à l'Ouest (Gibraltar) et de 1,0 cm/a à l'Est (au niveau de l'Arc Hellénique) où l'Afrique plonge sous la lithosphère égéenne. Cette subduction égéenne, initiée au Miocène supérieur, est mise en évidence par la distribution des séismes intermédiaires le long d'un arc situé en Mer de Crète, entre le Péloponnèse et Rhodes. Elle plonge vers l'intérieur de l'Arc . Un peu plus au Nord de l'arc sédimentaire, la montée du matériel chaud du manteau due à la subduction a créé une série de volcans qui forment l'arc volcanique des Cyclades. A l'extérieur de l'arc, le domaine égéen est bordé par les dépressions profondes (-5200 m) des fosses ioniennes à l'Ouest et des fosses de Pliny et Strabo au Sud. Située au Sud des fosses helléniques, la Ride Méditerranéenne joue le rôle d'un prisme d'accrétion . À l'intérieur de l'arc, l'extension domine. Ses directions témoignent d'un mécanisme d'écoulement gravitaire au dessus de la plaque africaine ce qui entraîne un amincissement important de la lithosphère égéenne . La particularité de l'Arc Hellénique par rapport aux arcs insulaires typiques est le fait que le taux de l'extension à l'intérieur de l'arc est beaucoup plus important que la vitesse de convergence entre l'Afrique et l'Europe. Le rapprochement entre l'Afrique et l'Europe se fait avec une vitesse approximative de 1 cm/a dans la direction N-S, alors que la vitesse de convergence de l'arc égéen par rapport à l'Afrique, à travers le fossé Hellénique, est de l'ordre de 5 à 10 cm/a vers le SW .