Tomographie de pente fondée sur l'état adjoint : un outil d'estimation de modèle de vitesse pour l'imagerie sismique

Adjoint slope tomography : a velocity macro-model building tool for seismic imaging

• Auteur / Autrice : Borhan Tavakolifaradonbeh
• Direction : Alessandra Ribodetti, Jean Virieux
• Type : Thèse de doctorat
• Discipline(s) : Sciences de la planète et de l'univers
• Date : Soutenance le 16/11/2017
• Etablissement(s) : Université Côte d'Azur (ComUE)
• Ecole(s) doctorale(s) : École doctorale Sciences fondamentales et appliquées (Nice ; 2000-....)
• Partenaire(s) de recherche : établissement de préparation : Université de Nice (1965-2019)
• Partenaire(s) de recherche : Laboratoire : Laboratoire Géoazur (Sophia Antipolis, Alpes-Maritimes) - Géoazur
• Jury : Président / Présidente : René-Edouard Plessix
• Jury : Examinateurs / Examinatrices : Alessandra Ribodetti, Jean Virieux, René-Edouard Plessix, Gilles Lambaré, Hervé Chauris, Mark Noble, Reda M'Hammed Baina
• Jury : Rapporteurs / Rapporteuses : Gilles Lambaré, Hervé Chauris

Mots clés

Mots clés contrôlés

Équation iconale

Inversion (géophysique)

Tomographie

Imagerie sismique

Mots clés libres

Imagerie sismique

Tomographie

Inversion

Eikonal

Méthode d'état adjoint

Seismic imaging

Tomographie

Inversion

Eikonal

Adjoint state method

Résumé

La construction du macro-modèle de vitesse est une étape cruciale de la chaîne d’imagerie sismique afin de produire un modèle adéquat pour la migration ou l'inversion de la formes d’onde complètes (Full Waveform Inversion). Parmi les approches possibles pour construire ce macro-modèle, la tomographie des pentes est fondée sur le pointé d’évènements localement cohérents caractérisés par leur temps de trajet et leurs pentes dans les collections de sismogrammes à source et réflecteur communs. Chaque évènement dans les observables est associé à un petit segment de réflecteur dans le sous-sol caractérisé par sa position et son pendage. Cette thèse propose une reformulation de la tomographie des pentes pour palier aux deux limitations sus-mentionnées. Les temps de trajet sont calculés à l’aide d’un solveur eikonal et la méthode de l’état adjoint est utilisée pour calculer efficacement le gradient de la fonction coût dans le contexte de méthodes d’optimisation locale. Le solveur eikonal est fondé sur une méthode aux différences finies pour la discrétisation des opérateurs différentiels. La méthode est implémentée pour des milieux 2D transverses isotropes avec un axe de symétrie dont l’orientation varie spatialement (milieux TTI). La méthode est évaluée avec différents examples synthétiques comme les cas du modèle isotrope complexe Marmousi et du modèle 2D TTI BP-salt. Dans le cas de milieux TTI, les couplages pouvant exister entre les paramètres de différente nature sont analysés avec un cas synthétique canonique. Cette thèse se conclut par une application à des données de sismique réflexion multitrace large bande (fournie par CGG) afin de reconstruire le modèle de vitesse vertical.

Velocity macro-model building is a crucial step in seismic imaging workflows as it provides the necessary background model for migration or full waveform inversion. Slope tomography, as a reliable alternative for conventional tomography, provides a tool to achieve this purpose where one picks the local coherent events rather than continuous events. These approaches are based on the slopes and traveltimes of the local coherent events which are tied to a reflecting/diffracting point (scatterer) in the subsurface. In this thesis, I introduce an anisotropic slope tomographic approach which aim at macro-model building for subsurface properties in 2D tilted transversely isotropic (TTI) media. In this method, I reformulate the stereotomography, as an slope tomographic tool, such that I replace the ray-based forward engine with a TTI eikonal solver and take advantage of the adjoint state method to calculate the gradients. In result, I can efficiently calculate the traveltimes for complex media and long offset acquisition on a regular grid of the subsurface and formulate a matrix-free framework for the inversion. Different synthetic examples including the isotropic Marmousi and 2D TTI BP-salt model are considered to assess the potential of the method in subsurface parameter estimations. Also, through a simple example, the footprint of parameter cross-talk is investigated for the Thomsen parametrization. As a real data application, the proposed method is applied on a 2D marine BroadSeis data set (provided by CGG) to retrieve the vertical velocity of the subsurface.