Brevet FR2798737A1 — publié le 20010323

L'invention concerne l'exploration sismique, et plus particulièrement la tomographie de réfractions. L'exemple spécifique vise à la correction d'anomalies proches de la surface. Un problème important dans l'exploration sismique est celui posé par des distorsions d'images sismiques dues à des anomalies de vitesse proches de la surface. L'analyse de données de réfraction pour une estimation des statiques afin de révéler des informations de la zone proche de la surface et, de façon correspondante, des statiques, est donc très importante dans le traitement sismique. Des- anomalies de vitesse proches de la surface produisent de fortes distorsions dans des images sismiques. Si l'on connaissait la structure de ces anomalies, le meilleur moyen de traiter ce problème serait d'effectuer un calage à la référence de l'équation d'ondes ou une migration profondeur de l'équation d'ondes à partir de la surface à travers la structure connue. Cependant, la migration profondeur et le calage sur une surface de référence avant sommation, en trois dimensions, sont des défis de calcul. Par conséquent, des applications des statiques supposant une propagation des rayons cohérente en surface à travers la proche surface sont restées l'outil principal pour tenir compte des anomalies proches de la surface. L'analyse des réponses de présommation de données de réflexion fonctionne mieux pour l'estimation de la partie à courte période de la statique (dite statique résiduelle), tandis qu'une estimation de la statique de longue période utilisant des données de réflexion seules est très instable et inefficace. Par conséquent, l'analyse de données de réfraction pour une estimation de la statique afin de révéler davantage d'information concernant la zone proche de la surface et, de façon correspondante, concernant la statique, est courante dans le traitement sismique classique. La méthode du retour inverse généralisé a été largement appliquée aux données à deux dimensions (2D) . Malheureusement, pour la sismique à trois dimensions (3D) , la méthode du retour inverse est difficile à appliquer du fait du manque de données de retour inverse. L'idée de temps de retard est également apparue utile pour des calculs des statiques de réfractions 3D. La méthode des temps de retard supposait un modèle de la proche surface constitué de couches localement plates sur l'échelle de la plage de déports. Les arrivées premières étaient supposées être le commencement d'ondes ou dromochroniques se propageant le long d'interfaces réfringentes de ces couches. Les temps des pointés des arrivées premières étaient décomposés en temps de retard et en vitesses dans les couches réfringentes. Les temps de retard étaient ensuite convertis en épaisseurs de couche et en vitesses en supposant un angle critique d'incidence sur les couches réfrigentes. Une technique connue sous le nom de "Réfraction MISER & , disponible auprès de la cessionnaire de la présente demande constitue un exemple de la méthode des temps de retard. La méthode des temps de retard ou délais est également utilisée dans d'autres approches, telles que l'inversion linéaire généralisée (ILG) et la tomographie de réfractions d'ondes dromochroniques. Dans ces méthodes, au lieu d'une inversion en deux temps par des délais, des temps de parcours étaient inversés directement pour des épaisseurs de couche et des vitesses. Cependant, ces méthodes avaient également à traiter du problème de l'ambiguïté vitesse/profondeur. Pour aborder ce point, il était de pratique courante de fixer la vitesse dans la couche altérée avant une estimation de la profondeur. L'incorporation de données de réflexions pour une inversion combinée avec des données de réfractions pouvait parfois réduire les ambiguïtés de vitesse/profondeur. Malheureusement, les réflexions proches de la surface étaient très difficiles à pointer dans des données de production normales. L'un des problèmes principaux posés par les méthodes utilisant les ondes dromochroniques était que dans des zones de géologie complexe et de terrain accidenté, le modèle simple habituellement utilisé n'était pas suffisant pour expliquer d'importantes particularités des données. De plus, pour adapter la courbure non linéaire observée d'arrivées premières, on devait soit limiter la plage de déports, soit incorporer davantage de couches dans le modèle. Les deux approches rendaient l'inversion encore plus instable et plus ambiguë. on pensait que les amplitudes théoriques d'ondes réfractées ou dromochroniques étaient inférieures de plusieurs ordres de grandeurs à celles des ondes de volume et que, par conséquent, les ondes réfractées devaient être pauvrement représentées dans les données. Par ailleurs, même un faible gradient positif de vitesse verticale pouvait produire une variété d'ondes ou de modes d'interférence qui comprenaient des ondes mirages (tournantes) ayant des amplitudes très supérieures à celles des ondes réfractées. La cinématique des ondes mirages n'est pas linéaire. Dans la plupart des cas, une modélisation des ondes mirages s'adaptait mieux aux courbures observées des premières arrivées qu'une modélisation des ondes réfractées. Les raisons indiquées ci-dessus ont constitué la tomographie d'ondes mirages, en supposant que les arrivées premières sont des arrivées d'ondes mirages, une alternative très intéressante aux méthodes à ondes réfractées pour les calculs des statiques. Dans cette approche, le milieu était pararrétré sous la forme d'un certain nombre de cellules, des rayons mirages étaient tracés à travers le modèle, et les tenps de parcours résiduels étaient rétroprojetés ou inversés pour des perturbations de lenteur dans chaque cellule traversée par des rayons. Etant donné que les ondes mirages s'adaptaient mieux aux courbures non linéaires des arrivées premières que les ondes réfractées, des plages de déports plus larges pouvaient être incluses pour le traitement, réduisant l'ambiguïté dans l'inversion. Cependant, la limitation de l'ouverture de l'angle d'incidence au niveau de la couche altérée rendait encore ambiguë l'inversion pour les paramètres de cette couche. Un autre problème de la tomographie était posé par le très grand nombre d'inconnues, en particulier pour des levées en 3D. Ceci nécessitait des efforts spéciaux pour rendre l'inversion possible et stable. Des approches tomographiques précédentes utilisaient des ondes réfractées et des ondes mirages. Les méthodes basées sur les ondes réfractées étaient en général solides, car les paramètres des modèles dépendent à peu près linéairement des temps de parcours observés. Par contre, les méthodes basées sur les ondes mirages comprenaient une plage très large de déports, et la relation entre les paramètres des modèles et les temps de parcours devenait non linéaire du fait d'une dépendance notable du traitement vis-à-vis des trajectoires sismiques supposées. Ceci rendait les résultats tomographiques sensibles au modèle initial. La solution tomographique pouvait être stabilisée par l'utilisation d'un paramétrage de Fourier, mais la solution restait sensible au modèle initial du fait de l'utilisation d'un lancer de rayons. La tomographie de réfractions serait plus souhaitable si l'on pouvait éviter la dépendance du modèle initial vis- à-vis du lancer de rayons. Il a été démontré qu'une inversion des temps de parcours possède une formulation linéaire qui lui est propre dans le domaine tau-p. Cependant, la mise en pratique de la transformation en tau- p pour des temps de parcours n'était pas directe, en particulier dans le cas des données de levée en 3D. D'autres approches ont été proposées pour éviter un lancer de rayons en tomographie. Ces méthodes utilisaient une décomposition des temps de parcours dans l'espace limitée par l'équation de lleikonal. La méthode de lleikonal était basée sur un profil sismique vertical et une géométrie de transmission puits à puits. Quoique bien adaptées au profilage sismique vertical et aux géométries puits à puits, ces techniques étaient difficiles à appliquer à des données de réfractions. On peut donc voir que le paramétrage classique de cellules utilisé dans la tomographie était, d'après ce que l'on sait, généralement inefficace pour une modélisation de structures géologiques proches de la surface. En outre, il introduisait un grand nombre d'inconnues. Dans la formation basée sur des cellules, le problème tomographique était mal posé et nécessitait une régularisation qui n'était pas directe. En outre, la qualité du modèle initial nécessitait un lancer de rayons et une inversion non linéaire reposant sur des connaissances d'analyste, et ceci pouvait conduire finalement à une erreur systématique dans la solution finale. En bref, l'invention effectue une tomographie de réfractions de données sismiques en utilisant des délais différentiels. Conformément à l'invention, des données sismiques obtenues au cours d'un arpentage sismique à partir de plusieurs stations sismiques sont traitées pour éliminer des données l'effet de la statique. Des fonctions temps de trajet/déport sont formées pour les stations sismiques à partir des pointés d'arrivées premières obtenus à partir des données. Les fonctions temps de parcours/déport sont transformées en fonctions vitesse/profondeur pour dériver un modèle de la zone proche de la surface. Des statiques de longue période dans les données sont ensuite déterminées à partir du modèle de la zone proche de la surface dérivé. Des statiques de courte période sont estimées par une décomposition compatible en surface des temps de parcours résiduels. Les statiques déterminées de longue période et de courte période sont ensuite éliminées des données sismiques, lesquelles peuvent alors être utilisées pour l'affichage ou un traitement subséquent à des fins d'analyse. Les techniques de l'invention estiment le modèle de vitesse proche de la surface et produisent des corrections statiques au moyen d'une tomographie de temps de parcours d'arrivées premières sans lancer de rayons. L'inversion de temps de parcours exécutée conformément à l'invention implique plusieurs étapes. Premièrement, on sélectionne les données devant être inversées à partir d'une plage de déports spécifiée et on les ajuste, si cela est nécessaire. Ensuite, on estime une fonction 3D à partir de pointés d'arrivées premières en utilisant une méthode combinant une modélisation 1D et une régression non paramétrique. Ensuite, des délais différentiels sont utilisés pour calculer des fonctions locales vitesse/profondeur qui constituent le modèle de la zone proche de la surface. Puis des statiques de longue période sont calculées en utilisant le modèle dérivé. Des statiques de courte période sont estimées par une décomposition compatible en surface des temps de parcours résiduels. Les statiques de courte période et de longue période peuvent ensuite être éliminées des données sismiques. L'invention combine une approche d'une décomposition spatiale sans contrainte des temps de parcours avec une inversion linéaire dans le domaine tau-p. L'invention n'exécute aucun lancer de rayons explicite. L'approche utilisée combine la solidité des méthodes à ondes réfractées (délais), car elle ne nécessite pas un modèle initial ou un tracé de rayon, et la souplesse de la tomographie d'ondes mirages, car elle inverse à la fois les ondes réfractées et les ondes mirages sur toute la plage de déports. 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 un diagramme schématique d'un arpentage sismique illustrant la présence de statiques dans les levées ; la figure 2 est un diagramme schématique du parcours d'ondes réfractées à travers le sol pendant un arpentage sismique ; la figure 3 est un diagramme schématique du parcours d'ondes mirages à travers le sol pendant un arpentage sismique ; la figure 4 est une vue en plan d'une partie de la terre illustrant des emplacements de sources et de récepteurs pendant un arpentage sismique la figure 5 est un graphique comparatif de résultats de données obtenus avec la présente invention en comparaison avec des résultats obtenus avec un lancer de rayons ; les figures 6A, 6B et 6C sont des graphiques de données sismiques illustrant des résultats avant un traitement et après un traitement avec deux techniques différentes ; et la figure 7 est un organigramme montrant une séquence d'étapes pour le traitement de données sismiques selon l'invention. Sur les dessins, une vue en coupe transversale d'une couche altérée L au-dessous d'une partie 10 de la surface du sol E est montrée au voisinage d'un levé sismique souterrain. De l'énergie sismique est introduite dans le sol en un point de tir S, habituellement à une certaine profondeur convenable dans la couche altérée L. L'énergie sismique s'éloigne du point de tir S à travers la couche altérée L et de celle-ci à travers une base 12 ou une interface entre la couche altérée L et des formations souterraines plus profondes. L'énergie sismique s'enfonce dans le sol E et des parties de l'énergie sont réfléchies par des interfaces, telles que celles indiquées en 14, entre des formations géologiques souterraines indiquées de façon générale par la référence numérique 16. L'énergie réfléchie s'élève depuis les réflecteurs ou miroirs à travers les formations souterraines 16 et à travers la couche altérée L jusqu'à un détecteur ou un dispositif de détecteurs indiqué schématiquement en D, placé à la surface du sol 10 ou à proximité de celle-ci. Les réponses des réflecteurs, que se soit sur un dispositif à deux dimensions (comme représenté) ou un dispositif à trois dimensions, sont ensuite captées et traitées. On doit également comprendre que l'invention peut être aussi bien utilisée avec des données de levées de sismique marine. Le but de l'arpentage par réflexion est d'identifier des formations ou des particularités souterraines auxquelles on s'intéresse. Cependant, la couche altérée L n'est pas d'une épaisseur ou hauteur uniforme et la composition des matériaux et leur densité à l'intérieur de la couche L varient. Ces effets sont généralement appelés statiques. En général, il y a deux catégories de statiques. Des statiques de longue période sont généralement provoquées par des variations d'épaisseur et de vitesse de la couche altérée L au-dessus d'un niveau de référence 18. Des statiques de courte période sont dues à des variations de la couche altérée L et également à un couplage de surface de la source S et du détecteur D avec le sol. Dans le traitement de données sismiques, il est souhaitable de corriger ces statiques et d'éliminer autant que possible les effets de variations présentes dans la couche altérée L et d'autres statiques sur les données sismiques. Il est très souhaitable que les données sismiques présentent une forme comme si elles résultaient d'un arpentage exécuté sur une référence ou un plan sensiblement plat 18, à une hauteur constante dans le sol. Pour compenser les statiques, il est nécessaire de déterminer la valeur du temps de retard ou délai introduit par le parcours de l'énergie sismique au-dessus du niveau de référence 18, puis d'éliminer les effets de ce délai des données sismiques. En fait, une correction des statiques compense des différences de hauteur de la source S et du détecteur D au-dessus du niveau de référence 18, des différences de couplage par surface de la source et des récepteurs et des variations de densité des matériaux terrestres au-dessus du niveau de référence 18. La tomographie de réfractions est une méthode pour trouver la répartition des vitesses à partir des temps de parcours d'arrivées premières. La tomographie de réfractions est la formation de modèles de particularités souterraines basés sur la réponse des formations de terrain à des ondes sismiques de réfractions. Des ondes de réfractions sont formées lorsqu'une partie de l'énergie sismique arrive à une interface entre deux formations, à ce qui est connu sous le nom d'angle critique comme c'est le cas en un emplacement 15. L'angle critique est basé sur les vitesses sismiques relatives dans les couches au-dessus et au-dessous de l'interface particulière et de parcours le long d'une telle interface. L'onde dénergie sismique de réfraction se déplace le long de cette interface jusqu'à un emplacement 17 où l'onde de réfraction sort de nouveau sous le même angle critique pour s'élever et être captée au récepteur R. La détection de l'arrivée d'une telle onde de réfraction est connue sous le nom de pointé d'arrivée première. Il y a deux classes principales de méthodes de tomographie de réfractions. La première classe traite d'un modèle proche de la surface formé de couches localement plates sur l'échelle de la plage des déports, ou distance entre les emplacements 15 et 17. La figure 2 illustre schématiquement le trajet de parcours 20 d'une onde réfractée telle que celle décrite en regard de la figure 1. L'onde réfractée 20 se déplace d'une source S jusqu'à un détecteur D le long d'une interface 22 entre des formations 24 et 26. Le trajet de réfraction 20 de l'onde réfractée est représenté sous la forme d'une ligne droite sur la figure 2, bien qu'il suive en réalité l'interface 22. Des arrivées premières sont supposées être le commencement d'ondes réfractées se propageant le long d'interfaces de réfractions entre de telles couches. La cinématique des ondes réfractées est linéaire. Les temps de pointés d'arrivées premières sont décomposés en temps de retard ou délais et en vitesse dans les couches réfringentes. Les délais sont ensuite convertis en épaisseurs et vitesses de couche en supposant une certaine valeur pour un angle critique d'incidence sur les douches réfringentes. Une technique de traitement connue sous le nom de Réfraction MISER@ disponible auprès de la cessionnaire de la présente demande constitue un exemple de cette approche. Il y a également d'autres modes connus de réfraction d'ondes réfractées dans l'art antérieur. La seconde classe de la tomographie de réfractions porte sur ce qui est connu sous le nom d'ondes mirages ou ondes tournantes, comme montré dans un modèle (figure 3) dans les mêmes emplacements que ceux de la figure 2. Des ondes mirages présentent un gradient de vitesse vertical. Les arrivées premières au détecteur D en provenance de la source S sont supposées être le commencement d'ondes mirages ou tournantes, comme montré schématiquement en 28. La cinématique des ondes mirages n'est pas linéaire. Le milieu est paramétré sous la forme d'un certain nombre de cellules, et des rayons sont tracés à travers le modèle. Des résidus entre des temps de parcours observés et calculés sont rétroprojetés ou inversés pour des perturbations de lenteur dans chaque cellule traversée par les rayons. Un certain nombre de mises en application différentes de la tomographie de réfraction d'ondes mirages sont décrites dans la littérature. L'un des problèmes principaux avec ces méthodes est que des formations terrestres réelles sont plus complexes que montrées sur les f igures 2 et 3. Dans des zones ayant une géologie complexe et des terrains accidentés, le modèle simple habituellement utilisé ne suffit pas à l'explication de particularités de données importantes. De plus, pour ajuster la courbure non linéaire observée d'arrivées premières, on doit soit limiter la plage de déports, soit incorporer davantage de couches dans le modèle. Les deux approches rendent l'inversion encore plus instable et ambiguë. Cependant, des méthodes basées sur les ondes réfractées sont en général solides, car les paramètres du modèle dépendent à peu près linéairement des temps de parcours observés. Par contre, pour les méthodes basées sur des ondes mirages qui incorporent une plage de déports plus large, la relation entre les paramètres du modèle et les temps de parcours devient non linéaire du fait de sa dépendance notable vis-à-vis des trajets des rayons. Ceci rend les résultats tomographiques sensibles au modèle initial. Pour tenir compte de ce qui précède, le procédé de l'invention ne fait pas appel à un paramétrage explicite du modèle et il évite donc les problèmes mentionnés ci-dessus. L'invention obtient une mesure de statiques par une inversion des temps de parcours qui implique trois étapes. Dans la première étape, des données de temps de parcours sont divisées en un ensemble de plages secondaires de déports. Les temps de parcours dans chaque plage secondaire sont décomposés dans l'espace afin qu'un pointé de temps de parcours observé t (R,, Rd) pour la source et le détecteur avec des coordonnées R., et Rd, respectivement, soit modélisé par : où le déport est X, et X R.-Rd Ici, T(R, X) est la fonction temps de parcours modélisé pour un emplacement de surface R, le vecteur des coordonnées définissant l'emplacement de la surface. Les coordonnées R, et R2 correspondent aux points centraux entre les coordonnées R,, pour la source S et le point milieu, et la coordonnée Rd pour le détecteur D et le point milieu, respectivement, comme on le voit sur la figure 4. La décomposition des temps de parcours est basée, dans une forme appréciée de réalisation, sur une régression de troisième ordre et un hybride de 1D et sur une régression non paramétrique. Premièrement, pour chaque plage secondaire de déports, une grille spatiale à deux dimensions est définie. Puis les restes entre les pointés modélisés ajustés à une dimension sont rétroprojetés sur la grille. En général, plusieurs itérations de ce processus sont utilisées. En variante, la décomposition des temps de parcours peut également être effectuée par plusieurs itérations d'une application de t(R#;, Rd) sur T(R, X) en utilisant un filtre de Kalman en temps que récurrence de rétroprojection pour une mise à jour du modèle. Les propriétés.du filtre de Kalman dans une application à la tomographie sont décrites, par exemple, dans le document "Traveltime Tomography Using Regularized Récursive Least Squaresli-, 1998, Geophysical J. Int. , Vol. 134, pages 545 à 553 de Yao et al. , dont la Demanderesse est co-auteur. Un lissage spatial dans la décomposition résulte de l'utilisation d'une pondération gaussienne similaire à la régularisation de diffusion introduite comme décrit dans le document 113D Nonlinear Inversion by Entropy of Image Contrast Opmtimization", 1995, Nonlinear Processes in Geophysics, Vol. 2, pages 228 à 240, de Ryzhikov et al. La décomposition a pour résultat que les valeurs estimées de T(R, X) pour toutes les plages secondaires de déports forment une fonction temps de parcours/déport variant dans l'espace. Dans la deuxième étape de l'invention, des vitesses et temps de retard apparents sont calculés d'après la fonction estimée temps de parcours/déport, où Xi est le déport central pour la iiè,e plage secondaire de déports. et T,R est le temps de retard défini sous la forme Par conséquent, une fonction tau-p en 3D a été estimée à partir de la décomposition des temps de parcours des arrivées premières. Dans la troisième étape, les vitesses apparentes dérivées peuvent être utilisées pour calculer Z (ViR) # la profondeur de pénétration de l'onde réfractée, pour toutes les stations au moyen de la formule suivante : où x,=o et Z (V.R) =0. Cette f ormul e suppose un accroissement monotone de ViR avec le déport et, de f açon correspondante, avec la profondeur, ce qui a pour résultat d'ignorer les inversions de vitesse. Une autre formulation, plus avantageuse, utilise des temps de retard ou délais différentiels et des vitesses apparentes dérivées pour calculer Z (V,R) par une estimation séquentielle de l'épaisseur de la couche du haut vers le bas du modèle : où le temps de retard ou délai différentiel est Pour chaque procédé, la vitesse apparente Vi R représente la vitesse de l'intervalle instantané à la profondeur Z (V,R) correspondant à la profondeur à laquelle l'onde miroir tourne ou à la profondeur du réfracteur pour une onde réfractée. Les valeurs de Z (Vi5) pour toutes les plages secondaires de déports forment une fonction vitesse/profondeur au-dessous d'un point donné de la surface. Les fonctions vitesse /profondeur pour un ensemble de points de commande composent une estimation du modèle de la zone proche de la surface. Par conséquent, avec la présente invention, l'utilisation des propriétés géométriques des ondes réfractées, associée au modèle de vitesse dérivé, donne une reconstruction tomographique sans lancer explicite de rayons. La résolution spatiale de V,R est limitée par le lissage imposé sur la décomposition des temps de parcours, et dans la limite correspond à Xi/2. Le degré optimal de lissé dépend du niveau de bruit des données et peut être choisi en fonction d'une approche de régularisation basée sur des données. Un exemple de cette technique est décrit dans un article dont la Demanderesse est co-auteur : "Data Based Regularization in Traveltime Tomographyll, 199G, Osypov et Scales, 66th Annual International Mtg., Society of Exploratory Geophysicists, pages 2003 à 2006. Une fois dérivé, le modèle de vitesse est utilisé pour un calcul de statiques. La statique dans une station particulière est calculée sous la forme de la différence entre deux termes. Le premier terme est le temps de parcours vertical à travers le modèle de vitesse depuis la hauteur de la station en descendant jusqu'à une référence intermédiaire. Le second terme est le temps de parcours vertical depuis la référence intermédiaire en remontant jusqu'à la référence finale en utilisant une vitesse de remplacement. Ces statiques sont de période longue car elles correspondent au modèle de vitesse lissé. Les statiques de période courte sont estimées par une décomposition itérative compatible en surface des temps de parcours résiduels entre les temps de parcours observés et les temps de parcours modélisés des arrivées premières. Le procédé de l'invention a été appliqué à un arpentage en 3D à la dynamite dans la région de ceinture du bassin Nord-Ouest de l'Argentine. Cet arpentage comprend une commande supplémentaire portant sur la zone proche de la surface sous la forme de mesure de temps vertical pour douze trous de 150 m de profondeur. Des variations topographiques sont d'environ 150 m dans les 100 kM2 de la prospection. A des fins d'étude, aucune information provenant de l'expérience de gueule de trou n'a été utilisée dans les procédés appliqués. Les mêmes temps de pointés ont été également inversés en utilisant un algorithme de tomographie d'ondes mirages à lancer de rayons. La figure 5 montre des temps verticaux observés pour douze trous de 150 m en fonction des temps de parcours unidirectionnels pour les modèles de vitesse estimés par les deux méthodes tomographiques. Des barres d'erreurs reflètent les incertitudes des pointés. On peut voir que les valeurs estimées (les carrés) obtenues avec l'invention peuvent présenter une meilleure corrélation avec les temps verticaux que les résultats du procédé à lancer de rayons de 1 ' art antérieur (losanges) . Ce procédé à lancer de rayons est décrit dans un article de la Demanderesse : A Comparative Study Between 3-D Diving Wave Tomography and Head Wave Refraction Methods", 1998, Osypov, G8th Ann. Internat. Mtg., Society of Exploratory Geophysicists, pages 1222 à 1225. Cet article montrait que la tomographie d'ondes mirages à lancer de rayons donnait des résultats mieux adaptés à la commande en gueule de trou que des solutions basées sur des ondes réfractées. La figure GA est une composition brute en ligne de données sismiques avec l'application de statiques classiques de terrain. La forte corruption par des anomalies statiques dans une zone indiquée en 60 correspond à la partie de la composition utilisée pour l'étude comparative basée sur un algorithme de tomographie d'ondes mirages à lancer de rayons conformément à l'article de 1998 de la Demanderesse indiqué ci-dessus. La figure 6B est une composition en ligne des mêmes données que celle de la figure GA après l'application de la tomographie d'ondes mirages à lancer de rayons. La figure 6C montre des résultats provenant de données traitées conformément à la tomographie de réfractions de la présente invention. La qualité globale de la composition de la présente invention est comparable à celle pour la visualisation de la figure 6B, sinon meilleure. De plus, les temps des événements sur la partie de gauche de la composition de la figure 6C diffèrent d'environ 50 ms par suite de différences portant sur les valeurs estimées des statiques de longue période. Ces différences de statiques correspondent partiellement aux différences des estimations temps-profondeur sur la figure 5. Ceci est un exemple de l'ambiguïté de vitesse qui avait une influence sur les résultats du fait de la contamination par le bruit des données de déports rapprochés. Enfin, l'invention a pour avantage de ne nécessiter que quelques cycles d'essai avec différentes contraintes de lissage, en comparaison avec le résultat de lancers de rayons de la f igure 6B qui a été choisi comme étant le meilleur parmi un large ensemble de solutions obtenues à partir de différents modèles initiaux. La f igure 7 est un organigramme de processus d'une séquence d'étapes selon l'invention. Le résultat de ce traitement permet d'obtenir des résultats analogues à ceux de la figure 6C à partir de données traitées en appliquant des statiques de terrain classiques, telles que celles montrées sur la f igure 6A, et faisant contraste avec les techniques de l'art antérieur à lancer de rayons, demandant beaucoup de temps et beaucoup de calculs, les résultats étant montrés sur la figure 6B. Au cours d'une étape 70, des données sismiques de terrain pour une levée sismique dans trois dimensions sont enregistrées. Les données sismiques peuvent être enregistrées pendant un arpentage sismique terrestre ou un arpentage sismique marin dans une masse d'eau. Les autres étapes de traitement présentées ci-dessous peuvent être exécutées dans une plate-forme de calcul appropriée, de base Unix. Pour une décomposition des délais, dans des arpentages dans l'eau, les données sont ajustées en fonction des stations en étant déplacées vers le fond de l'eau en utilisant la vitesse de l'eau. Pour des arpentages terrestres, les données sont ajustées sur la surface libre en utilisant les temps à la gueule du trou ou temps verticaux. Des statiques de longue période et de courte période sont cependant déterminées sur la base de leurs positions d'origine. Au cours d'une étape 72, des temps de parcours de réfraction à la fois pour des ondes réfractées et des ondes mirages sont sélectionnés ou pointés. Ensuite, durant une étape 74, les temps de parcours de réfraction pointés sont décomposés dans le domaine espace-déport. Les données de temps de'parcours sont divisées en un ensemble de plages secondaires de déports et décomposées spatialement comme décrit précédemment. Comme indiqué à l'étape 76, une valeur estimée résultante d'une fonction tau-p à trois dimensions est ensuite formée de la manière décrite précédemment. Au cours d'une étape 78, des délais différentiels sont alors calculés et les vitesses apparentes associées sont déterminées. A partir des résultats de vitesses apparentes dérivés obtenus au cours de l'étape 78, on obtient au cours d'une étape 80 les valeurs Z (V,R de fonctions locales vitesse/profondeur au-dessous de chacun des points de surface auxquels on s'intéresse. Ceci peut être réalisé en utilisant l'une des méthodes décrites précédemment. Les fonctions vitesse/profondeur Z (V,R) pour un ensemble de points de commande sont ensuite assemblées au cours d'une étape 82 afin que des statiques de période longue soient calculées. De plus, au cours d'une étape 84 après l'étape de décomposition 72, les statiques résiduelles résultantes ou restantes sont utilisées pour former une estimation de statiques de période courte. Les statiques de période courte sont avantageusement estimées par une décomposition itérative compatible en surface des temps de parcours entre les temps de parcours des arrivées premières observés et modélisés. on doit comprendre qu'on pourrait utiliser d'autres méthodes en plus d'une rétroprojection pondérée de ce type. D'autres techniques d'estimation de statiques a court terme telles que les approches Gauss-Zeidel ou les approches par gradients conjugués peuvent toutes aussi bien être utilisées. Après une détermination des statiques de période longue au cours de l'étape 82 et des statiques de période courte au cour de l'étape 84, les statiques ainsi déterminées sont alors appliquées aux données enregistrées au cours d'une étape 86. Les données corrigées par les statiques résultant de l'étape 86 peuvent ensuite être utilisées pour un traitement classique subséquent et un affichage à des fins d'analyse comme indiqué à une étape 88. Les affichages ainsi formés sont indiqués à titre d'exemple sur la figure 6C et peuvent être produits par tout système ou poste d'affichage de données sismiques convenable, disponible dans le commerce. La tomographie de réfractions avec la présente invention est donc un procédé stable et souple de modélisation de la zone proche de la surface et de production de correction de statiques. Ce procédé d'inversion des temps de parcours diffère des approches tomographiques classiques par le fait qu'il n'effectue aucun lancer de rayons explicite. L'approche combine la fiabilité des procédés à temps de retard ou délais, car elle ne nécessite aucun modèle initial ou lancer de rayons initial, à la souplesse de la tomographie, car elle inverse à la fois les ondes dromochroniques et mirages sur toute la plage de déports. L'inversion des données de réfraction sur la plage complète des déports minimise l'ambiguïté vitesse/profondeur. Cependant, des données de réfraction proches des déports sont souvent bruyantes ou inutilisables. Dans ce cas, une information a priori, telle que des temps à la gueule du trou, des données concernant le puits, des données de gravité et toute connaissance géologique peuvent être utilisées avec la présente invention pour contenir les vitesses de la partie la moins profonde du modèle. L'un des problèmes de la tomographie est le traitement d'inversions de vitesse ou de zones "aveugles". En général, une inversion pour des zones de basse vitesse est ambiguë. Néanmoins, avec la présente invention, une épaisseur de couche de basse vitesse peut être estimée en utilisant des délais différentiels de la manière décrite ci-dessus. Pour toutes les méthodes de réfraction, la qualité des résultats dépend notablement de la qualité des pointés. Le procédé de l'invention produit des estimations statistiques des erreurs des modèles sur la base d'analyses des erreurs des pointés. Des estimations de ce type peuvent être obtenues, par exemple, d'après le signal de sortie d'un capteur de pointés d'arrivées premières a réseau neuronal classique. Des points fondamentaux, propres à toutes les méthodes utilisant des statiques de réfraction, comprennent la différence de vitesse entre des ondes réfractées et réfléchies, et la violation de la cohérence de surface dans les zones contenant un permagel épais ou de l'eau profonde sur le dessus de coulées boueuses. Le premier point est de en partie à l'écart des temps pointés à partir des commencements d'arrivées "réelles", et en partie aux différences entre des vitesses de groupes d'ondes d'interférence réfractées et des vitesses d'ondes réfléchies. Le point de la violation de la cohérence de surface peut être résolu par l'application d'une correction spéciale à base de rayons exécutée conformément à l'invention ou par l'utilisation d'un calage sur référence ou d'une migration de l'équation d'ondes à partir de la surface par l'intermédiaire du modèle tomographique. On peut voir d'après ce qui précède que la décomposition des temps de parcours dans trois dimensions conformes à l'invention est réalisée avec une régression à la fois paramétrique et non paramétrique, plutôt que selon une approche complètement paramétrique. Avec la présente invention, la décomposition des temps de parcours est plus souple et demande moins de calculs. La présente invention est mise en #uvre sur des données de réfraction de surface et n'est pas limitée à des géométries particulières uniquement présentes dans un profilage sismique vertical et des géométries puits à puits. De plus, les contraintes imposées à l'approche par l'équation de lleikonal de l'art antérieur sont évitées. En outre, l'invention réalise une inversion à une dimension dans le domaine tau-p. Elle inverse des délais différentiels plutôt que des délais complets (tau/2). Ceci réduit les erreurs qui peuvent être provoquées par des structures à trois dimensions. L'invention est donc mieux adaptée à une tomographie dans trois dimensions, utilisant la fonction tau-p à trois dimensions obtenue à partir de la décomposition des temps de parcours. En outre, la rétroprojection régularisée, préférée, évite des approches de programmation linéaires ou quadratiques plus coiâteuses. Bien que l'invention ait été décrite sous la forme d'étapes spécifiques, il va de soi que les spécialistes de la technique peuvent substituer des étapes comparables, des additions ou des modifications sans sortir du cadre de l'invention. REVENDICATIONS 1. Procédé de traitement de données sismiques obtenues pendant un arpentage sismique à partir de plusieurs stations sismiques pour éliminer l'effet de statiques des données, caractérisé en ce qu'il comprend : la formation de fonctions temps de parcours/déport pour les stations sismiques à partir de pointés d'arrivées premières obtenus à partir des données ; la transformation des fonctions temps de parcours/déport en fonctions vitesse/profondeur pour dériver un modèle de la zone proche de la surface la détermination de statiques de période longue dans les données provenant du modèle dérivé de la zone proche de la surface ; et l'élimination des statiques déterminées de période longue des données sismiques. 2. Procédé selon la revendication 1, caractérisé en ce que les temps de parcours résiduels restent à la suite d'une détermination de statiques de période longue, et en ce qu'il comprend en outre : la formation d'une estimation des statiques de période courte dans les données sismiques ; et l'élimination des estimations formées des statiques de période courte des données sismiques. 3. Procédé selon la revendication 2, caractérisé en ce que la formation d'une estimation des statiques de période courte comprend l'exécution d'une décomposition compatible en surface des temps de parcours résiduels. 4. Procédé selon la revendication 1, caractérisé en ce que la formation des fonctions temps de parcours/déport comprend la division des données sismiques en un groupe de plages secondaires de déports. 5. Procédé selon la revendication 4, caractérisé en ce que la formation de fonctions temps de parcours/déport comprend en outre l'exécution d'une décomposition spatiale des temps de parcours pour chacun des groupes de plages secondaires de déports afin d'obtenir des pointés de temps de parcours modélisés pour chacune des plages secondaires de déports. 6. Procédé selon la revendication 5, caractérisé en ce que l'exécution d'une décomposition spatiale comprend la détermination d'un pointé de temps de parcours observé t(R,Rd) pour une source avec des coordonnées vectorielles R,, et pour un récepteur avec des coordonnées vectorielles Rd, conformément au modèle oÙ T(R, X) est la fonction de temps de parcours modélisé pour un emplacement de surface R, et RI correspond à un point central entre la source et le point milieu, et R2 correspond à un point central entre le récepteur et le point milieu. 7. Procédé selon la revendication 5, caractérisé en ce que la formation de la fonction temps de parcours estimé/déport comprend la détermination d'une vitesse apparente à partir du pointé de temps de parcours pour chacune des plages secondaires de déports. 8. Procédé selon la revendication 7, caractérisé en ce que la détermination d'une vitesse apparente comprend la détermination d'une vitesse apparente ViR conformément à une fonction temps de parcours estimé/déport : où Xi est le déport central dans la i'è` des plages secondaires de déports. 9. Procédé selon la revendication 7, caractérisé en ce qu'il comprend en outre la détermination d'un délai pour chacune des plages secondaires de déports. 10. Procédé selon la revendication 9, caractérisé en ce que la détermination d'un délai comprend la détermination d'un délai -C,R où 11. Procédé selon la revendication 9, caractérisé en ce qu'il comprend en outre la formation d'une fonction tau-p à trois dimensions à partir de la décomposition des temps de parcours d'arrivées premières. 12. Procédé selon la revendication 7, caractérisé en ce que la détermination de statiques de période longue comprend la détermination de la profondeur de pénétration de l'onde réfractée sur la base de la vitesse apparente déterminée. 13. Procédé selon la revendication 12, caractérisé en ce que la détermination de la profondeur de pénétration comprend la détermination de la pénétration conformément à la relation suivante : où la valeur de Z (ViR) pour toutes les plages secondaires de déports forme une fonction vitesse/profondeur au-dessous d'un point de surface donné et représente une estimation d'un modèle d'une zone proche de la surface. 14. Procédé selon la revendication 7, caractérisé en ce que la détermination de statiques de période longue comprend la formation d'une estimation séquentielle d'épaisseurs de couches dans le modèle. 15. Procédé selon la revendication 14, caractérisé en ce que la formation d'une estimation séquentielle d'épaisseurs de couches AZ (V,R) pour une couche i est conforme à l'équation suivante : où 16. Procédé de prospection sismique pour obtenir des données sismiques avec des effets réduits sur la statique dans les données, caractérisé en ce qu'il comprend : l'introduction d'énergie sismique dans le sol à partir d'au moins une position de source ; l'obtention de données sismiques à partir de la réponse de formations souterraines à l'énergie sismique la formation de fonctions temps de parcours/déport pour les stations sismiques à partir des pointés d'arrivées premières obtenus à partir des données ; la transformation des fonctions temps de parcours/déport en fonctions vitesse/profondeur pour dériver un modèle d'une zone proche de la surface ; la détermination de statiques de période longue dans les données à partir du modèle dérivé de la zone proche de la surface ; et l'élimination des statiques déterminées de période longue des données sismiques. 17. Procédé selon la revendication 16, caractérisé en ce que l'introduction d'énergie sismique comprend l'introduction d'énergie sismique dans le sol à partir d'au moins une position sonique afin que cette énergie se déplace au moins en partie sous la forme d'ondes sismiques de réfraction. 18. Procédé de traitement de données sismiques obtenues pendant une prospection sismique à partir de plusieurs stations sismiques pour éliminer l'effet de statiques des données, caractérisé en ce qu'il comprend l'obtention de pointés d'arrivées premières à partir des données sismiques ; la formation de fonctions temps de parcours/déport pour les stations sismiques à partir des pointés d'arrivées premières ; la transformation de fonctions temps de parcours/déport en fonctions vitesse/profondeur pour dériver un modèle d'une zone proche de la surface ; la détermination de statiques de période longue dans les données à partir du modèle dérivé de la zone proche de la surface ; et l'élimination des statiques déterminées de période longue des données sismiques. 19. Procédé selon la revendication 18, caractérisé en ce que l'introduction d'énergie sismique comprend l'introduction d'énergie sismique dans le sol à partir d'au moins une position sonique afin que cette énergie se déplace au moins en partie sous la forme d'ondes sismiques de réfraction. 20. Procédé selon la revendication 18, caractérisé en ce qu'il comprend en outre la formation d'un affichage des données sismiques après l'élimination des statiques déterminées de période longue. 21. Procédé de traitement de données sismiques obtenues au cours d'une prospection sismique à partir de plusieurs stations sismiques pour éliminer des données l'effet des statiques de période longue et de période courte, caractérisé en ce qu'il comprend : la formation de fonctions temps de parcours/déport pour les stations sismiques à partir des pointés d'arrivées premières obtenus à partir des données ; la transformation des fonctions temps de parcours/déport en fonctions vitesse/profondeur afin de dériver un modèle d'une zone proche de la surface ; et la détermination de statiques de période longue dans les données à partir du modèle dérivé de la zone proche de la surface - l'exécution d'une décomposition compatible en surface des temps de parcours résiduels restant après la détermination des statiques de période longue ; la formation d'une estimation de statiques de période courte à partir de la décomposition compatible en surface ; l'élimination des statiques déterminées de période longue des données sismiques ; et l'élimination des estimations formées des statiques de période courte des données sismiques. 22. Procédé selon la revendication 21, caractérisé en ce qu'il comprend en outre la formation d'un affichage des données sismiques après l'élimination des statiques de période longue et de période courte.