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. <br />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.<br />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.

Mieux appréhender la relation entre fluides, sismicité et déformation asismique est crucial en terme de risques et de ressources. Dans les zones d'injections de fluide une augmentation du taux de sismicité est observée, certains événements dépassant Mw=5. Quel est alors le rôle des fluides dans le déclenchement et le contrôle de la sismicité ? Une série d'injections de fluide à haute-pression a été réalisée dans les séries carbonatées du LSBB (Rustrel), dans la zone endommagée d'une faille inactive à 280 m de profondeur. Ces expériences in-situ ont permis l'étude des réponses sismiques et hydromécaniques (enregistrées par un large réseau de capteurs) de différentes structures géologiques à une stimulation hydraulique. Seuls certains tests ont été impactés par des séismes bien qu'une rupture ait été mesurée au point d'injection par un extensomètre. 215 séismes ont été détectés et se caractérisent par un contenu haute-fréquence (0.6 à 3 kHz) et de faible magnitude (-4.1 à -3.1). Leur localisation absolue et relative (précision de 1.5 m) a mis en évidence un manque de séismes à proximité du puits d'injection. En comparant le moment sismique cumulé et un moment équivalent de déformation, plus de 96 % de la déformation est asismique. Deux comportements sismiques distincts ont montré qu'au moins une partie de la sismicité était contrôlée par un transfert de contraintes. L'interprétation jointe des données géologiques, mécaniques, hydrogéologiques et sismiques a permis de reconstruire le mouvement des blocs de roches au point d'injection. Ainsi, la stimulation hydraulique de faille génèrerait un mouvement asismique, qui par transfert de contrainte, déclencherait la sismicité. / Better understanding how fluids pressure produce seismic or aseismic motion along faults is an important goal for seismic hazard assessment and for geological reservoir monitoring. Seismicity rate increase in fluid injection areas where some events may reach magnitude greater than 5. How fluids may induce and control seismicity? High-pressure fluid injections were performed in limestones, in the damaged zone of an inactive fault at 280m depth. These in-situ experiments allow to study the seismological and hydromechanical responses (recorded by 31 sensors) of different fracture types to a fluid perturbation. Only a few tests have generated seismicity even if ruptures are observed with a displacement sensor at the injection point. 215 earthquakes were detected and are characterized by high frequency content (0.6 to 3 KHz) and weak magnitude (-4.1 to -3.1). The relative and absolute locations (1.5m accuracy) indicate a lack of events in the vicinity of injection borehole. Comparing cumulated seismic moment with an equivalent deformation moment, more than 96% of the deformation is aseismic. Two distinct seismic behavior show that at least one part of the seismicity might be controlled by a stress transfer in the medium. Moreover, the joined interpretation geological, mechanical, hydrogeological and seismic data allow to build bloc motions at the injection point. Finally, our experiments showed that fluid injection mainly drives aseismic motion and the seismicity might be only an indirect effect related to stress transferred from the volume deformed by fluid pressurization.

