Ces dernières années ont vu le développement d'une nouvelle méthode d'imagerie des structures géologiques basée sur l'utilisation du bruit sismique continu. Dans ce travail nous avons utilisé cette approche dans deux problématiques différentes. La première consiste à réaliser le suivi temporel des vitesses des ondes sismiques dans la croûte lors de séismes lents qui ont eu lieu dans la région de Guerrero au Mexique. Les séismes lents sont des glissements asismiques et transitoires qui ont été découverts récemment dans la lacune sismique de Guerrero. Ils sont considérés comme une part importante de la relaxation des déformations dans le cycle sismique. Les séismes lents affectent le processus de chargement et déchargement de l'interface, il est donc important de comprendre le comportement mécanique de la subduction dans cette région pour mieux évaluer le risque sismique. Dans notre étude, nous avons mesuré les perturbations de vitesse des ondes associées à deux séismes lents en 2006 et 2009-2010 à partir des enregistrements continus du bruit sismique. Pour chacun des deux séismes lents on observe une chute de vitesse : elle s'élève à 0.2% pour celui de 2006 et à 0.8 % pour celui de 2009-2010. Au cours de ces séismes lents, les ondes de longues périodes (>10 s) sont perturbées. A courtes périodes, aucune variation de vitesse n'est observée ce qui suggère qu'un endommagement des couches superficielles de la croûte ne peut pas être à l'origine du changement de vitesse. Par ailleurs, la perturbation de vitesse est reliée au taux de dé- formation plutôt qu'à la déformation elle-même. Cette observation suggère que pendant de forts séismes lents, la croûte chevauchante présente un comportement mécanique non linéaire. Nous pouvons donc utiliser les variations de vitesse comme des marqueurs du taux de déformation du milieu. Enfin, une corrélation entre les trémors non volcaniques et les variations de vitesse suggère qu'une part importante de la déformation résultant des séismes lents est accommodée par la croute chevauchante. La deuxième problématique abordée dans cette thèse est l'imagerie de structures à fort contraste de vitesse et dans lesquelles la propagation des ondes de surface est complexe. Imager et comprendre la propagation des ondes dans la vallée de Mexico est crucial pour l'estimation du risque sismique à la capitale. Nous mesurons la dispersion des ondes de Rayleigh reconstruites à partir d'intercorrélations de bruit de fond sismique. Pour identifier les modes nous utilisons une mesure du rapport spectral des composantes horizontales sur la composante vertical (H/V) sur la coda des séismes que l'on compare avec le rapport H/V théorique. Grâce à cette identification des modes, nous pouvons retrouver le modèle de vitesse de la structure.

