FR 2876458 A1 20060414 FR 0410672 A 20041008 PERFECTIONNEMENT AUX TRAITEMENTS SISMIQUES POUR LA SUPPRESSION DES REFLEXIONS MULTIPLES La présente invention est relative aux procédés de traitement sismique. Plus particulièrement, elle concerne les procédés pour la suppression des "réflexions multiples ". DOMAINE TECHNIQUE GENERAL io On sait qu'une technique répandue de recherche de pétrole ou de gaz consiste à effectuer une prospection sismique du sous-sol. Pour imager la structure du sous-sol, le géophysicien utilise les techniques de sismique réflexion . Celles-ci consistent à émettre à la surface du sol des signaux acoustiques 15 et à les enregistrer après leurs réflexions successives sur les interfaces entre couches géologiques. En sismique terrestre, à partir de plusieurs points de la surface du sol appelés points de tir, on émet des vibrations sismiques (ondes de compression et de cisaillement), et l'on recueille en différents points de la surface, à l'aide de récepteurs sismiques (géophones), les ondes acoustiques réfléchies par les interfaces entre couches géologiques, appelées réflecteurs. Les récepteurs sismiques convertissent en signaux électriques les ondes réfléchies. Ces récepteurs sont disposés de façon telle et en nombre suffisant pour que les signaux enregistrés, appelés traces, constituent des données sismiques et permettent de reconstituer la configuration des couches géologiques. Dans la pratique, le signal enregistré par un récepteur sismique présente une amplitude qui varie constamment en fonction du temps et les pics enregistrés correspondent normalement à des réflecteurs entre couches. En réalité, l'interprétation des données d'un tir est beaucoup plus complexe. En effet, d'une part, les ondes peuvent traverser une première interface entre couches, pour se réfléchir sur l'interface suivante, puis sur la première interface et ainsi de suite avant de parvenir au géophone, et d'autre part, l'amplitude des enregistrements diminue très rapidement en fonction du temps. Les enregistrements comportent donc des pics correspondant à des réflexions multiples, ou réflexions parasites, qu'il convient d'éliminer pour reconstituer correctement la topographie du soussol. Une technique analogue est utilisée pour la prospection sismique marine, lo le tir étant effectué dans ce cas à quelques mètres sous la surface de l'eau et les récepteurs sismiques, ou hydrophones, étant eux-mêmes disposés à la même profondeur. Dans cette méthode de prospection marine, les inconvénients rappelés cidessus sont toutefois amplifiés par le fait que le fond marin est hautement 15 réflecteur, de même que l'interface eau/air. Un but général de la présente invention est donc de fournir un procédé de traitement sismique permettant d'éliminer au moins partiellement les signaux parasites correspondant aux réflexions multiples dans les enregistrements effectués par les récepteurs sismiques lors d'un tir de prospection. ETAT DE LA TECHNIQUE Il a déjà été proposé un certain nombre de procédés de traitement sismique permettant d'atténuer les réflexions multiples des enregistrements 25 sismiques. Ces procédés sont principalement basés sur deux techniques: la première consiste à utiliser les données pour obtenir la topographie du sous-sol, la seconde consiste à modéliser les réflexions multiples afin de les soustraire aux données sismiques et ainsi ne conserver que les données sismiques utiles, à savoir les réflexions primaires, pour l'obtention de la topographie du sous sol. Verschuur et Al., 1992, [1], ont présenté un procédé dans lequel les données sismiques sont utilisées pour modéliser les réflexions multiples dans s le domaine deux dimensions. Egalement van Dedem et Verschuur, 1998, [2], ont présenté une généralisation du procédé au domaine trois dimensions. Dans les deux cas, le procédé de modélisation de réflexions multiples basé sur les données sismiques requiert des convolutions entre les traces de données acquises et io les traces du point de tir situées aux positions de récepteur. Cependant, ce procédé présente un inconvénient du fait qu'il n'y a pas un point de tir acquis pour chaque emplacement de géophone. Wiggins, 1988, [3] a proposé un autre procédé utilisant quant à lui des techniques d'extrapolation d'onde. Dans son procédé, Wiggins propose de propager les données acquises en surface à travers la couche d'eau jusqu'au fond marin de telle manière que les ondes primaires descendantes (ondes incidentes) coïncident dans l'espace avec les réflexions multiples propagées en inverse (réflexions multiples reflétées). Les opérateurs adaptatifs de soustraction entre ces deux champs d'ondes représentent ainsi la réflectivité du fond marin. Cependant, ce procédé présente l'inconvénient de nécessiter la connaissance de la topographie du fond marin pour permettre l'adaptation à la bonne position des deux champs d'onde propagés. En outre, le procédé a été présenté exclusivement dans le domaine deux dimensions. Dans des travaux plus anciens, Berryhill et Kim, 1986, [4], ont également proposé d'employer l'extrapolation d'équation d'onde pour modéliser les peglegs. Leur procédé consiste à effectuer une extrapolation d'équation d'onde (redatuming) dans le domaine deux dimensions, en fait par une intégrale de Kirchhoff, à partir des données sismiques enregistrées jusqu'au fond marin sélectionné, et ensuite encore jusqu'à la surface. Ici aussi, seul la connaissance de la vitesse de propagation de l'onde dans la couche d'eau et la connaissance de la topographie du fond marin choisi sont nécessaires, en plus des points de tir d'entrée. Très récemment, tout comme Berryhill et Kim (1986), Lokshtanov, 2000, [5], a montré une implémentation pour la prédiction d'équation d'onde actuelle, dans le domaine de c-p cette fois. Avec ce procédé, on effectue ensuite une soustraction des peg-legs de la couche d'eau en utilisant la topographie de fond marin. Le long du fond marin, les coefficients de réflexion sont censés être lo égaux à 1 pour n'importe quel angle de propagation. Cependant, un inconvénient des deux procédés cités ci-dessus est la nécessité de la connaissance de la topographie du fond marin. Or, dans certains cas (exploration sismique marine pour des profondeurs de fond marin très importantes), cette topographie est inconnue des géologues, ou ne peut être déterminée assez précisément pour obtenir des résultats exploitables. Un but de la présente invention est de fournir un procédé permettant de palier au moins l'un des inconvénients précités. PRESENTATION DE L'INVENTION A cet effet on prévoit selon l'invention un procédé de traitement sismique, caractérisé en ce que, pour supprimer des réflexions multiples sur des données sismiques, on migre les données sismiques en temps ou en profondeur et on traite les données ainsi migrées pour déterminer une approximation des réflexions multiples que l'on soustrait des données sismiques. On notera que contrairement aux techniques de modélisation à partir des données sismiques, connues dans l'art antérieur, une telle modélisation des réflexions multiples superficielles à partir de la migration des données sismiques, peut être mise en oeuvre sans qu'il ne soit nécessaire de disposer de géométries d'acquisition denses. En outre, le procédé présenté peut être appliquée à l'acquisition au s moyen de câbles de fond d'océan, communément appelés OBC (abréviation de l'anglais Ocean Bottom Cable). Par ailleurs, il n'est aucunement besoin de sélectionner la ou les profondeur(s) du fond marin, car les formes d'ondes des événements migrés sont employées comme réflecteurs à leurs propres emplacements. io En outre, non seulement le fond marin, mais également tout réflecteur situé sous le fond marin concourt à la génération de réflexions multiples, aussi profondément que la bonne connaissance du champ de vitesse de propagation le permet. Enfin, le procédé de modélisation de réflexions multiples lui-même ne 15 requiert aucune manipulation des données sismiques. En effet, le procédé ne nécessite que les coordonnées des traces et la section migrée. La contribution des données sismiques n'intervient qu'à l'étape conventionnelle de migration. PRESENTATION DES FIGURES D'autres caractéristiques et avantages de l'invention ressortiront encore de la description qui suit, laquelle est purement illustrative et non limitative et doit être lue en regard des dessins annexés, sur lesquels: la figure 1 est une illustration schématique de l'acquisition de données sismiques marines montrant des exemples de réflexions à partir du fond marin et de réflecteurs sous-jacents, - la figure 2 est un organigramme du procédé de traitement sismique, la figure 3 est un graphique illustrant un champ des vitesses du sous-sol pour une onde, les figures 4 à 8 sont des schémas illustrant des étapes du procédé de traitement sismique sur une vue en coupe comprenant une couche d'eau, un fond marin et un sous-sol constitué d'un empilement de couches géologiques la figure 9 est un schéma illustrant une réflexion multiple de type peg-leg et des réflexions primaires. la figure 10 est un deuxième organigramme du procédé de traitement sismique, la figure 11 est un autre schéma illustrant des étapes du procédé de traitements sismiques. DESCRIPTION Généralités La figure 1 illustre schématiquement une prospection sismique marine. On a en particulier représenté une masse d'eau 10 s'étendant au dessus d'un fond marin 12, avec à la surface de l'eau (surface 11) un navire 13 de prospection sismique qui remorque un ou plusieurs câbles 14 incluant une pluralité de récepteurs sismiques 15 (hydrophones). Une ou plusieurs sources sismiques 16 destinées à émettre une onde acoustique dans l'eau sont déployées à partir du navire 13. La source 16 peut être un canon à air comprimé, un vibrateur marin ou d'autres types de sources connus de l'homme de l'art. La vibration sismique émise descend le long de plusieurs trajectoires et est réfléchie aux interfaces entre des matériaux d'impédance acoustique différente tels que le fond marin 12 de l'océan, et les réflexions étant captées par les récepteurs 15. L'onde se propage également dans le soussol et est 10 réfléchie à l'interface 19. Pour simplifier l'illustration, une seule de ces interfaces réfléchissantes est représentée: en réalité, il peut y avoir un grand nombre de ces interfaces réfléchissantes (réflecteurs). Les trajectoires des ondes représentées en 17 et 18 sont appelées ondes directes ou primaires, car la trajectoire ne comprend qu'une réflexion à partir d'une interface ou le fond marin 12 avant une détection par un récepteur 15. Une réflexion multiple de type multiple de fond marin est également représentée par la trajectoire sismique 20 sur la figure 1, qui est réfléchie deux io fois à partir du fond marin 12 et une fois à partir de la surface 11 de la masse d'eau. Les réflexions multiples constituent une difficulté importante pour le traitement des données sismiques car elles masquent des événements sismiques primaires réels qui ont lieu dans un temps comparable à celui que met l'onde soumise à une réflexion multiple pour atteindre le même hydrophone et ils constituent une source de données parasites. Traitement de modélisation/soustraction Un exemple de réalisation du procédé selon l'invention va maintenant être décrit en référence aux figures 2 à 11. Ce procédé se déroule en deux phases. La première phase consiste à calculer un modèle de champ de vitesses de propagation des ondes et à migrer les données sismiques enregistrées. La seconde phase consiste à calculer un modèle des réflexions multiples. Migration des données enregistrées Une première étape (100) consiste à calculer un champ des vitesses de propagation des ondes sismiques. En effet, une migration temps ou profondeur consiste à déplacer les points du volume d'acquisition sur la base d'un champ des vitesses de propagation. s Le champ des vitesses peut être calculé par tout type de procédé connu de l'homme de l'art comme un procédé d'inversion sismique. Le procédé d'inversion sismique le plus simple consiste à tester tous les signaux d'entrée possibles, et à observer les sorties qu'ils engendrent. Avec ce procédé d'inversion sismique, la réponse modélisée qui s'approche le plus de la io réponse observée est considérée comme étant la solution de l'inversion. On obtient ainsi un modèle des vitesses. Un exemple de modèle des vitesses 80 est illustré à la figure 3. Il s'agit du modèle correspondant à la zone de sous-sol illustrée à la figure 4. Ce modèle de vitesses est représenté dans un repère dont l'abscisse 25 et l'ordonnée 26 sont graduées en kilomètres. Ainsi, le modèle de vitesses illustré à la figure 3 modélise les vitesses sur une distance de 9 km et une profondeur de 3 km. Ce modèle de vitesses comprend différentes zones 21, 22, 23, 24 correspondant à quatre couches géologiques. A titre indicatif, la première couche 21 située en dessous du fond marin présente une vitesse de propagation de 1500 m/s la seconde couche 22 située en dessous de la première couche 21 présente une vitesse de propagation de 2500 m/s, la troisième couche 23 comprise dans la seconde couche 22 présente une vitesse de propagation de 2000 m/s, et la quatrième couche 24 située en dessous de la seconde couche 22 présente une vitesse de propagation de 2000 m/s. Une fois que l'on a obtenu le champ des vitesses, la deuxième étape (110) du procédé consiste à migrer les traces obtenues au niveau de chaque récepteur sismique, qui constituent les données sismiques. Cette opération de migration des traces vise à restituer les formes correctes des interfaces géologiques. Les traces peuvent être migrées en temps ou en profondeur, selon la complexité du modèle géologique de vitesse utilisé pour la migration. A vitesse s constante ou à vitesse variable verticalement, les deux types de migration sont équivalents. La migration dans un cas comme dans l'autre peut être effectuée avant ou après sommation. Une migration en profondeur avant sommation des données sismiques consiste à déterminer, sur une collection de traces, les événements qui io décrivent le sous-sol à l'aplomb d'une position de surface (x, y). Une telle collection de traces est généralement rangée par classes de distances source-récepteur (distances également appelées distances de déport ou d'offset) et également par classes d'orientation des paires source-récepteur. En imagerie avant sommation (prestack imaging), les traces sismiques associées à une paire source-récepteur unique sont traitées de manière individuelle. En imagerie après sommation (post-stack imaging), en revanche les traces issues d'une pluralité de paires source-récepteur sont combinées par sommation (stacking) avant d'être traitées. La technique post-stack imaging augmente le rapport signal sur bruit des données sismiques et réduit le nombre de traces sismiques devant être traitées. La section migrée obtenue en sortie de la deuxième étape (110) est ensuite employée comme modèle de réflectivité du sous-sol, incluant de ce fait le fond marin mais également tout réflecteur situé dans la zone du sous-sol sous le fond marin. L'utilisation de la section migrée obtenue comme modèle de réflectivité dispense d'avoir à étudier la topographie du fond marin, une telle étude pouvant cependant être justifiée pour réduire l'espace mémoire d'ordinateur exigé, lorsque la profondeur du fond marin change considérablement. i0 Un exemple de section migrée 40 comprenant des réflecteurs 32a, 32b est illustrée à la figure 4. La partie 50 de la section migrée 40 illustrée à la figure 4 correspond aux quatre zones 21, 22, 23, 24 du modèle de champ des vitesses de propagation d'onde illustré à la figure 3. Cette section migrée est représentée dans un repère dont l'axe des x représente une distance et l'axe des y représente une profondeur exprimée en mètres. Le trait horizontal situé à l'altitude 0 correspond à la surface 11 de la masse d'eau. On passe ensuite à la deuxième phase du procédé de traitement sismique consistant à modéliser les réflexions multiples. i0 Modélisation des réflexions multiples Cette deuxième phase utilise la technique d'extrapolation d'équation d'onde pour modéliser les réflexions multiples. Plusieurs approximations de l'équation d'ondes à sens unique peuvent être employées pour l'extrapolation d'onde exigée par ce procédé, selon le degré de complexité du champ des vitesses: déphasage pour une vitesse de propagation constante, déphasage plus interpolation lorsque l'on traite les variations 20 progressives de champ de vitesses, extrapolateurs implicites à sens unique pour des variations latérales plus rapides de vitesses, extrapolateurs explicites pour des variations latérales fortes de vitesses. La troisième étape (120) du procédé est illustrée à la figure 5. Dans cette troisième étape (120), une onde directe 31 est émise vers le bas (par extrapolation d'onde), à partir de la position donnée (x, y, z) d'une source 30, pour illuminer les réflecteurs sismiques 32a, 32b matérialisés par la section sismique migrée. Il Dans le cas simple d'une vitesse constante, l'extrapolation des ondes est émulée en traçant des rayons droits. Le délai correspondant au temps de parcours de l'onde depuis la source jusqu'au réflecteur 32a est appliqué plus tard (à la quatrième étape) dans la condition de réflecteur détonnant pour le 5 réflecteur considéré. Dans le cas d'une vitesse variable, une impulsion à largeur de bande limitée est extrapolée en employant l'équation d'onde à sens unique jusqu'à ce que les réflecteurs soient illuminés. La quatrième étape (130) du procédé est illustrée à la figure 6. Dans cette quatrième étape, chaque échantillon 33a, 33b, 33c, 33d de la section migrée 40 est illuminé par l'onde directe 31 en multipliant chaque échantillon migré 33a, 33b, 33c, 33d, dans le domaine x-y-z (si l'on a effectué une migration en temps NMO) ou le domaine de x-y-t (si l'on a effectué une migration en profondeur DMO), par les valeurs de l'onde directe 31 descendante ou par une exponentielle complexe pondérée contenant le temps de propagation de l'onde depuis l'emplacement de source 30. Une fois que la section migrée 40 a été illuminée, alors l'énergie reflétée primaire est construite en démigrant, c'est-à-dire en migrant vers le haut, la section migrée "illuminée" (multiplié). Cette démigration est faite par tranches de fréquence constantes soit dans le domaine f-kx-ky si l'on traite des champs de vitesses constantes, soit dans le domaine f-x-y si l'on traite des champs de vitesses variables. Un champ d'ondes migrées 37 est donc émis (extrapolé) vers le haut, chaque onde 37a, 37b, 37c, 37d du champ d'ondes migrées 37 étant issue d'un échantillon illuminé 33a, 33b, 33c, 33d de la section migrée. Par ailleurs, il est à noter que pour ne pas surcharger le schéma de la figure 6, chaque onde 37a, 37b, 37c, 37d du champ d'onde 37 a été représentée par une flèche unique dirigée vers le haut, étant entendu que les ondes 37a, 37b, 37c, 37d se propagent dans toutes les directions à l'instar de l'onde 31. La cinquième étape (140) du procédé est illustrée à la figure 7. Dans cette cinquième étape (140), lorsque le champ d'onde migré 37 atteint la surface 11, un coefficient de réflexion est appliqué (par exemple égale à -1), et les données sont extrapolées à nouveau vers le bas dans le modèle, en employant les mêmes techniques d'extrapolation que pour le premier chemin ascendant (quatrième étape (130) du procédé). Ce dernier champ d'onde 38 extrapolé vers le bas est employé pour "illuminer" (en multipliant dans le domaine x-y-z ou de x-y-t) la section migrée 50. La sixième étape (150) du procédé est illustrée à la figure 8. Dans cette sixième étape, la section migrée 50 illuminée par le champ d'onde descendant 38 est migrée, comme dans la quatrième étape (130), vers la surface. Le résultat obtenu en grille régulière selon les axes x et y, est ensuite interpolé (160) dans les dimensions x et y vers les positions réelles de récepteur des données de point de tir recueillies pour être modélisées. Enfin, ces données sont transformées (160) par transformée FFT du domaine fréquentiel au domaine temporel, permettant ainsi l'obtention des réflexions multiples de surface du premier ordre. Pour obtenir les réflexions multiples du second ordre, les cinquième et sixième étapes décrites ci-dessus sont répétées en employant le dernier champ d'onde extrapolé (réflexions multiples du premier ordre) pour illuminer à nouveau la section migrée initiale. Cette opération est répétée récursivement pour obtenir chaque réflexion multiple d'ordre supérieur. Dans la description ci-dessus, on vient de présenter la modélisation d'un multiple de type fond marin. Cependant, d'autres types de multiples peuvent être modélisés par le procédé de traitement sismique tel que les peg-legs. Une réflexion multiple de type peg-leg est illustrée à la figure 9. Comme on le voit sur la figure 9, une vibration sismique en un point de tir S situé près de la surface 1 de l'eau peut donner naissance à une onde 2 qui franchit le fond marin 3 pour pénétrer dans une formation géologique sous-jacente, puis se réfléchit sur l'interface 4 entre cette formation géologique et la couche suivante. L'onde réfléchie franchit le fond marin 3 pour revenir dans la couche d'eau. Elle se réfléchit à nouveau contre l'interface eau/air 1 puis contre le fond marin 3 avant de parvenir au récepteur sismique R, qui enregistre un signal correspondant au trajet parcouru par l'onde 2. Ce trajet a le profil représenté sur la figure 9, qualifié de "peg-leg" (en anglais jambe de bois) à cause de sa forme dissymétrique.". Le procédé de traitement sismique permettant la modélisation des 10 réflexions multiples de type peg-leg est illustré en référence aux figures 10 et 11. Les trois premières étapes (200, 210, 220) du procédé restent inchangées. La quatrième étape (230) consiste à illuminer chaque échantillon du 1s réflecteur 32a de la section migrée par l'onde directe 31. Un premier champ d'ondes 99 est émis vers le bas - c'est-à-dire extrapolé (ou migré) vers le bas - pour illuminer le deuxième réflecteur 32b de la section migrée, chaque onde du premier champ d'ondes 99 étant issue d'un échantillon illuminé de la section migrée 40. Lorsque le deuxième réflecteur 32b est illuminé, on passe à la cinquième étape (240) du procédé. Un deuxième champ d'onde 37' est démigré (extrapolé vers le haut) à partir de chaque échantillon du deuxième réflecteur 32b illuminé. Lorsque le deuxième champ d'ondes 37' atteint la surface 11 de la masse 25 d'eau, on passe à la sixième étape (250) du procédé sismique consistant à extrapoler vers le bas un troisième champ d'ondes 38'. Lorsque le troisième champ d'ondes intersecte le premier réflecteur 32a, on passe à la septième étape (260) du procédé. Un quatrième champ d'onde est extrapolé vers le haut. Une fois les données interpolées et transformées par FFT, on obtient une modélisation des multiples de type peg-leg. Tous les types de réflexions multiples sont ainsi modélisables au moyen du procédé de traitement sismique ci-présenté. Par ailleurs tous ces types de 5 réflexions multiples sont modélisables à tous les ordres. Une fois les réflexions multiples modélisées, celles-ci sont soustraites des données sismiques brutes. Le procédé de traitement sismique décrit ci-dessus permet de modéliser les réflexions multiples pour tout type d'offset résultant de la réverbération entre la surface libre et les réflecteurs du sous-sol dans les domaines deux dimensions (2D) et trois dimensions (3D). Il est à noter que la modélisation des réflexions multiples qui vient d'être décrite peut être réalisée sur tous les réflecteurs du sous-sol modélisés par la section migrée obtenue en sortie de la première phase. II est également à noter que dans le cas du présent procédé, l'opération de migration des traces sismiques pour restituer les formes correctes des interfaces géologiques est effectuée avant la suppression des données parasites que constituent les réflexions multiples, contrairement aux techniques traditionnelles dans lesquelles la migration des traces est effectuée après avoir éliminé les données parasites. Par ailleurs, la description ci-dessus de réflexions multiples dans le cas de la sismique marine n'est nullement limitative. Des réflexions multiples ont lieu aussi en sismique terrestre et la présente invention peut être utilisée tout aussi bien pour déterminer des temps de parcours de réflexions multiples dans ce contexte. REFERENCES [1] J. Verschuur, A.J. Berkhout, et C.P.A. Wapenaar, 1992, Adaptive surfacerelated multiple elimination . Geophysics, 57(9) : 1166-1177. [2] van Dedem, E.J. et Verschuur, D.J., 1998, 3-D surface-related multiple elimination and interpolation , 68th Ann. Internat. Mtg: Soc. of Expl. Geophys., 1321-1324. [3] Wiggins, J.W., 1988, Attenuation of complex water-bottom multiples by wave-equation-based prediction and subtraction . Geophysics, 53(12), 1527-1539. [4] Berryhill, J.R., et Kim, Y.C., 1986, Deep-water peg legs and io multiples . Geophysics, 51, 2177-2184. [5] Lokshtanov, D., 2000, Suppression of water-layer multiples--from deconvolution to wave-equation approach , 70th Ann. Internat. Mtg: Soc. of Expl. Geophys., 1981-1984. REVENDICATIONS 1. Procédé de traitement sismique, caractérisé en ce que, pour supprimer des réflexions multiples (2, 20) sur des données sismiques, on migre les s données sismiques en temps ou en profondeur (110) et on traite les données ainsi migrées pour déterminer une approximation des réflexions multiples que l'on soustrait des données sismiques (120 à 160). 2. Procédé selon la revendication 1, caractérisé en ce qu'avant la migration, lo on modélise (100, 200) un champ des vitesses de propagation des ondes (80). 3. Procédé selon la revendication 1, caractérisé en ce que pour déterminer une approximation des réflexions multiples, on extrapole (120, 220) une onde sismique (31) à partir d'une source sismique (30) jusqu'à un réflecteur (32a, 32b) de la section migrée (40). 4. Procédé selon la revendication 3, caractérisé en ce qu'après cette étape d'extrapolation: - on illumine (130, 230) chaque point du réflecteur (32a) de la section migrée (40) intersectant l'onde sismique (31) extrapolée vers le bas, - on extrapole (130) vers le haut un premier champ d'ondes sismiques (37) à partir des points du réflecteur (32a) de la section migrée illuminés par l'onde sismique extrapolée vers le bas. 5. Procédé selon les revendications 4, caractérisé en ce qu'à l'issue de ces étapes d'illumination et d'extrapolation: - on illumine (140) chaque point de l'interface (11) eau/air intersectant les ondes du premier champ d'ondes (37), - on extrapole (140) vers le bas un deuxième champ d'ondes sismiques (38) à partir des points de l'interface (11) eau/air illuminés par les ondes du premier champ d'ondes (37). 6. Procédé selon la revendication 5, caractérisé en ce qu'à l'issue de ces nouvelles étapes d'illumination et d'extrapolation on illumine (150) chaque point du réflecteur (32a) de la section migrée (50) intersectant les ondes du deuxième champ d'ondes (38), i0 on extrapole (150) vers le haut un troisième champ d'ondes sismiques (39) à partir des points du réflecteur de la section migrée illuminés par les ondes du deuxième champ d'ondes (38). 7. Procédé selon la revendication 6, caractérisé en ce qu'il comprend en outre l'étape consistant à interpoler vers les positions réelles de récepteurs, les données sismiques ainsi extrapolées et à transformer du domaine fréquentiel au domaine temporel les données sismiques interpolées afin d'obtenir un modèle des réflexions multiples du premier ordre. 8. Procédé selon l'une des revendications 6 ou 7, caractérisé en ce qu'il comprend en outre les étapes consistant à : illuminer chaque point de l'interface (11) eau/air intersectant les ondes du troisième champ d'ondes (39), - extrapoler vers le bas un quatrième champ d'ondes sismiques à partir de chacun des points de l'interface (11) eau/air illuminés par les ondes du troisième champ d'ondes (39), illuminer chaque point du réflecteur (32a) de la section migrée (50) intersectant les ondes du quatrième champ d'ondes, extrapoler vers le haut un cinquième champ d'ondes sismiques à partir des points du réflecteur (32a) illuminés par les ondes du quatrième champ d'ondes. interpoler vers les positions réelles de récepteurs, les données sismiques ainsi extrapolées et transformer du domaine fréquentiel au domaine temporel les données sismiques interpolées afin d'obtenir un modèle des réflexions multiples du deuxième ordre, répéter ces différentes étapes pour obtenir un modèle des réflexions multiples d'ordre supérieur.