Afin de mieux comprendre les interactions entre les modifications des contraintes induites par l'exploitation minière et la génération d'activité sismique, une zone profonde de la mine de Garpenberg (Suède) a été instrumentée par l’Ineris avec de sondes microsismiques et de cellules géotechniques. L’analyse spatio-temporelle des événements microsismiques enregistrés entre 2015 et 2016, ainsi que leurs paramètres à la source, ont mis en évidence deux types de réponses sismiques : une locale et courte dans le temps directement induite par les tirs de production, l’autre plus persistante et distante des excavations étant principalement contrôlée par des hétérogénéités géologiques. L’analyse des données géotechniques a montrée l’occurrence de déformations asismiques, ainsi que de phénomènes de fluage induits par l’exploitation. De plus l'activité sismique décroît proportionnellement au taux de diminution des déformations mesurées. Cette dernière observation implique que le fluage peut être un autre mécanisme menant à la sismicité, s’ajoutant au changement de contrainte immédiat induit par les tirs de production. Dans la dernière partie de cette thèse, un modèle géomécanique élasto-plastique 3D a été réalisé et ses résultats ont été comparés aux données géophysiques. Cette comparaison a montré que les modèles numériques à l'échelle de la mine peuvent être des outils puissants pour étudier la sismicité induite à grande échelle. Cependant, il y a certains aspects de la sismicité induite que le modèle ne peut expliquer entièrement. Cela est le cas pour la sismicité déclenchée à distance des excavations, alors que de meilleures corrélations sont trouvées lorsque l'on considère la sismicité à proximité des zones de production. Les résultats de cette thèse ont démontré qu'une approche combinée associent les données sismiques et géotechniques à la modélisation numérique peut améliorer considérablement notre compréhension de la réponse des massifs rocheux à l'exploitation minière. La combinaison de ces méthodologies dans une approche intégrée peut réduire considérablement leurs limites explicites qui sont évidents lorsque ces instruments sont considérés séparément. / With the aim of better understanding interactions between stress modifications induced by mining and the generation of seismic activity, a deep area of Garpenberg mine (Sweden) was instrumented by Ineris with microseismic probes and geotechnical cells. Spatiotemporal analysis of recorded seismicity between 2015 and 2016, as well as seismic source parameters, have highlighted two types of seismic rock mass responses: one local and temporally short directly induced by production blasts, the other long-lasting over time and remote from excavations being mainly controlled by geological heterogeneities. Geotechnical data analysis showed the occurrence of aseismic deformations, as well as creep phenomena induced by mining exploitation. In addition, seismic activity decays proportional to the decaying rate of measured strains. This latter observation implies that creep may be another mechanism driving seismicity, in addition to the immediate stress change induced by blasting. In the last part of this thesis, a 3D elasto-plastic geomechanical model has been realized and its results have been compared with geophysical data. This comparison showed that mine-wide numerical models can be suitable for the analysis of mining induced seismicity at large-scale. However, there are some aspects of the induced seismicity that the model cannot fully explain. This is particularly true for remote seismicity occurring at a distance from excavations, while better correlations are found when considering seismicity close to production areas. Results of this thesis demonstrated that a combined approach which associates seismic and geotechnical data with numerical modelling can significantly improve our understanding of the rock mass response to mining. The combination of these methodologies in an integrated approach can significantly reduce their straightforward limitations, which appears evident when these instruments are considered separately.

Le gisement de gaz de Lacq est exploité depuis 1957 et la sismicité induite a débuté en 1969. De magnitude modérée (M <= 4.4), cette sismicité est induite par l'extraction des fluides. Après nettoyage de la base de données disponible (1989-1997), un Pointé Automatique dont les paramètres sont optimisés par Recuit Simulé a permis la mise à jour de la base de données. Ajoutés aux données pointées entre 1974 et 1988 et grâce à un modèle de vitesse 3D, plus 1600 évènements ont été localisés, la plupart dans et au dessus du réservoir. La confrontation des localisations des hypocentres et des mécanismes aux foyers avec des informations apportées par les forages et la sismique 3D a permis de localiser les différents essaims de séismes au niveau de failles relativement importantes et orientées préférentiellement pour accommoder la perturbation de contrainte. Alors que le gisement présente une orientation WNW-ESE, la sismicité est localisée selon<br />un axe préférentiel SSW-NNE sur les failles principales d'orientation WNW-ESE. Les zones WNW et ESE du gisement sont quant à elles déformées de manière asismique, comme le montre les profils de nivellement. La confrontation de la sismicité avec l'historique et le modèle de réservoir a permis de mettre en avant (i) un essaim déclenché par une injection d'eau involontaire, (ii) une connexion importante entre le centre du réservoir et les réserves de gaz au nord-est et la déplétion associée, (ii) ainsi qu'une réinterprétation du modèle de réservoir dans certaines zones mal connues. La seconde voie d'étude a été portée sur l'analyse des formes d'ondes par décomposition atomique en chirplets, généralisation de la transformée en ondelettes et de la transformée de Gabor. Après une généralisation à 7D des paramètres de construction des chirplets pour une adaptation aux ondes sismologiques et l'utilisation d'algorithme de séparation d'ondes, il a été possible de décomposer les sismogrammes en un nombre limité d'atomes. Chacun des atomes est qualifié de 9-tuplet, caractérisé par des attributs morphologiques : (i) un temps d'arrivée, (ii) un fréquence centrale, (iii) une énergie, (iv) une phase, (v) une durée, (vi) un taux de modulation de fréquence, (vii) un type de modulation de fréquence, (viii) une asymétrie d'enveloppe et (ix) une plateau constant de l'enveloppe. Grâce à cette caractérisation poussée des ondes, il a été possible de pratiquer des classifications des sismogrammes. La classification et la représentation en dendrogramme ont montré une différenciation nette entre les séismes déclenchés par surpression de fluide des autres, révélant que la source des évènements déclenchés est plus courte et impulsive. En se basant uniquement sur les attributs morphologiques et en filtrant la composante de temps et fréquence des arrivées des ondes, nous avons montré qu'il était possible de comparer des séismes de localisations différentes.

