OA10264A - Procédé de réalisation d'un cube 3d en traces proches à partir de donnés acquises en sismiques marine réflexion - Google Patents

Description

Procédé de réalisation d’un cube 3D en traces proches à partir de données acquises en sismiques marine réflexion

La présente invention concerne un procédé deréalisation d'une image du sous-sol en trois dimensionsdite cube 3D en traces proches à partir de donnéesacquises en sismique réflexion marine et elle serapporte, plus particulièrement, à la génération d'unmini cube 3D réalisé à bord d'un navire d'exploration sedéplaçant dans la zone à explorer suivant uneconfiguration bien connue.

La sismique marine à trois dimensions ou 3D estutilisée de préférence à la sismique marine à deuxdimensions ou 2D, pour améliorer le positionnementlatéral des réflecteurs sur lesquels se réfléchissent desondes émises par une source d'émission, les ondes réflé-chies étant reçues et enregistrées sur des récepteurs-enregistreurs régulièrement espacés le long d'un câbleenregistreur ou flûte (streamer) remorqué dans ladirection de navigation du navire. La sismique marine 3Dpermet également d'améliorer la résolution latérale dansla direction perpendiculaire (cross-line) à la directiond'acquisition donnée par la direction de la flûte.

En sismique marine 3D, on remorque simultanémentune ou plusieurs sources et un groupe de flûtes sensi-blement parallèles, de manière à acquérir sur un trajetdu bateau (ligne naviguée) plusieurs lignes de donnéessismiques, dites lignes d'acquisition dont chacune estdéfinie par un couple source-flûte ; ces lignes permet-tent, après un traitement adéquat, d'obtenir une image entrois dimensions du sous-sol, généralement appelé cubesismique. En utilisant plusieurs sources d'émission et un . > 'i . t k » Λ * * v i b i, t) 4 plus grand nombre de flûtes, par exemple 2 ou 3 sources et de 4 à 6 flûtes, on couvre à moindre coût une plus grande zone d'exploration avec des données finales de cohérence accrue.

La nature des sources d'émission, l’écartementdesdites sources et des flûtes, la longueur de chaqueflûte qui varie entre 1000 et 5000 m, le nombre derécepteurs-enregistreurs prévus sur chaque flûte (desmilliers d'hydrophones sur des centaines de canauxindépendants) , la distance entre récepteurs-enregistreursconsécutifs ou inter-trace sont adaptés à 1'imageried'objectifs situés à des profondeurs de quelques kilo-mètres. Les sources et les récepteurs-enregistreurs sontdotés de moyens pour les localiser spatialement de tellesorte qu'on puisse traiter les données sismiquesrecueillies en position géographique exacte à moins d'unedizaine de mètres près.

Les points milieux géométriques de chaque couplesource-récepteur sont regroupés à l'intérieur d'unegrille de "bins" ou cellules rectangulaires formant unmaillage spatial régulier dans un plan X, Y, X étant depréférence la direction de référence de navigation et Yétant la direction perpendiculaire à la direction X. Depréférence, les cellules de la grille ont un pas égal àla moitié de l'inter-trace suivant la direction X et dansla direction Y, à la moitié de la plus petite distancelatérale soit entre sources lorsqu'il y en a plusieurs,soit entre flûtes dans leur configuration géométriquethéorique. La grille de bins utilisée dynamiquement esten réalité une grille géographique statique calculée unefois pour toutes avant la campagne d'exploration ensupposant une navigation du navire et une géométrie dudispositif remorqué (source(s) et flûtes) théoriques etparfaites. j ύ1G284

Pour l'exploitation finale des données sismiquesenregistrées, suivant les lignes de traitement en X ("online") ou en Y ("crossline") définies par le maillage, ilest nécessaire d'éviter que le remplissage en temps réelde la grille de tins (binning en anglais) d’une ligned'acquisition sur l'autre se fasse avec des trous ou unesurcouverture. Cependant, les coordonnées des positionsfinales des sources d'émission et des récepteurs sontseulement obtenues dans une phase ultérieure de post-calcul, de manière à filtrer les données de position-nement d'une façon plus élaborée qu'en temps réel. Lespositions qui se révèlent aberrantes en seconde analysesont ultérieurement recalculées afin de fournir en fin demission une qualité de données de positionnementpratiquement constante.

