Brevet FR2759785A1 — publié le 19980821

=2759785 Méthode de traitement de signaux sismiques pour réduire les effets de pendage dans l'analyse de courbures d'indicatrices La présente invention concerne une méthode de traitement de signaux sismiques pour réduire les effets de pendage dans l'analyse de courbures d'indicatrices. Elle se rapporte plus particulièrement à une méthode de correction du type de celle qui est désignée par les spécialistes o10 comme étant une correction DMO (DIP MOVE OUT) dont la propriété essentielle est de rendre les courbures des indicatrices indépendantes du pendage en les ramenant à celles des pendages nuls. Depuis 1978, de nombreuses méthodes dites de correction DMO ont été proposées, le document rédigé par DEREGOWSKI en 1986 et intitulé "What is DMO ?" (FIRST BREAK, vol. 04, July 1986, pages 07-24) constituant un résumé des opérateurs DMO utilisés à cette date. En effet, le traitement sismique en couverture multiple nécessite une sommation des traces enregistrées à l'aide d'un ensemble émetteurs- récepteurs, et groupées en une collection de traces en point milieu commun (PMC). Pour effectuer cette sommation, il faut corriger l'obliquité des traces avant de les sommer. Classiquement, on effectue d'abord une correction dynamique appelée NMO (NORMAL MOVE OUT) suivie d'une autre correction dynamique dite DMO. Pour s'en convaincre, il suffit de se reporter à l'étude faite par DEREGOWSKI, et de remarquer que l'écart temps At DMO a pour valeur une différence de trajet qui est dépendante du temps, du déport et de l'angle de pendage. C'est la raison pour laquelle la correction DMO est faite après la correction dynamique NMO. De plus, il faut noter que les corrections d'obliquité ne sont valides que pour des séries de pendages monotones et qu'en cas de conflits de pendages, comme par exemple en tectonique complexe, alors les opérateurs DMO connus ne peuvent s'appliquer. D'autres auteurs ont tenté d'étendre le DMO à des cas de champs de vitesses dont la valeur varie uniquement en fonction de la profondeur, ce qui revient à considérer que les courbes ou les surfaces de vitesse constante sont des droites ou des plans horizontaux. D'autres propositions ont été faites, notamment par FOREL et GARDNER, ("A three dimensional perspective on two dimensional move out", GEOPHYSICS vol. 53, 1988, pages 604-610), qui consistent à effectuer le DMO avant le NMO. Une telle proposition n'est pas satisfaisante car elle prend en considération des données approximatives relatives au déport h et au temps t de manière à compenser la différence non linéaire qui existe entre le cheminement NMO (h) suivi par DMO (h) et le cheminement DMO (h) suivi par NMO (h). Pour atteindre un même point par les deux chemins, FOREL et GARDNER définissent un pseudo-déport h' et un temps de remplacement t1 qui est substitué au temps to du DMO classique. Cette procédure permet de faire des analyses de vitesses après application du DMO, les analyses de vitesses étant effectuées habituellement tous les 500 mètres environ. Une interpolation bilinéaire (en x et t) est effectuée ensuite entre ces analyses de vitesses. I1 est aisé de comprendre que l'énorme pas d'échantillonnage spatial choisi (environ 500 mètres) a pour conséquence que la vitesse est toujours spatialement aliasée et que le redressement des traces distales est alors imparfait entre les points de mesure. o De plus, toutes les méthodes précédentes nécessitent une connaissance a priori du champ de vitesses de sommation (stack) corrigé de DMO ce qui n'est jamais le cas et oblige les techniciens à procéder par itérations successives. Par ailleurs, pour l'analyse de courbures d'indicatrices, Eric de BAZELAIRE a introduit en 1988 la notion de temps (tp) de profondeur de mise au point (NMO revisited: inhomogeneous media and curved inter- faces: GEOPHYSICS, vol. 53, n 2, pages 143-157) qui présente l'avantage d'être indépendant du temps to. La méthode d'analyse de vitesses développée par Eric de BAZELAIRE est utilisée pour la réalisation de sections sommes améliorées, appelées sections "POLYSTACK" et elle consiste à construire un document appelé "bap" qui est associé à une collection en point milieu commun de traces sismiques enregistrées. La "bap" associée à une collection de traces en point milieu commun (PMC) est une autre collection de traces dont chacune est le résultat de la sommation des traces de la collection PMC, après application à 3 2759785 chacune de ces traces d'une correction de type statique (indépendante du temps) et différente d'une trace à l'autre. Une définition plus détaillée du procédé classique POLYSTACK et du calcul des tp est par exemple donnée dans une autre demande de brevet de la demanderesse déposée le 19 octobre 1994, publiée sous le n 2 726 091 et concernant une autre méthode originale appelée DELTASTACK. La présente invention a pour but de proposer l'application d'un nouvel opérateur DMO dans laquelle le DMO utilisé est une correction statique faisant intervenir notamment le temps de profondeur de mise au point tp et surtout une méthode dans laquelle les deux opérateurs NMO et DMO spécifiques sont statiques et commutables. Un autre avantage de la présente invention réside dans le fait qu'elle ne demande aucune connaissance a priori du champ de vitesses de sommation. Un autre avantage de la présente invention est qu'elle peut être appliquée tant en tectonique calme qu'en tectonique complexe, ce qui n'avait pu être obtenu jusqu'à présent. La méthode selon l'invention est caractérisée en ce qq'elle consiste à effectuer les étapes suivantes: a - émettre à partir d'au moins une source d'émission des ondes acoustiques dans un milieu, b - recevoir sur des récepteurs des ondes réfléchies dans ledit milieu et enregistrer lesdites ondes sous forme de traces sismiques, chaque trace sismique étant au moins repérée par un déport (h) et un point milieu (x0) qui sont liés au couple source-récepteur qui a généré ladite trace sismique, deux points milieux consécutifs x0, pour un même déport h, étant séparés par un écart constant Ax, c - définir une suite de N indices entiers j compris entre deux limites (jmin, jmax) et représentatifs des courbures des indicatrices d'au 3 0 moins une collection de traces prédéterminée, d - calculer un temps de profondeur de mise au point tp pour chaque valeur de j, e - définir une suite de P nombres entiers u déterminant un intervalle (P- 1)Ax, f recopier, pour chaque valeur de j, chaque trace identifiée par ses paramètres x0 et h, autant de fois qu'il y a de nombres entiers u dans l'intervalle (P-1)Ax, g - appliquer sur chacune desdites traces recopiées au moins une correction statique (dt) dépendant au moins des paramètres j, h et u Ax, h - recommencer les étapes f et g pour toutes les valeurs de x0 et dej. Selon une caractéristique de l'invention, la méthode est caractérisée en ce que la collection de traces prédéterminée est constituée par un ensemble de collections en point milieu commun (PMC), et en ce qu'elle consiste en outre à: - ordonner les collections en point milieu commun (PMC) suivant les valeurs croissantes de x0 et les traces de chaque collection en point milieu commun (PMC) suivant les valeurs croissantes de h, puis à appliquer une autre correction statique (At), dite NMO PSCAN, sur chacune des traces de chaque collection en point milieu commun, avant l'application de ladite correction statique (dt). Selon une caractéristique de l'invention, la méthode est caractérisée en ce que la collection de traces prédéterminée est constituée par un ensemble de collections (H) à déport constant h, les traces de chaque collection (H) à déport constant étant ordonnées suivant les valeurs croissantes de x0, les collections (H) étant ordonnées suivant les valeurs croissantes de h, et en ce qu'elle consiste à effectuer les étapes f et g pour toutes les traces, lorsqu'elles existent, de chaque collection (H). Selon une caractéristique de l'invention, la méthode est caractérisée en ce qu'elle consiste en outre à réaliser pour chaque valeur de j, un bloc comprenant les collections en point milieu commun (PMC) ordonnées suivant les valeurs croissantes de x0 et les collections (H) à déport constant ordonnées suivant les valeurs croissantes du déport h, ladite correction statique NMO PSCAN At est appliquée pour toutes les valeurs de j et à toutes les traces de toutes les collections en point milieu commun. Selon une caractéristique de l'invention, la méthode est caractérisée en ce qu'elle consiste en outre et dans chaque collection (H) à déport constant, à: - sélectionner dans l'intervalle (P-1)Ax une position d'abscisse xu pour laquelle il existe une trace sismique, - positionner la trace x0 au point xu, - appliquer à ladite trace positionnée en xu ladite correction (dt), - sommer la trace corrigée de statique (dt), échantillon par échantillon, à la trace existant au point d'abscisse xu. Selon une caractéristique de l'invention, la méthode est caractérisée en ce que les opérations effectuées sur les traces x0 et xu sont répétées pour toutes les valeurs de u associées à l'intervalle (P-1)Ax, puis pour toutes les valeurs x0 de chaque collection (H) et pour toutes les collections (H) de chaque bloc j. Selon une caractéristique de l'invention, la méthode est caractérisée en ce qu'elle consiste en outre à réaliser des plans ITP dont chacun est caractérisé par un indice j et composé de traces localisées aux différentes valeurs de x0 pour lesquelles il existe une trace sismique, lesdites traces étant le résultat des deux corrections statiques (At et dt), de façon à obtenir un plan ITP corrigé de migration à déport nul (MZO). Selon une caractéristique de l'invention, la méthode est caractérisée en ce que les plans ITP corrigés de MZO sont rassemblés dans un bloc d'axes x, t etj pour constituer un bloc ITP corrigé de MZO. Selon une caractéristique de l'invention, la méthode est caractérisée en ce que la correction statique NMO PSCAN (At) est effectuée avant la correction statique (dt) dite DMO PSCAN. Selon une caractéristique de l'invention, la méthode est caractérisée en ce que la correction statique (dt) dite DMO PSCAN est effectuée avant la correction NMO PSCAN (At). Selon une caractéristique de l'invention, la méthode est caractérisée en ce que la collection de traces prédéterminée est une collection point de tir de traces élémentaires et en ce qu'elle consiste à effectuer les 3 0 étapes suivantes: a) identifier chaque trace élémentaire par l'abscisse x0 du point milieu correspondant, le déport h, la valeur u dans l'intervalle (P-l)Ax et un indice j, 6 2759785 b) réaliser, pour chaque valeur de j, un plan ITP(j) composé de traces nulles localisées à chaque valeur de x0 pour laquelle il existe une trace sismique, c) sélectionner une trace élémentaire d'abscisse x0 et la recopier en un point d'abscisse XU pour lequel il existe une trace sismique dans le plan ITP(j), d) appliquer une correction statique MZO (A't) sur ladite trace recopiée, ladite correction statique MZO étant la somme de deux corrections statiques NMO PSCAN (At) et DMO PSCAN (dt), la correction statique NMO PSCAN (At) dépendant seulement des paramètres j et h et la correc- tion statique DMO PSCAN (dt) dépendant des paramètres j, h, u et Ax, e) sommer ladite trace recopiée et corrigée de MZO, à la trace existant au point d'abscisse xu, f) recommencer les étapes précédentes pour toutes les valeurs de u, toutes les valeurs de h et toutes les valeurs de j, g) effectuer les étapes b) à f) pour chacune des traces de ladite collection point de tir. Selon une caractéristique de l'invention, la méthode est caractérisée en ce que la correction statique NMO PSCAN (At) est donnée 2 o par la formule At(j,h) = - tp (j) + tp (j)2 + (h / Vo)2 dans laquelle V0 est la vitesse d'entrée dans le milieu. Selon une caractéristique de l'invention, la méthode est caractérisée en ce que la correction statique DMO PSCAN (dt) est donnée par la formule: dt(j,h,u)=-tp(j) 1 - 1 h 2 Selon une caractéristique de l'invention, la méthode est caractérisée en ce que la correction statique globale MZO (A't) est donnée par la formule A't(j,h,u)= t2(j)±- -t (j) 1- (h) V, h' dans laquelle V0 est la vitesse d'entrée dans le milieu. D'autres avantages et caractéristiques ressortiront mieux à la lecture de plusieurs modes de réalisation de l'invention, ainsi que des dessins annexés sur lesquels: - la figure 1 est une représentation schématique d'un bloc de plans, chaque plan représentant une collection de traces en point milieu commun (PMC). - la figure 2 est une représentation de blocs pour diverses valeurs de j telles que Jmin, j et jmax, chaque bloc comprenant des plans représentant des collections de traces un PMC et des plans H représentant des collections de traces triées à déport constant ou iso-déport. - la figure 3 est une représentation schématique d'un plan H avec un opérateur DMO PSCAN. - la figure 4 est une représentation schématique d'un plan ITP pour un j donné. - la figure 5 est une représentation schématique d'un bloc de plans bap et ITP. Selon une première variante de la méthode de la présente invention, on constitue un ensemble de traces sismiques qui sont obtenues à l'aide d'un dispositif comportant une ou plusieurs sources d'émission, désignées comme étant des points de tir, et des récepteurs-enregistreurs, les sources et les récepteurs étant disposés à la surface d'un milieu qu'on veut explorer. On dispose donc d'un ensemble de traces sismiques qui sont échantillonnées en temps avec un pas d'échantillonnage ôt et groupées en plusieurs collections point de tir. Chaque collection point de tir correspond à l'ensemble des traces enregistrées pour une émission (tir) donnée. Chaque trace sismique est identifiée par plusieurs paramètres tels que les coordonnées des positions spatiales du couple point de tir ou source d'émission (xs, Ys) et du récepteur (xr, Yr), ayant généré ladite trace sismique. 8 2759785 L'identification comprend également les valeurs calculées (xo, Yo) des coordonnées du point milieu (PM) et du déport (h) associés audit couple, à partir des formules suivantes bien connues: Xs+ Xr Ys +Yr x 2 Y 2 (1) h = (Xr,-xs)2 + (yr-y)2 En sismique à deux dimensions (2D), on admet que Yo reste constant tout le long de la ligne sismique et que deux valeurs consécutives de x0 sont séparées d'une valeur constante Ax appelée distance inter PMC pour des collections en point milieu commun. Dans une première étape, on extrait des valeurs précédentes la valeur du déport maximum hmax. On mesure à l'aide des méthodes connues la bande passante moyenne B des signaux enregistrés sur les traces sismiques. On définit une suite de N indices entiers j à partir desquels on calcule les valeurs de tp qui vont être utilisées pour constituer le bloc de baps. Ces valeurs sont obtenues en appliquant la formule t (j)! (hm'B kj (2) 2 V02 kj B Les indices j (indices tp) sont des entiers positifs pour les valeurs de tp positives et négatifs pour des valeurs de tp négatives. L'indice-tp nul correspond au tp infini. Le nombre d'indices-tp est limité par l'opérateur. L'indice jmax correspond à la courbure la plus élevée des indicatrices qui sont balayées et il est choisi de manière à couvrir toutes les courbures des indicatrices supposées dans le milieu à explorer. L'indice Jmin est au plus égal à -jmax. Le paramètre k est une constante définissant le pas d'échantillonnage en courbure; il a classiquement une valeur comprise entre 0,5 et 0,7. Enfin V0 est une constante représentant la vitesse en m/s du milieu dans lequel on réalise la mesure des courbures en termes de temps de profondeur de mise-au-point (la profondeur de mise-au-point p est égale à l 5 v, t 9 2759785 Dans une deuxième étape et selon un premier mode de réalisation, on groupe les collections de traces point de tir suivant une première série de collections PMC qui sont ordonnées suivant les valeurs croissantes des points-milieux x0 sur l'axe x, chaque collection PMC étant située dans un plan d'abscisse x0, x1....., un seul plan PMC 1 d'abscisse x0 étant représenté sur la figure 2. Les traces de chaque plan PMC sont ordonnées suivant les valeurs croissantes de h. Dans une troisième étape, on applique aux traces des collections PMC une correction statique (NMO PSCAN) de la manière suivante a) on duplique N fois chaque collection PMC, b) aux traces d'une collection PMC dupliquée correspondant à une valeur de j, on applique une correction statique At, dite NMO PSCAN, définie par la formule: At(jh) = - tp(j)+ /tp(j)2 +(h/V0)2 c) on répète la sous-étape b) pour toutes les valeurs j et pour toutes les collections PMC. Pour chaque valeur de j, on dispose d'un bloc j de collections 2 o PMC dont les traces sont corrigées de NMO PSCAN. Comme en général le calcul de la correction statique At(j, h) ne donne pas un nombre entier de fois le pas d'échantillonnage ât, les échantillons sont interpolés linéairement. Dans une quatrième étape, et pour chaque bloc j, on trie les traces à déport constant de manière à obtenir, pour chaque valeur de h, une collection iso-déport située dans un plan H, parallèle au plan xot. Les plans H sont ordonnés suivant les valeurs croissantes de h et les traces de chaque plan H sont ordonnées suivant les valeurs croissantes de x0. Sur la figure 2, on a représenté trois blocs 2, 3, 4 correspondant aux valeurs Jmin, j et jmax, chacun des blocs comprenant un nombre n de plans H dont quelques uns sont représentés et un nombre q de plans PMC dont un seul est représenté. Dans le plan PMC du bloc 4 correspondant à Jmax, on a représenté la correction statique NMO PSCAN pour la valeur Jmax et pour une valeur h. Le résultat de la correction est schématisé par une indicatrice sur-corrigée 7. Dans le bloc 3 correspondant à l'indice j, le résultat de la correction statique NMO PSCAN est schématisé par une indicatrice corrigée horizontale 6 tandis que dans le bloc 2 correspondant à imin, le résultat de la correction statique NMO PSCAN est schématisé par une indicatrice sous-corrigée 5. Une comparaison des indicatrices corrigées 5, 6, 7 avec l'indicatrice initiale 8 du bloc de départ 9 permet de constater que la bonne correction est celle représentée sur le bloc 3 (indicatrice horizontale). Pour chaque plan élémentaire H tel que représenté sur la figure 3 caractérisé par son déport h, on calcule une valeur limite xmax, à l'aide de la formule h2 Xmax 2 V t p (j)2 + h(4) Chaque plan élémentaire H contient toutes les traces positionnées aux différentes valeurs x0 définissant les points milieux PM qui sont ordonnées suivant les valeurs croissantes de x0. On détermine pour chacune desdites valeurs x0 définissant une position PMC dudit plan élémentaire H, une suite de P nombres entiers u tels que uAx soit compris entre (x0 - xmax) et (x0 + xmax). On ne prend pas en considération le résidu inférieur à Ax aux deux extrémités. Cela donne la formule suivante: xo- x,., trace du plan H existant à ladite abscisse xu. iii) Cette opération est répétée pour chaque valeur x0 du plan H, ce qui implique de déplacer le gabarit 8 et de le centrer sur la nouvelle valeur de x0. Dans une étape suivante et pour chaque valeur de j, on réalise un plan ITP(j) caractérisé par son indice j et composé de traces nulles localisées aux différentes valeurs de x0 pour lesquelles il existe une trace sismique. Le résultat de la sommation pour toutes les valeurs de x0 des étapes ii et iii, qui correspond à autant de traces sommes qu'il y a de valeurs x0 dans le plan H, est reporté dans le plan ITP(j) aux points d'abscisses correspondantes. Comme en général le calcul de la correction statique dt(j, h, u) ne donne pas un nombre entier de fois le pas d'échantillonnage 5t, les échantillons sont interpolés linéairement avant d'être sommés. De la sorte, on obtient des traces somme corrigées de DMO PSCAN. Dans le plan ITP(j) de la figure 4, on a représenté plusieurs positions du gabarit 8 d'un même plan H, par exemple celui de la figure 3, chacune des positions étant centrée sur une valeur différente de x0, mais passant toutefois par un même point de coordonnées x0, to. On recommence les opérations précédentes pour tous les plans H dublocj. Le résultat de cette série d'opérations est l'obtention d'un plan ITP caractérisé par son indice j et corrigé de migration à déport nul MZO (Migration Zero Offset), puisque chaque trace a subi deux corrections statiques, l'une à l'aide de l'opérateur NMO PSCAN At(j, h) et l'autre à l'aide de l'opérateur DMO PSCAN dt(j, h, u). Comme les deux corrections définies par les formules (3) et (6) sont indépendantes du temps, nous sommes bien en présence de deux corrections statiques. Dans des étapes suivantes, on réalise des plans ITP corrigés de MZO, pour toutes les valeurs de j. L'ensemble de plans ITP corrigés chacun de MZO constitue alors un bloc bap corrigé de MZO, dont les axes sont x, t et j (figure 5). Une section de ce bloc par un plan à x0 constant constitue une bap élémentaire que l'on peut introduire dans un traitement du type POLYSTACK à la place des baps classiques. Dans un deuxième mode de réalisation, on effectue un tri des traces enregistrées à déport constant h de manière à obtenir n collections de traces iso-déport, chaque collection étant située dans un plan H parallèle au plan xot. Les opérations de correction statique DMO PSCAN dt de la formule (6) décrites précédemment sont effectuées sur les traces des plans H de chaque bloc j. Les traces de chaque bloc j, corrigées de DMO PSCAN, sont réorganisées en collections PMC dans chacune desquelles les traces PMC sont ordonnées suivant les valeurs croissantes de h. Dans chaque collection PMC on applique sur chacune des traces une correction NMO PSCAN (At) telle que décrite précédemment; on somme les traces corrigées de NMO PSCAN pour obtenir au point d'abscisse x0 une trace somme corrigée de MZO. La trace somme est positionnée dans le plan ITP de même indice j que celui de la collection traitée de sorte qu'on obtienne un plan ITP () corrigé de MZO. On répète ces opérations pour toutes les valeurs de j de manière à 2 o obtenir le bloc de la figure 5. Suivant un troisième mode de réalisation plus simplifié et donc plus facile à mettre en oeuvre et qui nécessite moins de temps de calcul et moins d'opérations distinctes, la méthode selon l'invention consiste dans une étape initiale, à l'aide des formules précédentes, à calculer les coordonnées (Xo, yo) du point milieu commun, du déport h (formule 1), la valeur du déport maximum hmax, les valeurs de tp(j) (formule 2), les valeurs limites Xmax (formule 4), l'intervalle (P-l)Ax et la suite des P nombres entiers u. Dans une autre étape suivant ce troisième mode de réalisation, on crée directement et pour chaque valeur de j un plan nommé ITP(j) composé de traces nulles localisées à chaque valeur de x0 pour laquelle on dispose d'une trace sismique. On trie ces N plans suivant les valeurs croissantes de j pour obtenir ainsi un bloc bap dont les axes sont x, t et j. Dans une étape suivante, on prend une trace élémentaire issue des collections point de tir. Cette trace est caractérisée par ses paramètres x0 et h. Pour chaque valeur de j, pour chaque valeur de h et pour chaque valeur de u définies ci-dessus, on renouvelle l'opération suivante: on recopie ladite trace élémentaire, on la place au point d'abscisse xu dans le plan ITP(j), on lui fait subir une correction statique de valeur A't définie par l'équation ci-dessous et on la somme échantillon par échantillon à la trace déjà présente dans le plan ITP(j) et d'abscisse xu. L'équation de cette correction statique dite MZO PSCAN est la suivante , h 5A't(j,h,u)= p(j)±-o-t p(j) 21 h2 (7) VO h Comme en général le calcul de la correction statique S't(xo, j, h, u) ne donne pas un nombre entier de fois le pas d'échantillonnage Et, les échantillons sont interpolés linéairement avant d'être sommés. L'opération est renouvelée pour toutes les P valeurs de u, pour toutes les n valeurs de h et pour toutes les N valeurs de j. Une fois que ceci est accompli, on passe à la trace suivante. Le traitement est donc monocanal. Quand on a analysé toutes les traces incidentes, on obtient un bloc bap corrigé de MZO dont les axes sont x, t et j (figure 5). Une section de ce bloc par un plan à x0 constant constitue une bap élémentaire que l'on peut introduire dans un traitement de type POLYSTACK à la place des baps classiques. De ce qui précède, il est remarquable de noter qu'on peut effectuer une correction statique DMO PSCAN avant ou après la correction statique dite NMO PSCAN ou encore conjuguer les deux corrections statiques pour ne faire qu'une correction statique directement sur les traces de départ ou d'origine de la collection point de tir. Une fois le bloc de baps corrigé de DMO PSCAN obtenu, l'analyse automatique de courbures peut être effectuée qui permettra ensuite 3 0 de définir le champ de vitesses corrigé de DMO PSCAN. REVENDICATIONS 1. Méthode de traitement des signaux sismiques pour l'analyse de courbures des indicatrices caractérisée en ce qu'elle consiste à effectuer les étapes suivantes: a - émettre à partir d'au moins une source d'émission des ondes acoustiques dans un milieu, b - recevoir sur des récepteurs des ondes réfléchies dans ledit milieu et enregistrer lesdites ondes sous forme de traces sismiques, chaque trace sismique étant au moins repérée par un déport (h) et un point milieu (x0) qui sont liés au couple source-récepteur qui a généré ladite trace sismique, deux points milieux consécutifs x0 pour un même déport h, étant séparés par un écart constant Ax, c - définir une suite de N indices entiers j compris entre deux limites (jmin, jmax) et représentatifs des courbures des indicatrices d'au moins une collection de traces prédéterminée, d - calculer un temps de profondeur de mise au point tp pour chaque valeur de j, e - définir une suite de P nombres entiers u déterminant un intervalle (P- 1)Ax, f recopier, pour chaque valeur de j, chaque trace identifiée par ses paramètres x0 et h, autant de fois qu'il y a de nombres entiers u dans l'intervalle (P-l)Ax, g - appliquer sur chacune desdites traces recopiées au moins une correction statique (dt) dépendant au moins des paramètres j, h et u Ax, h - recommencer les étapes f et g pour toutes les valeurs de x0 et dej. o 2. Méthode selon la revendication 1, caractérisée en ce que la collection de traces prédéterminée est constituée par un ensemble de collections en point milieu commun (PMC), et en ce qu'elle consiste en outre àa: - ordonner les collections en point milieu commun (PMC) suivant les valeurs croissantes de x0 et les traces de chaque collection en point milieu commun (PMC) suivant les valeurs croissantes de h, puis à appliquer une autre correction statique (At), dite NMO PSCAN, sur chacune des traces de chaque collection en point milieu commun, avant l'application de ladite correction statique (dt). 3. Méthode selon la revendication 1, caractérisée en ce que la collection de traces prédéterminée est constituée par un ensemble de collections (H) à déport constant h, les traces de chaque collection (H) à déport constant étant ordonnées suivant les valeurs croissantes de x0, les collections (H) étant ordonnées suivant les valeurs croissantes de h, et en ce qu'elle consiste à effectuer les étapes f et g pour toutes les traces, lorsqu'elles existent, de chaque collection (H). 4. Méthode selon les revendications 1 à 3, caractérisée en ce qu'elle consiste en outre à réaliser pour chaque valeur de j, un bloc comprenant les collections en point milieu commun (PMC) ordonnées suivant les valeurs croissantes de xo et les collections (H) à déport constant ordonnées suivant les valeurs croissantes du déport h, ladite correction statique NMO PSCAN At est appliquée pour toutes les valeurs de j et à toutes les traces de toutes les collections en point milieu commun. 5. Méthode selon les revendications 3 et 4, caractérisée en ce qu'elle consiste en outre et dans chaque collection (H) à déport constant, à: sélectionner dans l'intervalle (P-l)Ax une position d'abscisse xu pour laquelle il existe une trace sismique, - positionner la trace x0 au point xu, - appliquer à ladite trace positionnée en xu ladite correction (dt), sommer la trace corrigée de statique (dt), échantillon par échantillon, à la trace existant au point d'abscisse xu. 6. Méthode selon la revendication 5, caractérisée en ce que les opérations effectuées sur les traces x0 et xu sont répétées pour toutes les valeurs entières de u associées à l'intervalle (P-l)Ax, puis pour toutes les valeurs x0 de chaque collection (H) et pour toutes les collections (H) de chaque bloc j. 7. Méthode selon les revendications 1 à 6, caractérisée en ce qu'elle consiste en outre à réaliser des plans ITP dont chacun est caractérisé par un indice j et composé de traces localisées aux différentes valeurs de x0 pour lesquelles il existe une trace sismique, lesdites traces étant le résultat des deux corrections statiques (At et dt), de façon à obtenir un plan ITP corrigé de migration à déport nul (MZO). 8. Méthode selon la revendication 7, caractérisée en ce que les plans ITP corrigés de MZO sont rassemblés dans un bloc d'axes x, t et j pour constituer un bloc ITP corrigé de MZO. 9. Méthode selon l'une des revendications précédentes, caractérisée en ce que la correction statique NMO PSCAN (At) est effectuée avant la correction statique (dt) dite DMO PSCAN. 10. Méthode selon l'une des revendications 1 à 8, caractérisée en ce que la correction statique (dt) dite DMO PSCAN est effectuée avant la correction NMO PSCAN (At). 1l. Méthode selon la revendication 1, caractérisée en ce que la collection de traces prédéterminée est une collection point de tir de traces élémentaires et en ce qu'elle consiste à effectuer les étapes suivantes: a) identifier chaque trace élémentaire par l'abscisse xo du point milieu correspondant, le déport h, la valeur u dans l'intervalle (P-l)Ax et un indice j, b) réaliser, pour chaque valeur de j, un plan ITP(j) composé de traces nulles localisées à chaque valeur de x0 pour laquelle il existe une trace 2 o sismique, c) sélectionner une trace élémentaire d'abscisse x0 et la recopier en un point d'abscisse xu pour lequel il existe une trace sismique dans le plan ITP(j), d) appliquer une correction statique MZO (A't) sur ladite trace recopiée, ladite correction statique MZO étant la somme de deux corrections statiques NMO PSCAN (At) et DMO PSCAN (dt), la correction statique NMO PSCAN (At) dépendant seulement des paramètres j et h et la correc- tion statique DMO PSCAN (dt) dépendant des paramètres j, h, u et Ax, e) sommer ladite trace recopiée et corrigée de MZO, à la trace existant au point d'abscisse xu,, f) recommencer les étapes précédentes pour toutes les valeurs de u, toutes les valeurs de h et toutes les valeurs de j, g) effectuer les étapes b) à f) pour chacune des traces de ladite collection point de tir. 12. Méthode selon l'une des revendications 2, 4, 7, 9, 10 et 11, caractérisée en ce que la correction statique NMO PSCAN (At) est donnée par la formule: At(j,h)= -tp(j)+ /tp(j)2 +(h/Vo)2 dans laquelle Vo est la vitesse d'entrée dans le milieu. 13. Méthode selon l'une des revendications 2, 3, 5, 7, 9, 10 et 11, caractérisée en ce que la correction statique DMO PSCAN (dt) est donnée par la formule: dt(j,h,u)=-to(j) V h2 14. Méthode selon la revendication 11, caractérisée en ce que la correction statique globale MZO (A't) est donnée par la formule: A't(j,h,u)=t2()+ G-tp(j) - 2 dans laquelle V0 est la vitesse d'entrée dans le milieu.