Ces dernières années ont vu le développement d'une nouvelle méthode d'imagerie des structures géologiques basée sur l'utilisation du bruit sismique continu. Dans ce travail nous avons utilisé cette approche dans deux problématiques différentes. La première consiste à réaliser le suivi temporel des vitesses des ondes sismiques dans la croûte lors de séismes lents qui ont eu lieu dans la région de Guerrero au Mexique. Les séismes lents sont des glissements asismiques et transitoires qui ont été découverts récemment dans la lacune sismique de Guerrero. Ils sont considérés comme une part importante de la relaxation des déformations dans le cycle sismique. Les séismes lents affectent le processus de chargement et déchargement de l'interface, il est donc important de comprendre le comportement mécanique de la subduction dans cette région pour mieux évaluer le risque sismique. Dans notre étude, nous avons mesuré les perturbations de vitesse des ondes associées à deux séismes lents en 2006 et 2009-2010 à partir des enregistrements continus du bruit sismique. Pour chacun des deux séismes lents on observe une chute de vitesse : elle s'élève à 0.2% pour celui de 2006 et à 0.8 % pour celui de 2009-2010. Au cours de ces séismes lents, les ondes de longues périodes (>10 s) sont perturbées. A courtes périodes, aucune variation de vitesse n'est observée ce qui suggère qu'un endommagement des couches superficielles de la croûte ne peut pas être à l'origine du changement de vitesse. Par ailleurs, la perturbation de vitesse est reliée au taux de dé- formation plutôt qu'à la déformation elle-même. Cette observation suggère que pendant de forts séismes lents, la croûte chevauchante présente un comportement mécanique non linéaire. Nous pouvons donc utiliser les variations de vitesse comme des marqueurs du taux de déformation du milieu. Enfin, une corrélation entre les trémors non volcaniques et les variations de vitesse suggère qu'une part importante de la déformation résultant des séismes lents est accommodée par la croute chevauchante. La deuxième problématique abordée dans cette thèse est l'imagerie de structures à fort contraste de vitesse et dans lesquelles la propagation des ondes de surface est complexe. Imager et comprendre la propagation des ondes dans la vallée de Mexico est crucial pour l'estimation du risque sismique à la capitale. Nous mesurons la dispersion des ondes de Rayleigh reconstruites à partir d'intercorrélations de bruit de fond sismique. Pour identifier les modes nous utilisons une mesure du rapport spectral des composantes horizontales sur la composante vertical (H/V) sur la coda des séismes que l'on compare avec le rapport H/V théorique. Grâce à cette identification des modes, nous pouvons retrouver le modèle de vitesse de la structure. / Recent years have seen the development of a new method for imaging geological structures based on continuous seismic noise. In this work we used this approach in two different problems. The first is to monitor seismic waves velocity in the crust during slow slip events that occurred in the region of Guerrero in Mexico. These slow slip events are aseismic transients that were observed recently in the seismic gap of Guerrero. They are considered an important part of the strain relaxation in the seismic cycle. Since slow slip events affect the process of loading and unloading of the interface, it is important to understand the mechanical behavior of the subduction in this region to better assess the seismic risk. In our study, we measured wave velocity perturbations associated with two slow slip events in 2006 and 2009-2010 from continuous recordings of seismic noise. For both events we observed a drop in wave speed : it reaches 0.2 % in 2006 and 0.8 % in 2009-2010. During these slow slip events waves of long periods (> 10 s) are disturbed. At short period, no velocity variation is observed suggesting that damage of the superficial layers of the crust cannot produce such velocity perturbation. Moreover, the wave speed change is related to the strain rate rather than the deformation itself. This observation suggests that during strong slow slip events the overriding crust presents a nonlinear mechanical behavior. We can therefore use the velocity variations as a proxy of the strain rate of the medium. Finally, a correlation between non-volcanic tremor and changes in waves speed suggests that part of the deformation resulting from the slow slip events is accommodated by the overriding crust. The second issue addressed in this thesis is imaging geological structures with high velocity contrast in which the propagation of surface waves is complex. Charac- terizing and understanding wave propagation in the Valley of Mexico is crucial for the estimation of seismic risk in Mexico City. We measure the dispersion of Rayleigh waves reconstructed from cross-correlations of seismic noise. To identify the modes of Rayleigh waves we use a measure of the spectral ratio of the horizontal components to the vertical component (H / V) in the coda of earthquakes which are compared with the theoretical H / V. With this identification method, we can find the velocity model of the structure.

La corrélation de bruit sismique est aujourd'hui une technique largement utilise en sismologie. L'aspect attrayante de cette technique est lie a la possibilité d'obtenir une estimation de la fonction Green partout sur la surface de la Terre, aussi en absence de source sismique comme le tremblement de terre. Les ondes des surface ils sont facilement obtenu par corrélation de bruit, et elles sont très utiliser pour l'imagerie de la croute et du manteau superieur. Dans la possible observation des ondes de volume par corrélation de bruit, il peut avoir le potentielle pour résoudre aussi les structures de la Terre profonde.Dans la première partie de la thèse, on utilise les ondes de surface estime par corrélation de bruit sismique au fin de produire en modèle tomographique de la croute en Finlande.Dans la deuxième partie, on se focalise sur la nouvelle possibilité d'obtenir des ondes de volume par corrélation de bruit. Suite a la démonstration aue des ondes des volume son bien observe par corrélation de bruit a toutes les échelles, on utilise ces ondes pour l'imagerie de la Terre profonde. Ces nouveaux et prometteur résultats ils permettent d’incrémenter la connaissance de la structure profonde de la Terre, avec un incroyable augmentation de la résolution de la méthode d'imagerie. / Ambient seismic noise correlation technique is today widely applied in seismology. The attractive aspect of this method relies in the possibility of obtaining an estimated Green's function everywhere in the world, also in absence of explosions or earthquakes. As surface waves dominate the estimated Green's fucntion, they can be used for high-resolution imaging of the shallow Earth. Similar observation of body waves would provide the required resolution to solve the deeper Earth structures.In the first part of our work we focus on resolving the crustal structure of northern Finland. Using surface waves reconstructed from noise correlation, we reconstructed a threedimensional S waves model that gives acces to the ancient structures of the crust.In the second part we analyse how body waves can emerge from seismic noise correlation. We show that body waves travelling from regional to teleseismic distances are well reconstructed from noise correlation. As these waves contain information about the Earth structure they represent a new and original dataset to improve the knowledge of our planet.