La précision finale (de l'ordre de quelquesmètres), demandée par certaines étapes du traitement desdonnées sismiques, approche les limites des méthodesmises en oeuvre, au moins en ce qui concerne les posi-tions latérales des sources et des récepteurs, souventmoins bien contraintes par les mesures de positionnementque les positions longitudinales, et sujettes à certainsbiais systématiques que les données de positionnementseules ne permettent pas forcément de détecter et decorriger. C'est principalement pour détecter et corrigerces erreurs et ces biais que la réalisation à bord d'uncube 3D de données sismiques sommées en couverturemultiple réduite à été adoptée, dès qu'elle a étépossible (depuis 1992) . La couverture multiple réduiteutilisée pour l'obtention du cube traces proches selon l'art antérieur consiste à travailler sur les traces de déports courts après correction dynamique.

En produisant, à partir d'un cube de données traitées conventionnellement (mais en couverture réduite), des coupes à temps constant (figure 1) et des sections "crossline" c'est-à-dire perpendiculaires à la direction X (figure 2), on peut effectivement mettre en évidence des erreurs de positionnement, qui se traduisent par des déphasages temps entre les traces associées à une 5 même ligne naviguée ou entre les traces associées à deslignes naviguées adjacentes.

On observe que parmi les cubes en couvertureréduite ce sont ceux en couverture unité utilisant laseule trace la plus proche de chaque ligne de traitement 10 qui offre la plus grande sensibilité aux erreurs depositionnement. Dans ce cas, on observe sur l'image dufond de l'eau ou des horizons ou réflecteurs les plusproches dudit fond, par exemple situés à 200 ms sousledit fond, les deux phénomènes suivants (figure 1). 15 Le premier phénomène est lié à la géologie des premières couches du sous-sol ; il s'agit de l'image parexemple de paléochenaux quaternaires remplis par dessédiments d'impédance acoustique plus forte que le fondde l'eau et désignés par la référence 1. Une impédance 20 acoustique plus forte engendre un coefficient deréflexion plus important c'est-à-dire un effet de fonddur.

Le deuxième phénomène est dû aux artefacts liésaux imperfections de l'acquisition et du traitement, et 2 5 il est matérialisé sur la figure par un motif approxima-tivement périodique sous forme de rubans orientés dans ladirection des trajets de navigation et désignés par laréférence 2.

Lorsque les coordonnées utilisées dans le calcul 30 d'un cube 3D en traces proches sont entachées d'erreurs,ces dernières engendrent des décalages en temps qui sontégalement orientés sous forme de rubans, éventuellementlocalement périodiques s'il s'agit de biais systéma-tiques . G 1 G 2 6 4

La partie principale du motif périodique observéest un artefact du traitement et de la mise en oeuvre ducube 3D en traces proches selon l'art antérieur. Unecause sinon principale, du moins importante, des arte-facts est la distorsion fréquentielle différentiellerésultant de la correction dynamique entre des déports(offsets) courts variables d'une ligne de traitement oud'acquisition à l'autre. Or, sur les horizons ou réflec-teurs proches du fond de l'eau, la distorsionfréquentielle est très importante et varie rapidement enfonction du déport, le déport étant la distance séparantune source d'un récepteur donné. L'effet de distorsion fréquentielle s'exprimepour un milieu homogène par la formule : F' (V.TO) dans laquelle : F' est la fréquence après correction dynamique, F est la fréquence avant correction dynamique, X est le déport pour un couple source-récepteurdonné, V est la vitesse de correction dynamique,

To est le temps double pour un déport nul (sourceet récepteur confondus).

Cet effet de distorsion fréquentielle estd'autant plus important que les horizons sont super-ficiels ou très proches du fond de l'eau et que le déportminimum varie d'une ligne de traitement à l'autre sur unemême ligne naviguée. Dans le cas de la sismique 3Dmarine, cette variation de déport est en fait due à unevariation de la composante latérale Y de déport et elle 6 augmente en particulier avec le nombre de lignes d'acquisition ou de flûtes, à taille de bin égale.

Lorsque ces effets sont constatés et que leur élimination est jugée nécessaire, il est proposé généra- 5 lement un traitement consistant à appliquer sur lessections à temps constants un filtre spatial (en Kx, Ky) ,transformée de Fourier spatiale dans le domaine descoordonnées x et y, respectivement dans la direction desflûtes et dans la direction perpendiculaire auxdites 10 flûtes, de manière à effectuer ainsi une déconvolutionspatiale globale du phénomène.

Cependant, un tel traitement réduit l'effet etnon la cause et il risque d'éliminer d'autres distorsionsassociées à des erreurs de positionnement. En effet, avec 15 deux sources et un groupe de quatre flûtes, un biaissystématique de position d'une source ou d'une flûte setraduit par un décalage respectivement toutes les deux outoutes les quatre lignes de traitements Or, lapériodicité normale des déports proches dans cette 20 configuration est de huit lignes avec des harmoniquesspatiaux à quatre et deux lignes. Un filtre spatialoptimisé pour éliminer ces artefacts risquerait donc deréduire la capacité de détection de biais de positionnement des sources ou des flûtes. 25 De plus, comme toute convolution nécessitant une transformée de Fourier, il se produit des effets de bordsgênants lorsqu'on utilise un cube incomplet. Or, c'est encours de mission lorsqu'on dispose seulement d'un cubetrès incomplet, qu'on a besoin d'avoir des informations 30 sur les premières couches du sous-sol marin.

On pourrait envisager la réalisation à bord du navire d'un cube 3D final en couverture complète.Toutefois, cela nécessiterait des moyens embarqués prohi-bitifs (informatique très lourde, personnel supplémen- 35 taire et très qualifié), ce qui, en général, n'est pas υ 1 g 2 6/- justifié par le gain de temps escompté dans la fourniture du résultat final.

La présente invention a pour but de proposer unprocédé simple qui atténue grandement le phénomène desartefacts précédemment décrits dus au traitement des don-nées sismiques, au moins pour les premières couches dusous-sol marin, et qui permet ainsi d'identifier deserreurs de positionnement qui autrement auraient étémasquées par lesdits artefacts.

La présente invention a pour objet un procédé deréalisation d'un cube 3D en traces proches à partir dedonnées acquises en sismique marine réflexion, du typedans lequel : a) on procède à l'acquisition desdites données endéplaçant dans une direction de référence D, i) au moins une source, ii) plusieurs flûtes sismiques espacées de lasource, distantes latéralement les- unes desautres et comportant chacune des récepteurs-enregistreurs régulièrement espacés le long dechaque flûte, et iii) des moyens de localisation spatiale de chaquesource et de chaque flûte, de manière à obtenir simultanément, pour chaque lignenaviguée, d'une part des données sismiques sous la formede traces, fournies par les récepteurs-enregistreurs, lelong d'autant de lignes d'acquisition, de préférencesensiblement parallèles à la direction de référence D,que le produit du nombre de flûtes par le nombre desources utilisés et d'autre part, des données de locali-sation spatiales de chaque source et des récepteurs-enregistreurs pour fournir la distance ou déport séparantchaque source de chaque récepteur ainsi que les coordon-nées du point milieu du couple source/récepteur-enregistreur associé à chaque trace, Ο b) on réalise une image du sous-sol en trois dimensionsdans un système d'axes (X, Y, T) sous la forme dite cube3D en attribuant à chaque cellule (bin) d'un ensemble decellules formant un maillage spatial régulier dans unplan X, Y, représentatif du plan d'acquisition desdonnées sismiques, une trace centrale élaborée à partirdes traces enregistrées et représentée au centre deladite cellule suivant l'axe temporel (T), ledit procédé se caractérisant en ce que pour chaqueligne naviguée on constitue la trace centrale de chaquecellule à partir d'au moins une première trace enregis-trée dont le point milieu associé est situé à l'intérieurde ladite cellule et dont le déport est le plus prochepossible d'un déport constant égal au plus petit déportcommun à au moins l'ensemble des traces de la lignenaviguée.

