Brevet FR2794249A1 — publié le 20001201

L'invention a trait au domaine du traitement des données sismiques. L'invention concerne plus particulière- ment un procédé de prédiction des temps d'arrivée d'événements de réflexion multiple qui masquent des événements sismiques réels dans des enregistrements d'exploration sismique à partir de mesures de temps de parcours d'événements de réflexion primaire. Dans la sismologie de réflexion, des ondes acoustiques sont introduites dans le sol, généralement par l'activation d'une source sismique. Des récepteurs acoustiques détectent les ondes acoustiques après leur réflexion à partir de strates et d'interfaces souterraines. Une analyse des ondes acoustiques est utilisée, en association avec la position connue de la source et du récepteur, pour produire une image du sous-sol. Des réflexions multiples de surface apparaissent dans des données sismiques lorsque des ondes acoustiques se déplaçant d'une source sismique à un récepteur sismique sont réfléchies vers le haut par des changements souterrains de l'impédance acoustique au moins deux fois et vers le bas à partir de la surface au moins une fois. Par ailleurs, on entend par réflexions primaires des événements sismiques qui ne subissent q'une réflexion à partir d'un changement d'impédance acoustique. En général, des réflexions multiples de surface dans un ensemble de données sismiques sont considérés comme des bruits indésirables qui interfèrent avec le signal souhaité et le compliquent. On a consacré beaucoup d'efforts à la conception de l'acquisi- tion de données sismiques et au traitement de données sismiques pour limiter l'impact des réflexions multiples sur les profils sismiques traités et obtenus en sortie. Même ainsi, dans de nombreuses zones, la qualité des données sismiques est abaissée, parfois sensiblement, par la présence de réflexions multiples. Divers procédés ont été essayés dans l'art antérieur pour éliminer les réflexions multiples de surface de traces enregistrées. Pour des situations géométriques simples, les temps d'arrivée de réflexions multiples de surface peuvent être calculés par voie analytique. Par exemple, lorsque le fond de l'eau et les horizons réfléchissants, ou miroirs, sont plats, le temps d'arrivée pour une réflexion primaire est donné par la relation bien connue 2 =t20 + 2 (1) V2 o t est le temps d'arrivée avec un déport horizontal source-récepteur de x, to est le temps d'arrivée pour un déport source-récepteur nul et Ve est une vitesse effective. L'équation (1) est valable pour des valeurs du déport source-récepteur qui sont communément utilisées dans l'acquisition de données sismiques. La réflexion multiple à partir du même miroir possède un temps de parcours à déport nul qui est double de celui de la réflexion primaire, et le temps d'arrivée de la réflexion multiple est donné par 2 = 4 t20 + x2 (2) V2 Lorsque les réflecteurs souterrains ne sont pas uni- formément plats et dans une géologie plus complexe, on utilise des principes de traçage de rayons complexes qui génèrent des multiples synthétiques et les soustrait de l'enregistrement sismique réel pour obtenir un enregistre- ment qui est supposé sans multiple. Cependant, ces procédés sont très incommodes par le fait qu'ils nécessitent d'importantes connaissances de la structure sous-marine ainsi que de la configuration du fond de l'océan avant que les multiples synthétiques puissent être générés. Des multiples synthétiques similaires peuvent être générés en utilisant des procédés plus précis n'impliquant pas directement un lancer de rayons, par exemple des techniques de propagation de champs d'ondes, mais ils exigent également une connaissance détaillée au moins du fond de l'océan, ainsi que de la forme des interfaces sous-marines, et ils ne sont donc pas aussi commodes qu'on le souhaite- rait. La publication de brevet WO 97/1130 décrit un procédé pour atténuer des multiples liés à l'eau dans des données sismiques. Le procédé repose sur un processus efficace de lancer de rayons pour déterminer des temps de parcours pour chaque multiple. Les caractéristiques des données sismiques sont utilisées pour estimer la forme d'onde du multiple devant être atténué. Le temps de parcours et la forme d'onde estimée sont ensuite utilisés pour atténuer le multiple par rapport aux données d'origine. Le procédé est exécuté dans un système d'équations non linéaires qui sont résolues en utilisant des techniques bien connues et une fonction supposée ou mesurée caractérisant le fond marin. Dans une forme appréciée de réalisation, après avoir été estimés, des multiples d'ordre inférieur sont utilisés pour estimer les multiples d'ordre supérieur, réduisant ainsi le temps de calcul. Le brevet des États Unis d'Amérique n 5 587 965 délivré à Dragoset et Jericevic décrit un procédé basé sur l'équation de l'onde acoustique et l'extrapolation du champ d'onde sismique à partir d'observations enregistrées. Le procédé décrit ici peut demander un long temps de calcul pour des levers sismiques importants. Le brevet des États Unis d'Amérique n 4 887 243 délivré à Pann décrit un procédé de génération des multiples dans une trace en combinant des paires de traces réelles, lesquelles paires ont un point extrême en commun et leurs autres extrémités identiques aux positions de source et de récepteur de la trace à laquelle on s'intéresse. L'étape de combinaison est exécutée par une convolution des traces appariées, comme cela est bien connu dans la technique. Toutes les paires possibles de multiples qui, ainsi, "s'additionnent" à la trace à laquelle on s'intéresse sont examinées conformément au principe de Huygens, selon lequel l'une de plusieurs ondes possibles dont le temps de parcours est le plus court dans tout système réel est celle qui apparaît réellement. En conséquence, la paire combinée de traces réelles avec l'onde ayant le temps total de parcours le plus court équivaut à la trajectoire d'onde réellement prise par le multiple. Dans une forme appréciée de réalisation, les traces synthétiques sont "sommées" pour générer une trace synthétique du temps de parcours minimal. Les traces soumises à la convolution qui sont proches du point de réflexion réel s'additionnent de façon cohérente tandis que celles s'éloignant davantage s'additionnent de façon incohérente, donnant une réalisation approximative d'une trace multiple. La trace multiple synthétique ainsi générée peut être soustraite de la trace soumise à une analyse pour donner une trace sans multiple. La génération d'une trace unique ayant des multiples exige donc la convolution d'un grand nombre de traces séparées puis la sommation de leurs sorties respectives. La raison pour laquelle un grand nombre de traces séparées doivent être mises en convolution est d'obtenir un nombre suffisant de traces soumises à une convolution pour que la sommation soit efficace. Ce procédé demande donc un calcul intense. Dans de nombreuses applications, il suffit d'identi- fier des réflexions multiples uniquement d'après leur cinématique (c'est à dire leurs temps de parcours). Une fois que l'identification a été réalisée, des principes de filtrage appropriés peuvent être utilisés pour filtrer les multiples. En variante, un interprète examinant une section sismique contenant une réflexion multiple identifiée ignore simplement la réflexion multiple identifiée à des fins d'interprétation. Pour de telles applications, il est souhaitable de disposer d'un procédé de calcul rapide pour prédire la cinématique d'événements de réflexions multiples. L'invention satisfait à cette nécessité. L'invention concerne un procédé pour calculer les temps de parcours d'événements sismiques réfléchis de façon multiple dans un profil sismique ou un modèle pourvu que des temps de parcours pour les événements primaires soient disponibles. Les temps de parcours pour les événements primaires peuvent être obtenus par un lancer de rayons à travers un modèle de vitesse du sous-sol ou par la prise de temps de parcours à partir d'horizons de réflexion primaire, ou miroirs primaires, présents dans un ensemble de données sismiques réelles. Les temps de parcours pour les événements primaires définissent une table de temps de parcours pour diverses combinaisons de positions de sources et de récepteurs. Le temps de parcours pour un multiple correspondant à des positions de sources et de récepteurs spécifiées est ensuite déterminé par sommation des temps de parcours pour des événements primaires à partir de la source jusqu'à un point intermédiaire et à partir du point intermédiaire jusqu'au récepteur: la position correcte du point intermédiaire correspondant à une trajectoire sismique de réflexion multiple réelle à partir de la source jusqu'au récepteur est celle qui, conformément au principe de Fermat, possède la valeur minimale du temps de parcours sommé. Ce procédé est applicable à tous les événements primaires, en sorte qu'on peut déterminer les temps de parcours de multiples souterrains aussi bien que de multiples du fond de l'eau. L'invention sera décrite plus en détail en regard des dessins annexés à titre d'exemples nullement limitatifs et sur lesquels: la figure 1 est une illustration schématique de l'ac- quisition de données sismiques marines montrant des exemples de réflexions à partir du fond de l'eau et de miroirs sous-jacents; la figure 2 illustre les trajectoires sismiques impli- quées dans une réflexion multiple; la figure 3 illustre l'utilisation du principe de réciprocité pour remplir les tables de temps de parcours; la figure 4 est un organigramme des étapes principales de la présente invention; la figure 5 est un organigramme détaillé montrant la façon de déterminer des temps de parcours pour des réflexions multiples; et la figure 6 illustre l'utilisation d'une ouverture limitée de recherche pour le point de réflexion intermé- diaire dans un arpentage sismique dans trois dimensions. En référence à présent à la figure 1, celle-ci montre une illustration schématique d'une prospection sismique marine. Une masse d'eau 5 ayant une surface supérieure 10 et un fond 20 est représentée. Un navire 12 de prospection sismique entraîne dans la masse d'eau un câble 16 d'une flûte qui porte des détecteurs sismiques 18. Pour une acquisition dans un environnement marin, les détecteurs communément utilisés sont des hydrophones qui sont sensibles à des pulsations de pression dans l'eau. Une source sismique 14 destinée à émettre une onde acoustique dans l'eau est déployée à partir du navire 12. La source peut être un canon à air comprimé, un vibrateur marin ou d'autres types de sources connus des spécialistes de la technique. On peut également utiliser des dispositifs de ces sources. L'onde acoustique émise descend le long de plusieurs trajectoires et est réfléchie à partir des interfaces entre des matériaux d'impédance acoustique variable tels que le fond 20 de l'océan, et elle est réfléchie en remontant pour être captée par les détecteurs 18. L'onde s'enfonce également dans le fond de la mer et est réfléchie à l'interface 22. Pour simplifier l'illustra- tion, 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 (miroirs). Les trajectoires des ondes représentées en 34 et 32 sont appelées ondes directes ou primaires, car la trajectoire ne comprend qu'une réflexion à partir d'une interface dans l'océan ou le fond de la mer avant une détection par un détecteur. Une réflexion "multiple" est également représentée par la trajectoire sismique 36 sur la figure 1, qui est réfléchie deux fois à partir du fond 20 de l'océan et une fois à partir de la surface 10 de l'océan. Il existe également d'autres types de multiples dans lesquels les trajectoires correspondent à des réflexions à partir de miroirs tels que 22 ou 24 et également à partir du fond 20 de l'eau, ainsi que des réflexions à partir de plus d'un miroir tel que 22 et 24. Il s'agit de trajectoires sismiques selon lesquelles l'onde arrive à la surface plus d'une fois. Comme décrit précédemment, ces multiples constituent une source d'aggravation importante du traitement des données sismiques car ils masquent des événement sismiques primaires réels qui ont lieu dans un temps comparable à celui mis par l'onde soumise à une réflexion multiple pour atteindre le même hydrophone et ils constituent une source de données parasites. En consé- quence, il est souhaitable de disposer d'un procédé pour identifier ces événements sismiques multiples dans des enregistrements sismiques. La description ci-dessus de réflexions multiples dans un arpentage sismique marin entend n'être qu'un exemple de l'apparition de réflexions multiples dans une prospection sismique. Des réflexions multiples ont lieu aussi bien dans des prospections sismiques terrestres et la présente invention peut être utilisée tout aussi bien pour déterminer des temps de parcours de réflexions multiples dans des arpentages sismiques terrestres. Un ensemble de données sismiques dans deux dimensions dans le stade de pré-sommation consiste en un grand nombre de traces qui sont caractérisées chacune par la position de la source qui a généré une impulsion acoustique et la position du récepteur (ou du détecteur) qui a détecté et enregistré la réponse du sol à cette impulsion. Les positions de la source et du récepteur le long d'une ligne sismique sont spécifiées au moyen d'une seule coordonnée, appelée le numéro de poste. Ainsi, une trace peut être spécifiée par une paire ordonnée de nombres (Si, Rj), o Si signifie que la source correspondant à la trace était placée au numéro de poste i, et Rj signifie que le récepteur était placé au numéro de poste j. Il convient de noter que le cas spécial selon lequel i = j signifie une trace dite à déport nul, c'est à dire que la source et le récepteur sont dans la même position. Dans la présente invention, les traces dans un ensem- ble de données sismiques sont analysées pour identifier les temps de parcours d'un ou plusieurs événements de réflexion primaire, tels que celui correspondant à la trajectoire sismique 32 décrite ci-dessus en référence à la figure 1. Les données de temps de parcours pour la n-ième réflexion primaire peuvent être représentées par une table de temps de parcours tpn (Si, Rj). Ces données sont simplement des temps de parcours agencés dans une table ou une matrice à deux dimensions. Chaque rangée de la table est associée à un numéro de poste de source particulier et chaque colonne à un numéro de poste de récepteur particulier. Le tableau 1 donne un exemple d'une telle table de temps de parcours. Tableau 1: Table des temps de parcours pour une réflexion primaire n 1 Enreg. - 1 2 3 4 5 6 7 8 Source I... 1 1,5 1,503 1,513 1,530 1,552 1,581 1,616 1,655 2 1,503 1,5 1,503 1,513 1,530 1,552 1,581 1,616 3 1,513 1,503 1,5 1,503 1,513 1,530 1,552 1,581 4 1,530 1,513 1,503 1,5 1,503 1,513 1,530 1,552 1,552 1,530 1,513 1,503 1,5 1,503 1,513 1,530 6 1,581 1,552 1,530 1,513 1,503 1,5 1,503 1,513 7 1,616 1,581 1,552 1,530 1,513 1,503 1,5 1,503 8 1,655 1,616 1,581 1,552 1,530 1,513 1,503 1,5 L'invention sera le mieux expliquée en considérant un exemple spécifique. La figure 2 montre un multiple générique de surface marine du premier ordre dont la trajectoire allant de la source Si en 114 jusqu'au récepteur Rj en 118 est SiP1WkP2Rj o Pl et P2 sont des points de réflexion vers le haut à partir des miroirs 132 et 134, et Wk est le point de réflexion vers le bas à la surface de l'eau. Sans lancer de rayons, les positions réelles des points P1, Wk et P2 ne sont pas connues. Cependant, on connaît deux facteurs essentiels concernant cet événement multiple. Premièrement, conformément au principe de Fermat, la trajectoire réelle du multiple, quel qu'il puisse être, est la trajectoire qui, parmi toutes les trajectoires possibles, possède le temps de parcours minimal. Deuxièmement, il ressort clairement de la figure 1 que le temps de parcours du multiple est la somme des temps de parcours des deux réflexions primaires qui constituent le multiple. Autrement dit, tml(Si,Rj) = tpl(Si,Wk)+tp2(Wk,Rj) o tm1 désigne le temps de parcours du multiple. Par conséquent, si la position du numéro de poste k était connue, le temps de parcours du multiple pourrait être obtenu simplement par la sommation des temps de parcours des deux événements primaires qui constituent le multiple. On obtient les temps de parcours pour les deux événements primaires simplement en regardant dans les tables de temps de parcours pour les deux horizons réfléchissants ou miroirs. L'invention détermine le temps de parcours du multiple en, premièrement, déterminant et sauvegardant les temps de parcours multiples pour chaque position possible du point Wk et, deuxièmement, sélectionnant le plus petit de ces temps de parcours déterminés. Conformément au principe de Fermat, ce résultat doit être le temps de parcours réel du multiple et la position correspondante de Wk est le point réel de réflexion à partir de la surface pour le multiple. Les spécialistes de la technique reconnaîtront que l'énoncé du principe de Fermat donné ci-dessus couvre la grande majorité des temps de parcours de réflexion, et la situation dans laquelle le temps de parcours est un maximum n'est habituellement pas importante. Pour la construction des tables de temps de parcours, l'invention comprend une étape additionnelle qui peut être nécessaire, utilisant le principe de réciprocité. Ceci peut être le mieux compris en se référant à la figure 3 sur laquelle une réflexion multiple à partir de la source 214 vers la position 218 du récepteur, avec une réflexion double sur le miroir 232, est montrée. Le point de réflexion en surface, dans ce cas, est une position 222 et il est clair que, dans une acquisition sismique marine dans laquelle la source et les récepteurs sont remorqués derrière un navire (comme dans la direction indiquée sur la figure 1), il n'y a pas de données ou de temps de parcours de réflexion primaire correspondant à une source dans la position 214 et à un récepteur dans la position 222. Dans un tel cas, le principe de réciprocité indique que le temps de parcours de réflexion pour une source dans la position 214 et un récepteur en 222 est le même que pour une source dans la position 222 et un récepteur dans la position 214. Lorsque le navire navigue de la droite vers la gauche, ces données sont enregistrées à la suite des données pour lesquelles la source est dans la position 214. Selon un aspect de l'invention, le principe de réciprocité est utilisé pour remplir des lacunes dans la table des temps de parcours qui résultent simplement de la géométrie d'acquisition. Un second aspect de l'invention implique la situation dans laquelle, pour certaines combinaisons de positions de sources et de récepteurs, on ne dispose d'aucune donnée pour le temps de parcours d'une réflexion primaire. Dans un tel cas, la table des temps de parcours est remplie au moyen d'une valeur importante qui est supérieure à toutes valeurs auxquelles on s'attend. Étant donné que des temps de parcours multiples sont toujours pris en choisissant la plus petite des sommes des temps de parcours d'événements primaires, ceci assure que les temps de parcours manquants ne contribuent aucunement aux résultats obtenus finalement. Les spécialistes de la technique reconnaîtront qu'un temps de parcours pourrait manquer du fait que, dans une géologie compliquée, certains miroirs peuvent ne pas posséder de trajectoire connectant chaque position source-récepteur. En référence à présent à la figure 4, on décrira les étapes de la présente invention. L'objectif de l'invention, comme décrit précédemment, est de déterminer des temps de parcours pour des réflexions multiples entre plusieurs positions de sources et de récepteurs lorsqu'on dispose de temps de parcours pour des événements de réflexion primaire. En 300, une table de temps de parcours pour des événements de réflexion primaire est formée par compila- tion. Un exemple d'une telle table de temps de parcours a été décrit cidessus en référence au tableau 1. Les temps de parcours qui entrent dans cette table pourraient être obtenus par de nombreux procédés, incluant l'exécution de prélèvements sur une section sismique ou par un procédé de modélisation utilisant un modèle de vitesse. De tels procédés sont bien connus des spécialistes de la technique. Ensuite, la table des temps de parcours est soumise à une extrapolation 320 lorsque cela est nécessaire en utilisant le principe de réciprocité. Celui-ci a été décrit ci-dessus en référence à l'acquisition de données sismiques marines o les positions des récepteurs sont nécessairement sur un côté de la source. Le point jusqu'auquel on effectue l'extrapolation dépend de la géométrie souterraine; les spécialistes de la technique reconnaîtront qu'avec un fond de l'eau en pente, le point de réflexion de l'un des événements primaires pourrait se trouver du côté de la source opposé à celui de la position du récepteur. Dans certains cas, l'espacement entre les positions des sources dans un arpentage sismique peut être différent de l'espacement entre les positions des récepteurs. Pour utiliser le procédé de l'invention, il est souhaitable que les deux espacements (c'est à dire l'espacement des sources et l'espacement des récepteurs) soient égaux. Si les s espacements sont différents, la table des temps de parcours est soumise à une interpolation 320 afin que les deux espacements soient égaux. Selon l'invention, des procédés d'interpolation tels qu'une interpolation linéaire ou une interpolation par polynôme peuvent être utilisés. Un cas spécial d'une interpolation par polynôme qui peut être utilisée est un spline cubique. De tels procédés d'interpolation sont connus des spécialistes de la technique et ne seront pas davantage décrits. Comme indiqué précédemment, dans certains cas, des valeurs manquantes des temps de parcours peuvent devoir être complétées. Ceci est réalisé, si cela est nécessaire, et une table complète de temps de parcours est construite, en 340, pour des événements primaires pour un ensemble choisi de positions de tir. En référence à présent à la figure 5, on décrira la détermination des temps de parcours pour des multiples. Des positions de tir et de récepteur sont sélectionnées en 350. Pour les positions sélectionnées de tir et de récepteur, un point de réflexion intermédiaire est sélectionné en 360 comme décrit ci-dessus au moyen de la référence 122 sur la figure 2. Pour le point de réflexion intermédiaire, la somme des temps de parcours pour les événements primaires à partir du point de réflexion intermédiaire jusqu'à la source et jusqu'au récepteur est déterminée en 370. Cette étape implique l'utilisation d'une table de temps de parcours primaires, dont un exemple est donné ci-dessus dans le tableau 1. Lorsque l'interface de réflexion souterraine pour le multiple est la même pour les deux branches du multiple, on peut alors voir qu'une seule table, telle que le tableau 1, suffit. Cependant, si les interfaces de réflexion sont différentes, par exemple si l'une est une réflexion au fond de l'eau et l'autre est une interface de réflexion souterraine, des données provenant de deux tables sont nécessaires. Les temps de parcours sommés sont délivrés en sortie en 385, en même temps qu'une information de position concernant le point de réflexion intermédiaire, et un contrôle est réalisé pour déterminer si toutes les valeurs du point de réflexion intermédiaire ont été évaluées. Si tel n'est pas le cas, un autre point de réflexion intermédiaire est sélectionné en 390 et le processus revient en 360 pour la détermination d'un temps de parcours obtenu par sommation. Une fois que toutes les valeurs raisonnables du point de réflexion intermédiaire ont été évaluées, le temps de parcours minimal pour les valeurs de sortie 385 est déterminé et ce temps est le temps de parcours pour les positions de tir et de récepteur ayant l'espacement sélectionné. Cette valeur est sauvegardée dans une table de temps de parcours de sortie et un contrôle est réalisé pour déterminer si toutes les combinaisons sources- récepteurs ont été traitées. Si la réponse est non, une nouvelle paire tir-récepteur est sélectionnée en 420 et le processus est répété à partir de 350. Une fois que toutes les paires tir-récepteur ont été traitées, une table complète de temps de parcours est délivrée en sortie en 425. Le processus est répété pour d'autres réflexions primaires jusqu'à ce que les multiples pour les événements de réflexion primaire souhaités aient été évalués. Ces multiples sont appelés multiples de premier ordre. Les spécialistes de la technique reconnaîtront que le processus entier décrit en référence aux figures 4 et 5 peut être exécuté en utilisant, au lieu de seulement des tables de temps de parcours primaire, une combinaison de tables de temps de parcours primaire et de multiples du premier ordre pour donner des tables de temps de parcours pour des multiples du second ordre. Ce processus peut en outre être répété si cela est souhaité en utilisant les temps de parcours pour les multiples primaires, du premier ordre et du second ordre dans une combinaison quelconque. Il est prévu que toutes ces applications entrent dans le cadre de la présente demande. On donnera maintenant un exemple de la façon dont le minimum est déterminé sur la figure 5. À cet effet, les données du tableau 1 sont utilisées en tant que première table de temps de parcours primaires et les données du10 tableau 2 ci-dessous sont utilisées en tant que seconde table de temps de parcours primaires. Tableau 2: Table de temps de parcours pour une réflexion primaire n 2 Enreg. 1 2 3 4 5 6 7 8 Source _ _ _ 1 2,5 2,502 2,507 2,518 2,532 2,550 2,571 2,596 2 2,502 2,5 2,502 2,507 2,518 2,532 2,550 2,571 3 2,507 2,502 2,5 2,502 2,507 2,518 2,532 2,550 4 2,518 2,507 2,502 2,5 2,502 2,507 2,518 2,532 4 2,532 2,518 2,507 2,502 2,5 2,502 2,507 2,518 6 2,550 2,532 2,518 2,507 2,502 2,5 2,502 2,507 7 2,571 2,550 2,532 2,518 2,507 2,502 2,5 2,502 8 2,596 2,571 2,550 2,532 2,518 2,507 2,502 2,5 Le tableau 3 donne la somme de temps de parcours primaires depuis une source dans la position 1 jusqu'à une position de réflexion intermédiaire et depuis la réflexion de position intermédiaire jusqu'à un récepteur dans la position 8. Ceci est obtenu en utilisant les données provenant des tableaux 1 et 2. Tableau 3: Somme des temps de parcours de la source dans la position 1 jusqu'à une position intermédiaire et de la position intermédiaire jusqu'à un récepteur dans la position 8: Position 1 2 3 4 5 6 7 8 intermé- diaire Temps 4,096 4,074 4,063 4,062 4,07 4,088 4,118 4,155 Dans le tableau 3, les temps de parcours obtenus par sommation correspondant aux positions intermédiaires différentes sont énumérés. On voit que le temps de parcours minimal est de 4,062 et qu'il apparaît lorsque le point de réflexion intermédiaire est dans la position 4. L'exemple montré ci-dessus est donné uniquement à titre illustratif: dans l'invention proprement dite, des positions de sources et de récepteurs beaucoup plus nombreuses seraient présentes et des positions de réflexions intermédiaires beaucoup plus nombreuses seraient évaluées. En outre, les calculs pourraient être effectués avec un degré de précision plus élevé (dans les tableaux 1 à 3, les données ont été arrondies à la milliseconde la plus proche, c'est à dire 0,001 seconde). Facultativement, une interpolation des temps de parcours sommés déterminés peut être effectuée en utilisant des procédés connus pour obtenir le temps de parcours sommé minimal. Un autre aspect de l'invention réside dans l'utilisa- tion d'une ouverture limitée pour la recherche du point de réflexion intermédiaire. En réduisant le nombre de consultations de table nécessaires, ceci réduit le temps global de calcul pour la détermination des temps de parcours multiples. Dans un cas à deux dimensions, tel que celui décrit ci-dessus, la réduction du temps global de calcul peut ne pas être notable. Cependant, dans l'acquisition de données sismiques en trois dimensions, ceci peut entraîner une réduction importante du temps de calcul. Ceci est illustré en référence à la figure 6 o une vue en plan d'un lever sismique en trois dimensions est représenté. Plusieurs lignes de récepteurs 501ae 501b,... 501n, ayant chacune des récepteurs sismiques dans des positions 520, sont représentées. Un exemple d'une ligne 510 de sources est représenté, la ligne de sources ayant des positions de sources 511a, 511b,... 511n. La mise en place des positions des sources et des récepteurs suivant des lignes droites n'entend pas être limitative de la présente invention, non plus que l'orientation des lignes de sources par rapport aux lignes de récepteurs. Les spécialistes de la technique reconnaîtront qu'avec un grand nombre de positions de sources et de récepteurs, la construction et l'accès aux tables de temps de parcours lors d'une consultation de table peuvent devenir longs. Pour limiter le temps de calcul, dans une forme de réalisation de la présente demande, seule une partie spatiale limitée de la matrice de positions sismiques dans deux dimensions est utilisée. Par exemple, lorsqu'on souhaite déterminer un temps de parcours multiple à partir d'une position de source spécifique telle que la position 511a jusqu'à un récepteur spécifique tel que le récepteur 530, la recherche pour le point de réflexion intermédiaire est limitée à une région telle que celle indiquée en 550. La sélection de cette région de recherche, ou ouverture pour le point de réflexion intermédiaire, est basée sur la connaissance de l'orientation des interfaces réfléchissan- tes souterraines. Une connaissance précise de l'orientation des interfaces réfléchissantes souterraines n'est pas nécessaire. Sur la figure 6, cette région de recherche est représentée uniquement à titre illustratif sous la forme d'un rectangle; d'autres formes arbitraires pourraient être utilisées. Facultativement, dans un lever à trois dimensions, une interpolation à deux dimensions des temps de parcours sommés peut être effectuée pour l'obtention du temps de parcours minimal. L'interpolation à deux dimensions est un spline bicubique dans une forme appréciée de réalisation. Cette technique d'interpolation est connue des spécialistes de la technique et ne sera pas davantage décrite. Il va de soi que de nombreuses modifications peuvent être apportées au procédé décrit et représenté sans sortir 0 du cadre de l'invention. REVENDICATIONS 1. Procédé de détermination d'au moins un temps de parcours multiple sur un lever sismique à partir d'une position spécifiée d'une source (14) jusqu'à une position spécifiée d'un récepteur (18), le lever sismique portant plusieurs réflexions primaires, une position de source, une position de récepteur, une interface réfléchissante (22) et un temps de parcours étant associés à chacune de ces réflexions primaires, le procédé étant caractérisé en ce qu'il comprend: (a) la détermination d'un premier temps de parcours pour une réflexion primaire depuis une source si- tuée dans la position spécifiée de la source jus- qu'à un récepteur situé dans une position inter- médiaire; (b) la détermination d'un second temps de parcours pour une réflexion primaire depuis une source si- tuée dans la position intermédiaire jusqu'à un récepteur situé dans la position spécifiée du ré- cepteur; (c) l'obtention d'un temps de parcours sommé en tant que somme du premier temps de parcours et du se- cond temps de parcours; (d) la répétition des étapes (a) à (c) pour plusieurs positions intermédiaires afin de produire plu- sieurs temps de parcours sommés; et (e) la détermination du temps de parcours pour la réflexion multiple en tant que minimum de la plu- ralité des temps de parcours sommés. 2. Procédé selon la revendication 1, caractérisé en ce que les positions spécifiées de la source et du récepteur se trouvent dans une masse d'eau (5). 3. Procédé selon la revendication 1, caractérisé en ce que la pluralité des temps de parcours de réflexion primaire et les temps de parcours associés sont stockés dans une table de temps de parcours, et la détermination des premiers et seconds temps de parcours comprend en outre une consultation de table. 4. Procédé selon la revendication 2, caractérisé en ce qu'au moins l'une des interfaces réfléchissantes associées au premier ou au second temps de parcours primaire est le fond (20) de la masse d'eau. 5. Procédé selon la revendication 1, caractérisé en ce qu'il comprend en outre l'utilisation du principe de réciprocité pour obtenir le temps de parcours associé à une position d'une source et à une position d'un récepteur. 6. Procédé selon la revendication 1, caractérisé en ce que le lever sismique est un lever à trois dimensions et le point de réflexion intermédiaire est à l'intérieur d'une ouverture spatiale limitée du lever à trois dimensions. 7. Procédé selon la revendication 6, caractérisé en ce que la détermination du minimum des temps de parcours sommés comprend en outre l'utilisation d'une technique d'interpolation sélectionnée entre (i) une interpolation par spline bicubique et (ii) une interpolation par polynôme. 8. Procédé selon la revendication 1, caractérisé en ce que la détermination du temps de parcours associé à une position d'une source et à une position d'un récepteur comprend en outre l'exécution d'une interpolation entre des temps de parcours disponibles proches de la source et du récepteur. 9. Procédé selon la revendication 8, caractérisé en ce que l'interpolation est choisie entre (i) une interpolation linéaire et (ii) une interpolation par polynôme. 10. Procédé selon la revendication 1, caractérisé en ce que la détermination du minimum de la pluralité des temps de parcours sommés comprend en outre l'utilisation d'une technique d'interpolation choisie entre (i) une interpolation par spline, et (ii) une interpolation par polynôme. 11. Procédé selon la revendication 1, caractérisé en ce que le temps de parcours multiple, au nombre d'au moins un, comprend en outre plusieurs temps de parcours multiples, le procédé comprenant en outre la répétition des étapes (a) à (e) pour donner plusieurs temps de parcours multiples déterminés du premier ordre. 12. Procédé selon la revendication 11, caractérisé en ce qu'il comprend en outre l'utilisation de la pluralité des temps de parcours multiples déterminés du premier ordre pour déterminer plusieurs temps de parcours multiples du second ordre.