L'inversion des formes d'ondes (FWI) est une méthode d'imagerie sismique quantitative multiparamètre actuellement en plein essor. La prise en compte de l'ensemble des phénomènes de propagation enregistrés, en particulier des ondes de surface, en font une technique prometteuse pour des applications de la subsurface sur des problématiques géotechniques ou environnementales. Cependant, l'adaptation au contexte de la subsurface d'une méthode sismique validée numériquement ou sur des applications d'imagerie profonde est une tâche difficile. Nous proposons ici une approche par modélisation physique à échelle réduite à l'aide d'un laboratoire de Mesures laser-Ultrasonore Sans Contact (MUSC) mis en place dans le cadre de cette thèse. Les différents éléments du laboratoire ont été caractérisés pour la modélisation sismique à échelle réduite et une attention particulière a été portée à la simulation de la source. Les données acquises ont été validées par confrontation avec des données obtenues par modélisation numérique à l'aide des codes fournis par Géoazur dans le cadre du projet ANR SEISCOPE. Dans un deuxième temps, le potentiel de la FWI pour des applications de la subsurface a été évalué par l'approche conjointe de la modélisation numérique et de la modélisation expérimentale. Nous avons ainsi pu montrer par cette double approche que la FWI permet d'obtenir une information sur la morphologie de cavités souterraines grâce à la prise en compte des ondes de surface. Par ailleurs, sur les différents modèles multicouches étudiés, l'information sur les zones les plus superficielles fournie par les ondes de surface dispersives a pu être exploitée par la FWI.

Les expériences de sismique active sont couramment utilisées pour estimer la valeur d'un modèle de vitesse de propagation desondes P dans le sous-sol. Les méthodes dites d'« analyse de vitesse par migration » ont pour but la détermination d'un macro-modèle de vitesse, lisse, et responsable de la cinématique de propagation des ondes. Dans une première étape de « migration », une image de réflectivité est obtenue à partir des données enregistrées en utilisant une première estimation du macro-modèle. Cette image dépend d’un paramètre additionnel permettant dans un second temps d’estimer la qualité du macro-modèle puis de l'améliorer. Les images de réflectivité obtenues par les techniques de migration classiques sont cependant contaminées par des artefacts, altérant la qualité de la remise à jour du macro-modèle. En particulier, elles ne prennent pas en compte les réflexions multiples, habituellement retirées des données avant traitement. Cette étape reste cependant délicate et on se prive alors de l'information supplémentaire contenue dans les multiples.Nous proposons dans cette étude une stratégie d’optimisation imbriquée en itérant l'étape de migration avant de remettre à jour le macro-modèle. La migration itérative produit des images de réflectivité satisfaisantes pour l'analyse de vitesse et s’étend naturellement aux réflexions multiples. Un désavantage de la méthode est son coût de calcul. Un pseudo-inverse de l'opérateur de modélisation est alors utilisé comme préconditionneur pour limiter le nombre d’itérations dans la boucle interne. Une autre difficulté est l'instabilité de la remise à jour du modèle de vitesse calculée pour des modèles de réflectivité successifs proches les uns des autres. Une nouvelle approche plus robustesse est proposée, valide aussi dans le cas de multiples. Son efficacité est testée sur des jeux de données synthétiques 2D. / Active seismic experiments are commonly used to recover a model of the P-wave propagation velocity in the subsurface. “Migration Velocity Analysis” techniques aim at deriving a smooth background velocity model controlling the kinematics of wave propagation. First, a reflectivity image is obtained by “migration” of observed data using a first estimate of the background velocity. This image depends on an additional “subsurface-offset” parameter allowing to assess the quality of the background velocity model with a focusing criterion and to correct it. However classical migration techniques do not provide a sufficiently accurate reflectivity image, leading to inconsistent velocity updates. In particular they do not take into account multiple reflections, usually regarded as noise and removed from the data before processing. Multiple removal is however a difficult step, and additional information contained in multiples is discarded.In this thesis, we propose to determine the reflectivity model by iterative migration before subsequent velocity analysis, leading to a nested optimisation procedure. Iterative migration yields accurate reflectivity image and extends naturally to the case of multiples. One of its disadvantages is the associated increased computational cost. To limit the number of iterations in the innerloop, a preconditioner based on a pseudo-inverse of the modelling operator is introduced. Another difficulty is the instability of the velocity update obtained with very close successive reflectivity models. We propose a modified approach, valid in the presence of multiples, and discussed through applications on 2D synthetic data sets.

