FR 2735244 A1 19961213 FR 9506887 A 19950606 La présente invention concerne une méthode de traitement de données de prospection sismique relatives à une zone souterraine, permettant de synthétiser des données sismiques correspondant au déport nul à partir de données sismiques obtenues avec des déports multiples. Elle trouve son application par exemple pour la mise en oeuvre de la technique appelée par les spécialistes migration profondeur après sommation ("poststack depth migration" en anglais). I) Etat de la technique L'imagerie sismique de structures géologiques complexes notamment celles mettant en jeu des variations latérales de vitesse importantes requiert l'utilisation de techniques spécifiques: par exemple, des techniques de migration profondeur en 3D, idéalement de migration 3D avant sommation ("3D prestack depth migration" en anglais). Ces techniques ne pouvant être pleinement valorisées tant que le modèle de vitesse n'est pas suffisamment précis, une phase importante de l'élaboration du modèle de vitesse repose sur l'utilisation de l'imagerie profondeur 3D après sommation, dite "poststack", par exemple pour la mise en oeuvre d'une technique de "mise en place successive des couches" telle que décrite par exemple dans: D.G. Rockliff, et al.; "3-D depth migration in the Diapir province of the central North Sea", 62th Ann. Internat. Mtg., Soc. Expl. Geophys., Expanded Abstracts, 925-926. Différentes démarches (cf. Fig. 1) sont utilisées classiquement pour l'élaboration du modèle de vitesse par utilisation de l'imagerie profondeur 3D après sommation. On peut distinguer deux familles de procédures pour une telle élaboration. 1) Les techniques cinématiques utilisant donc les temps 3 0 d'arrivée des réflexions: par exemple la migration de cartes, l'inversion des temps d'arrivée (tomographie de réflexion), etc., l'interprétation (le pointé) des réflexions dans le "bloc 3D de données sismiques après sommation" pouvant être facilité par l'utilisation de la migration temps 3D après sommation (poststack) suivi d'une démigration. On peut considérer ces techniques cinématiques comme un cas particulier des techniques d'imagerie profondeur 3D après sommation puisque permettant de construire, à partir du "bloc 3D de données sismiques apres sommation", un modèle 3D de la structure géologique: le modèle de vitesse dans le domaine profondeur. 1 0 2) Les techniques reposant sur l'utilisation du résultat de la migration profondeur 3D après sommation. Dans ce cas, l'interprétation est faite sur le "bloc 3D de données sismiques migrées profondeur" et différentes techniques peuvent être utilisées pour l'élaboration ou la mise à jour du modèle de vitesse telles que, par exemple, les techniques de "mise en place successive des couches" (cf. D.G. Rockliff, et al.; "3-D depth migration in the Diapir province of the central North Sea", 62th Ann. Internat. Mtg., Soc. Expl. Geophys., Expanded Abstracts, 925- 926 déjà cité), les techniques de balayage ("scan") de vitesse, etc. L'utilisation effective de ces techniques demande non seulement un bon logiciel d'imagerie 3D post-sommation (poststack) mais également de bonnes données ("bloc 3D de données sismiques après sommation"), c'est-à-dire de données résultant d'une sommation convenable, ce qui est difficile à réaliser sur des structures géologiques complexes. Le but des techniques de sommation classiques est d'améliorer le rapport signal sur bruit en ramenant à un déport (ou offset) nul (correspondant au cas o, on le rappelle, la source sismique et le récepteur d'ondes sont confondus) les données 3 0 sismiques obtenues avec différents déports non nuls. Les données sismiques après sommation (poststack) sont donc souvent considérées comme des données sismiques à déport nul et l'imagerie profondeur (par exemple la migration profondeur) de données après sommation (poststack) est assimilée à l'imagerie profondeur de données sismiques à déport nul. Cette assimilation concerne en général les données sismiques pour lesquelles le trajet aller (avant réflexion) coïncide avec le trajet retour (après réflexion), ce qui correspond ainsi à une expérience dite de "réflecteurs explosifs". Cette technique dite des "réflecteurs explosifs" o la réflexivité ("reflectivity" en anglais) est considérée comme source sismique excitant le milieu au voisinage du temps t = 0 et o les vitesses de propagation sont divisées par 1 0 deux de façon à ce qu'un trajet simple corresponde à un trajet aller- retour, est décrite par exemple dans: D. Loewenthal, et al; 1976, "The wave equation applied to migration", Geophysical Prospecting, 24, 380-399. Dans la suite de la description, on assimilera donc migration profondeur de données à déport (offset) nul et migration profondeur après sommation (migration "poststack"). Les techniques classiques de sommation reposent sur des corrections dynamiques d'obliquité dites corrections NMO, définies par des vitesses de sommation résultant elles-mêmes de techniques d'analyses de vitesse. Ces techniques sont justifiées si l'on suppose que les indicatrices (courbes de temps d'arrivée) sont des fonctions hyperboliques du déport (offset) pour un point milieu (CMP) donné, et si l'on dispose des bonnes vitesses de sommation. L'existence de difficultés résultant de pendages des réflecteurs, a amené à introduire, préalablement à la sommation classique, des corrections de pendage dites corrections DMO. Il est connu que le choix des algorithmes de corrections DMO est assez critique. Il en existe qui permettent de sommer des événements sismiques associés à de forts pendages, sous réserve 3 0 cependant d'hypothèses restrictives liées à l'utilisation de ces corrections de pendage. L'existence d'importantes variations latérales de vitesse introduit une difficulté supplémentaire, la sommation traditionnelle (même après les corrections de pendage DMO) ayant peu de chance de produire une sommation constructive. Dans ce cas, le bloc 3D de données sismiques migrées profondeur est difficilement exploitable de par sa mauvaise qualité (Fig. 2) Il est alors naturel de chercher à refaire la sommation, par exemple en la limitant aux déports faibles. Une telle opération n'est pas si simple, car, encore faut-il pouvoir détecter les courbures des hyperboles sensées représenter les indicatrices. De plus, dans cette opération, on risque de se trouver handicapé par le mauvais rapport signal sur bruit, contrepartie de la limitation aux seuls déports proches. Les réflexions multiples (importantes dans l'exemple d'illustration) et la qualité moyenne de la restitution des événements sismiques profonds à l'aplomb du dôme, risquent de gêner gravement l'interprétation requise pour l'élaboration du modèle de vitesse. L'image obtenue par migration 2D avant sommation appliquée à un profil d'acquisition (Fig. 4) montre une nette 2 0 amélioration, notamment pour la restitution des événements sismiques profonds à l'aplomb du dôme. Cependant l'image n'est pas complètement satisfaisante puisqu'elle ne prend pas en compte les aspects tridimensionnels de la propagation contrairement au résultat de la Fig. 2. La méthode selon l'invention vise à exploiter les qualités des images obtenues par migration avant sommation pour la mise en oeuvre de l'imagerie profondeur 3D après sommation, évitant par là les inconvénients des méthodes antérieures. La méthode selon l'invention a pour objet de synthétiser des 3 0 données sismiques à déport nul, à partir de données acquises dans le cadre d'opérations de prospection sismique-réflexion (notamment à trois dimensions) et à déports multiples, d'un 2735244 milieu souterrain. Elle est caractérisée en ce qu'elle comporte les étapes suivantes: a) on sélectionne une distribution de vitesse suivant chacune des dimensions du milieu exploré (modèle de vitesse dans le domaine profondeur); b) on applique aux données sismiques une migration profondeur avant sommation (migration "prestack"), réalisant ainsi une sommation des événements sismiques dans le domaine profondeur; et c) on synthétise les données sismiques à déport nul dans le domaine temps en appliquant une technique de modélisation de propagation d'ondes à déport nul, aux résultats obtenus par la dite sommation dans le domaine profondeur, le modèle de vitesse utilisé pour cette modélisation étant le même que celui utilisé 1 5 pour la migration avant sommation. La méthode peut comporter en outre une migration profondeur des données sismiques après sommation suivie éventuellement d'une élaboration ou d'une mise à jour, à partir de l'interprétation des données sismiques résultantes, du modèle de vitesses sélectionné en utilisant par exemple une technique de "mise en place successive des couches" ou bien encore une technique dite de balayage (scan) de vitesse de migration. Elle peut aussi comporter une élaboration ou une mise à jour du modèle de vitesse à partir de l'interprétation dans le domaine temps des données sismiques synthétisées par la méthode, que l'on fasse ou non appel à des techniques de migration temps pour réaliser cette interprétation. La méthode présente de nombreux avantages. * Le principal est d'offrir la possibilité d'extrapoler à déport 3 0 nul les informations sismiques obtenues pour des déports non 6 2735244 nuls, même dans le cas de cinématiques compliquées (c'est-à-dire d'indicatrices non hyperboliques sur des collections point milieu (CMP)) résultant par exemple d'importantes variations latérales de vitesse dans le milieu. Elle offre par ailleurs les avantages complémentaires suivants: a Elle permet une atténuation des bruits (bruits aléatoires, bruits organisés tels que réflexions multiples, ondes converties, ondes de surface, ondes guidées, etc.) grâce à la sommation des événements dans le domaine profondeur. La possibilité de sommation le long d'indicatrices non hyperboliques d'une part, et le très grand degré de couverture réalisé lors de la sommation dans le domaine profondeur d'autre part, ont pour effet que l'atténuation des bruits est beaucoup plus efficace que dans une 1 5 sommation traditionnelle (basée sur des corrections dynamiques NMO) o la couverture est beaucoup plus faible. En corollaire, la méthode permet une atténuation très raisonnable des bruits à partir d'un sous-ensemble très restreint des données enregistrées: une mise en oeuvre de la méthode sur des données décimées (décimation des tirs, décimation ou réduction de la gamme de déports) permet une réduction très sensible des coûts de calcul. * Contrairement aux techniques de sommation traditionnelle (sommation NMO et éventuellement DMO) qui reposent sur une analyse de vitesse, la méthode selon l'invention est entièrement automatique, une fois défini un modèle de vitesse (dans le domaine profondeur). La qualité de ce modèle (c'est à dire la précision avec laquelle il représente la vraie distribution de vitesse de propagation), influe sur la qualité de la sommation c'est-à-dire sur le renforcement du signal et sur l'atténuation des bruits. Mais c'est là la seule incidence du modèle de vitesse. Même si la sommation s'effectue dans le domaine profondeur, la cinématique des événements à déport nul synthétisés par la méthode n'est pas affectée par le modèle de vitesse choisi: elle est intrinsèque. Notons que l'on peut éventuellement contrôler la qualité de la sommation réalisée, par exemple par examen de l'amplitude normalisée des échantillons sismiques sommés dans le domaine profondeur (Fig. 6). D'autres caractéristiques et avantages de la méthode selon l'invention, apparaîtront à la lecture de la description ci-après d'un exemple de mise en oeuvre non limitatif, en se référant aux dessins annexés o: 1 0 - la Fig. 1 montre comment la démarche suivie par la méthode selon l'invention se démarque des méthodes antérieures; -la Fig. 2 montre une section sismique suivant deux axes X,Z, représentant respectivement la position horizontale le long d'un profil sismique et la profondeur, extraite d'un bloc de données 3D migrées profondeur après sommation, obtenue dans le cadre de l'exploration sismique d'une structure salifère de mer du Nord, le modèle de vitesse utilisé pour mettre en oeuvre cette migration (représenté sur les Fig. 3a, 3b) ayant été obtenu par la technique de migration de cartes; -les Fig. 3a et 3b montrent le modèle de vitesse 3D suivant trois axes X, Y, Z représentant respectivement la position horizontale le long d'un profil sismique, la position horizontale transverse au profil sismique et la profondeur; qui est utilisé tant pour l'obtention des résultats présentés sur la Fig. 2 que pour la mise en oeuvre de la méthode selon l'invention (pour simplifier la représentation graphique, ont été simplement représentées, sur la Fig.3a, les bases des couches 1 et 2 et, sur la Fig.3b, les bases des couches 4 et 6, alors que le modèle comprend en réalité 7 couches); les vitesses varient entre 1800 m/s (Ire couche) et 4500 m/s (7e couche); - la Fig. 4 montre la section sismique obtenue par migration profondeur 2D avant sommation ( migration "prestack"), des données de la ligne d'acquisition coïncidant avec la section sismique de la Fig. 2; - la Fig. 5 montre le modèle de vitesse utilisé pour l'obtention du résultat de la figure 4, ce modèle étant une tranche verticale passant par le profil d'une version lissée du modèle de vitesse de la Fig. 3, les vitesses variant entre 1627 m/s (parties les plus claires) et 4500 m/s (parties les plus sombres) -la Fig. 6 montre la section sismique à déport (T représentant le temps) nul synthétisée par la méthode selon l'invention, appliquée au profil d'acquisition considéré; - la Fig. 7 montre, pour le profil d'acquisition considéré, un résultat de contrôle de qualité de la sommation (stack) dans le domaine profondeur par examen de "l'amplitude normalisée des échantillons sismiques sommés"; - la Fig. 8 montre la section, associée au profil considéré, extraite du bloc sismique 3D après sommation (poststack) non migré traditionnel -les Fig. 9 et 10 montrent une comparaison des résultats de migration profondeur 2D après sommation (poststack) appliquée aux données respectivement des Fig. 6 et 8 avec le modèle de vitesse de la Fig. 5. II) Description de la méthode La démarche suivie on le rappelle, à sélectionner une distribution de vitesses (ou modèle de vitesse) suivant chacune des dimensions du milieu exploré, à appliquer aux données sismiques une migration profondeur avant sommation réalisant ainsi une sommation des événements sismiques dans le domaine 3 0 profondeur, et à synthétiser les données sismiques à déport nul dans le domaine temps en appliquant une technique de modélisation de propagation d'ondes à déport nul aux résultats obtenus par ladite sommation dans le domaine profondeur, le modèle de vitesse utilisé pour cette modélisation étant le même que celui utilisé pour la mise en oeuvre de la migration avant sommation. Les événements migrés avant sommation étant toujours stationnaires au voisinage du déport nul, la sommation des données sismiques qui est effectuée dans le domaine profondeur, est toujours cohérente dans ce voisinage et le résultat de cette sommation dans le domaine profondeur est donc l'image migrée à déport nul (avec la distribution de vitesse choisie) des données. L'appellation dite ici de sommation dans le domaine profondeur se justifie puisque l'on réalise une migration profondeur avant 1 5 sommation (migration "prestack") des événements sismiques. Comme on l'a déjà vu, la sommation est toujours cohérente au voisinage du déport nul quel que soit le modèle de vitesse sélectionné. Sur un plan pratique, la disponibilité de l'information sismique pour les faibles déports est essentielle pour garantir l'aspect constructif de la sommation et, en conséquence, l'amélioration du rapport signal sur bruit sera d'autant plus grande que le modèle de vitesse est précis. Des contrôles de qualité (par exemple la représentation de "l'amplitude normalisée des échantillons sismiques sommés"), permettent d'analyser la 2 5 qualité de la sommation réalisée dans le domaine profondeur. Comme défini par Yilmaz: "Seismic data processing", Investigation in geophysics N 2, Society of Exploration Geophysicists, Tulsa, 1987, on appelle "amplitude normalisée des échantillons sismiques sommés" ("normalized stacked amplitude") la quantité, toujours comprise entre -1 et 1, (,ei)/(z(|ei,) o i est l'indice de sommation des échantillons e1. o 2735244 1 0 La sommation dans le domaine profondeur ayant donc fourni comme résultat l'image sismique migrée à déport nul, la synthèse des données sismiques à déport (offset) nul (c'est-à-dire le retour dans le domaine temps) est obtenue en appliquant à ce résultat (considéré alors comme modèle de réflexivité ("reflectivity" en anglais) du sous- sol) une technique de modélisation de propagation d'ondes à déport nul, le modèle de vitesse utilisé pour cette modélisation étant le même que celui utilisé pour la migration avant sommation (prestack). Pour la modélisation, on peut utiliser la technique dite des "réflecteurs explosifs" citée plus haut. Dans ce cas, la migration avant sommation ayant fourni l'image migrée des données à déport nul pour le modèle de vitesse choisi, le résultat de cette modélisation sera effectivement les données sismiques à déport 1 5 (offset) nul puisque la modélisation dite par réflecteurs explosifs et la migration après sommation (poststack) sont deux opérations inverses l'une de l'autre. Ce résultat reste vrai sur le plan cinématique même si l'on n'utilise pas, pour la modélisation, la simplification par réflecteurs explosifs. La cinématique des événements sismiques à déport nul synthétisés est donc indépendante du modèle de vitesse choisi: elle est intrinsèque. On peut envisager différentes stratégies de mise en oeuvre de la méthode et, pour chaque stratégie, différents algorithmes de calcul. II 1 Différentes stratégies de mise en oeuvre 1. Utilisation de la migration profondeur 3D avant sommation Cette stratégie permet d'aborder des problèmes très difficiles o la prise en compte des effets 3D est particulièrement importante. Le coût du calcul handicape une telle stratégie, à 3 0 l'heure actuelle. Du fait de la diminution attendue de la tarification des calculs, et aussi de l'évolution de l'acquisition sismique (qui pourra rendre moins lourds les calculs requis pour la migration 3D avant sommation), la méthode présente un très grand intérêt pour l'amélioration du modèle de vitesse par toute technique connue telle celle de "mise en place successive des couches", utilisant des résultats de migration 3D post-sommation, car il coûte moins cher d'itérer sur les résultats d'une migration profondeur 3D après sommation que sur ceux d'une migration profondeur 3D avant sommation. 2. Utilisation de la migration profondeur 2D avant sommation, pour des mises en oeuvre en 3D. Dans le cas o l'on dispose de données sismiques par profils sériés, la méthode peut être mise en oeuvre en remplaçant la migration 3D avant sommation mentionnée au II 1 1 ci-dessus, par une série de migrations 2D avant sommation. La mise en oeuvre de cette stratégie demande alors la réalisation des étapes suivantes: * l'organisation des données 3D en profils rectilignes (ou quasi rectilignes) par des rabattements appropriés ou toute autre technique, * le découpage du modèle de vitesse 3D en tranches verticales associées à ces profils, * Pour chaque profil, il faut - obtenir, à partir des données du profil, une section migrée avant sommation avec un algorithme quelconque de migration 2D avant sommation (prestack) en utilisant comme modèle de vitesse la tranche verticale associée au profil (ou une approximation de ce modèle), - obtenir les données sismiques à déport nul associées au profil à l'aide d'une modélisation de propagation d'ondes à 3 0 déport nul, par "réflecteur explosif" ou par une autre technique, utilisant comme données de base, la section migrée 2D avant sommation et le même modèle de vitesse que celui ayant servi à l'obtention de cette section migrée. organiser les données à déport nul synthétisées pour chaque profil en bloc de données 3D. Cette opération est triviale si les profils ont été choisis rectilignes, parallèles et suffisamment serrés; sinon, il faut avoir recours à des techniques d'interpolation classiques. Le bloc de données 3D synthétisé est prêt à servir de données exploitables pour une interprétation ou pour servir d'entrée à des algorithmes de migration profondeur 3D après sommation. Une telle mise en oeuvre est particulièrement adaptée au cas d'acquisitions 3D de sismique marine lorsque les profils de 1 5 capteurs ne sont pas trop éloignés des profils de tirs ce qui rend accessible l'information sismique à faible déport dont nous avons souligné l'aspect essentiel. L'intérêt de cette stratégie de mise en oeuvre par rapport à la stratégie mentionnée au II.1 1, réside dans la diminution considérable (du moins avec les acquisitions sismiques actuelles) du coût de calcul résultant du remplacement de la migration 3D avant sommation par une série de migrations 2D avant sommation. La contrepartie de cette simplification est que la méthode ne fournira de résultats vraiment satisfaisants que si l'on peut raisonnablement considérer que les ondes enregistrées sur chaque profil se sont propagées suivant des trajectoires (rayons) qui ne s'écartent pas trop du plan vertical passant par le profil. Il faut noter cependant la robustesse de la méthode tant vis à vis 3 0 de cette hypothèse simplificatrice que vis à vis de l'exactitude du modèle de vitesse choisi: la méthode procède par sommation d'événements qui, de toute façon, sont stationnaires au voisinage du déport nul; le seul enjeu est le rapport signal sur bruit et la cinématique des données à déport nul (zéro offset) synthétisées sera quant à elle toujours correcte. La méthode a été appliquée sur la ligne d'acquisition coïncidant avec la section de la Fig. 2. La section à déport nul (zéro offset) obtenue par mise en oeuvre de la méthode, est donnée sur la Fig. 6. On note qu'elle est sensiblement différente de la section correspondante, extraite du bloc sismique 3D post- sommation non migré traditionnel (Fig. 8), et qu'elle est plus appropriée que celle-ci à l'imagerie de la structure, comme le montre la comparaison des résultats de migration profondeur 2D post-sommation obtenus à partir de ces sections (Fig. 9 et 10), ainsi que la comparaison avec le résultat de migration profondeur 3D post-sommation de la Fig. 2. La mise en oeuvre de la méthode sur les différents profils d'acquisition permet l'obtention d'un bloc sismique 3D à déport nul (zéro offset) prêt à servir de données exploitables pour une interprétation ou pour servir d'entrée à des algorithmes de migration 3D post-sommation. Il. 2 Différents algorithmes de calcul Un modèle de vitesse ayant été choisi, la méthode demande la mise en oeuvre d'algorithmes: - pour la migration avant sommation - pour la modélisation de propagation d'ondes. 1. Migration avant sommation Qu'il s'agisse d'une mise en oeuvre basée sur l'utilisation de la migration 3D avant sommation (II.1 1) ou sur l'utilisation de la migration 2D avant sommation (II.1 2), différents algorithmes de migration avant sommation peuvent être envisagés correspondant à: - différents assemblages des données sismiques (algorithmes de migration par déport (offset), par point de tir, par ondes planes, par récepteur, s-g, etc.), tels que présentés par exemple par: N.O. Whitmore et al.; 1993, "Common offset and common angle prestack depth migration applied to North Sea data", 63rd Ann. Internat. Mtg., Soc. Expl. Geophys., Expanded Abstracts, 886- 887. ou par A.J. Berkhout; 1984, "Seismic migration: B. practical aspects", Elsevier. - différentes techniques de modélisation de propagations d'ondes (modélisation numérique de forme d'ondes, calcul de fonctions de Green par des techniques cinématiques, calcul de fonctions de Green par formes d'ondes); - différents principes d'imagerie présentés par exemple dans J.F. Claerbout; 1971, "Towards a unified theory of reflector mapping", Geophysics, 36, 467-481, ou bien dans A.J. Berkhout; 1984, "Seismic migration: A. theoretical aspects", Elsevier. Même si les calculs sont organisés parfois d'une façon différente, tous ces algorithmes ont en commun le même principe (d'o le terme migration avant sommation): - procéder d'abord à une remise en place des événements (migration) par modélisation de propagations d'ondes utilisant le modèle de vitesse choisi et par application d'un principe d'imagerie; - procéder ensuite à la sommation des événements sismiques remis en place. 2. Modélisation de propagations d'ondes à déport nul Qu'il s'agisse de la mise en oeuvre présentée au II.1 1 ou bien de celle présentée au II. 1 2, cette étape consiste à calculer la réponse sismique à déport nul à partir: - du résultat de la migration avant sommation considéré comme réflectivité à offset nul du milieu; - du modèle de vitesse utilisé pour la mise en oeuvre de cette migration avant sommation. Pour la mise en oeuvre présentée au II.1 1 il faut modéliser des propagations d'ondes à 3D, et pour celles présentées au II.1 2 il faut procéder à une série de modélisations de propagations d'ondes à 2D. Une grande variété d'algorithmes est envisageable aussi pour la réalisation de cette étape: par exemple, les techniques numériques de modélisation de formes d'ondes (basées sur la résolution numérique d'une équation d'ondes ou d'une de ses approximations telle que celles reposant sur l'approximation paraxiale) ou bien la technique dite des rayons. Une approche particulièrement attrayante pour cette modélisation de propagation d'ondes à déport nul est celle qui repose sur la simplification des réflecteurs explosifs telle que décrite par exemple par: D. Loewenthal, L. Lu, R. Robertson, J. Sherwood; 1976, "The wave equation applied to migration", Geophysical Prospecting, 24, 280-399. Si l'on suppose que les trajets aller (avant réflexion) et retour (après réflexion) sont identiques, tout revient, du moins au niveau de la cinématique, à considérer, on le rappelle, la réflexivité ("reflectivity") du milieu comme terme "source" agissant au voisinage du temps t = O, pour peu que l'on ait divisé la vitesse par 2 (de façon à ce qu'un trajet simple corresponde à un trajet aller-retour). Si le principe d'imagerie utilisé pour la migration profondeur à déport nul et le terme "source" utilisé pour la modélisation par réflecteurs explosifs sont convenablement définis, alors on peut montrer dans le cas o les propagations d'ondes sont modélisées 1 0 par une approximation paraxiale de l'équation d'ondes, que ces deux opérations sont inverses l'une de l'autre, ce qui justifie que les données sismiques à déport nul synthétisées par la méthode ne dépendent pas du modèle de vitesse et ont donc un caractère intrinsèque (comme mentionné au début de la section II). Dans le 1 5 cas o l'on ne prend pas de telles précautions pour la définition du principe d'imagerie pour la migration profondeur à déport nul et du terme "source" pour la modélisation par réflecteurs explosifs, la méthode reste applicablepuisque le résultat demeure vrai quant à la cinématique des données qui constituent l'information sismique utile pour l'imagerie des structures géologiques. III Applications de la méthode La méthode selon l'invention a des applications intéressantes notamment pour l'imagerie sismique 3D de structures géologiques présentant de fortes variations latérales de vitesse, c'est-à-dire celles o le recours à l'imagerie profondeur est importante pour la restitution des structures. La qualité du modèle de vitesse ayant une incidence sur la qualité du résultat, la méthode trouve son intérêt au stade o un premier modèle de 3 0 vitesse dans le domaine profondeur a été créé. Un tel modèle a pu être obtenu à partir des données post-sommation traditionnelles par "migration de cartes", opération bien connue des spécialistes, ou par mise en oeuvre de la technique de "mise en place successive des couches" définie précédemment, ou encore par d'autres techniques (tomographie, balayage (scan) de vitesse de migration, etc.). Les données à déport nul synthétisées par la méthode sont, de par leur qualité, particulièrement adaptées pour une interprétation ou comme données d'entrée pour un processus de migration profondeur 3D postsommation: - pour l'obtention finale d'un cube de données sismiques 3D migrées en profondeur lorsqu'un modèle de vitesse a été définitivement choisi; - pour l'élaboration d'un modèle de vitesse, par mise en oeuvre de la technique de "mise en place successive des couches" évoquée précédemment, la qualité des données d'entrée exploitées facilitant largement le pointé des limites de couches sur le résultat de la migration profondeur 3D post-sommation, ou bien encore - pour l'élaboration du modèle de vitesse par toute méthode utilisant les résultats de la migration profondeur 3D post- sommation. REVENDICATIONS 1) Méthode pour synthétiser des données sismiques à déport nul, à partir de données acquises dans le cadre d'opérations de prospection sismiqueréflexion à deux ou trois dimensions à déports multiples, d'un milieu souterrain, caractérisée en ce qu'elle comporte les étapes suivantes: a) on sélectionne une distribution de vitesse suivant chacune des dimensions du milieu exploré; b) on applique aux données sismiques une migration profondeur avant sommation, réalisant ainsi une sommation des événements sismiques dans le domaine profondeur; et c) on synthétise les données sismiques à déport nul dans le domaine temps en appliquant une technique de modélisation de propagation d'ondes à déport nul, aux résultats obtenus par la dite sommation dans le domaine profondeur, le modèle de vitesse utilisé pour cette modélisation étant le même que celui utilisé pour la mise en oeuvre de la migration avant sommation. 2) Méthode selon la revendication 1, caractérisée en ce que l'on synthétise les données à déport nul dans le domaine temps par une technique de modélisation dite à réflecteurs explosifs. 3) Méthode selon la revendication 1 ou 2, caractérisée en ce que, dans le cadre d'opérations de prospection sismique-réflexion en trois dimensions ou 3D, elle comporte en outre une migration profondeur 3D après sommation des données sismiques synthétisées. 4) Méthode selon la revendication 1 ou 2, caractérisée en ce que, dans le cadre d'opérations de prospection sismique-réflexion en deux dimensions ou 2D, elle comporte en outre une migration profondeur 2D après sommation des données sismiques synthétisées. ) Méthode selon la revendication 3 ou 4, caractérisée en ce qu'elle comporte en outre, à partir des données résultantes interprétées, une élaboration ou une mise à jour dudit modèle de vitesses. 6) Méthode selon la revendication 5, caractérisée en ce que l'on réalise l'élaboration ou la mise à jour du modèle de vitesses par une technique de "mise en place successive des couches". 7) Méthode selon la revendication 5, caractérisée en ce que l'on réalise l'élaboration ou la mise à jour du modèle de vitesses par une technique dite de balayage de vitesse de migration. 8) Méthode selon la revendication 1 ou 2, caractérisée en ce qu'elle comporte une interprétation des données sismiques synthétisées. 9) Méthode selon la revendication 8, caractérisée en ce que l'interprétation des données sismiques synthétisées comporte l'utilisation d'une technique de migration temps. ) Méthode selon la revendication 8 ou 9, caractérisée en ce qu'elle comporte en outre une élaboration ou une mise à jour du modèle de vitesse à partir des données interprétées.