Dans le cadre d’un programme de recherche mené par le GISOS (Groupement d’Intérêt Scientifique de Recherche sur l’Impact et la Sécurité des Ouvrages Souterrains), le site pilote de Cerville-Buissoncourt (Lorraine, France) a fait l’objet d’une importante instrumentation géophysique et géotechnique pour assurer la surveillance d’une cavité saline à 200 m de profondeur, depuis son état stationnaire jusqu’à l’effondrement des terrains du recouvrement. Les objectifs principaux de cette thèse consistaient à : 1) valider la technique de surveillance basée sur l’écoute microsismique dans un contexte salin, et 2) modéliser numériquement le comportement mécanique complexe du recouvrement, particulièrement l’initiation des microfissures et leur propagation. L’analyse de la sismicité induite enregistrée a permis de caractériser l’état initial de la cavité en terme d’activité microsismique. Deux types d’événements ont été identifiés : les événements isolés correspondant aux ruptures localisées, et les événements en rafale, d’une dizaine de secondes de durée. D’après les résultats de localisation d’hypocentres, la totalité de la sismicité enregistrée est générée au niveau de la cavité dans le gisement de sel, ou bien dans les faciès marneux qui composent le toit immédiat de la cavité actuelle. Les déclenchements en rafale seraient liés à des phénomènes de délitement puis de décrochement de blocs de marne, suivis des chutes de blocs dans la cavité remplie de saumure. Le travail de modélisation numérique a été focalisé sur la possibilité de rendre compte de l’endommagement dans les couches fragiles du recouvrement. Nous avons mis en oeuvre un modèle géomécanique hybride à l’échelle du site pilote qui intègre les différentes formations géologiques présentes dans le recouvrement, ainsi que l’initiation, la propagation et la coalescence des microfissures dans le banc raide, à l’aide des logiciels FLAC et PFC2D. La calibration du modèle discret PFC2D pour reproduire le comportement en traction du banc raide a été vérifiée numériquement à l’échelle du site pilote. Cette vérification a été basée sur la comparaison, en termes de la réponse élastique et d’apparition des ruptures dans le banc raide, entre l’approche hybride FLACPFC 2D et la modélisation purement continue avec FLAC. Le modèle hybride ainsi défini pourra être utilisé dans le cadre d’une retro-analyse une fois que les mesures in-situ, notamment les enregistrements microsismiques et les données de déformation, seront disponibles à Cerville-Buissoncourt / Within the framework of a research program carried out by the GISOS (Scientific Grouping of Research Interest on the Impact and Safety of Underground Works), the pilot site of Cerville-Buissoncourt (Lorraine, France) was the subject of a large geophysical and geotechnical instrumentation to ensure the monitoring of a salt cavity at 200 m depth, from its stationary state to the final overburden collapse. The main objectives of this work consisted on : 1) the validation of the microseismic monitoring technique in a salt mine environment, and 2) the numerical modelling of the mechanical behavior of the overburden, particularly the initiation and the propagation of microcracks. The analysis of the recorded induced seismicity allowed to characterize the initial state of the cavity in terms of microseismic activity. Two types of events were identified : isolated events corresponding to localized ruptures, and swarms of events, of tens of seconds of duration. According to hypocenter location results, the totality of the recorded seismicity is generated either in the cavity surroundings within the salt layer, or in the marly facies of the current cavity roof. Swarms would be related to delamination of clayley marls in the immediate roof, followed by rock debris falling in the brine filled cavity. The numerical modelling was focused on the possibility of accounting for the damage in the fragile layers of the overburden. We implemented a hybrid geomechanical model of the pilot site which integrates the various geological formations present in the overburden, as well as the initiation, the propagation and the coalescence of microcracks in the stiff layer, using FLAC and PFC2D softwares. The calibration of the discrete PFC2D model to reproduce the tensile behaviour of the stiff layer was numerically checked on the site scale. The validation was based on the comparison, in terms of the elastic response and the damage onset in the stiff layer, between the hybrid approach FLAC-PFC2D and the purely continuous modelling using FLAC. The hybrid model thus defined would be used for back-analysis studies once in-situ measurements, in particular microseismic recordings and deformation data, will be available at Cerville-Buissoncourt