Un avantage de la présente invention réside dansle fait qu'outre une réduction substantielle du phénomènedes artefacts, elle permet de restituer avec une plusgrande précision les premiers milieux par le cube tracesproches.

Un autre avantage est la possibilité de détecterplus facilement les biais de positionnement des récep-teurs de tête des flûtes et/ou des sources et defaciliter l’application de traitements ultérieurs plussophistiqués tels que le DMO (Dip Move Out) 3D, lesanalyses de vitesse fines pour la migration profondeuravant somme. D'autres avantages et caractéristiques ressorti-ront mieux à la lecture du procédé selon 1'inventionainsi que des dessins annexés sur lesquels : - la figure 1 est une section à temps constant (92 ms) obtenue à partir de l'exploitation d'un cube traces proches en couverture un selon l'art antérieur, 9 G 1 02 64 sur laquelle sont représentés les deux phénomènes qu'on observe pour un fond marin situé à environ 90 ms, la fiqure 2 est une section temps dans ladirection perpendiculaire (cross-line) au trajet dunavire obtenue selon l'art antérieur, à partir del'exploitation du même cube traces proches utilisé pourla figure 1. - la figure 3 est une vue schématique de dessusd'un dispositif d’acquisition de données pour la mise enoeuvre du procédé selon l'invention, en sismique marine, - les figures 4 à 7 sont des courbes illustrantles phénomènes de distorsion fréquentielle due à lacorrection dynamique en fonction du déport, - la figure 8 est une section à temps constant obtenue selon 1'invention pour un fond marin situé à environ 140 ms, la figure 9 est une section temps dans la direction perpendiculaire (cross-line) au ' trajet du navire obtenue avec le procédé selon l’invention.

