FR 2801387 A1 20010525 FR 9914709 A 19991123 PERFECTIONNEMENTS AUX PROCEDES DE TRAITEMENT SISMIQUE METTANT EN #UVRE UNE DETERMINATION DES CORRECTIONS DYNAMIQUES RESIDUELLES La présente invention est relative à des perfectionnements aux procédés de traitement sismique mettant en #uvre une détermination des corrections dynamiques résiduelles. DOMAINE TECHNIQUE GENERAL ET PRESENTATION DE L'ETAT DE LA TECHNIQUE On rappelle que le principe de la prospection sismique consiste à provoquer, à l'aide d'une source sismique, un ébranlement dans le sous-sol et à enregistrer, à l'aide de capteurs, des données sismiques générées par cet ébranlement. Ces données sont ensuite traitées pour en tirer une information sur la géologie du sous-sol et notamment pour détecter la présence d'hydrocarbure. Classiquement, le traitement sismique comporte une phase de détermination du champ de vitesse et l'application aux traces sismiques d'une correction dynamique (NMO ou "normal move out" selon la terminologie anglo-saxonne généralement utilisée par l'homme du métier) destinée à compenser l'effet d'obliquité des trajets. Le but le plus évident de l'étape de détermination des corrections résiduelles de NMO est d'améliorer la qualité des traces-sommes (ou "stack" selon la terminologie anglosaxonne couramment employée). De nombreux traitements de corrections dynamiques résiduelles ont déjà été proposés. On connaît notamment des traitements permettant de les déterminer en s'appuyant sur des analyses ponctuelles de vitesses suivant des grilles à pas grossier économisant une partie du temps nécessaire à l'interprétation des spectres de vitesses (de Bazelaire, 1988 ; Doicin et al, 1995). Généralement, ces traitements consistent à améliorer l'énergie des traces-sommes à partir d'un balayage des résiduelles NMO fondé sur une sélection automatique des événements faisant appel à un critère basé sur la semblance. Le développement plus récent des études d'amplitude en fonction de l'offset (AVO ou "Amplitude Versus Offset" selon la terminologie anglosaxonne généralement utilisée) a rendu les déterminations de corrections dynamiques d'autant plus cruciales que les erreurs sur les déterminations de temps entraînent une mauvaise estimation des attributs d'AVO, sans compter d'autres difficultés liées à l'anisotropie, l'étirement des traces et la résonance des bancs minces. Le but de l'invention est de proposer un procédé de traitement sismique mettant en #uvre une détermination des corrections dynamiques résiduelles qui allie une résolution et une précision supérieures aux procédés courants. Il est courant en sismique d'utiliser la transformée de Radon parabolique pour atténuer les réflexions multiples. Dans le domaine espace-temps, une réflexion résiduelle après réalisation d'une correction dynamique préalable approchée peut être assimilée à la parabole : t = # + q x2 où t et x désignent respectivement les variables temporelle et spatiale, où q est un paramètre de courbure et où # est un paramètre temporel d'intercept parabolique à l'origine. La transformée de Radon parabolique directe consiste à effectuer la sommation des points qui se trouvent sur un jeu de paraboles de différentes courbures q pour des événements d'amplitude constante. On modélise les événements enregistrés au moyen de cette approximation et on résout un problème inverse pour adapter les paramètres du modèle à la réalité. Cette résolution est réalisée couramment par la méthode des moindres carrés (dite encore de norme L2) mais on peut la généraliser à des ordres de norme différents (Ln), notamment la norme L1 qui met en jeu la minimisation des différences prises en valeur absolue, ce qui rend le résultat moins sensible à la présence des queues de distribution. La transformée de Radon parabolique inverse consiste à repasser de l'espace (# - q) à l'espace (t - x) en modélisant les indicatrices à amplitude constante dérivées de cet espace (i# - q). Le premier traitement impliquant une transformée de Radon parabolique initialement proposé était réalisé dans le domaine temporel (Thorson and Claerbout, 1985). Il s'avéra lourd et gourmand en tempsmachine. Une solution dans le domaine fréquentiel a ensuite été développée (Hampson, 1986). Elle s'avéra rapidement très compétitive sur le plan industriel. Après une correction dynamique NMO initiale liée aux événements réels, chaque collection de traces (point miroir commun ou CDP "Common Depth Point" en 2D, pavé ou "bin" en 3D) est transformée dans l'espace parabolique de Radon pour fournir des événements idéalement focalisés selon des paquets distincts grâce à une inversion de norme Ln (L2 étant la norme standard). (Cette focalisation est toutefois souvent loin d'être optimale à cause des variations d'amplitude le long des événements et des effets de bord de collection). La séparation entre énergies primaire et multiple est fondée sur un critère lié à la courbure. Les événements voisins ou en deçà de la courbure nulle sont considérés comme des réflexions réelles et mutées. La transformée de Radon parabolique inverse produit un modèle de réflexions multiples à soustraire de la collection d'origine. La séparation entre réflexions réelles et multiples a d'abord été réalisée suivant une courbure constante jouant le rôle de limite de séparation. Ce traitement a été perfectionné par l'emploi d'un mute (mise à zéro des échantillons entre le début de la trace et une limite donnée) variable en temps # et courbure q pour éviter de briser un certain étalement des événements réels superficiels. Un tel traitement permet une meilleure séparation des événements superficiels, mais entraîne des allers et retours entre le domaine temporel et le domaine fréquentiel. Quoiqu' indispensable, il est de ce fait relativement coûteux en temps de calcul. De fait, on doit s'efforcer de comparer au mieux les focalisations obtenues avec l'approche classique d'énergie maximum dérivée d'une simple sommation le long des paraboles et sa contrepartie obtenue par inversion (ou déconvolution) suivant la norme usuelle L2. Plus récemment, des améliorations additionnelles importantes ont été obtenues en utilisant la norme L1 soit dans le domaine temporel (Hugonnet and Canadas, 1997), soit dans le domaine fréquentiel (Sacchi and Ulrych, 1995, Herrmann et al, 1999). Ces développements avaient surtout pour but d'une part d'éviter d'utiliser des programmes d'interpolation pour résoudre les problèmes de repliement spectral ("aliasing" selon la terminologie anglosaxonne) liés aux fortes courbures des multiples lents et d'autre part de garantir une meilleure préservation des amplitudes en réduisant notablement les artefacts liés aux fortes amplitudes et aux effets de bord. PRESENTATION DE L'INVENTION Le traitement proposé par l'invention permet de déterminer, à partir d'enregistrements en 2D ou 3D, des corrections dynamiques résiduelles en utilisant la transformée de Radon parabolique dans le cadre d'une inversion de norme Ln, au lieu d'avoir recours aux techniques mettant en #uvre un traitement sur les énergies des traces. Canadas et Dunand (1995) ont déjà proposé une technique mettant en oeuvre une inversion de norme Ln utilisant la transformée de Radon parabolique afin d'optimiser la correction de NMO résiduelle en vue des applications liées à l'AVO. Cette méthode est limitée à la prise en compte d'événements isolés dont une technique basée sur des corrélations peut affiner la détermination. L'invention propose quant à elle un procédé mettant en jeu une transformation morphologique globale de l'espace (# - q) après inversion mettant en oeuvre une transformée de Radon parabolique directe suivi d'une transformée de Radon parabolique inverse de l'espace (# - q) vers l'espace (t- x). Les phases peuvent se décrire comme suit : - on réalise une inversion mettant en #uvre une transformée de Radon parabolique, - on calcule une matrice d'énergie sur un corridor sélectionné dans l'espace de Radon, - on détermine les points de ladite matrice qui correspondent à des maxima discrets d'énergie, - on met en oeuvre un traitement de sélection sur ladite matrice pour ne garder, pour le corridor sélectionné, qu'un point d'énergie maximale pour chaque temps de l'espace de Radon, lorsqu'un tel point existe, - après transposition de la matrice d'énergie, on met en oeuvre un traitement d'interpolation et de corrections statiques suivant l'axe des courbures, qui concentre les points d'énergie maximale sur la courbure nulle, - on met en #uvre sur les données ainsi obtenues une transformée de Radon parabolique inverse. L'inversion de départ est une inversion de norme Ln, n étant un nombre entier ou fractionnaire. Notamment, on tire avantage globalement de l'amélioration de la focalisation des événements en recourant à une inversion de norme L2, voire de norme L1 pour un meilleur résultat, afin de déterminer de façon précise la position des événements significatifs réels de l'espace de Radon avant de les aligner suivant la courbure nulle. Un tel procédé est en outre avantageusement complété par les différentes caractéristiques suivantes : - le calcul de la matrice d'énergie met en #uvre une transformée de Hilbert ; - le corridor sélectionné dans l'espace de Radon est délimité par une courbure minimale et une courbure maximale ; - le traitement de sélection se fait en fonction de l'amplitude des points d'énergie maximale ou d'un paramètre, tel que la kurtose, calculé pour chaque point d'énergie maximale en fonction de son amplitude et de celle de points qui l'entoure ; - le traitement d'interpolation met en #uvre à partir des maxima sélectionnés, une loi de corrections statiques continue en #, appliquée suivant le contenu de l'espace de Radon ; - le traitement d'interpolation utilise un interpolateur polynômial à 9 coefficients, et détermine les points de l'espace de Radon pour lequel ce polynôme est maximal ; - le traitement d'interpolation met en oeuvre une interpolation linéaire ou une interpolation de convolution cubique des corrections statiques selon l'axe des #; - les corrections statiques sont appliquées le long de l'axe des courbures (#-q) à l'espace (-q) ; - on appliquera au même espace la transformation parabolique de Radon inverse ; - pour mettre en #uvre un traitement d'antimultiples simultané, on isole la partie de l'espace (i-q) correspondant aux multiples (fortes valeurs de courbure q) et on lui fait subir séparément la transformation de Radon parabolique inverse. On soustrait le résultat obtenu au résultat de transformée de Radon parabolique inverse global obtenu précédemment ; la partie de l'espace qui correspond aux multiples (courbures élevées) se détermine par rapport à la séparation entre événements réels et multiples définie par une ligne de "mute". On notera par ailleurs que lorsque, pour un temps de l'espace de Radon donné, il n'existe pas de point d'énergie maximale, on détermine une valeur de courbure qui correspondrait à un tel point par interpolation calculée à partir des maxima discrets d'énergie autour de cette valeur. PRESENTATION DES FIGURES D'autres caractéristiques et avantages de l'invention ressortiront encore de la description qui suit. Cette description est purement illustrative et non limitative et doit être lue en regard des dessins annexés sur lesquels - la figure 1a est un exemple d'enregistrement sismique synthétique dans le domaine (t, x) sur lequel on a porté une collection de traces sismiques appartenant à un CDP incluant des réflexions primaires mal corrigées avec un faible écart de NMO et quelques réflexions multiples à NMO plus marqué ; - la figure 1 b correspond dans le plan de Radon au résultat obtenu à partir de l'enregistrement de la figure 1a par une transformée de Radon parabolique correspondant à une sommation simple suivant une série de paraboles ; - la figure 1 c correspond dans le plan de Radon au résultat obtenu à partir de l'enregistrement de la figure 1a en mettant en #uvre une inversion de norme L2 impliquant une transformée de Radon parabolique ; - la figure 1d correspond dans le plan de Radon au résultat obtenu à partir de l'enregistrement de la figure 3a en mettant en #uvre un traitement du type de celui proposé par l'invention, ; - la figure 1e correspond au résultat obtenu dans le domaine (t, x) après transformation de Radon parabolique inverse du résultat de la figure 1d ; - la figure 2 illustre dans l'espace de Radon une transformée de Hilbert sur un corridor choisi dans ledit espace ; - la figure 3 illustre dans l'espace de Radon l'alignement des événements réels obtenus après mise en #uvre du traitement par interpolation d'une loi de corrections statiques et application de celle-ci ; - la figure 4 représente un enregistrement sismique réel dans le domaine (t, x) (collection de bin - partie gauche de la figure 4) et le résultat obtenu après mise en #uvre sur ledit enregistrement d'un traitement conforme à l'invention (partie droite de la figure 4). DESCRIPTION DETAILLEE D'UN MODE DE MISE EN #UVRE POSSIBLE POUR L'INVENTION Rappels théoriques liés à la transformée de Radon parabolique Dans le domaine temps-espace ( t - x ), après un traitement NMO classique dérivé d'une analyse de vitesses réalisée sur les noeuds d'une grille peu dense, une réflexion résiduelle peut être au premier ordre assimilée à une parabole : (1) t = # + q x2 Cette parabole correspond idéalement dans l'espace de Radon parabolique à un point ( # , q). Inversement, un point du temps-espace correspond à une droite pentée dans l'espace de Radon, ce qui explique que toute forte anomalie locale d'amplitude peut générer des artefacts linéaires. Soient Vr la vitesse réelle optimale pour la sommation (ou vitesse rms dans le cas de la migration) et Vc la vitesse obtenue par analyse préliminaire, le temps résiduel après un premier NMO de calcul est : (2) t' = t - At avec : (3) At = t - (t2 - x2 / VC2)112 D'où: (4) t' = (t2 - x2 / Vc2)1/2 = (i + x2 l Vr2 - x2 VC2)1/2 Ce qui, par approximation donne : (5) t' =#(1 + 1/2 (x2 / #2). (1 / Vr 2 1 /Vc2)) et fournit l'expression de la courbure (6) q = 1/2(1/#).(1/Vr2-1/Vc2) Présentation des différentes étapes d'un traitement conforme à un mode de mise en oeuvre possible pour l'invention On met initialement en #uvre sur un enregistrement sismique donné, ayant éventuellement subi des traitements antérieurs et notamment des traitements de filtrage, une inversion impliquant une transformée de Radon parabolique. Cette inversion est de norme Ln où n est un nombre entier ou fractionnaire, L2 étant la norme standard, L1 conduisant à de meilleurs résultats avec un temps de calcul supérieur. Comme l'illustrent en effet les figures 1a à 1 c, une sommation simple suivant une série de paraboles (figure 1 b) donne une focalisation de qualité moindre des événements dans le plan de Radon, tandis qu'avec une inversion de norme Ln (en l'occurrence de norme L2), les événements significatifs sont beaucoup mieux focalisés (figure 1c). Une inversion de norme L1 permet une focalisation encore meilleure comme cela a été vérifié par un certain nombre d'essais sur des données synthétiques et réelles. Une fois cette inversion de norme Ln effectuée, on met en #uvre les différentes étapes a) à f) suivantes. a) Dans un premier temps, on détermine un corridor de l'espace (# - q ) pour séparer les énergies primaire et multiple. Ce corridor est par exemple limité par des valeurs de qmin et qmax minimale et maximale. La valeur qmin correspond à la limite des réflexions primaires surcorrigées, tandis que la valeur de qmax correspond à la limite des réflexions sous-corrigées non considérées comme multiples. C'est ce qu'illustre l'encadré de la figure 1c, qui a été repris à une échelle plus importante sur la partie gauche de la figure 2. En l'occurrence, cet encadré a été choisi de façon à retenir deux ensembles d'échantillons de part et d'autre de la courbure nulle (voir figure 3). Cette détermination peut encore être raffinée au moyen de corridors variables en # et q pour améliorer la restitution des événements superficiels. b) Dans une deuxième étape, le corridor sélectionné subit un traitement destiné à le convertir en une matrice d'énergie discrète, dans l'espace (#, q). Ce traitement est effectué, ainsi que l'illustre la figure 2, par une transformation de Hilbert, dans le domaine des fréquences, le long de l'axe des #. c) On repère ensuite les maxima locaux du corridor transformé, en considérant qu'ils sont entourés par 8 à des échantillons de moindre amplitude. A cet effet, chaque maximum retenu de la matrice résultant de la transformation de Hilbert le long de l'axe des Test dans un premier temps remplacé par un maximum obtenu par interpolation. Cette interpolation utilise 9 coefficients à calculer (a à i) d'un polynôme défini comme suit : (7) A = a x2 y2 + b x2 y + c x y2 + d x2 + e x y + f y2 + g x + hy + i où x et y jouent les rôles de # et de q et où i correspond au centre d'amplitude maximale de chaque groupe de 9 points sélectionnés. Les 9 coefficients a à i sont déterminés en résolvant un système d'équation linéaire obtenu en affectant à l'amplitude de chacun des 9 points ses coordonnées x et y. L'interpolation consiste ensuite à résoudre le système d'équations résultant de l'annulation des dérivées partielles : 'DA/'Dx = 2 a x y2 + 2 b x y + c l + 2 d x + e y + g (8) ôA/ôy = 2ax2y + bx2 + 2cxy + ex + 2fy + h la solution de ce système fournissant en effet la position optimale du maximum. Cette interpolation peut être initialisée en résolvant le système linéaire : 2dx+ey+g=0 ex+2fy+h=0 (9) (On notera que ce choix initial se justifie du fait que la taille canonique (ou normalisée) de la moitié du côté du carré à 9 points est comprise entre 0 et 1). Si le système est mal conditionné, on peut revenir à une double interpolation parabolique effectuée à partir de trois points centraux le long de l'axe des t et trois points centraux le long de l'axe des q. Sinon, les premières estimations sont entrées alternativement dans les dérivées partielles pour déterminer des solutions du premier degré en x pour la première équation, en y pour la seconde. Au cours des itérations, quelques coordonnées peuvent s'aventurer au delà des limites du carré à 9 points d'origine. Elles sont alors remplacées par les valeurs de l'itération précédente. En cas de non-convergence après un nombre préfixé d'itérations, on a à nouveau recours à la double interpolation parabolique en x et y . L'interpolation qui vient d'être décrite dans le cas d'un traitement 2D (en particulier, pour des collections en point miroir résultat d'une acquisition en 3D) a été facilement étendue en 3D pour déterminer des maxima de spectres d'énergie. Dans ce cas, on utilise 27 points et les dérivées partielles sont calculées par rapport aux coordonnées x, y et z. On part de la même façon d'un système linéaire et on transpose les premières estimations dans les équations aux dérivées partielles. d) On met ensuite en #uvre sur les maxima locaux ainsi retenus un traitement de sélection, destiné à permettre de ne garder qu'un échantillon maximal pour chaque temps #, si tel est le cas. Cette sélection est fonction de l'amplitude d'un point maximal ou d'une valeur dérivée d'une double interpolation parabolique en # et q, ou encore de la kurtose de 9 points considérés pour chaque maximum (l'échantillon d'énergie maximale et les huit échantillons qui l'entourent), ou toute combinaison appropriée de ces paramètres. On rappelle que la fonction de kurtose s'exprime: (10) K = # h4/(#h2)2-3 où h sont les amplitudes des échantillons. e) On met ensuite en #uvre un traitement d'interpolation et d'application d'une loi de corrections statiques (transformation morphologique de l'espace de Radon). 1) Après avoir obtenu une série écrémée de maxima interpolés, on détermine une loi de corrections statiques. Celle-ci est interpolée selon l'axe des # à partir des maxima optimaux en mettant en oeuvre une interpolation linéaire ou une interpolation de convolution cubique traitant des points à espacements inégaux (Keys et Pann, 1993). 2) Des corrections statiques sont ensuite appliquées dans la direction q pour chaque échantillon interpolé en #. Ces corrections statiques permettent de concentrer les points d'énergie maximale sur la courbure nulle. 3) On peut dégager un champ de vitesses modifiées à partir des maxima initiaux qui peut être lissé au moyen d'une régression mobile locale en faisant glisser un élément de plan de moindres carrés (Grubb et Walden, 1995, procédé LOESS, et Adler, 1998), procédé complété par des techniques d'écrémage qui augmentent sa robustesse. Le résultat du traitement décrit jusqu'ici est illustré par la figure 3 dont la partie gauche représente l'espace de Radon (# - q) initial et la partie droite le résultat obtenu après application de la loi de corrections statiques. Comme on peut le voir sur cette partie droite de la figure 3, les événements réels sont alignés sur la courbure nulle. f) A l'issue du traitement par corrections statiques, on réintègre le corridor matriciel dans l'espace de Radon (figure 1d) et on met en #uvre une transformée de Radon parabolique inverse de l'espace de Radon parabolique globale. Comme on peut le constater sur la figure 1e, les événements réels sont alors horizontalisés. On a ainsi appliqué de façon optimale un jeu de résiduelles de NMO. On notera que sur la figure 1e, l'élément le plus superficiel n'a pas été corrigé en raison d'une courbure trop forte dépassant les limites du corridor sélectionné. La figure 4 montre des résultats obtenus sur un bin d'acquisition 3D réelle en mettant en #uvre un traitement conforme à la mise en #uvre de l'invention. On constate sur cette figure 4 que les événements réels sont horizontalisés (partie gauche avant traitement ; partie droite après traitement). On notera que le champ lissé des vitesses déterminé lors du traitement de corrections statiques peut être utilisé par lui-même pour la construction du modèle des vitesses des données avant migration. Par ailleurs, un contrôle de qualité peut être mis au point en comparant les champs lissés et non lissés pour isoler les zones douteuses, notamment dans le cas d'une application de station de travail interactive dédiée à l'analyse des vitesses. De plus, on peut associer ces expertises avec des sélections réalisées sur la base de pointé automatique d'événement quelconque ou d'horizons continus (Baixas et al, 1980 ; de Bazelaire, 1988 : Doicin et al , 1995) en 2D ou 3D. On notera en outre que le traitement qui vient d'être décrit peut être utilisé pour mettre en #uvre un traitement antimultiples. A cet effet, on sélectionne dans le domaine de Radon un sous-corridor qui correspond au domaine des réflexions multiples, on met en #uvre le procédé décrit par l'invention pour déterminer la transformation morphologique de l'espace de Radon parabolique, puis on effectue les transformées de Radon paraboliques inverses : - de l'espace global - de l'ensemble des réflexions multiples de cet espace caractérisé par les courbures élevées définies par rapport à la ligne de "mute", et on soustrait le second résultat du premier. BIBLIOGRAPHIE - Adler, F., 1998, Robust estimation of dense 3D stacking velocities from automated picking. 60th Eur. Assoc. Expl. Geophys., Extended abstracts, Session :P139. - Baixas F., Leflaive R. and Paturet D., 1980, Automatic residual statics in 3D seismic. 50th Annual Internat. Mtg., Soc. Expl. Geophys, Expanded Abstracts, Session: G107. - Canadas G. and Dunand J-P., 1995, AVO study in the parabolic Radon transform domain. 65th Annual Internat. Mtg., Soc. Expl. Geophys, Expanded Abstracts, Session : SP2. 2. - de Bazelaire E., 1988, NMO revisited ; inhomogeneousmedia and curved interfaces. Geophysics 53,2, 143-157. - Doicin D., Johnson C. and Hargreaves N., 1995, Improved 3D processing turnaround via automatic velocity picking. 57th Eur. Assoc. Expl. Geophys., Extended abstracts, Session : B033. - Grubb H. and Walden A., 1995, Smoothing seismically derived velocities. Geophysical Prospecting, 43, 1061-1062. - Hampson, D., 1986, Inverse velocity stacking for multiple attenuation. 56th Annual Internat. Mtg., Soc . Expl. Geophys, Expanded Abstracts, Session : S6. 7. - Herrmann P., Mojesky T., Magesan M., 1999, Amplitude preserving Radon demultiple: beyond sampling rate and aperture limitations. CSEG meeting abstracts , 73-74. - Hugonnet P. and Canadas G., 1997. Regridding of irregular data using 3D Radon decomposition. 67th Annual Internat. Mtg., Soc . Expl. Geophys, Expanded Abstracts, Session : SP3. 5. - Keys R. and Pann K., 1993, Geophysical applications of cubic convolution interpolation. 63th Annual Internat. Mtg., Soc . Expl. Geophys, Expanded Abstracts, Session : CT1. 5. - Sacchi M. and Ulrych T., 1995, Model re-weighted least-squares Radon operators. 65th Annual Internat. Mtg., Soc . Expl. Geophys, Expanded Abstracts, Session : PP4. 5 . - Thorson J. and Claerbout J , 1985, Velocity-stack and slant stochastic inversion . Geophysics, 50,12, 2727-2741. REVENDICATIONS 1. Procédé de traitement sismique selon lequel on met en #uvre sur des enregistrements de traces sismiques une détermination des corrections dynamiques résiduelles, caractérisé par les différentes étapes suivantes : - on réalise une inversion mettant en #uvre une transformée de Radon parabolique, - on calcule une matrice d'énergie sur un corridor sélectionné dans l'espace de Radon, - on détermine les points de ladite matrice qui correspondent à des maxima discrets d'énergie, - on met en oeuvre un traitement de sélection sur ladite matrice pour ne garder, pour le corridor sélectionné, qu'un point d'énergie maximale pour chaque temps de l'espace de Radon, lorsqu'un tel point existe, - on met en #uvre un traitement d'interpolation suivant l'axe des intercepts paraboliques (#) et de corrections statiques suivant l'axe des courbures, qui concentre les points d'énergie maximale sur la courbure nulle, - on met en #uvre sur les données ainsi obtenues une transformée de Radon parabolique inverse. 2. Procédé selon la revendication 1, caractérisé en ce que l'inversion est une inversion de norme Ln, où n est un nombre entier ou fractionnaire. 3. Procédé selon la revendication 2, caractérisé en ce que l'inversion est de norme L2. 4. Procédé selon la revendication 2, caractérisé en ce que l'inversion est de norme L1. 5. Procédé selon l'une des revendications précédentes, caractérisé en ce que le calcul de la matrice d'énergie met en #uvre une transformée de Hilbert. 6. Procédé selon l'une des revendications précédentes, caractérisé en ce que le corridor sélectionné dans l'espace de Radon parabolique est délimité par une courbure minimale et une courbure maximale. 7. Procédé selon l'une des revendications précédentes, caractérisé en ce que le traitement de sélection se fait en fonction de l'amplitude des points d'énergie maximale ou d'un paramètre, tel que la kurtose, calculé pour chaque point d'énergie maximale en fonction de son amplitude et de celle de points qui l'entoure, des tels points d'énergie maximale. 8. Procédé selon l'une des revendications précédentes, caractérisé en ce que le traitement d'interpolation met en oeuvre à partir des maxima sélectionnés, une loi de corrections statiques continue en #, appliquée suivant le contenu de l'espace de Radon. 9. Procédé selon la revendication 8, caractérisé en ce que le traitement d'interpolation utilise un interpolateur polynomial à 9 coefficients, et détermine les points de l'espace de Radon parabolique pour lequel ce polynôme est maximal. 10. Procédé selon l'une des revendications précédentes, caractérisé en ce que le traitement d'interpolation met en #uvre une interpolation linéaire des corrections statiques selon l'axe des #. 11. Procédé selon l'une des revendications 1 à 9, caractérisé en ce que le traitement d'interpolation met en #uvre une interpolation de convolution cubique des corrections statiques selon l'axe des T . 12. Procédé selon l'une des revendications précédentes, caractérisé en ce que pour mettre en #uvre un traitement antimultiples, on isole la partie de l'espace de Radon qui correspond aux multiples, on lui fait subir séparément une transformation de Radon parabolique inverse et on soustrait le résultat obtenu du résultat de transformée de Radon parabolique inverse global obtenu précédemment.