FR 2729766 A1 19960726 FR 9500730 A 19950123 Méthode d'exploration d'un milieu à tectonique complexe En exploration géophysique, les spécialistes cherchent à définir la géométrie notamment des pièges stratigraphiques et/ou des pièges structuraux. Dans le premier cas, le problème consiste à accroître la précision des mesures de temps et de courbure des réflecteurs sismiques. Dans le deuxième cas, le problème consiste surtout à définir un champ de vitesse géologique de tranche et une sommation temps permettant l'obtention d'une image en profondeur du sous-sol à explorer. Pour résoudre ces problèmes de mesure de temps et de vitesse, plusieurs techniques d'acquisition peuvent être mises en oeuvre en sismique réflexion. La sismique réflexion est la méthode usuelle pour obtenir une image sismique du sous-sol. Dans cette méthode, à partir de sources d'énergie appropriées, appelées émetteurs, on émet des ondes acoustiques qui se propagent dans le sous-sol à explorer et se réfléchissent sur les différents réflecteurs qu'il renferme. Les ondes réfléchies sont enregistrées en fonction du temps sur des récepteurs adaptés disposés à la surface du sol ou dans l'eau. Chaque enregistrement ou trace fourni par un récepteur est ensuite affecté à la position du point ou d'une zone (bin) qui est situé au milieu du segment reliant l'émetteur au récepteur considérés. Cette opération constitue le tri en point milieu commun (PMC). Une technique de prospection sismique devenue incontournable est la couverture multiple. Dans une telle technique, les émetteurs et les récepteurs sont disposés à la surface du milieu de telle sorte qu'un même point milieu regroupe plusieurs enregistrements. La série d'enregistrements associés au même point milieu forme ce qu'il est convenu d'appeler une collection d'enregistrements ou de traces en point milieu commun. L'ensemble des collections d'enregistrements est associé à une suite de points milieux différents situés de préférence, le long d'une même ligne de traitement. Pour l'obtention de telles collections, il convient de répartir à la surface du milieu les émetteurs et les récepteurs suivant une organisation prédéterminée. Actuellement deux types de répartition prédominent. Dans le premier type dit "acquisition 2D", les émetteurs et récepteurs sont alignés sur une seule et même ligne. Ainsi les enregistrements sont tous associés à des trajectoires d'onde situées dans un même plan, au moins dans le cas théorique d'un sous-sol structuré en couches planes et parallèles. Ce type d'acquisition est utilisé couramment, en particulier, pour les acquisitions à grande échelle et pour étudier des sous-sols présentant une tectonique calme. Dans le second type dit "acquisition 3D", les émetteurs et récepteurs sont répartis plus ou moins régulièrement sur une surface de sorte que les enregistrements associés à un même PMC, en l'occurrence appelé "bin" car on attribue en fait un même PMC à des positions milieux appartenant à des surfaces élémentaires ou bin, correspondent à des trajectoires n'étant pas toutes situées dans le même plan vertical. Ce dernier type d'acquisition est plus particulièrement utilisé à ce jour dans le cas où l'on recherche une grande résolution (étude réservoir par exemple) ou dans le cas de tectonique complexe. Ces deux types d'acquisition et les traitements sismiques classiques associés permettent à partir des collections en PMC d'obtenir une image sismique dans le plan vertical passant par l'ensemble de ces points milieux si l'on fait l'hypothèse classique d'un sous-sol homogène et isotrope, en couches planes ou pentées et parallèles ; dans ce cas, les réflexions associées à chacun des réflecteurs du sous-sol s'alignent le long d'indicatrices qui sont théoriquement des hyperboles centrées à la verticale du point-milieu ; dans le cas de couches planes parallèles pentées, une correction d'obliquité est applicable, soit après sommation par une migration , une mesure du pendage, soit avant sommation par D.M.O. rendant la formule de DIX applicable; - si les enregistrements sont parallèles aux directions des génératrices d'une série cylindrique monoaxiale. On rappelle que la formule de DIX, qui s'applique dans le cas de réflexions dans un sous-sol formé de couches planes et parallèles et pour de faibles déports entre émetteurs et récepteurs associés, exprime la vitesse de tranche Vn dans une couche n du sous-sol limitée par des réflecteurs n et n-l, par la relation: dans laquelle VRMS(n) et VRMS(n-l) sont les vitesses moyennes quadratiques pour les trajets à déport nul depuis la surface respectivement aux réflecteurs n et n-i, tn et ton l étant les temps desdits trajets. Dès que l'on s'éloigne trop de ces cas, les traitements classiques, tous basés sur l'hypothèse de couches planes et parallèles (conditions d'application de la formule de DIX), ne permettent plus l'obtention d'une image correcte du sous-sol. Quand la formule de DIX n'est pas applicable le passage de l'information temps à l'information profondeur par les traitements classiques conduit à autant d'images profondeur que de traitements appliqués. De plus, les images obtenues contiennent des défauts irréversibles tels que superpositions locales ou absences locales d'images. La présence de trous ou d'informations redondantes ne peut être traitée par aucune des techniques de traitement conventionnel permettant de prendre en compte les effets de pendage ou de passer des temps aux profondeurs, par exemple, inversion, migration, tomographie... La présente invention a pour but une méthode d'acquisition et de traitement de données sismiques particulièrement adaptée à l'étude de tectoniques complexes permettant de remédier aux inconvénients précités, quand on n'est plus dans le champ d'application de la formule de DIX et que les effets de pendage ne peuvent plus être compensés par des traitements correctifs tels que DMO. A titre d'exemple, des déformations du sous-sol se décrivant soit par une structure continue dont les axes tournent le long de son contour, soit par un empilement de structures continues dont les axes sont tournés les uns par rapport aux autres ou encore des déformations se décrivant par des structures discontinues du fait de la présence de failles perpendiculaires à l'axe de la ligne de traitement ou parallèles à lui, seront à cataloguer en tant que tectoniques complexes. De même, toute tectonique contenant des coins de vitesse devra être considérée comme une tectonique complexe, les indicatrices présentant dans ce cas des variations de courbures (effet de coin). De même, une zone calme située sous une zone de tectonique complexe sera à considérer comme complexe. Un autre but de la présente invention est de permettre l'obtention des paramètres caractérisant au moins le champ de vitesse de sommation associé à des éléments réflecteurs correspondant à un point donné de la surface dudit milieu. La présente invention a pour objet une méthode dans laquelle on enregistre quatre sismiques 2D en particulier du type sérié, (style 3D), ces quatre sismiques 2D étant orientées dans quatre directions tournées les unes par rapport aux autres d'un angle qui est de préférence égal à 45" Plus précisément, la méthode consiste à réaliser une collection de traces sismiques en point milieu commun dont le point milieu commun est ledit point donné, à soumettre les traces de cette collection à un traitement de corrections dynamiques et/ou statiques, et elle se caractérise par le fait qu'on constitue ladite collection de traces en point milieu commun associées audit point donné en enregistrant quatre sous-collections de traces en point milieu commun associées audit point donné suivant quatre directions angulaires passant par ledit point donné et qu'on applique à chacune des quatre souscollections de traces en point milieu commun des corrections dynamiques et/ou statiques de façon à déterminer pour chaque sous-collection de traces une valeur de la courbure et/ou de la vitesse optimisant I'énergie de la somme des traces de la sous-collection corrigée considérée, les quatre valeurs ainsi obtenues représentant les composantes du champ de vitesse recherché pour chacun des éléments réflecteurs correspondant au point donné, chaque élément réflecteur étant caractérisé par un temps vertical to de propagation avec réflexion pour un déport nul. Un avantage de la présente invention réside dans le fait qu'on peut réaliser une sismique 3D avec de la sismique 2D pour des milieux qui ne peuvent être explorés qu'en trois dimensions. Un autre avantage réside dans le fait qu'on peut obtenir des informations relatives à des zones ou trous situées entre des zones pour lesquelles on avait des informations ; il en résulte qu'on a en tectonique complexe la même continuité d'informations qu'on obtenait par la méthode 3D classique en tectonique calme ou non complexe. Selon une autre caractéristique de l'invention, les directions angulaires passant par le point donné, sont telles que deux directions consécutives forment un angle aigu compris entre 40 et 500 et de préférence d'environ 45". La méthode peut être appliquée en une pluralité de points donnés constituant des points milieux communs, situés aux noeuds d'une grille à mailles régulières, de préférence en forme de parallélogramme. L'acquisition des sous-collections de traces est réalisée avantageusement selon des lignes coïncidant avec les côtés ou les diagonales des mailles en forme de parallélogramme formant la grille. Pour la constitution d'un bloc 3D de traces sommes, en chaque point milieu formant un noeud de la grille on constitue une trace somme en juxtaposant des portions de traces centrées sur les valeurs des to associées aux réflecteurs à la verticale dudit point milieu, chaque portion de trace centrée sur une valeur to associée à un réflecteur étant obtenue en sommant séparément chacune des quatre sous-collections de traces 2D associées au point milieu et en sélectionnant la portion centrée sur to de celle des traces sommes correspondant au maximum d'énergie. D'autres avantages et caractéristiques ressortiront mieux à la lecture de la description d'une méthode selon la présente invention, ainsi que des dessins annexés sur lesquels - la figure 1 est une représentation schématique d'un plan de position des lignes d'acquisition des données sismiques selon l'invention, le cercle représentant les déports maximum du dispositif en point milieu commun qui est constitué par le centre de l'étoile. - la figure 2 est une représentation schématique de trois types de collections de traces en point milieu commun, - la figure 3 est une représentation synoptique d'un milieu à tectonique complexe, - la figure 4 est une représentation en vue de dessus des lignes sur lesquelles sont définies les composantes du champ de vitesses et les paramètres définissant la forme d'onde en surface, - la figure 5 est une représentation schématique d'une quadri-BAP représentative de trois événements synthétiques, - la figure 6 est une représentation schématique d'une quadri-BAP où les événements de la quadri-BAP de la figure 5 sont pointés, - la figure 7 est une représentation de la trace somme des maxima d'énergie des événements pointés de la figure 6, - la figure 8 est une représentation schématique d'un milieu à tectonique complexe monoaxial. Sur la surface d'un milieu à explorer, on dispose un dispositif d'émission, de réception et d'enregistrement de données sismiques issues du sous-sol du milieu, ledit dispositif étant du type utilisé en sismique réflexion. Les lignes d'acquisition (enregistrement) c'est-à-dire les lignes de séries d'émetteurs et de récepteurs-enregistreurs, sont représentées sur la figure 1 et elles sont dirigées dans au moins quatre directions 1 à 4 qui sont telles que l'angle entre deux directions consécutives (1,2 ; 2,3 ;; 3,4) est de préférence égal à 45". Avec la disposition de la figure 1, on réalise des enregistrements quatre fois 2D et en particulier des enregistrements 4 x 2D sérié, c'est-à-dire un ensemble de lignes d'acquisition parallèles, chaque série de traits parallèles représentant des enregistrements en 2D de type marine ou de type linéaire avec tir au centre, mais avec un pas spatial qui est, par exemple, de 400 m pour les lignes dirigées suivant les directions 1 et 3 et d'environ 283 m pour les deux autres séries de lignes à 45" et dirigées suivant les directions 2 et 4. Bien évidemment, ont peut diminuer le pas spatial jusqu'à l'inter-trace qui est généralement de 25 m.Chaque ligne est enregistrée avec des paramètres classiques de 2D, par exemple 96 traces espacées de 25 m et une distance inter-tir de 12,5 m (stack array). Les enregistrements des quatre séries de lignes sont groupés en point milieu commun (PMC), les groupements en PMC pouvant être classés en trois catégories (figure 2). La première catégorie classique, représentée par une ligne unique 20, correspond aux déports entre émetteur-récepteur le point central 201 représentant le point milieu commun. On utilisera ce type de PMC pour l'imagerie. La deuxième catégorie est représentée par une croix 21. La troisième catégorie est représentée par une étoile 22 regroupant quatre lignes 221 à 224 à 45" les unes des autres. Le point milieu commun est centré sur un point donné A de la surface. Cette troisième catégorie permet de réaliser quatre sous-collections de traces PMC centrées sur le point donné A. Chaque sous-collection de traces PMC est corrigée dynamiquement et/ou statiquement de manière connue en soi.Sur les PMC de catégorie trois, on effectue une analyse de courbure par BAP comme indiqué pour la méthode POLYSTACK de E. de BAZELAIRE (GEOPHYSICS, vol. 53, n 2, Fev. 88, p. 143-157), les analyses de courbure ou BAP des quatre lignes d'un même PMC sont rassemblées en un seul tableau à quatre plans ou quadri-BAP, une BAP étant une collection de traces dont chacune est issue d'une même collection en point milieu commun après application d'une correction hyperbolique de type statique, les n traces d'une BAP correspondant alors à une investigation selon n courbures différentes. Après quoi, on détecte les maxima d'énergie comme pour le POLYSTACK en considérant le tableau dans son ensemble. Une fois les to pointés, c'est-à-dire une fois que les temps de trajet to pour un déport nul. correspondant à une distance nulle entre émetteur et récepteur, ont été déterminés, on pointe les maxima dans chaque tableau individuel pour les quatre valeurs tpl(to) tp2(to), t 2 tp3(to) et tp4(to) et ceci pour chaque to déterminé. Ces divers tp correspondent à la somme optimale de chacune des quatre sous-collections de traces au voisinage de la valeur to considérée dans une fenêtre élémentaire de largeur égale à l/B, B étant la bande passante des signaux enregistrés. La figure 3 représente schématiquement un sous-sol d'un milieu dont la surface est matérialisée par le plan S et qui comprend deux couches C1 et C2, de vitesses différentes V1 et V2 et séparées par un réflecteur présentant la forme d'une surface torique ST, les épaisseurs à l'apex des couches C1 et C2 étant respectivement égales à el et e2. Sur la surface S, on a représenté le point donné A, point milieu commun de chaque sous-collection de traces, une direction de référence X dans le plan de symétrie et une direction 4 d'enregistrement orientée d'un angle a par rapport à l'axe X de référence. La surface torique ST a pour rayons principaux orthogonaux R1 et R2, R1 étant mesuré dans le plan de symétrie. L'équation du rayon de courbure R(a) de la surface torique ST autour du point O est donnée par la formule I cos2a sin2a + R(a) R1 R2 Dans le cas de l'approximation paraxiale, c'est-à-dire le cas pour lequel tous les rayons utilisés forment avec l'axe de symétrie de la famille de rayons des angles petits pour lesquels l'angle en radians est égal à son sinus ou à sa tangente (angles inférieurs à environ 15o), on peut admettre que tous les rayons appartenant au point milieu commun (PMC) sont contenus dans un même plan, celui défini par la ligne d'enregistrement 4 et la droite AOM où A est le point donné ou PMC, O l'apex de l'interface ou surface torique ST et M le centre de courbure du front d'onde incident. Le centre de courbure M'(a) du front d'onde de sortie dans la première couche C1 est l'image optique du point source M à travers l'interface ST d'apex O. Le point M' se déduit donc de la position du point M et des paramètres du sous-sol par la formule de conjugaison optique (SHAH et LEVIN, 1973 et E. de BAZELAIRE).Cette formule, donnée ciaprès sous une forme simplifiée, exprime la stationnarité au deuxième ordre des temps de trajet entre deux points quelconques de l'espace situés de part et d'autre de l'interface et elle est conforme au principe de Fermat. 1 1 1 1 1 1 # # # # x - = x - (1) V2 OM(&alpha;) R(&alpha;) V1 OM'(&alpha;) R(&alpha;) Le temps to pour un déport nul. c'est-à-dire pour un couple émetteur-récepteur confondu au point donné A, est donné par la formule: 2e1 2e2 to = + (2) V1 V2 Dans ce cas on peut utiliser l'approximation plane pour les nappes de rayons tournées de &alpha; La théorie PSCAN développée par E. de BAZELAIRE, défiait la profondeur de mise au point P(a) telle que : 2P(a)=Vltp(a) où tp est le temps de profondeur de mise au point. Dans le cas de la figure 3, P(a) représente la distance AM'(a). Le calcul de tp(a) donne 2@2V2 1 tp(&alpha;) = 2@1/V1+ V1 1 + (V@ - V@)@@[sin2&alpha;/R@ + cos2&alpha;/R1] (3) Si on considère les définitions suivantes et si on reporte ces termes dans l'équation (3) donnant tp(a), on obtient: #(&alpha;) = #mcos2&alpha;+#Msin2&alpha; (5) L'équation (5) est celle d'un tenseur elliptique qui est exprimé en courbures translatées. La nouvelle équation de tp est alors donnée par la formule: tP(0) tP(#/2)-t01(tP(0)cos2&alpha;+tP(#/2)sin2&alpha;) tp(&alpha;= tp(0)sin2&alpha;+tP(#/2)cos2&alpha;-t01 (6) Cette équation a une forme intrinsèque: elle ne dépend plus du modèle utilisé. Elle est représentative de l'ensemble des fronts d'onde que l'on peut rencontrer en sismique 3D quand les interfaces entre les milieux traversés sont des surfaces continues et à dérivées continues. On voit sur la formule (6) que, pour chaque valeur de temps pointée, il y a quatre paramètres indépendants à rechercher qui sont tp(0), tp(7t/2), tlO et l'angle a ou azimut. Pour les retrouver, on doit faire alors obligatoirement au moins quatre mesures indépendantes t1p, t2p, t3p, t4p pour une même position temps to. Les axes de référence étant défiais par le dispositif de mesure, I'axe de référence angulaire est donc OX. Par rapport à cet axe, on définit l'angle a avec l'axe du "Parapluie" de sommation de plus grand tp, c'est-àdire de plus faible courbure de la représentation temps du front d'onde 3D. Un "Parapluie" de sommation est la représentation temps en fonction des déports et des azimuts, du front d'onde tridimensionnel reçu en surface par les récepteurs. 11 s'agit d'une extension à trois dimensions de l'indicatrice hyperbolique représentative en temps d'un front d'onde circulaire qu'on recevrait dans un espace à deux dimensions. La représentation temps cidessus est constituée en projection sur le plan de la surface par une ellipse E. On appelle t#= 1/## (7) p cette valeur de référence. L'axe perpendiculaire à cette grandeur est celui sur lequel la valeur du tp est minimale, c'est-à-dire de courbure la plus forte. On appelle cette valeur I i - @/# (8) La valeur #3 de l'axe à 45 degrés se déduit des deux valeurs précédentes et du temps de transit T.Les paramètres qui définissent le front d'onde sont donc: fP# tP(0) #tP# correspondant aux # tP(#/2) # t01 &alpha; &alpha; (9) Les paramètres de droite ont été définis lors de l'analyse de l'équation donnant tp(a). Les paramètres de gauche sont liés aux mesures faites par le système d'équations suivant: 11 faut inverser ce système pour trouver les paramètres: a (azimut), @ (temps de transit), tp\\ (courbure d'indicatrice minimum), tpl (courbure d' indicatrice maximum) en fonction des mesures. tP1 tP2 tP3 tP4 (11) Ce calcul est tout-à-fait classique. On commence par rechercher la valeur de @ par une suite de combinaisons linéaires entre les quatre éqttons du système ci-dessus. On trouve une équation du second degré: il y a deux valeurs de T qui sont solutions. Cela veut dire qu'il existe deux fronts d'ondes seulement le long des prolongations des rayons dans le milieu de sortie pour lesquels le parapluie de sommation est représentable par un ellipsoïde de courbures. C'est sur ces deux parapluies de sommation et seulement sur eux que l'on a le droit de faire l'interpolation spatiale bidimensionnelle. Les deux valeurs de T sont: Pour les deux valeurs de @ la valeur de l'angle a est formellement la même, ce qui signifie que les rayons calculés sont restés coplanaires.On obtient les deux valeurs de a modulo # par la formule suivante 2 4 i 3 (@P-@P)(@P-@)(@P-@) tan2&alpha; = (tP3-tP1)(tP4-#)(tP2-#) (13) Les valeurs de tp(parallèle) et tp(perpendiculaire) se déduisent des relations suivantes 1 1 tP# = #+ tP# = #+ ## ## (14) avec Muni de ces quatre valeurs, il est alors possible de reconstituer l'ensemble de toutes les intersections du parapluie de sommation (stack) par un plan vertical quelconque passant par son axe central. Dans ce qui précède et notamment en se référant à la figure 1, on voit que l'invention permet de déterminer les quatre valeurs nécessaires pour la reconstitution de l'ensemble de toutes les intersections du parapluie de sommation par un plan vertical quelconque passant par son axe central. De plus et selon l'invention, il est préférable de sommer les quatre lignes séparément et choisir après la meilleure sommation parmi ces quatre sommes. En conséquence, l'imagerie doit donc être dissociée entièrement de l'estimation du champ de vitesse, lesdites vitesses étant obtenues uniquement en les points de croisement du système d'enregistrement (4 x 2D) pour lesquels il existe les quatre lignes en même temps. On peut ainsi estimer de façon complète les quatre valeurs nécessaires et en déduire les vitesses de tranche et les courbures principales du sous-sol. Sur la figure 5, on a représenté une quadri-BAP relative à un exemple synthétisé de trois événements complets déterminés suivant les quatre lignes à 45" les unes des autres. Chaque BAP est relative à une souscollection de traces suivant l'une des directions d'acquisition 0=0, 0=n/4, 8=il/2, et 8=3x/4 ; cette BAP optimise l'énergie de la somme des traces et elle est représentative des événements apparaissant aux temps tl à t3 en fonction de l'indice tp (profondeur de mise au point). Au temps tl, c'est l'événement de la ligne 0=0 et d'indice il des tp qui présente le maximum d'énergie. Au temps t2, c'est l'événement de la ligne 0=X/4 et d'indice i2 qui présente le maximum d'énergie, alors qu'au temps t3, c'est l'événement de la ligne située à 37r/4 et d'indice i3 qui présente le maximum d'énergie. C'est ce qui est représenté sur la figure 6 qui est un exemple de pointé des trois événements. La figure 7 est un exemple de construction de la trace somme des maxima d'énergie où on retrouve les trois événements aux divers temps tl, t2, t3. Dans le cas où la structure complexe est composée d'une série cylindrique (plis uni-axiaux) c'est-à-dire que tous les cylindres ont des génératrices parallèles entre elles et ne tournent pas les uns par rapport aux autres et lorsqu'on connaît la direction des axes, on peut s'affranchir des enregistrements suivant les axes à 45" pour se limiter aux directions DIP et STRIKE c'est-à-dire respectivement de plus fort et de plus faible pendage. Dans ces conditions, on combine les informations DIP et les informations STRIKE alors que jusqu'à maintenant on utilisait les informations DIP ou les informations STRIKE. En se reportant à la figure 8, on voit que nous sommes en présence d'un sous-sol de structure cylindrique mono-axiale c'est-à-dire que le sous-sol est assimilé à un cylindre d'axe AX dont tous les rayons de courbure parallèles à l'axe AX sont infinis. Lorsqu'on connaît ladite direction AX, L'équation (5) du tenseur elliptique dégénère et il suffit de mesurer les tp parallèles et perpendiculaires à l'axe AX pour connaître la totalité du tenseur qui devient symétrique. En fait, tous les termes en a disparaissent et on est ramené à une acquisition de données sismiques de type 2 x 2D sérié dans les directions AX (STRIKE) et AY (DIP). Il est évident qu'on peut choisir d'autres valeurs d'angle entre les lignes d'acquisition pour autant qu'on puisse se ramener à une base orthogonale dans un espace vectoriel à quatre dimensions. Une base est orthogonale si lorsqu'on fait varier un seul des quatre paramètres à mesurer, les trois autres ne varient pas. Lorsque la base n est pas orthogonale c'est-à-dire lorsque les directions des lignes d'acquisition ne sont pas exactement à 45" les unes des autres, on définit une autre grille ou étoile dont les axes sont à 45" et on effectue une projection des mesures des paramètres sur ces directions de manière à se ramener à une base orthogonale. Une telle projection n'est possible que si les directions des lignes d'acquisition de données sismiques sont suffisamment éloignées les unes des autres. Dans certains cas, en sismique marine par exemple, il peut être difficile de faire passer exactement les quatre lignes d'acquisition par un même point. On peut utiliser alors le centre de gravité de la zone formée par les intersections desdites quatre lignes d'acquisition comme centre de l'étoile et on projette les mesures des paramètres acquises avec lesdites lignes d'acquisition sur la base orthogonale correspondant à ce centre de gravité pour déterminer ensuite les paramètres du front d'onde. REVENDICATIONS 1. Méthode d'exploration d'un milieu en sismique réflexion, particulièrement adaptée à un milieu à tectonique complexe, pour l'obtention des paramètres caractérisant au moins le champ de vitesse de sommation associé à des éléments réflecteurs situés à la verticale d'un point donné (A) de la surface du milieu (S), ladite méthode consistant à réaliser une collection de traces sismiques en point milieu commun dont le point PMC est ledit point donné, à soumettre les traces de cette collection à un traitement de corrections dynamiques et/ou statiques, et se caractérisant en ce que l'on constitue la collection de traces associée au point donné (A) en enregistrant quatre sous-collections de traces en PMC associées audit point donné suivant quatre directions angulaires (1 à 4) passant par ledit point donné et l'on applique à chacune des quatre sous-collections de traces en PMC des corrections dynamiques et/ou statiques de façon à déterminer pour chaque sous-collection une valeur de la courbure et/ou de la vitesse optimisant l'énergie de la somme des traces de la sous-collection corrigée considérée, les quatre valeurs ainsi obtenues représentant les composantes du champ de vitesse recherché pour chacun des éléments réflecteurs correspondant au point donné et caractérisé par un temps vertical to de propagation avec réflexion pour un déport nul. 2. Méthode selon la revendication 1, caractérisée en ce qu'à partir des quatre composantes du champ de vitesse associé à un desdits éléments réflecteurs, on détermine des paramètres tels que a, T, tp\\ et tpl caractérisant la géométrie dudit élément réflecteur, a représentant l'azimut de l'élément réflecteur par rapport à une direction de référence prédéterminée, tp\\ et tpl désignant respectivement les valeurs de plus faible pendage et de plus fort pendage dudit élément réflecteur et t désignant le temps de transit. 3. Méthode selon la revendication I ou 2, caractérisée en ce que les directions angulaires passant par le point donné, sont telles que deux directions consécutives forment un angle aigu compris entre 40 et 50 et de préférence d'environ 45". 4. Méthode selon l'une des revendications I à 3, caractérisée en ce qu'elle est appliquée en une pluralité de points donnés constituant des points milieux communs, situés aux noeuds d'une grille à mailles régulières, de préférence en forme de parallélogramme. 5. Méthode selon la revendication 4, caractérisée en ce que l'acquisition des sous-collections de traces est réalisée selon des lignes coïncidant avec les côtés ou les diagonales des mailles en forme de parallélogramme formant la grille. 6. Méthode selon la revendication 5, pour la constitution d'un bloc 3D de traces sommes, caractérisée en ce qu'en chaque point milieu formant un noeud de la grille on constitue une trace somme en juxtaposant des portions de traces centrées sur les valeurs des t0 associées aux réflecteurs à la verticale dudit point milieu, chaque portion de trace centrée sur une valeur to associée à un réflecteur étant obtenue en sommant séparément chacune des quatre sous-collections de traces 2D associées au point milieu et en sélectionnant la portion centrée sur to de celle des traces sommes correspondant au maximum d'énergie. 7. Méthode selon la revendication 1, caractérisée en ce que lorsque la structure du sous-sol à explorer est une structure cylindrique monbaxiale dont la direction de l'axe longitudinal (AX) est déterminée, on réalise une acquisition de données sismiques du type 2 x 2D sérié uniquement dans les directions parallèles (STRIKE) et perpendiculaires (DIP) audit axe (AX).