Restauration par migration de signaux 2D sous-échantillonnés

Restauration par migration de signaux 2D sous-échantillonnés

Rebuilt by Migration of Undersempled 2D Signals

Alain Herique Wlodek Kofman 

CEPHAG URA 346 ENSIEG, BP 46 F-38402 Saint Martin d'Hères

Page: 
53-69
|
Received: 
11 October 1995
| |
Accepted: 
N/A
| | Citation

OPEN ACCESS

Abstract: 

Active geophysical prospecting techniques (seismics, ground penetrating radar) are based on the measurement of a reflected signal. This measurement is done for different points in the space and gives an image of reflected signal amplitude versus time (range) and versus one direction in space (x) . We use this data to rebuild an image of the subsurface. But the fact that directivity of transmitter and receiver is low, induces that the image versus x and t is not equivalent to the subsurface image (x, z). Migration techniques are developed to derive the subsurface image.

In this paper, we use the f - k migration technique which gives a solution of the wave equation in the frequency-wavenumber domain and we study the effect of undersampling in the time domain on this method . We describe a method to characterize this effect and estimate the real wavelet . We calculate the analytic result of migration of undersampled signal for a simple reflector. We find that the correctly sampled part of the signal is focused by migration and the aliased part is spread out from this calculation.

We deduce a method to rebuild the well-sampled 2D signal . The weakness of this technique is a decrease of the signal to noise ratio and the generation of spurious images . We discuss all aspects of this new technique on synthetic data .

The wavelet estimation and the frequency—wavenumber migration are applied to the ground penetrating data acquired in Antarctica by the radar developed for Mars'98 mission.

Résumé

Les techniques de prospection géophysique active (sismique, radar géophysique) reposent sur la mesure d'un signal réfléchi. Cette mesure est effectuée en différents points de l'espace, elle fournit une image de l'amplitude du signal réfléchi en fonction du temps et d'une dimension de l'espace (x) . A partir de ces données, nous cherchons à obtenir une coupe du sous-sol. Mais, du fait de la faible directivité des capteurs et des sources utilisées, l'image en temps et en espace n'est pas immédiatement transposable en une image fonction de deux dimensions de l'espace (x, z), d'où l'utilisation des techniques de migration pour la reconstruire.

Nous nous intéresserons ici à la migration f k qui donne une solution de l'équation de propagation dans l'espace des fréquences et des nombres d'onde, et à son comportement face au sous-échantillonnage temporel. Pour cela, nous développerons une méthode de caractérisation du sous-échantillonnage et d'estimation de I'ondelette émise. Ensuite, nous calculerons analytiquement le résultat de la migration d'un signal sous-échantillonné dans le cas d'un réflecteur simple. Ce calcul montre que la partie correctement échantillonnée du signal est focalisée par la migration, tandis que la partie repliée est étalée . Nous en déduirons une méthode de compensation du repliement fondée sur la cohérence spatiale du signal. Cette compensation — partielle — se fait au détriment du rapport signal sur bruit et peut générer des fantômes. Nous discuterons les différents aspects de cette nouvelle méthode : des simulations permettront de caractériser les performances de la méthode. Nous appliquerons l'estimation de I'ondelette et la migration aux signaux réels acquis en Antarctique par le radar géophysique développé pour la mission spatiale Mars'98.

Keywords: 

Sampling, Aliasing, 2D signals, f - k Migration, Geophysics, Seismics, Ground Penetrating Radar, Glaciology

Mots clés

Echantillonnage, Repliement, Signaux 2D, Migration f - k, Géophysique, Sismique, Radar géophysique, Glaciologie

1. Introduction
2. Migration
3. Le Sous-Échantillonnage
4. Migration Et Sous-Échantillonnage
5. Simulations
6. Données Réelles
7. Conclusion
  References

[Stolt] R .H . Stolt, « Migration by Fourier transform », Geophysics, Vol . 43, n°1 , 1978, pp . 23—48 .

[French] W.S .K . French, «Two dimensional and three dimensional migration of model—experiment reflector profils », Geophysical Prospecting, Vol . 39, n°3, 1974, pp . 269-277 .

[Max] J . Max et alü, Méthodes et techniques de traitement du signal et applications aux mesures physiques, Tome 1, 1985, Masson, Paris, 351 p .

[Barbin] Y. Barbin, W. Kofman, M . Elkine, M . Finkelstein, V. Glotov, V. Zolotarev , «Mars'96 subsurface radar », Int. Proc . Symp. on Radars and Lidars in Planetary Sciences, 2—4 September 1991, Cannes, France, pp . 51—58 .

[Herique] A. Herique, Reconstruction et Interprétation d'Images Géoradar, Application aux Tests en Antarctique du Radar de la mission Mars'98, 6 février 1995, Thèse de Docteur de 1'INPG, Grenoble, 267 p .

[Claerbout] J .B . Claerbout, Earth Sounding Analysis, Processing versus inversion , 1992, Boston, Blackwell Scientific Publications, pp. 186—197 .

[Spitz] S . Spitz, « Seismic trace interpolation in the F —X domain », Geophysics , Vol. 56, n°5, 1991, pp . 785—794 .

[Yilmaz] O . Yilmaz, Seismic Data Processing, 1987, Tulsa, Society of Exploration Geophysicist, 523 p .