L'objet de cette thèse est l'étude des possibles mécanismes élastiques expliquant l'amplification du bruit ambiant au droit de certains réservoirs multi-phasiques. Trois jeux de données sont traités. La signature spectrale observée d'un réservoir de vapeur géothermique est différente de celle d'un stockage de gaz. Dans une approche empirique, un algorithme de classification permet d'extraire et de cartographier les anomalies que l'on présume liées au réservoir. Un travail de modélisation est effectué pour tenter d'expliquer les anomalies mesurées. Dans les données réelles, une forte présence de modes supérieurs d'ondes de Rayleigh est détectée. On modélise numériquement en 2D la propagation de ces modes à travers un réservoir placé au sein d'une structure géologique réaliste. La réponse simulée du réservoir se révèle trop faible par rapport aux observations de terrain. Néanmoins, on parvient à inverser les faibles perturbations d'amplitude synthétiques pour la position du réservoir, dans des modèles de référence simples. Cette méthode pourrait être utilisable pour l'imagerie à partir de faibles variations d'amplitudes dans le cadre du monitoring. Pour ce qui est de fortes anomalies observées sur le terrain , il est à noter que les effets visco-élastiques, les effets 3D, et les effets liés à un éventuel champ incident diffus n'ont pas été pris en compte dans la modélisation. Ainsi ce travail n'exclut pas la possibilité de telles anomalies liées à la présence d'un réservoir. / This PhD work investigates the possible elastic mechanisms behind the ambient noise amplification above multi-phase fluid reservoirs. Three datasets are analysed above different reservoirs. The observed spectral signature is different in the gas storage and geothermal contexts. A non-supervised algorithm for amplitude spectrum classification is developed, allowing to extract and map the relevant attributes of a multi-phase fluid presence. As a first modelling step, a wavefield characterisation methodology is applied to determine the composition of the ambient noise. It reveals the presence of strong Rayleigh overtones. Numerical 2D elastic modelling is used to simulate the propagation of overtones across a reservoir within a realistic geological structure. The modelled reservoir response is too small compared to the real data. However, the small amplitude perturbations arising in the numerical simulations are successfully inverted for the position of the reservoir, in simple background models. The developed method could in theory be used for imaging small time-lapse amplitude variations (monitoring), despite the obstacles remaining to be overcome before a real-data application. Neither visco-elastic nor 3D effects are adressed. Thus this work does not exclude the possibility of strong reservoir-specific spectral anomalies.

