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.

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.