La donnée sismique depuis sa conception (modélisation d’acquisitions sismiques), dans sa phase de traitement (prétraitement et migration) et jusqu’à son exploitation pour en extraire les informations géologiques pertinentes nécessaires à l’identification et l’exploitation optimale des réservoirs d’hydrocarbures (interprétation), génère un volume important de calculs. Nous montrons dans ce travail de thèse qu’à chacune de ces étapes l’utilisation de technologies accélératrices de type GPGPU permet de réduire radicalement les temps de calcul tout en restant dans une enveloppe de consommation électrique raisonnable. Nous présentons et analysons les éléments sous-jacents à ces performances. L’importance de l’utilisation de motifs d’accès mémoire adéquats est particulièrement mise en exergue étant donné que l’accès à la mémoire représente le principal goulot d’étranglement pour les algorithmes abordés. Nous reportons des facteurs d’accélération de l’ordre de 40 pour la modélisation sismique par résolution de l’équation d’onde par différences finies (brique de base pour la modélisation et l’imagerie sismique) et entre 8 et 113 pour le calcul d’attributs sismiques. Nous démontrons que l’utilisation d’accélérateurs matériels élargit considérablement le champ du possible, aussi bien en imagerie sismique (modélisation de nouveaux types d’acquisitions à grande échelle) qu’en interprétation (calcul d’attributs complexes sur station de travail, paramétrage interactif des calculs, etc.). / During the seismic imaging workflow, from seismic modeling to interpretation, processingseismic data requires a massive amount of computation. We show in this work that, at eachstage of this workflow, hardware accelerators such as GPUs may help reducing the time requiredto process seismic data while staying at reasonable energy consumption levels.In this work, the key programming considerations needed to achieve good performance are describedand discussed. The importance of adapted in-memory data access patterns is particularlyemphasised since data access is the main bottleneck for the considered algorithms. When usingGPUs, speedup ratios of 40× are achieved for FDTD seismic modeling, and 8× up to 113× forseismic attribute computation compared to CPUs.

Le but de ce travail est une étude pluridisciplinaire de modélisation sismique 2D et d'inversion d'une coupe géologique de la bordure méridionale du Vercors. L'originalité de ce travail réside dans l'emploi de techniques qui mènent à la caractérisation sismique de réservoir à partir d'une coupe géologique. Le modèle géologique est construit à partir de l'étude sur le terrain des ameurements de la bordure de la plateforme urgonienne (Barrémien - Aptien inférieur). Les corrélations cartographiques et séquentielles de coupes lithologiques de détail permettent d'établir une coupe géologique de 16 kilomètres de long pour 800 mètres d'épaisseur environ. Les calcaires hémipélagiques du Barrémien inférieur et les calcaires bioclastiques du prisme de bas niveau et de l'intervalle transgressif de la séquence BA3 montrent une diminution d'épaisseur en direction du Nord. Cette géométrie résulte d'une érosion des couches du Barrémien inférieur et~Ùne disposition en biseau d'aggradation ("onlaps") des paraséquences du prisme de bas niveau et de l'intervalle transgressif de la séquence BA3 sur la limite SbB3. Ces calcaires bioclastiques ont pu être datés de la base du Barrémien supérieur. Les Calcaires urgoniens montrent une position rétrogradante par rapport aux calcaires bioclastiques de Glandasse, indiquant qu'ils appartiennent à une formation transgressive dont l'âge est Barrémien supérieur. L'échantillonnage des différents Iithofaciès identifiés permet de réaliser la mesure en laboratoire des propriétés pétrophysiques et acoustiques (Vp, Vs et densité) de ces roches. Des relations linéaires entre les vitesses et la porosité ou l'argilosité sont déterminées. Ces. relations sont utilisées pour construire les modèles d'impédance, à partir des valeurs de porosité d'un réservoir réel. La construction des modèles d'impédance est réalisée en accord avec les hypothèses de la stratigraphie séquentielle, toutes les limites des paraséquences sont des réflect,eurs potentiels. La modélisation sismique a été réalisée par deux approches différentes, tracé de rayons à déport nul et différences fmies "multioffset". Les signaux sources employés ont un contenu fréquentiel comparable à ceux généralement extraits des profils sismiques de prospection pétrolière. L'interprétation des profils sismiques est comparée à la coupe géologique. Cette comparaison montre qu'il est difficile de retrouver le découpage en séquences de dépôt "il partir de l'interprétation des relations géométriques et des terminaisons des biseaux des événements sismiques. La dernière partie de ce travail est consacrée à la mise en oeuye de deux techniques de caractérisation sismique de réservoirs sur les profils simulés. Dans un premier temps, une inversion stratigraphique a été réalisée sur les données modélisées par tracé de rayons. Les sections d'impédance acoustiques obtenues sont caractérisées par des valeurs d'impédance proches du modèle d'origine. Dans un second temps, les attributs sismiques après sommation sont employés pour réaliser une interprétation lithologique des données sismiques. L'emploi de l'impédance sismique comme attribut permet une première discrimination lithologique. Un second attribut (attribut instantané) permet de discriminer les classes lithologiques présentant la même impédance sismique. Les profils lithologiques obtenus après prédiction sont comparables au modèle géologique.