Le dispositif d'acquisition de données pour la mise en oeuvre du procédé selon l'invention (figure 3)comprend de manière conventionnelle un navire 3 remor-quant deux sources d'émission S]_ et S2, ainsi qu'ungroupe de quatre flûtes 4 à 7 par exemple. Les flûtes 4 à7 sont censées se déplacer parallèlement entre elles etsont disposées symétriquement par rapport à un axe longi-tudinal de symétrie 8 qui est sensiblement la directionde déplacement D ou de trajet du navire 3 (lignenaviguée). Les sources et S2 sont également disposéessymétriquement par rapport à l'axe longitudinal 8. Chaqueflûte (streamer en anglais) 4 à 7 comprend un grandnombre de récepteurs-enregistreurs qui reçoivent etenregistrent, chacun, des signaux ou traces qui sontreprésentatives des ondes émises par les sources d'émis-sion Slz S2 et réfléchies par les réflecteurs du sous-sol 10 marin. Une ligne d'acquisition de données sismiques estdéterminée par la ligne de déplacement d'une flûte surlaquelle sont régulièrement espacés les récepteurs. Entête de chaque flûte, il y a un récepteur de tête R14, 5 κΛς, R16, R3^, l'indice en exposant indiquant le rang durécepteur et les indices 4, 5, 6, 7 indiquant la flûte considérée. Dans l'exemple représenté, le récepteur detête R14 de la flûte 4 est séparé de la source Sx par unedistance ou déport égale à 225 m, le récepteur de tête 10 R3^ étant à 170 m de Slf les autres récepteurs de tête R16 et étant respectivement séparés de la source de 151 m et 190 m. En raison de la symétrie dudispositif, les récepteurs de tête R14, R1^ R16, R17 sont séparés de la source S2 respectivement de 190 m, 151 15 m, 170 m et 225 m. L'écartement des sources S]_ et S2, tout commel'écartement entre les flûtes 4 à 7, est obtenu à l'aidede panneaux divergents angulés appropriés de manière àobtenir une poussée hydrodynamique latérale adéquate. Le 20 positionnement des sources et des récepteurs est obtenud'abord en relatif par rapport au navire puis en absoluen utilisant la position géographique de ce dernier,déterminée par exemple par un système de positionnementpar satellites GPS (Global Positionning Systems), en mode 25 différentiel. Tous ces moyens ainsi que ceux nécessairespour une détermination correcte de la localisationspatiale des flûtes, récepteurs et sources sont bienconnus des spécialistes et ne seront donc pas décrits.

Lorsqu'on utilise les techniques classiques de 30 traitement sismique marine avec éventuellement une som-mation de traces, bien connues et non décrites ici, onobtient une section temps dans la direction cross-lineanalogue à celle représentée sur la figure 2 sur laquelleon constate que les traces enregistrées par exemple à 92 35 ms, 124 ms et 172 ms sont peu distinctes et, en 11 conséquence, il est difficile de déterminer l'origine des artefacts dont une des causes, comme précisé ci-dessus, est la distorsion fréquentielle.

Dans la réalisation d'un cube 3D en tracesproches dans un système d’axes (X, Y, T), on attribue àchaque cellule ou bin d'un ensemble de cellules du planX, Y, représentatif du plan d'acquisition des donnéessismiques, une trace centrale représentée au centre de lacellule considérée suivant l'axe temporel T.