Le développement et l'exploitation de réservoirs géothermiques profonds s'accompagnent généralement d'une sismicité induite - un effet secondaire indésirable. Cette recherche est axées sur l'utilisation de simulations numériques pour étudier la propagation des fractures hydrauliques et la réactivation de failles préexistantes lors de la stimulation hydraulique des réservoirs afin de mieux comprendre le comportement du réservoir fracturé et de réduire le risque potentiel de sismicité induite.La sismicité induite est d'abord étudiée du point de vue de l'utilisation de la loi de conservation de l'énergie afin d'expliquer le mécanisme de génération d'ondes élastiques à partir d'une rupture de roche. Ensuite, une approche approximative est proposée pour calculer les accélérations de pointe (PGA) induites par le glissement de faille. Les PGA calculés à la surface du sol servent à évaluer la perception humaine des ondes sismiques et le potentiel de dégradation des structures. / The development and the exploitation of deep geothermal reservoirs are usually accompanied with induced seismicity - an unwanted side effect. This research is focused on using numerical simulations to investigate the propagation of hydraulic fractures and the reactivation of pre-existing faults during the hydraulic stimulation of the reservoirs in an effort to better understand the fractured reservoir behavior and to reduce the potential risk of induced seismicity.The induced seismicity is studied first from the standpoint of using the law of energy conservation in order to explain the mechanism of generating elastic waves from rock failure. Then an approximate approach is proposed to calculate the peak ground accelerations (PGAs) that are induced by the fault slip. The computed PGAs on ground surface are used to assess the human perception of the seismic waves and the damage potential to structures.

Sur la base des faibles variations de contraintes poroélastiques qui sont à l'origine de la séismicité déclenchée localement par la baisse de pression du gisement de Lacq, on tente dans le chapitre 3 d'évaluer les échelles spatiales des connections fluides dans la croûte supérieure à l'aide de la migration des fronts de pression induits. Les séismes associés à la mise en eau d'une retenue artificielle (M 4-5, Monteynard, Vercors, France) permettent par exemple de cartographier des failles sismiques potentielles au sud de l'agglomération grenobloise, les séismes induits jouant alors le rôle de jauges de contrainte. Dans la région de Lacq, à l'aide de modélisations analytiques des transferts de contraintes visco-élastiques et poro-élastiques, on montre qu'on ne peut rejeter a-priori des interactions entre l'extraction du champ de Lacq, les séismes majeurs de la faille Nord-Pyrénéenne distants d'une trentaine de kilomètres, et les séismes locaux à l'aplomb du champ d'hydrocarbure. Le chapitre 4 constitue une synthèse des mécanismes d'instabilités sismiques déclenchées par les exploitations d'hydrocarbures. On isole trois classes de mécanismes de déclenchements des séismes correspondant à des types d'exploitations et à des réponses sismiques (taille-temps-espaces) bien différenciés: augmentation de pression due à l'injection de fluide; baisse de pression due à extraction de fluide, déficit de masse lors d'extraction massive. Les mécanismes de ruptures sont en accord avec les lois de contraintes effectives, et s'expliquent soit par des transferts poroélastiques soit par des compensations isostasiques. Si la sismicité induite par les extractions de fluides reste marginale par rapport au nombre de gisements exploités, la surveillance sismique fine effectuée sur les gisements montre que de nombreux petits séismes (M < 3) sont présents et permettent un suivi in-situ du comportement des réservoirs au cours de l'exploitation. Notre approche de la séismicité à l'aide des séismes induits par des activités humaines montre les limites des lois de la mécanique classique qui, si elles permettent de comprendre et d'évaluer les seuils critiques qui déclenchent les instabilités sismiques, ne peuvent expliquer ni la durée des phénomènes d'instabilités qui sont entretenus durant de nombreuses années ni la taille des instabilités sismiques (Mmax 7). Sur la base des cas mondiaux de séismicité induite on propose la notion de Systèmes Critiques Auto Organisés Induits où l'on peut observer la genèse, la pérennité et la disparition de systèmes critiques dans un contexte, proposé par ailleurs, d'état critique auto-organisé pour l'ensemble de la croûte supérieure.

