Projet : ESTIME

Précédent : Méthodes numériques pour la migration Remonter : Résultats nouveaux Suivant : Algorithmes de points intérieurs non

   
Inversion électromagnétique

Participants : Michel Kern, Mohammed Ghemires.

Mots clés : Équations de Maxwell, radar géologique .

 Nous avons poursuivi l'application des méthode issues de l'inversion sismique au cas de données électromagnétiques provenant d'un radar géologique. On cherche à retrouver la permittivité diélectrique et la conductivité, du proche sous-sol (quelques mètres), pour détecter des structures enterrées. Le modèle est fourni par les équations de Maxwell et nous nous sommes restreint au système TE (transverse électrique) qui conduit à un modèle 2D pour un champ scalaire.

Nous avons abordé la résolution du problème inverse, à partir du logiciel de G. Vigo permettant de calculer la fonction objectif et son gradient, couplé au logiciel M1QN3 (Modulopt). Les premières expériences ont porté sur des modèles stratifiés. Nous avons ainsi inversé avec succès des modèles à deux et quatre couches. Nous mettons en particulier en évidence l'importance de disposer de mesures en des points différents de la source (mesures à « offset non-nul » dans le vocabulaire géophysique. Des modèles 2D comprenant des hétérogénéités (tuyaux de transport de gaz) dans un milieu homogène ont également été retrouvés de façon satisfaisante.

  

Figure 2: Inversion par radar géologique d'un sous-sol contenant trois conduites

Ces résultats sont illustrés sur la figure 2. Le milieu de départ est un sous-sol homogène contenant trois conduites. La figure représente le résultat de l'inversion, en partant du milieu homogène. Les trois conduites ont été localisées correctement.

Précédent : Méthodes numériques pour la migration Remonter : Résultats nouveaux Suivant : Algorithmes de points intérieurs non