Les modèles géologiques sont couramment utilisés pour estimer les ressources souterraines, pour faire des simulations numériques, et pour évaluer les risques naturels ; il est donc important que les modèles géologiques représentent la géométrie des objets géologiques de façon précise. La première étape pour construire un modèle géologique consiste souvent à interpréter des surfaces structurales, telles que les failles et horizons, à partir d'une image sismique ; les objets géologiques identifiés sont ensuite utilisés pour construire le modèle géologique par des méthodes d'interpolation. Les modèles géologiques construits de cette façon héritent donc les incertitudes d'interprétation car une image sismique peut souvent supporter plusieurs interprétations structurales. Dans ce manuscrit, j'étudie le problème de réduire les incertitudes d'interprétation à l'aide des données sismiques. Particulièrement, j'étudie le problème de déterminer, à l'aide des données sismiques, quels modèles sont plus probables que d'autres dans un ensemble des modèles géologiques cohérents. Ce problème sera connu par la suite comme "le problème d'évaluation des modèles géologiques par données sismiques". J'introduis et formalise ce problème. Je propose de le résoudre par génération des données sismiques synthétiques pour chaque interprétation structurale dans un premier temps, ensuite d'utiliser ces données synthétiques pour calculer la fonction-objectif pour chaque interprétation ; cela permet de classer les différentes interprétations structurales. La difficulté majeure d'évaluer les modèles structuraux à l'aide des données sismiques consiste à proposer des fonctions-objectifs adéquates. Je propose un ensemble de conditions qui doivent être satisfaites par la fonction-objectif pour une évaluation réussie des modèles structuraux à l'aide des données sismiques. Ces conditions imposées à la fonction-objectif peuvent, en principe, être satisfaites en utilisant les données sismiques de surface (« surface seismic data »). Cependant, en pratique il reste tout de même difficile de proposer et de calculer des fonctions-objectifs qui satisfassent ces conditions. Je termine le manuscrit en illustrant les difficultés rencontrées en pratique lorsque nous cherchons à évaluer les interprétations structurales à l'aide des données sismiques de surface. Je propose une fonction-objectif générale faite de deux composants principaux : (1) un opérateur de résidus qui calcule les résidus des données, et (2) un opérateur de projection qui projette les résidus de données depuis l'espace de données vers l'espace physique (le sous-sol). Cette fonction-objectif est donc localisée dans l'espace car elle génère des valeurs en fonction de l'espace. Cependant, je ne suis toujours pas en mesure de proposer une implémentation pratique de cette fonction-objectif qui satisfasse les conditions imposées pour une évaluation réussie des interprétations structurales ; cela reste un sujet de recherche. / Subsurface structural models are routinely used for resource estimation, numerical simulations, and risk management; it is therefore important that subsurface models represent the geometry of geological objects accurately. The first step in building a subsurface model is usually to interpret structural features, such as faults and horizons, from a seismic image; the identified structural features are then used to build a subsurface model using interpolation methods. Subsurface models built this way therefore inherit interpretation uncertainties since a single seismic image often supports multiple structural interpretations. In this manuscript, I study the problem of reducing interpretation uncertainties using seismic data. In particular, I study the problem of using seismic data to determine which structural models are more likely than others in an ensemble of geologically plausible structural models. I refer to this problem as "appraising structural models using seismic data". I introduce and formalize the problem of appraising structural interpretations using seismic data. I propose to solve the problem by generating synthetic data for each structural interpretation and then to compute misfit values for each interpretation; this allows us to rank the different structural interpretations. The main challenge of appraising structural models using seismic data is to propose appropriate data misfit functions. I derive a set of conditions that have to be satisfied by the data misfit function for a successful appraisal of structural models. I argue that since it is not possible to satisfy these conditions using vertical seismic profile (VSP) data, it is not possible to appraise structural interpretations using VSP data in the most general case. The conditions imposed on the data misfit function can in principle be satisfied for surface seismic data. In practice, however, it remains a challenge to propose and compute data misfit functions that satisfy those conditions. I conclude the manuscript by highlighting practical issues of appraising structural interpretations using surface seismic data. I propose a general data misfit function that is made of two main components: (1) a residual operator that computes data residuals, and (2) a projection operator that projects the data residuals from the data-space into the image-domain. This misfit function is therefore localized in space, as it outputs data misfit values in the image-domain. However, I am still unable to propose a practical implementation of this misfit function that satisfies the conditions imposed for a successful appraisal of structural interpretations; this is a subject for further research.