Depuis les séismes dévastateurs de 1999 d’Izmit et de Duzce, la partie immergée de la Faille Nord Anatolienne (FNA)en Mer de Marmara fait l’objet d’une intense surveillance. Malgré cela, la micro-sismicité demeure mal connue. Par ailleurs, alors que la connexion avec le système pétrolier du Bassin de Thrace est établie, le rôle du gaz sur la sismicité n’a pas été identifié.Dans ce travail, nous avons analysé des données d’OBS (Ocean Bottom Seismometers) acquises dans la partie ouest de la Mer de Marmara (en avril-juillet 2011 et septembre-novembre 2014), à partir de méthodes non-linéaires –NonLinLocet d’un modèle 3D de vitesses. Une grande partie de la sismicité se produit à des profondeurs inférieures à 6 km environ : le long de failles secondaires, héritées de l’histoire complexe de la FNA ; ou dans des couches de sédiments superficiels (< 1 km) riches en gaz. Cette sismicité superficielle semble être associée à des processus liés au gaz, déclenchés par les séismes profonds de magnitude M1 > 4.5 qui se produisent régulièrement le long de la MMF.Par ailleurs, 2 familles de signaux de courte durée (<1s), dits ≪ SDE ≫ (pour Short Duration Event) apparaissent sur les enregistrements : 1) les SDE se produisant à raison de quelques dizaines de SDE/jour, en réponse à des causes locales (i.e. bioturbation, activité biologique, micro-bullage de fond de mer, mouvements à l’interface eau/sédiment), etc ; 2) lesSDE se produisant par ≪paquets≫, dont certains sont enregistrés sur les 4 composantes (y compris l’hydrophone) et apparaissent de manière périodique, toutes les 1.8 s environ, en réponse à diverses causes qui restent à déterminer (parmi lesquelles : les mammifères marins ; l’activité humaine ; la sismicité ; le dégazage ; les ≪trémors≫ sismiques ; etc). / Since the devastating earthquakes of 1999, east of Istanbul, the submerged section of the North Anatolian Fault (NAF), in the Sea of Marmara (SoM) has been intensively monitored, mainly using land stations. Still, the micro-seismicity remains poorly understood. In addition, although the connection of the SoM with the hydrocarbon gas system from the Thrace Basin is now well established, along with the presence of widespread gas within the sedimentary layers, the role of gas on seismicity is still not recognized.Here, we have analyzed Ocean Bottom Seismometer (OBS) data from two deployments (April-July 2011 and September-November 2014) in the western SoM. Based on a high-resolution, 3D-velocity model, and on non-linear methods (NonLinLoc), our location results show that a large part of the micro-seismicity occurs at shallow depths (< 6 a 8 km): along secondary faults, inherited from the complex history of the North-Anatolian shear zone; or within the uppermost (< 1 km), gas-rich, sediment layers. Part of this ultra-shallow seismicity is likely triggered by the deep earthquakes of intermediate magnitude (Ml > 4.5) that frequently occur along the western segments of the MMF.In addition, OBSs also record at least two families of short duration (<1 sec) events (SDEs): 1) “background SDEs” occurring on a permanent, at a rate of a few tens of SDEs/day, resulting from many possible, local causes, e. g.: degassing from the seafloor, biological activity near the seabed, bioturbation, etc; 2) “swarmed SDEs”, among which some are recorded also on the hydrophone, and characterized by a periodicity of ~ 1.8 seconds. The causes of these SDEs still remain to be determined (among which: anthropogenic causes, marine mammals, gas emissions, regional seismicity, tremors from the MMF, etc).

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)

Ce travail détermine l'influence du milieu (en terme de géologie et de mécanique) sur les paramètres à )a source de séismes induits par l'exploitation d'un gisement d'hydrocarbures. Tout d'abord, un rappel des études antérieures est proposé: la relocalisation par inversion des séismes et la bonne connaissance du site (70 forages) ont permis d'attribuer à chaque événement une appartenance géologique. D'autre part, des essais de mécanique des roches ont été effectués sur les niveaux les plus sismiques. Une chaîne de numérisation a été mise au point afin d'obtenir des signaux digitaux. Ces signaux sont ensuite traités avec une méthode classique d'analyse spectrale. On remonte alors aux paramètres des sources à l'aide d'un modèle dynamique de faille plane circulaire. Ces paramètres sont : le moment sismique, le rayon de la source, la chute de contrainte, la dislocation et l'énergie sismique libérée. Tous les facteurs intervenant dans le calcul sont évalués afin de réduire au maximum les incertitudes sur les mesures. L'interprétation des résultats est divisée en 2 parties. La première, dite classique, montre des anomalies dans les distributions temporelle et en profondeur des différents paramètres; mais ces répartitions ne peuvent être expliquées. C'est pourquoi la deuxième partie intègre la connaissance du site. Une analyse par niveau des paramètres est détaillée. On observe une grande dépendance des rayons avec l'épaisseur des bancs et/ou la préfracturation du milieu. Les chutes de contrainte sont elles influencées par le comportement mécanique du milieu. Plus le milieu sera compétent, plus les chutes de contrainte maximales seront élevées. De plus, le réservoir, compacté par drainage, présente des chutes de contrainte minimales plus élevées que celles de toutes les autres couches. Ces résultats démontrent l'importance du milieu sur les lois d'échelle en particulier. L'analyse de l'énergie sismique libérée en fonction de la production du gisement semble montrer une corrélation qui augmente avec le temps, mais cette tendance demande à être confirmée par une période d'étude plus large.