Selon l'invention, on atténue les artefacts enchoisissant en entrée du traitement sismique en vue del'obtention du cube traces proches des traces à déportconstant. Pour cela, on prend une première source, parexemple, et on choisit comme déport constant le déport de225 m, distance qui sépare la source S^ du récepteur detête R14 et qui représente le plus petit déport communentre la source et les divers récepteurs des flûtes 4à 7 ou encore le plus grand des déports courts. Si ontrace un arc de cercle 12 ayant comme centre et S1R14comme rayon, l'arc de cercle 12 coupe théoriquement la flûte 5 en un point qui est proche ou confondu avec lerécepteur R55, la flûte 6 au voisinage du récepteur R66et la flûte 7 au voisinage du récepteur R47. Dans ce cas,l'obtention de l'image du cube traces proches selonl'invention passera par exemple par la sélection desenregistrements associés aux couples S1-R14, S1-R55, S^-

R66 et S1-R47 au fur et à mesure de l'acquisition, puispar le traitement de ces données, les bins remplis à cestade seront ceux contenant les milieux des segmentsjoignant les couples Si_-récepteur considéré R14, R55, R66et R47 par exemple les milieux PM et PMj_ des segments S-^-R17 et S^-R47 et les lignes de traitement passant par cesbins, telles que la ligne 20 passant par milieux PM etPMV 12 .j d ,1 ί» ο y υ ι υ ζ ο k

On recommence l'opération avec la source S2 et le déport de 225 m déterminé par la distance et représentant le plus petit déport commun entre la sourceS2 et les divers récepteurs des flûtes 4 à 7, ce quidonne un arc de cercle 13 qui intersecte la flûte 4 auvoisinage du récepteur R44, la flûte 5 au voisinage durécepteur R^g, la flûte 6 au voisinage du récepteur R^get la flûte 7 au récepteur R3^. Sur la trace uniquesélectionnée pour chaque couple source-flûte, on effectueun traitement sismique qui se ramène essentiellement àune correction dynamique.

Avec le procédé selon l'invention, on constateque la variation de la distorsion fréquentielle estdiminuée d'un facteur 6 environ. En effet, la distanceentre deux récepteurs consécutifs sur chaque flûte estgénéralement de 12,5 m, qu'on désigne également parinter-trace. Comme la variation du déport qui est àl'origine de 75 m (225-150 m) est ramenée au maximum àl'inter-trace, soit 12,5 m, il est aisé de constater quele gain est de 75/12,5 = 6.

Une comparaison des figures 8 et 9 qui sont dessections du cube 3D traces proches selon l'invention,avec les figures 1 et 2 de l'art antérieur, permet devérifier que les bandes 15 des artefacts (figure 8) sontconsidérablement élargies et que les premières couches dusous-sol sont mieux restituées sur le cube 3D tracesproches obtenu selon l'invention même si les sectionsreprésentées sur les figures ne concernent pas les mêmesmilieux explorés et sont à des temps constants dif-férents. La figure 9 met en évidence que les traces sontplus lisibles et plus nettes les unes des autres, ce quipermet éventuellement d'identifier plus facilement deserreurs de positionnement.

Sur les figures 4 à 7, on a représenté des courbes de distorsion fréquentielle, dont la variation 13 est sensiblement linéaire en fonction du déport pour la gamme de valeurs (déport, temps, vitesse) considérée.

La figure 4 représente une courbe F'/F enfonction du déport avec To = 92 ms, un déport maximum de250 m et une vitesse de correction dynamique NMO (NormalMove Out) de 1500 m/s. La variation du rapport F'/F estcomprise entre 0,525 et 0,675 soit 0,15 pour unevariation du déport comprise entre 150 et 225. Le procédéselon l'invention permet de réduire la variation durapport F'/F à une valeur sensiblement égale à 0,025 soitun gain d'environ 6.

La figure 5 représente la courbe F'/F en fonctiondu déport avec To = 92 ms, vitesse NMO — 1600 m/s et undéport maximum de 300 m, alors que la figure 6 représenteune courbe similaire avec To = 124 ms, vitesse NMO : 1600m/s et un déport maximum de 3 00 m ; la courbe de lafigure 7 est réalisée avec To = 172 ms, vitesse NMO =1700 m/s et un déport maximum de 300 m. Dans tous les caset dans la partie linéaire de chaque courbe F’/F, onconstate un gain d'environ 6.