La péninsule de la Gaspésie (Québec, Canada) présente un potentiel pétrolier dans les roches datées de l'Ordovicien tardif au Dévonien inférieur. Dans l'est de la ceinture de Gaspé, des puits producteurs de gaz prouvent la présence de réservoirs pétroliers. Des nouvelles études structurales et des études de bassin dans la ceinture de Gaspé permettent de mieux évaluer le potentiel pétrolier de cette région. Nous présentons ici une nouvelle méthodologie pour la prospection d'hydrocarbures dans les régions de piémont de type « fold and thrust belt » comme la ceinture de Gaspé. Nous avons développé cette méthodologie en la testant au niveau de la partie centrale de la ceinture de Gaspé. Nous avons intégré les données géologiques et géophysiques disponibles pour l'étape d'imagerie sismique 2D afin de construire une image sismique directement en profondeur, ce qui a permis d'améliorer l'interprétation structurale, notamment la caractérisation des structures profondes et des failles majeures. Ce travail est suivi d'une modélisation de bassin afin d'évaluer le potentiel pétrolier. Cette étude s'effectue en plusieurs étapes : 1) La construction du modèle structural : L'intégration des données géologiques dans l'étape de la migration en profondeur avant sommation permet d'améliorer le rendu des images sismiques. Ces nouvelles images sont plus fiables et, étant migrées directement en profondeur, rendent les interprétations plus proches de la géometrie réelle du sous-sol. Ces informations permettent de construire un modèle géologique plus complexe et de mieux contraindre le modèle structural de la ceinture acadienne. Les nouvelles interprétations ont permis en particulier, de mieux comprendre la relation entre les ceintures acadienne et taconienne. 2) L'évolution du modèle structural : il a été possible de valider la cohérence de la géométrie structurale grâce aux techniques de restauration. Cependant ce procédé n'a été appliqué qu'au niveau du Synclinal du Lac des Huit-Miles sur les successions stratigraphiques siluro-dévoniennes de la ceinture acadienne : les formations cambroordoviciennes ont été déformées par les orogenèses taconiennes et acadiennes ce qui rend impossible leur restauration. Ce scénario cinématique a été utilisé pour comprendre l'évolution géodynamique de la ceinture de Gaspé et ainsi permettre de proposer une nouvelle géométrie plus favorable à la production et à la migration des hydrocarbures. 3) Évaluation du système pétrolier : Suite à l'étape de restauration, la modélisation de bassin avec le logiciel Temis2D® a été appliquée à la succession stratigraphique Silurien-Dévonien du synclinal du Lac des Huit-Milles et à l'anticlinal de Causapscal. Temis2D® a permis de prédire l'évolution de la roche mère et le degré de maturation ainsi que la génération et l'expulsion des hydrocarbures, en utilisant le modèle structural et les données géochimiques des puits de la ceinture de Gaspé.

Nous nous intéressons à la déconvolution aveugle de signaux bruités et plus précisément de signaux d'imagerie sismique. Pour réaliser l'inversion, nous souhaitons effectuer une sélection des statistiques d'ordre supérieur adaptées à la distribution du signal à déconvoluer. Pour cela, nous nous appuyons sur l'hypothèse de blancheur ou de non-gaussianité. Nous proposons une approche avec le taux d'information mutuelle comme mesure de blancheur et une autre basée sur la non-gaussianité du signal de sortie mesurée par la néguentropie. Après le développement d'algorithmes dans le domaine temporel et fréquentiel, nous caractérisons l'influence sur les critères du bruit additif présent sur les données. Nous démontrons que l'hypothèse de non-gaussianité est plus robuste à la présence d'un bruit additif sur les données blanc et gaussien. Cette approche permet pour des données synthétiques et réelles un très bon compromis entre la qualité de la déconvolution et l'amplification du bruit.

Cette thèse porte sur l'auscultation de fondations de pylônes électriques par contrôle non destructif. Une procédure d'imagerie par ondes sismiques est utilisée, ce qui consiste à générer une onde sismique au voisinage de la fondation avec une source vibratoire et à mesurer la vitesse de l'onde résultante à l'aide d'un réseau de géophones. Le travail présenté s'intéresse spécifiquement à la résolution du problème inverse. L'objectif est de concevoir un algorithme permettant de déterminer la géométrie de la fondation auscultée à partir des données enregistrées par les capteurs. Il s'agit d'un problème de tomographie de diffraction. Son caractère mal posé ainsi que la forte non-linéarité du modèle direct rendent sa résolution particulièrement difficile. Deux méthodes d'inversion ont été développées. Elles consistent toutes deux à reconstruire plusieurs cartes du sous-sol, chaque carte étant associée à une caractéristique physique du sous-sol. Dans les deux cas, un critère des moindres carrés régularisé est minimisé de façon itérative. Pour la première approche, peu d'a priori sont pris en compte. Elle se base sur une régularisation markovienne et la différence entre pixels voisins est pénalisée de manière à favoriser la reconstruction de régions lisses séparées par des contours francs. La seconde approche tient compte d'informations a priori plus nombreuses et plus précises. Elles sont liées à la structure particulière des fondations et aux caractéristiques de la terre et du béton. La méthode se base sur une segmentation du milieu à imager en deux régions. Des tests ont été réalisés sur des données simulées afin d'évaluer et de comparer les performances de ces deux méthodes.