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.

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.

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.