Le procédé selon l'invention est également appli-cable en couverture multiple réduite. A cet effet, àchaque bin est affecté outre une première trace telle quecelle associée aux points milieux PM ou PM^, au moins unedeuxième trace supplémentaire permettant de constituerladite couverture multiple réduite. Chaque trace supplé-mentaire peut être soit une trace choisie comme lapremière trace mais avec un déport constant différent,soit une trace consécutive à ladite première trace. Dansce cas, les traces ainsi sélectionnées y compris lapremière trace sont corrigées de dynamique puis somméespour obtenir une trace somme qui est utilisée comme tracecentrale de 1a. cellule (bin) dans laquelle sont situésles points milieux associés aux traces sélectionnées. 14 » La tendance actuelle est à l'augmentation du L nombre de flûtes ; cela augmente en conséquence les | différences entre déports courts d'une même ligne navi- | guée et aggrave potentiellement la distorsion fréquen- 5 tielle différentielle affectant les traces à déportconstant. Avec le procédé selon l'invention et quel quesoit le mode de couverture simple ou multiple, l'augmen-tation du nombre de flûtes n'aggrave pratiquement pas ladistorsion fréquentielle différentielle du fait qu'on 10 travaille avec un déport constant ou quasi constant encouverture simple et une gamme de déports constants encouverture multiple.

Claims (4)

15 REVENDICATIONS

1. Procédé de réalisation d'un cube 3D en traces proches à partir de données acquises en sismique marineréflexion, du type dans lequel : a) on procède à l'acquisition desdites données endéplaçant dans une direction de référence D, i) au moins une source (S-j_, S2), ii) plusieurs flûtes sismiques (4 à 7) espacées dela source, distantes latéralement les unes desautres et comportant chacune des récepteurs-enregistreurs (1Ο4, 1Ο5, R1^, R^7 · · · ) régulière-ment espacés le long de chaque flûte, et iii) des moyens de localisation spatiale de chaquesource et de chaque flûte, de manière à obtenir simultanément, pour chaque lignenaviguée, d'une part des données sismiques sous la formede traces, fournies par les récepteurs-enregistreurs, lelong d'autant de lignes d'acquisition, de préférencesensiblement parallèles à la direction de référence D,que le produit du nombre de flûtes par le nombre desources utilisés et d'autre part, des données de locali-sation spatiales de chaque source et des récepteurs-enregistreurs pour fournir la distance ou déport séparantchaque source de chaque récepteur ainsi que les coordon-nées du point milieu du couple source/récepteur-enregistreur associé à chaque trace, b) on réalise une image du sous-sol en trois dimensionsdans un système d'axes (X, Y, T) sous la forme dite cube3D en attribuant à chaque cellule (bin) d'un ensemble decellules formant un maillage spatial régulier dans unplan X, Y, représentatif du plan d'acquisition desdonnées sismiques, une trace centrale élaborée à partir 16 des traces enregistrées et représentée au centre deladite cellule suivant l'axe temporel (T), ledit procédé se caractérisant en ce que pour chaqueligne naviguée on constitue la trace centrale de chaque 5 cellule à partir d'au moins une première traceenregistrée dont le point milieu associé est situé àl'intérieur de ladite cellule et dont le déport est leplus proche possible d'un déport constant égal au pluspetit déport commun à au moins l'ensemble des traces de 10 la ligne naviguée.

2. Procédé selon la revendication 1, caractérisé en ce que le plus petit déport commun est égal au plus grand(S-[_ R17 ou S1-R14) des déports courts pour une ligne naviguée.

3. Procédé selon la revendication 1 ou 2, carac-térisé en ce que ladite trace centrale consiste uniquement en ladite première trace de manière à obtenirun cube 3D traces proches en couverture unité.'

4. Procédé selon la revendication 1 ou 2, carac- 20 térisé en ce que, pour réaliser une couverture multipleréduite, on affecte à chaque cellule, en plus de laditepremière trace, au moins une trace supplémentaire qui estsoit une trace dont le déport est constant mais différentde celui de ladite première trace, soit une trace 25 consécutive à ladite première trace, on corrige de dyna-mique les traces ainsi sélectionnées y compris lapremière, puis on somme les traces corrigées pour obtenirune trace somme, qui est utilisée comme trace centrale dela cellule (bin) dans laquelle sont situés les points 30 milieux associés aux traces sélectionnées.