Brevet FR2544869A1 — publié le 19841026

L'invention concerne l'exploration géophysique et plus particulièrement la collecte simultanée de données d'ondes primaires ou secondaires classiques, converties, suivie d'un traitement des traces résultantes dans lequel les longueurs des trajets non symétriques des ondes primaires incidentes et des ondes converties réfléchies sont prises en compte avant la sommation des traces pour transformer une section de déport fini résultante en une section de déport zéro réel quelle que soit la profondeur ou l'inclinaison du réflecteur L'utilisation de récepteurs multicomposantes à des fins de collecte est prévue. La collecte et le traitement classiques de données de sismique réflexion commence avec la génération séparée d'ondes classiques de pression (ondes primaires ou ondes P) ou d'ondes transversales classiques (ondes secondaires ou ondes S), suivie de leur enregistrement séparé sur des récepteurs de composante unique, c'est-à- dire des récepteurs qui comportent des éléments actifs réagissant à des mouvements des ondes réfléchies dans une seule direction Si l'on suppose une source sismique orientée verticalement, des ondes P classiques descendent dans la terre et sont réfléchies par une de récepteurs, dont les éléments actifs réagissent unique- ment à un mouvement d'ondes élastiques orienté verticale- ment, enregistre les ondes P De façon similaire, dans le cas d'une exploration par ondes transversales ou de cisaillement, les ondes S produites par une source sis- mique orientée horizontalement sont réfléchies par des réflecteurs similaires sous forme d'ondes S et sont en- registrées par le dispositif de récepteurs, d'une manière similaire, sauf que les éléments actifs des récepteurs réagissent exclusivement au mouvement d'ondes orientées horizontalement. Le traitement des données des ondes P et/ou des ondes S est en outre compliqué par le fait que la collecte est généralement réalisée à l'aide de procédés de "progression" à point milieu commun (procédés CMP). De tels procédés utilisent des dispositifs superposés de récepteurs en combinaison avec des sources "progressant vers l'avant" suivant une ligne de relevé pour générer des nombres substantiels de traces smiques "redondantes". Autrement dit, ces traces sont redondantes par le fait qu'un certain nombre de traces peuvent être associées au même point central commun situé à mi-distance entre plu- sieurs paires sources-récepteurs respectives qui produisent les traces dans le premier emplacement Après l'applica- tion de décalages dans le temps à ces traces (appelés corrections statiques et dynamiques), un groupement à point milieu commun (CMP) est créé Puis les traces associées de ce groupement font l'objet d'une sommation pour donner une caractéristique signal/bruit améliorée, (compte tem de l'importance de la compréhension de la relation entre les coordonnées de collecte, dans lesquelles des traces sont identifiées soit par des coordonnées de positions-sources (s) et d'emplacements-récepteurs (g) suivant la ligne de relevé, soit par des coordonnées associées à la distance de déport (f) des stations source- à-récepteur, et la position du point milieu (y) entre la source et le récepteur des paires respectives, on peut se référer, en détail, à l'ouvrage de John F Claerbout "FUNDAMENTALS OF GEOPHYSICAL DATA PROCESSING-, McGraw- Hill, 1976, pages 228 et suivantes). Quand bien même le groupement des traces, obtenu par sommation, est amélioré (en raison de la sommation), des interprétations peuvent encore être rendues difficiles par le fait qu'aux limites entre différents types de roches, une conversion partielle se produit entre un type d'ondes et un autre, si l'on suppose que l'angle d'incidence-de l'onde est supérieur à zéro Par exemple, une onde incidente P peut être partiellement convertie en une onde réfléchie Sv En variante, une onde incidente Sv peut être partiellement convertie en un signal réfléchi d'une onde P. Bien que les équations de Zoeppritz déterminent les amplitudes des ondes réfléchies et converties, elles ont été rarement utilisées par les déchiffreurs de données géophysiques malgré le fait que les procédés modernes de collecte en sismique réflexion, tels lue les procédés CMP, utilisent de longs déports et font appel à un angle d'incidence important La raison en est que dans le cas de réflecteurs à plus grande profondeur, les angles d'incidence sont relativement faibles et les contrastes de vitesse et de densité entre les couches sont supposés être faibles Voir, par exemple, pages 21 et suivantes, de l'ouvrage de Kenneth H Waters, intitulé "A TOOL FOR ENERGY RESOURCE EXPLORATION", John Wiley and Sons, 1978. De plus, la complexité de l'application de ces équations aux nombreux niveaux d'énergie différents associés aux diverses ondes réfléchies à tous les angles d'incidence et aux divers contrastes des matières présentes dans le champ peut être à l'origine d'un nombre suffisant de données pour submerger simplement l'interprète ou le déchiffreur Ce dernier peut trouver trop difficile d'appliquer les équations de Zoeppritz de façon systéma- tique, en particulier lorsque les données de champ sont collectées par des procédés modernes à point milieu commun (CMP) A cet égard, quand bien même les points centraux/points de réflexion peuvent ne pas être alignés verticalement, l'interprète ignore habituellement le fait, c'est-à-dire qu'il ignore les différences de longueur des trajets des ondes P converties en ondes Sv sur des projections verticales passant par des points centraux situés à mi-distance entre les paires respectives sources- récepteurs. Autrement dit, avec des ondes incidentes et réfléchies classiques, les points de réflexion de réflec- teurs plats et horizontaux sont situés directement au- dessous des projections verticales des points milieu de paires respectives soures-récepteursassociées aux traces intéressantes Par conséquent, les traces associées aux points communs de réflexion (ou points de profondeur) sur des réflecteurs plats, bien que provenant de diffé- rentes paires sources-récepteurs, peuvent être additionnées (traitées par sommation), après corrections appropriées pour aligner les traces Cependant, avec des ondes converties dans les mêmes conditions, les points de réflexion ne sont pas situés au-dessous des projections partant des points milieux des paires respectives sources- récepteurs, mais ils sont décalés d'une certaine distance par rapport à ces projections. L'art antérieur le plus proche décrit le problème des longueurs de trajets non symétriques et se trouve dans l'ouvrage "DIGITAL PROCESSING OF TRANSFORMED REF Ii ECTED XAVES", SOVIET GEOLOGY AND GOEPHYSICS, V 21, N O 4, pages 51-59. T.T Nefedkina et collaborateurs y décrivent l'utilisation d'ondes P converties en ondes Sv dans des régions de la Sibérie et autres régions analogues présentant un pergélisol Un processus de sommation pour ces ondes converties donne l'avantage de faire varier le point de sommation des groupements conformément-à une série de valeurs de normalisation associées à un code numérique de traitement soviétique spécial appelé "Kondakova alpha-langage" Mais, étant donné que le processus utilise une terminologie de traitement particulière, des groupes de données inférieures et des modèles simplistes, une utilisa- tion classique de ce travail dans le contexte de méthodes modernes d'exploration, en particulier lorsque l'on prévoit d'utiliser des réflecteurs à immersion, est impossible. Conformément à l'invention, les longueurs de trajets non symétriques d'ondes primaires incidentes et d'ondes réfléchies converties,que l'onde incidente soit ou non une onde P ou une onde Sv, sont prises en compte de manière que des traces de l'onde Sv ou P convertie puissent être correctement collectées en groupes dans lesquels les traces associées à chaque groupe échantil- lonnent essentiellement le même point de réflexion sur un réflecteur cible commun Autrement dit, les traces d'ondes Sv ou P converties peuvent être identifiées par des coordonnées de points de réflexion communs (CRP) qui ont été correctement transformées à partir de coordonnées d'un point source (SP) et d'un poste de détection (D) afin de tenir compte des longueurs de trajets non symétriques des ondes incidentes et converties. La présente invention est basée en partie sur le fait que les longueurs de trajets non symétriques sont une fonction du type de la source, du rapport de la vitesse des ondes P à la vitesse des ondes S, c'est-à- dire le rapport Vp/Vs des terrains de couverture au-dessus du réflecteur cible, ainsi que l'angle d'inclinaison et la profondeur du réflecteur Pour mettre en relation les variables précitées, la présente invention prévoit la collecte simultanée de données d'ondes P ou S classiques et converties, sur place, par la mise en oeuvre de méthodes classiques de collecte à points milieux communs (CMP). Ces techniques consistent à actionner séquentiellement une source sismique classique placée en une série de positions de sources-points (SP) suivant une ligne de levée-ou de prospection et à collecter de façon redondante les réflexions à la fois des ondes classiques et de leurs phases converties au moyen d'une série de détecteurs multicomposantes Ces détecteurs sont positionnés le long de la ligne de prospection dans plusieurs stations (D) de détecteurs En conséquence, chacune des traces classiques et converties peut être associée à une paire sourcepoint/station de détecteur dont les coordonnées de l'emplacement sont connues (SP par rapport à D) Ensuite, les traces classiques subissent une transformation de coordonnées et un accroissement ou une amplification afin de tenir compte (a) des différents trajets de parcours de l'onde classique dans le terrain de couverture au-dessus du réflecteur cible, et (b) de la profondeur et de l'in- clinaison de chaque réflecteur cible Après que la pro- fondeur et l'inclinaison du réflecteur cible ont été ainsi déterminées, l'invention règle la pente de lignes droites imaginaires de groupement sur un diagramme de sommation de coordonnées (SP en fonction de D) qui identifie des traces converties communes Cet ajustement tient compte des différences de longueurs de trajet et des angles d'incidence et de réflexion des ondes incidentes et réfléchies, ainsi que de l'inclinaison et de la profondeur des réflecteurs cibles La sommation des traces remises en ordre a lieu alors de façon à former une section de décalage zéro réel des traces converties. Dans le cas de réflecteurs plats, les coordon- nées du point commun de réflexion (CRP) pour un groupement de traces classiques projetées sur un plan horizontal de référence, sont en relation avec les coordonnées de la source-point (SP) et de la station du détecteur (D) conformément à: CRP = (D+SP)/2. La transformation ci-dessus peut être consi- dérée comme étant un processus pour déterminer les coordon- nées d'un point de réflexion de la cible, qui ont été projetées sur le plan de référence par multiplication d'une constante (k) qui tient également compte du rapport de vitesses (des ondes incidentes et des ondes réfléchies dans le terrain de couverture), multiplié par les coordon- nées (SP) et (D) de paires communes respectives sources- récepteurs constituant le groupement de traces, conformé- ment à une équation de transformation de la forme CRP = k D + ( 1-k>SP, o k est égal à 0,5. Conformément à la présente invention, les coordonnées du point commun de réflexion (CRP) pour un groupement de traces converties peuvent être de la même façon mises en relation avec les coordonnées de la source- point (SP) et de la station de détecteur (D), c'est-à- dire conformément à des équations similaires de transfor- mation, de la forme: Conversion P-Sv: CRP k D + ( 1-k)SP Conversion Sv-P: CRP ( 1-k)D + KSP o SP et D sont des coordonnées de source et de détecteur respectivement et k est une constante qui tient compte du rapport des vitesses du terrain de couverture (Vp/Vs), de l'angle d'inclinaison a et de la profondeur (h) du réflecteur cible, ainsi que de la distance de départ source-récepteur X Vp et Vs sont les vitesses des ondes P et Sv, respectivement, du terrain de couverture. Une formule approchée pour k, qui englobe tous ces effets, est: k (Vp/Vs) k=p l(Vp/Vs) + {l ±(X/h)sin x}(l+f)-1/2 l o f= s) l+X/h)2 cos 2 x lV/Vs) + 1 Vp, Vs ont trait au rapport des vitesses du terrain de couverture; l'angle dp (a) et la profon- deur (h), ainsi que le déport X source-récepteur sont tels que définis précédemment Cette formule s'est révélée précise pour des inclinaisons atteignant;a = 30 , et des déports aussi importants que le double de la profondeur du réflecteur, c'est-à-dire 2 h, et elle peut également être appliquée à des ondes classiques P-P (ondes longi- tudinales) si Vp/Vs est égal à un. Dans de nombreuses applications de déport (o le déport X est très inférieur à la profondeur h et l'inclinaison a est nulle) k est donné approximativement par l'expression plus simple: k = (V /Vs) (Vp/Vs) + 1 Dans ce cas, avec des réflecteurs cibles plats placés au-dessous d'une couverture pour laquelle Vp/vs ='2,4, l'évaluation de k à l'aide des équations ci-dessus se réduit à: P-Sv: CRP = 0,73 D + 0,27 SP - Sv-P: CRP = 0,27 D + 0,73 SP DEFINITIONS Dans la présente invention, certains termes fondamentaux ayant trait à la collecte et au traitement de données-sismiques à points multiples seront utilisés et sont définis ci-dessous. On suppose que chaque trace collectée CMP est décrite par la fonction W (SP, D) et que les coordon- nées de position de la source (SP) et les coordonnées de la position de récepteur (D)-sont des variables indé- pendantes. En réalité, la position de la source (SP) et l'emplacement du récepteur (D) n'ont pas une distribution continue suivant la ligne (ou l'axe) de prospection définie par des coordonnées x, mais ils sont habituellement suffi- samment rapprochés pour qu'il soit possible, par simple interpolation, de trouver W pour toutes coordonnées (SP) et (D) De plus, le long de l'axe x se trouvent les coordon- nées (f> de déport de distance de la source au récepteur et les coordonnées (CMP) des emplacements des points milieux communs, qui sont orthogonales entre elles, mais qui intersectent le plan (SP, D) sous un angle donné dépen- dant-des paramètres de collecte sur site Si la distance de "progression" incrémentielle est ASP=AD, l'angle d'intersection est alors de 450 et les coordonnées de déport et du point milieu sont en relation avec les coordon- nées de la source et du récepteur conformément à f = D SP CMP = (D+ SP)/2. 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 vue en plan d'un système de collecte à points milieux communs (CMP) illustrant comment des données CMP sont collectées de façon classique sur place, à l'aide, par exemple d'une ligne de détecteurs Dl Dm associés à des sources situées en des sources- points SP 1 SP, l'activation de la source en SP, per- mettant l'enregistrement de données par les détecteurs D 1 Dm et l'activation de la source en SP 2 permettant l'enregistrement de données dans les positions D 2 Dm+ 1 les figures 2 A à 2 D sont des coupes verti- cales d'une formation de terrain ayant fait l'objet d'un levé à l'aide du système de collecte CMP de la figure 1, ces coupes illustrant en détail comment des réflexions classiques sont enregistrées pour une source, et montrant également le fait qu'en l'absence de conversion de mode au réflecteur concerné, les longueurs des trajets de l'onde incidente et de l'onde réfléchie sont symétriques par rapport à des points de réflexion associés sur un réflecteur plat de sorte que les traces associées à des points centraux de points milieux communs entre des posi- tions de source et des emplacements de récepteur respectifs coïncident, que l'onde générée par la source soit ou non (i) une onde P (figure 2 A), (ii) une onde Sh (figure 2 B) ou (iii) une onde Sv (figure 2 C), pourvu que le récepteur associé possède une capacité de correspondance de réponse aux composantes (figure 2 D); les figures 3 et 4 sont des coupes verticales d'une formation de terrain soumise à un levé à l'aide du système de collecte CMP montré sur la figure 1, ces vues illustrant en détail le changement d'emplacement du point de réflexion en fonction du mode de propagation d'onde élastique initiale (c'est-à-dire suivant qu'il s'agit d'une onde P ou d'une onde Sv), o l'inclinaison de la cible est égale à O et la profondeur, les coordonnées Dl Dn de la station du détecteur et les coordonnées S Pl de l'emplacement de la source-point restent constantes; les figures 5 et 6 sont des coupes d'une formation de terrain montrant en détail que des traces associées à un groupement de points milieux communs coïncident (ou ne co Incident pas) avec un point de profon- deur commune associée sur un réflecteur plat, suivant qu'une conversion d'onde ou non s'est produite; la figure 7 est un diagramme de sommation dans lequel les coordonnées sources-récepteurs associées à des traces produites par des étapes de collecte CMP sont superposées sur un système de coordonnées de points milieux communs en fonction du déport, afin de mieux illustrer les processus associés au procédé de la présente invention; les figures 8 à 11 sont des coupes d'une formation de terrain illustrant en détail des techniques de modélisation de traces de rayon utilisées conformément à l'invention pour déterminer le degré de variation des points de réflexion le long de réflecteurs plats '(figure 8) et le long de réflecteurs inclinés (figures 9 à 11) en fonction du déport, lorsqu'une conversion d'onde élas- tique se produit; la figure 12 est un organigramme d'un processus mettant en oeuvre le procédé de l'invention en utilisant un calculateur numérique convenablement programmé, dans lequel des traces converties peuvent être groupées suivant des lignes de groupement appropriées pour tenir compte de trajets de parcours non svm triqes des rayons incidents et réfléchis à l'aide d'une équation choisie de transfor- mation de coordonnées; la figure 13 est une vue de détail à éèhelle agrandie du diagramme de sommation de la figure 7, sur lequel les lignes de points communs de réflexion (CRP) le long desquelles des traces converties peuvent être groupées, sont établies par une équation choisie de trans- formation de coordonnées; les figures 14 et 15 sont des organigrammes "'un traitement de traces classiques et converties, collectées simultanément, permettant la mise en oeuvre du procédé de l'invention; les figures 16 à 18 sont des coupes-films et des parties de coupes illustrant la capacité de diagnos- tic du procédé de l'invention à résoudre des structures complexes dans des exemples sur site réel; et les figures 19 et 20 sont des courbes de vitesses de sommation Vp, Vs, ainsi que de leurs moyennes géométriques, en fonction de l'emplacement du point de tir, d'événements séparés dans les coupes-films des figures 16 à 18. Avant de décrire en détail différentes formes de réalisation de l'invention, on procédera briève- ment à une description des concepts permettant de mieux comprendre l'invention. Dans le cas d'ondes sismiques classiques pro- duites et reçues, telles que des ondes P produites par une source d'ondes P et enregistrées sur des récepteurs orientés verticalement, les points dommuns de réflexion (CRP) de réflecteurs plats sont des projections verticales des points milieux communs compris entre les paires res- pectives sources-récepteurs associées aux traces enregis- trées Cette supposition est valable dans des zones d'exploration de structure ou d'inclinaison ou pendage limité et elle permet aux responsables de l'exploration de grouper des données sismiques commodément et efficace- ment Cependant, dans le cas d'ondes de conversion, les angles d'incidence et de réflexion sur des réflecteurs se trouvant au-dessous de la surface sont inégaux Le résultat est que les points de réflexion ne sont pas des projections verticales des points milieux entre les paires sourcesrécepteurs. La figure 1 est une vue en plan d'un système de collecte sismique illustrant comment sont établies les expressions intéressantes du présent mémoire. Par exemple, les expressions "point central" et "point milieu" correspondent à un emplacement géogra- phiqu se trouvant à mi-distance entre une coordonnée de source, par exemple la position de source S Pl d'une série de positions-de source SP 1 SP 4 d'un système de collecte sur site, et une série de positions de récepteurs, par exemple Dl d'une série de position sde récepteurs Dl Dm, à un horizon de référence proche de la surface de la terre, coïncidant avec la ligne de prospection ou de levé 11 Les points centraux sont désignés Cl C 4 et chaque point central apparaît comme pouvant être associé de façon adressable à une paire sources-détecteur choisie qui produit une trace donnée. Dans une collecte de points milieux communs (CMP), différents jeux de dispositifs de détecteurs "progressent" vers l'avant, dans le sens d'une flèche 12, en association avec le positionnement incrémentiel analogue, vers l'avant, d'une source dans de nouvelles positions suivant la ligne de levé 13 Des excitations dans les positions de la série peuvent alors s'effectuer suivant une séquence Autrement dit, la source est excitée de façon séquentielle dans les positions de source SP 2 * SSP 4. Le résultat est que des traces enregistrées aux emplacements des différents récepteurs peuvent être mises en relation avec des points milieux choisis, situés dans la zone souterraine prospectée Par exemple, si l'interface de réflexion est un horizon plat, le point de réflexion ou une réflexion se produit définit des lignes verticales qui passent pas les points centraux C 1 C 4 intéressants. L'application de corrections statiques et dynamiques à une trace obtenue sur site équivaut (dans les conditions indiquées ci-dessus) à placer la source au point central, à activer cette source, puis à la rem placer par un détecteur au mrêe emplacement et à enregis- trer la trace Si toutes les traces associées à un point milieu commun sont recomposées sur une base côte à côte, le jeu de traces résultant est appelé un groupement CMP. Si les traces sont additionnées, la trace résultante est un groupement CMP de sommation Par conséquent, l'inter- prétation sismique demande au traitement CMP classique au moins deux critères: (i) l'amélioration du rapport signal/bruit (SNR) des réflexions sur le groupement de sommation résultant, par sommation des traces associées à un point souterrain commun; (ii) les projections des points de réflexion souterrains passent par les points milieux des paires sources-détecteurs (SP,D) d'emplacements de coordonnées connues; par conséquent, les emplacements des réflexions de structure sont connus si les longueurs des trajets d'incidence et de réflexion sont sensiblement égales. Les figures 2 A à 2 C illustrent des phénomènes de réflexion d'un modèle à trois couches typique d'une section géologique peu profonde, jeune 20, constituée d'un grès 21 compris entre les couches 22 et 23, tel qu'on en trouve dans le golfe du Mexique, montrant que, même si des sources séparées d'énergie 28, 29 ou 30, produisant des ondes P, des ondes Sh ou des ondes Sv sont placées dans des positions SP 1, SP 2, SP 3, comme montré respectivement sur les figures 2 A, 2 B et 2 C, puis excitées, les longueurs des trajets d'incidence et de réflexion de l'énergie associée à chaque paire sources- récepteurs sont sensiblement égales Ceci suppose que les récepteurs 31, 32 et 33 dans les positions Dl, D 2 et D 3 soient préparés de façon à ne recevoir que l'énergie dominante de l'onde produite Des projections verticales 34, 35 et 36 des points de réflexion classiques communs CMP 1, CMP 2 ou CMP 3 coupent les points centraux (points médians)C 1, C 2, C 3 des coordonnées des paires respectives position de source-position de récepteur, comme montré. A cet égard, le terme "classique" est utilisé pour qualifier l'énergie reçue par les récepteurs 31, 32 et 33, énergie dans laquelle le mouvement dominant des particules correspond à celui de l'onde produite, que la source produise de l'énergie sous forme d'une onde P, d'une onde Sh ou d'une onde Sv dans son mode principal d'activation. Il convient de noter sur la figure 2 A que la source 28, située à l'emplacement ponctuel SP 1, produit une onde incidente qui se déplace vers l'extérieur de l'emplacement de cette source, sous la forme d'une série de fronts d'onde Chaque front d'onde définit une surface sphérique commune qui joint des points de la zone souter- raine o un mouvement est sur le point de commencer Si le milieu de propagation présente des propriétés indé- pendantes de la position et de la direction du parcours, les fronts d'onde forment un groupe de sphères concentri- ques centrées sur l'emplacement de la source Il est commode de "suivre" de tels fronts à l'aide de procédés de traçage de rayons classiques dans la technique, procédés dans lesquels on conçoit que l'énergie des fronts se déplace vers le bas dans la terre suivant un grand nombre de pyramides de section infinitésimale, et la ligne centrale de l'une quelconque de ces pyramides, dans une directicnchoisie, est considérée comme un rayon qui trace des trajets 37, 38, 39 Ces trajets passent à travers les strates 21 et 22 de la coupe géologique 20. Il convient de noter qu'au réflecteur 23, l'angle d'incidence de l'onde incidente est égal à l'angle de réflexion de l'onde réfléchie Par conséquent, la longueur incrémentielle du trajet de l'onde incidente et de l'onde réfléchie dans le terrain de couverture, c'est-à- dire dans les strates 21 et 22, est la même. La source 28 d'ondes P sur la figure 2 A est généralement une charge de dynamite enterrée ou un vibrateur monté sur des camions et vibrant verticalement sur la surface du sol La plupart des sources communes de cisaillement (onde S des figures 2 B et 2 C) comprennent un vibrateur qui transmet au sol des secousses orientées latéralement plutôt que de haut en bas Sur la figure 2 B, le vibrateur Sh 29 produit des secousses horizontales perpendiculairement à la direction de la flèche 12 de collecte CMP Si le vibrateur est tourné horizontalement de 90 de manière que le mouvement suive la direction de la flèche 12, comme indiqué sur la figure 2 C, la source est appelée source de mouvement de cisaillement ou mouve- ment transversal, de type Sv Il existe une autre diffé- rence entre l'énergie des ondes Sh et celle des ondes Sv. Sur les figures 2 A, 2 B et 2 C, on note également dans la strate 21, que la direction dumouvement des particules de l'énergie incidente et de l'énergie réfléchie, telle qu'indilquée par des paires de flèches 40 a, 40 b; 41 a, 41 b; et 42 a et 42 b, peut ou noni changer lorsqu'une réflexion se produit sur la strate réfléchissante 23. Sur la figure 2 A, par exemple, la flèche 40 a associée à l'onde incidente est dirigée vers le bas suivant le trajet 37, tandis que l'onde réfléchie associéeà la flèche 40 b est dirigée vers le haut De façon similaire, sur la figure 2 C, la flèche 42 a associée à l'onde incidente est dirigée vers le haut et s'éloigne du trajet 39 {per- pendiculairement à lui), tandis que l'onde réfléchie est dirigée vers le bas par rapport au trajet du rayon. Par ailleurs, sur la figure 2 B, la flèche 41 a de l'onde incidente définit les mouvements de particules et est perpendiculaire au plan du dessin Etant donné que le mouvement des particules est parallèle à la surface réflé- chissante, les ondes Sh ne sont pas affectées par une conversion de mode par réflexion ou réfraction sur le réflecteur cible Autrement dit, la source 29 d'ondes Sh produit des rayons d'énergie qui, à la suite d'une réflexion sur des lits plats, ne donnent que des ondes Sh qui, lorsqu'elles sont enregistrées par le récepteur 32, exigent seulement que ce récepteur 32 soit orienté per- pendiculairement à la flèche 12 indiquant la direction de la prospection selon laquelle s'effectue la collecte. Par ailleurs, des ondes P et des ondes Sv arrivant sur un réflecteur ne donnent pas seulement des ondes classiques, de type analogue, mais produisent égale- ment des ondes de conversion Lorsque les deux types d'ondes arrivent à la surface, des réflexions Sv sont enregistrées sur le segment de mouvement horizontal radial des récepteurs 31 et 33, alors que des ondes P sont enregistrées principalement sur le segment de réponse verticale des détecteurs. La figure 2 D illustre comment les directions de réponse dominantes des récepteurs 31, 32 et 33 peuvent être accordées pour réagir à un mouvement de particules de l'énergie d'ondes incidentes P ou S. Comme représenté, si la direction de la réponse des récepteurs 31, 32 et 33 est verticale par rapport au champ de gravité terrestre, à savoir suivant la flèche 43, l'énergie de toute onde P arrivant est détectée si la réponse du récepteur est horizontale, dans la direction de la flèche 44, l'énergie de toute onde Sh arrivant est détectée; de façon similaire, si la direction de la réponse est horizontale, dans la direction de la flèche 45, l'énergie de toute onde Sv arrivant est de la même manière détectée. Les figures 3 et 4 montrent que le point de réflexion sur un réflecteur plat (pour des ondes classi- ques) n'est pas une projection verticale du point milieu entre des paires sources-récepteurs respectives. Sur la figure 3, la source 49 située en une source-point S Pl produit une onde Sv possédant des fronts qui suivent des trajets rayonnants incidents 50, 51 et des trajets de réflexion 52, 53, respectivement, dans la couverture 54 Des réflexions sur le réflecteur 55 se produisent en des points CRP 1, CRP 2 CRP'n En raison de la production de réflexions d'ondes P de conversion sur le réflecteur 55, l'angle de réflexion (r) de l'onde P de conversion semble être supérieur à l'angle d'inci- dence (i) Il convient également de noter qu'à des angles d'incidence importants, non seulement l'amplitude de l'onde P de conversion augmente, mais, en outre, le point de réflexion CRP' sur le réflecteur 55 n'est pas aligné n avec une projection verticale passant par le point milieu de la paire source-récepteur respective, c'est-à-dire le point situé à midistance entre la source-point S Pl et la position Dn du récepteur. La figure 4 illustre le même principe, mais d'une manière réciproque. Comme représenté, une source 60 d'ondes P, située en une source-point SP 1, apparaît comme produisant des ondes P dont des fronts suivent des trajets rayonnants incidents 61, 62 et des trajets de réflexion 63, 64 dans un terrain de couverture 65 Cependant, les pentes de ces trajets apparaissent comme étant inversées par rapport à celles montrées sur la figure 3 En raison de la pro- duction d'ondes Sv de conversion aux points de réflexion CE? 1,1 CRP' et CRP'n, les angles de réflexion (r) sont inférieurs aux angles d'incidence (i) Mais il convient également de noter que le degré d'inégalité (entre les angles d'incidence et de réflexion) augmente avec le déport L'amplitude de l'onde de conversion augmente de façon similaire. Si l'on revient à la figure 1, il convient à présent de se rappeler que les détecteurs situés dans les stations Dl, D 2 Dm et les sources situées aux sources- points SP 1, SP 2 SP 4 sont utilisés de façon redondante, de sorte que les coordonnées de sources et de récepteurs similaires sont liées à un certain nombre de traces produites En partant de l'activation de la source située à la source-point SP 1, de l'énergie est détectée aux posi- tions des récepteurs Dl Dm o un mouvement du sol est enregistré pendant une période de temps spécifique, souvent 6 secondes Une telle période laisse suffisamment de temps à l'énergie pour descendre et être réfléchie vers le haut sur des réflecteurs, jusqu'au détecteur placé dans les stations D 1 Dm Ensuite, la source est 'avancée" jusqu'à la source-point SP 2 et activée Dans le même temps, les détecteurs sont positionnés dans les stations D 2, D 3 Dm+ 1 pour enregistrer les réflexions Lorsque la séquence décrite ci-dessus est répétée dans le sens indiqué par la flèche 12 le long de la ligne de prospection 13, on obtient, comme résultat, une série d'enregistrements de traces se superposant, pouvant être identifiés avec des coorddnnées redondantes de sources et de récepteurs et des points de réflexion échantillonnés de façon simi- laire, comme décrit précédemment. La figure 5 illustre comment des traces classiques, non converties, associées à des groupes différents de paires sources-récepteurs, échantillonnent le même point de réflexion sur un réflecteur cible. Comme représenté, on peut noter qu'un point commun 68 de réflexion est situé sur un réflecteur plat 69, en alignement vertical avec une projection imaginaire 70 qui passe par le point central Co entre toutes les paires possibles sources-détecteurs Le point commun 68 de réflexion du groupe illustré de paires sources- détecteurs est évidemment obtenu par traçage du rayon de l'onde à partir de sa source-point SPI, S Pl-1, SP 1-2 o. vers le bas, jusqu'au réflecteur 69, puis vers le haut jusqu'à la station particulière de détecteur DJ, Dj+ 1, Dj+ 2 Ainsi, les coordonnées du point milieu commun sont égales à toutes les coordonnées possibles des paires de positions de sourceset de détecteurs, avec la même valeur moyenne Autrement dit: CM Pi,j = CM Pi-l,j+l = CM Pi-2,j+ 2 (I) pour autant de paires qui échantillonnent le môme point de réflexion 68 D'autres observations peuvent être relevées sur la figure 5. Par exemple, on peut noter que, quand bien même les longueurs des trajets des rayons associés aux differentes paires sources-récepteurs sont sensiblement différentes, les longueurs des trajets des ondes incidentes et de réflexion d'une paire quelconque sont identiques. De plus, pour cette même paire, l'angle d'incidence est égal à l'angle de réflexion sur le réflecteur 69. Lorsque les traces associées à ces paires sources-détecteurs sont transformées en groupements de points milieux communs, les traces résultant par sommation sont dites échantillonner le même point de réflexion sur tout réflecteur cible plat, la relation des coordonnées de la réflexion avec les coordonnées de la source-point (SP) et de la station de détecteur (D) étant conforme à une équation de transfor- mation de coordonnées de la forme:- CM Pi,j = l(SP)i + (D)jl/2 (II) La figure 6 montre le fait que des ondes de conversion, enregistrées dans des stations Dj,Dj+i - ne donnent pas des points de réflexion qui sont alignés verticalement avec les coordonnées de points milieux des paires respectives sourcesdétecteurs desquelles les traces sont dérivées. 44869 Comme représenté, les points d réflexion 72, 73, 74, 75, 76 et 77 sur le réflecteur plat 78 ne peuvent pas être alignés avec la projection verticale 79 qui passe par le point central Co de toutes les paires sources-détecteurs se trouvant à la surface de la terre. Les raisons pour lesquelles il en est ainsi seront décrites brièvement ci-dessous. Dans le cas d'ondes de conversion P en Sv, l'angle de réflexion r du rayon Sv de conversion n'est pas égal à l'angle d'incidence i du rayon P, par suite de l'application de la loi de Snell qui s'applique à des rayons en optique, en acoustique ainsi que dans le cas de la propagation d'ondes élastiques dans la terre. La loi de Snell stipule que les angles (i) et (r) sont en relation par la vitesse de propagation des ondes incidentes et réfléchies Dans ce cas, la relation est Sin(r)/sin(i) = Vs/Vp (III) o Vp et Vs sont les vitesses des ondes de compression et de cisaillement, respectivement, dans le terrain de couver- ture situé au-dessus du point de réflexion. Dans des matières solides telles que des roches sédimentaires, la vitesse de l'onde P, à savoir Vp, est toujours Supérieure à la vitesse Vs de l'onde de cisaillement, souvent d'un facteur d'environ 2 Ceci a pour effet que l'angle (r) est toujours inférieur à l'angle d'incidence (i) Par conséquent, les points réels de réflexion 73-77 ne peuvent pas être alignés verticalement avec les coordonnées du point milieu entre les paires sources-récepteurs. Autrement dit: les points réels de réflexion de l'onde de conversion P en Sv tendent à être éloignés de la position du point milieu commun CMP 1, sur le réflecteur 78, l'éloignement correspondant à des valeurs choisies dans la direction des emplacements des détecteurs. Si les données d'ondes de conversion enregistrées par ces paires de sources et de détecteurs sont classiquement décalées dans le temp S pour être amenées en alignement sur la projection 79, il apparaît alors une perte de résolution latérale en raison d'un étalement de ces points sur le réflecteur cible. Conformément à l'invention, l'asymétrie des trajets des rayons incidents et réfléchis est compensée de manière que les points communs de réflexion (CRP) des réflecteurs cibles correspondent exactement à des coordonnées connues des sources-points (SP) et des stations de détecteurs (D) du système de collecte CMP le long de la ligne de-prospection. Brièvement décrit, le procédé de l'invention comprend les deux étapes fondamentales suivantes: La première étape est un calcul du tracé des rayons dans le cas d'ondes classiques (non converties) passant à travers une structure stratifiée afin de donner une section de déport de zéro classique Par itération numérique, une série de réflecteurs cibles est positionnée par exemple à certaines inclinaisons et profondeurs et constitue des événements situés le long de l'axe des temps de chaque trace Les points de réflexion des ondes classi- ques sont calculés par des étapes de transformation de coordonnées de points milieux, les valeurs obtenues étant augmentées pour tenir compte de la profondeur et de l'in- clinaison des réflecteurs Grâce à un artifice géométrique, le trajet du rayon de départ de l'itération est estimé à l'aide de la vitesse efficace de la formule du mouvement hyberbolique. La seconde étape consiste à utiliser l'infor- * mation d'inclinaison et de profondeur du réflecteur de la première étape, décrite ci-dessus, puis à choisir la combinaison appropriée de coordonnées de sources et de détecteurs afin que des groupements de traces échantillonnent le même point de réflexion dans la surface souterraine pour des ondes de conversion Ceci s'avère être, de façon heureuse, une ligne droite inclinée de groupements sur le diagramme de sommation (SP en fonction de D) de la figure 13, tournée par rapport aux lignes de sommation de phase classiques L'exécution de ce groupement de données exige une transformation de coordonnées relative- ment directe, un réadressage du programme qui trie les données par des points de réflexion communs (CRP) plutôt que par des coordonnées de points milieux. Avant la description détaillée du procédé de l'invention, une brève description de deux ensembles différents de systèmes de coordonnées de site, utilisés avec la présente invention, s'avère avantageuse pour la compréhension de certains aspects de l'invention et est donnée ci-dessous en regard de la figure 7. Le haut de la figure 7 est une vue en plan d'un système de collecte CMP similaire à celui illustré sur la figure 1, sauf que les positions des sources-points SP 1, SP 2 précèdent l'avance du dispositif de détecteurs Dl, D 2 Dm, au lieu de suivre le dispositif comme montré précédemment La direction de progression est dans le sens indiqué par la flèche 90 suivant la ligne de levé ou de prospection 91 Par suite d'avances incrémentielles similaires entre le dispositif et les positions des sources, c'est-à-dire chacun est avancé d'une position incrémen- tielle par cycle de collecte, de manière que ASP =AD, des traces peuvent être associées non seulement aux posi- tions respectives de source (SP) et de récepteur (D) au moyen d'axes orthogonaux 92 et 93, respectivement, sur le diagramme 94 de sommation, mais elles peuvent égale- ment être identifiées avec les coordonnées de point milieu commun (CMP) et de déport (f) qui s'étendent suivant des axes orthogonaux 96, 97, conformément aux équations de transformation de coordonnées données précédemment, à savoir CMP = (D + SP)/>, f = D-SP. Etant donné que la transformation de coordon- nées pour des ondes de conversion varie sensiblement par rapport à ces formules, une brève description de la base théorique pour exécuter la transformation de coordon- nées de traces de conversion conformément à l'invention sera à présent donnée. Brièvement décrites, à cet égard, les équations de transformation impliquant des réflecteurs plats seront développées initialement, suivies d'une description dé- taillée de la dérivation des équations de transformation en relation avec l'inclinaison des réflecteurs. Le traçage de rayons d'ondes de conversion dans le cas de réflecteurs plats sera à présent décrit. Dans le cas de phases classiques telles que l'onde P générée puis enregistrée sur des récepteurs réagissant verticalement, les points de réflexion communs (CRP) sont supposés être les points milieux communs (CMP) entre les paires sources-récepteurs associées Pour une coordonnée de déport donné, par exemple pour un déport donné le long de l'axe 97 de la figure 7, des traces associées à un point milieu donné et la coordonnée 96 peuvent être groupées simultanément Autrement dit, dans des zones sans structure ou sans inclinaison, des données CMP, peuvent être convenablement regroupées confor- mément aux coordonnées de points milieux de chaque paire source-récepteur, puis traitées par correction d'augmen- tation du temps due à la distance du point de tir à la trace considérée (correction NMO), par statique et par sommation pour donner des coupes-films finales précises. Cependant, dans le cas d'ondes de conversion, étant donné que les angles d'incidence et de réflexion sont inégaux, les points de réflexion associés à chaque cible ne peuvent être alignés avec les positions des points milieux précités Mais ils peuvent être déterminés avec précision, puis éliminés, conformément aux étapes de la présente invention. La figure 8 illustre des aspects du traçage de rayons de conversion de la présente invention pour l'établissement d'un alignement. Cette figure représente un empilage de couches désignées globalement en 100,à savoir des couches 1, 2, 3, n, dont une partie représente la couverture. La ligne ioi en traits pleins indique le trajet des rayons d'une onde P incidente produite par une source à des coordonnées SP à la surface 102 de la terre, et une onde de conversion P en Sv qui est réfléchie sur l'interface 103 d'une couche N en un point 104 de réflexion La réflexion est enregistrée par un récepteur placé dans la position D. Pour une distance X donnée récepteur-déport entre les coordonnées de la source (SP) et les coordonnées du récepteur (D), les coordonnées du rayon doivent être déterminées et la distance de déport de réflexion Xr peut être déterminée Le problème de conversion Sv-P est une image réfléchie de celui montré sur la figure 8, mais la solution générale est la même pour les angles d'incidence dans chaque couche. En bref, pour trouver les angles d'incidence afin de tracer la carte des rayons, l'invention transforme le traçage des rayons en passant de la résolution d'un problème de réflexion à celle d'un problème de transmis- sion Ceci peut être effectué, comme montré sur la figure 8, par construction d'une image correspondant à la symétrie des rayons Sv réfléchis par rapport à une interface 103 De cette manière, le trajet total des rayons semble être équivalent à la transmission vers le bas à travers deux groupes de couches, c'est-à-dire les couches 1, N et les couches n+ 1, n+ 2, 2 N avant qu'une réflexion se produise au point 105 de réflexion, c'est-à-dire le long des segments descendants de la ligne 101 en traits pleins et de la ligne 106 en traits pointillés, comme représenté. Il convient de noter que le groupe supérieur de couches est associé uniquement à des vitesses d'ondes P dans la couverture, tandis que le groupe inférieur de couches n'est associé qu'à des vitesses d'ondes Sv. Si les couches sont numérotées de haut en bas, comme montré, et qu'on leur attribue les valeurs correctes de vitesses, d'angles et d'épaisseurs, la distance totale X de déport et le temps T de-parcours du rayon sont 2544 a 69 2 n X = E hi tg O i i=l 2 N -l T = Z hîcos ei/vi ( 2) i= 1 Pour trouver les angles d'incidence 81, e 2 i dans les couches, on minimise le temps T de parcours du trajet en se contraignant à un déport X fixe Ceci détermine le paramètre p des rayons à partir duquel tous les angles peuvent être calculés par la loi de Snell. A cet égard, Taner et Koehler ( 1966) décrivent une technique pour tracer uniquement des réflexions clas- siques (mais qui ne convient pas à des ondes de conversion et qui n'utilise pas des trajets de transmission), cette technique étant intéressante ("VELOCITY SPECTRA-DIGITAL COMPUTER DERIVATION AND APPLICATIONS OF VELOCITY FUNCTIONS", Geophysics, Vol 34, pages 859-881). La fonction (F) sera à présent définie. F = T p X ( 3) 2 n F = Z (hivi l{(%tg 2 Oi + 1)1/2 p vit 4 i} ( 4) i=l o p est un paramètre variationnel. Pour minimiser T avec un déport fixe X, on considère la différentielle de F et on lui donne une valeur égale à 0, c'est-à-dire: 2 N 2 &F = Z (hi/vi) {etg ei/(tg 28 i + 1)12 p vi} 6 tg oi i=l Ceci est possible pour des différentiels d'angles arbi- traires seulement si P = 1 i ( 1 + cos 2 Oi)/2 = sin i( 5) c'est-à-dire, la loi de Snell s'applique, o p est le paramètre du rayon, avec des unités d'inertie horizontale apparente Ceci peut également être écrit sous la forme sin ai = p Vi = Vi/C ( 6) o c est la vitesse horizontale inconnue du rayon En substituant ceci aux angles dans l'équation 1, on obtient une équation simple pour c 2 n X=E h 1 vi /(c 2 V)/ i=l 1 i V ( 7) Cette équation n'a pas de solution algébrique pour c, mais elle peut être résolue numériquement par itération, exigeant une évaluation initiale du trajet des rayons, ce qui permet au procédé d'aboutir à une solution par itération Conformément à l'invention, la formule de courbure hyperbolique donne des valeurs de départ étonnamment bonnes pour c, de la manière décrite ci-dessous. En se référant de nouveau à la figure 8, on voit que le trajet de transmission totale à travers les deux groupes de couches est simplement égal à la moitié d'un trajet bidirectionnel classique de réflexion descendant jusqu'à l'interface 2 N et remontant à la sur- face La réflexion remontante arrive avec un déport 2 X. Les formules normalisées pour l'augmentation du temps due à la distance du point de tir à la trace considérée permettent d'estimer les temps des arrivées et les inerties en termes de vitesses efficaces en descente jusqu'à l'in- terface 2 n Des temps unidirectionnels jusqu'à l'interface 2 N avec un déport X sont obtenus à partir de 1 amême formule de courbure en divisant simplement par deux les variables de distance et de temps. Ainsi, dans le problème considéré de convertie, la formule approximative est: T 2 = T 2 + (X/Vp)2 ps s réflexion ( 8) o T est le temps 'e déport nul bidirectionnel cour les réflexions converties, égal à 2 n Tps = E hi/vi i= 1 ( 9) et V est la vitesse efficace égale à: ps 2 N 2 n Vps = l(E V? ti)/(E ti)l/2 i=l 1 i= 1 ( 10) Dans l'expression ( 10), les termes ti sont les temps de couche unidirectionnels donns par: les temps de couche unidirectionnels donnés par ti = hi/Vi ( 11) o les temps P sont obtenus pour i de 1 Sv sont obtenus pour i de n+ 1 à 2 n. A partir de l'équation ( 8), initiale de co peut être donnée pour la apparente du rayon Soit: à N et les temps une évaluation vitesse horizontale co = d X/d T = Vps lT 2 + (X/Vps)2 l/2/ (X/Vps) ps ( 12) Il convient de noter que co est infini à un déport zéro et diminue de façon asymptotique vers Vps lorsque X tend vers l'infini. Pour trouver la valeur exacte de c pour le modèle de stratification, c O est substitué à c dans l'équation ( 7) et des itérations sont effectuées vers une solution ayant la précision demandée. c étant déterminé, le déport Xr du point de réflexion peut être calculé à l'aide de l'équation ( 7), mais par sommation de i seulement 4,smu Ct à n, on obtient n Xr = E hi vi / (c 2-v 2 y/2 f 138 La vitesse efficace de l'onde de conversion peut être mise en relation avec les vitesses efficaces des ondes séparées P et Sv mesurées sur les réflexions séparées non converties Si Tp et Ts sont les temps de déport zéro, bidirectionnels, pour les événements non convertis, la vitesse efficace de l'onde de conversion dans l'équation ( 10) peut être écrite sous la forme V Ps = Ps = ps o V et Vs sont les vitesses efficaces des ondes P et p Sv, respectivement, en descendant vers l'interface n, obtenues à partir de carottes de puits adjacents ou par d'autres moyens classiques. L'expression ( 14) montre que Vps est une moyenne Ps 2 2 pondérée des vitesses moyennes quadratiques V et V de sorte que, en principe, des corrélations entre les phases classiques P et Sv et les vitesses et les temps de conversion à des fins d'identification peuvent être effectuées. Il convient également de noter que le problème des ondes de conversion a à présent été remplacé par un problème de courbure classique équivalente Les varia- tions de vitesses sont tout à fait extrêmes, car les vitesses P et Sv sont une fonction composée. Il est intéressant de noter que dans le cas d'une couche unique, l'équation ( 14) peut être davantage réduite Etant donné que les vitesses efficaces pour de petits déports sont égales aux vitesses réelles P et Sv, on a h 1 = Tp V = T Vs En introduisant ces éealités pp S S dans l'équation ( 14) Vps varie conformément à: ps V = (V { 1/2A 15) ps p S / pour la vitesse efficace de l'onde de conversion à utiliser dans l'équation ( 12). Les coordonnées des points de réflexion peuvent à présent être estimées La translation latérale à partir du point milieu est de préférence utilisée pour une couche unique à l'aide de la loi de Snell. Si le déport horizontal X des coordonnées de la source (D) à la position du détecteur (SP) est défini comme: X = D -SP ( 16) le déport réel Xr du point de réflexion peut alors être défini sous la forme d'un rapport de ces deux quantités: R = Xr/(X/2) = 2 tg ep/(tg ep + tg e SV) ( 17) o ep est l'angle d'incidence et e Sv est l'angle de réflexion La multiplication des membres supérieur et inférieur par cos ep et l'utilisation de la loi de Snell donnent: R = 2 (Vp / Vs)/(Vp / Vs + cos Op/cos esv) ( 18) Par un développement en séries de puissance pour chaque terme cosinus, par l'utilisation de la loi de Snell et par simplification, on obtient: R = 2 (Vp/Vs)/{(Vp/Vs + 1 (ap 2/2)( 1 (Vp/Vs2)} ( 19) o ep est l'angle d'incidence, en radians. P Le rapport approximatif du déport Xr du point de réflexion au déport X/2 du point milieu est d'une préci- sion s'élevant à ep = /4 = 45 degrés Pour des déports P X inférieurs à la profondeur du réflecteur, le terme en ep peut être ignoré, de sorte que: R = 2 (Vp/Vs) / (Vp/Vs + 1) ( 20) Il convient de noter que le rapport R dépend uniquement du rapport (VP/Vs) pour tous les réflecteurs dont la profondeur est supérieure à la distance de déport. L'expression ( 20) peut donc être utilisée pour grouper des traces associées à dos paires choisies sources-récepteurs présentant des événements, qui échan- tillonnent le même point sur une interface, de la manière suivante. Il convient de se rappeler que, dans des ondes classiques, les angles 8 p et esv sont égaux, de sorte que le point de réflexion est: CR Pi, = S Pl + 0, 5 (Dj S Pi) ( 21) = 0,5 (Dj + S Pi) Pour des réflexions de conversion P en Sv, le déport du point de réflexion est: CR Pi,j = S Pl + 0,5 R (Dj SP i) La transformation ci-dessus peut être considérée comme un processus pour déterminer les coordonnées du point de réflexion projetées sur le plan de référence horizontal par multiplication d'une constante (k) qui tient compte du rapport de vitesse de la couverture multi- plié par les coordonnées (SP) et (D) des paires sources- récepteurs respectives, conformément à: = k Dj + ( 1-k)S Pi( 22) o k = R/2 = (Vp/Vs)/(Vp/Vs + 1) ( 23) Il ressort de l'équation ( 22) que les points réels de réflexion pour des ondes de conversion sont égaux à une moyenne pondérée des coordonnées des positions des sources et des détecteurs, les pondérations s'additionnant jusqu'à l'unité Pour des réflexions classiques, les pondérations sont chacune de 0,5 et leur somme atteint 1. Du fait des expressions ( 22) et ( 23), les traces qui sont associées à des paires de coordonnées de sources et de détecteurs peuvent être groupées pour une position connue CRP d'un point de réflexion. Il convient de noter, sur les expressions ( 22) et ( 23), que k n'est qu'une fonction du rapport Vp/Vs pour le terrain se trouvant au-dessus du réflecteur, c'est-à-dire le rapport Vp/Vs pour la couverture, si le réflecteur est plat Cependant, si le réflecteur est incliné d'un angle a par rapport à une ligne horizontale perpendiculaire au champ de gravité terrestre (au lieu d'être plat), les formules correspondant aux coordonnées CRP i du point de réflexion sont plus complexes. Le traçage des rayons d'ondes classiques et converties avec des réflecteurs inclinés sera à présent décrit. La figure 9 montre la relation géométrique des points de réflexion CRP pour des ondes de conversion provenant d'une couche inclinée 110 Le paramètre k dépend dans ce cas, de deux facteurs, à savoir l'angle d'incli- naison a et le rapport du déport à la profondeur X/h Une formule exacte pour k est la suivante: k = 1/2 l 1 + (X/2 h) sin ( 24) l'expression ( 24) précédente étant développée de la manière suivante. Une source présente des coordonnées de posi- tion (SP) et tire vers le bas, sous une certaine incli- naison, dans un récepteur placé dans une position (D). Alors que pour des réflexions classiques, les angles d'incidence et de réflexion de l'interface inclinée 110 sont égaux, ils ne le sont pas dans le cas d'ondes de conversion Le but des expressions suivantes est de trouver une relation entre le déport réel Xr' du point de réflexion le long de l'interface 110 et le déport X'/2 du point milieu le long de la même interface, tels que mesurés à partir de coordonnées s' et d' qui sont, sur cette interface 110, les points des images projetées des points de surface SP et d, respectivement. R peut être défini comme le rapport du déport du point de réflexion au déport du point milieu à la surface de la terre, indiquée en 112, réellement enregistré. Par géométrie, on détermine que le mmie rapport des coordonnées d'inclinaison le long de l'interface 110 est égal à: R = Xr/(X/2) = Xr'/(X'/2) ( 25) Si l'angle d'inclinaison a était zéro, Xr' X'/2 et R serait égal à 1, comme prévu D'après la géo- métrie de la figure 9, étant donné que Xr' = h tg ap X' = ( 2 h + X sin a) tg p ( 26) une substitution de ces expressions donne k, comme suit: k = R/2 = 1/211 + (X/2 h) sin'al ( 27) comme indiqué précédemment dans l'équation ( 24). Pour une inclinaison nulle, l'expression de l'équation ( 27) se réduit à k = 1/2, ce qui signifie que le point réel de réfiexion correspond à la position du point milieu entre la source et le récepteur Pour une inclinaison non nulle, cette équation est commode à utiliser pour l'estimation de la position du point réel de réflexion, en tant que fonction combinée de l'inclinaison, du déport et de la profondeur. Ainsi, étant donné que k = R/2, conformément à l'expression ( 22), les points réels de réflexion CR Pij pour des ondes classiques peuvent être établis comme: CR Pi,j = k D + ( 1-k)Sp; Cependant, l'expression ( 27) ne traite que des réflexions classiques à l'interface 110 et elle n'applique pas les réflexions de conversion. La figure 10 illustre la relation géométrique des trajets des rayons de conversion par rapport à une interface inclinée. Dans ce cas, le paramètre k dépend de trois _ facteurs, à savoir le rapport Vp/Vs, l'angle d'inclinaison a et le rapport du déport à la profondeur X/h Une formule exacte pour k est donc: (Vp/Vs) {(VK/Vs) + l 1 + (X/h) sin el l 1 + fl-1/2} ( 28) o l de la manière suivante. Il convient d'abord de noter, sur la figure 10, que l'angle d'incidence (ep) et l'angle de réflexion (esv) ne sont plus égaux Par conséquent: Xr' = h tg op ( 30) et X' = h tg op + (h + X sin a) tg Osv ( 31) X' = h tg Op ( 1 + ( 1 + (X/h) sin a) tg Osv/tg Op) En les substituant dans l'équation ( 25), on obtient: R = 2/tl + li + (X/k) sinr al tcr sv/tg api ( 32) Hormis le rapport des tangentes, cette formule est similaire à celle établie pour des réflexions classi- ques. L'équation ( 32) est exacte, mais le rapport des tangentes ne peut être déterminé avec précision sans un traçage itératif des rayons Cependant, une simple approximation peut être réalisée et donne des réponses très proches. Si le point de réflexion est déplacé de CRP' à CRP" comme indiqué sur la figure 11, les angles d'incidence et de réflexion augmentent tous deux de ep P à ep,, et de e Sv à esv, respectivement Cependant, le rapport des angles des tangentes change très légèrement, de sorte qu'on obtient une bonne approximation: tg asv/tg ap tg)Sv"/tg p ( 33) Lorsque la profondeur du réflecteur augmente, cette approximation devient plus précise Ainsi, d'après la figure 11, on a: ( 34) cos 8 p" = h/(h 2 + Xr" 2) 1/2 sin op" = Xr"/(h 2 + y 2) 1/2 ( 35) o (pour des couches profondes) Xr" = (Vp/Vs)X cos a/(Vp/Vs + 1) ( 36) La loi de Snell indique que les angles d'in- cidence et de réflexion sont dans la relation: sin e 5 sv" = (Vs/Vp) sin 3 p" ( 37) En utilisant deux fois la loi de Snell, on peut calculer le rapport des tangentes pour les angles perturbés, conformément à: tg esv"/ tg ep" = (Vs/Vp) cos Op"/co S Osv' (Vs/Vp) cos Op" l 1-(V 5/Vp)2 sin 2 opl 11/2 ( 38) En substituant ( 34) et ( 35) dans cette équation pour cos e p" et sin e ", après simplification algébrique, P P on obtient: tg sv"/ tg vp = (Vs/Vp) l 1 + fl-1/2 ( 39) o f = (Xr-/h)2 ( 1-(Vs/Vp)2) ( 40) Pour exprimer f par les termes des distances X de déport de surface, on substitue ( 36) dans cette équa- tion, ce qui donne: =s(vp S j X/h)2 cos 2 a l(Vp/Vs) + 1 l(X/h)2 cs 2 ( 41) Enfin, le paramètre k = R/2 est obtenu par substitution du rapport des tangentes de ( 39) dans ( 32), donnant ainsi: (Vp/Vs) {(Vp/Vs) + l 1 + (X/h)sin alIll + fl- 11/2} ( 42) On peut observer, ainsi qu'il ressort de ce calcul, que les équations ( 41) et ( 42) vérifient les for- mules établies initialement dans les équations ( 29) et ( 28). Bien que la formule donnant k soit complexe, elle est également tout à fait générale et se réduit à des formes plus simples à la fois pour les ondes de conversion et pour les ondes classiques, avec et sans réflecteurs inclinés. En donnant à (Vp/Vs) une valeur égale à 1, on obtient la formule des ondes classiques pour k dans le cas de couchies inclinées (car les ondes d'incidence et de réflexion sont de môme type) On peut voir, d'après l'équation ( 41), que f tend vers O et que, par conséquent, k, dans ( 42), se réduit à: k = 1/2 l 11 + (X/2 h) sin al ( 43) ce qui est en conformité avec ( 24), calculé précédemment pour des ondes classiques. Dans le cas d'ondes de conversion, en établis- sant l'angle d'inclinaison a = 0, on réduit le paramètre k à la formule: (vp/vs) {(Vp/Vs) + 1 +f 11//} ( 44) o l(vp/vs) 1 l f = (X/h)2 ( 45) I(Vp/V 5) + il L'effet de f dans cette expression est de corriger les variations de déport affectant X, qui devien- nent importantes lorsque X atteint une valeur aussi grande que la profondeur du réflecteur h. Si l'utilisateur est intéressé uniquement dans des conditions de déport proche, o le rapport X/h est très inférieur à 1, f peut alors être établi comme étant égal à zéro, entrainant une simplification supplémentaire de k, de la manière suivante: k = Vp/Vs/(Vp/Vs + 1) ( 46) qui est l'approximation la moins complexe dans le procédé de é 'nvention, l'aproximation étant utile pour de faibles déports dans des zones à inclinaison nulle Cette expression pour k a été développée précédemment pour l'équation ( 234. Il ressort également du développement précédent que le paramètre k peut être calculé pour des ondes classi- ques et des ondes de conversion, englobant les effets du rapport du déport à la profondeur (X/h), de l'angle d'inclinaison a, et du rapport de vitesse Vp/Vs Ainsi, k est une fonction générale de ces trois paramèrres: k = k (X/h, Vp/Vs, a) ( 47) qui peut être calculée pour tout cas intéressant Comme décrit précédemment, k est nécessaire pour le calcul des points communs de réflexion (CRP) pour chaque type d'onde. La présente invention sera à présent décrite -en détail Sur la base théorique solide indiquée ci-dessus, on peut décrire les opérations mettant en oeuvre le procédé de l'invention, en regard des figures 12, 13, 14 et 15. La figure 12 montre comment s'effectue la transformation d'une série de traces des terrains confor- mément à l'invention Elle représente un point de vue global. Comme montré, l'organigramme 199 détermine la séquence souhaitée d'étapes commandant le fonctionne- ment d'un calculateur numérique, tel que le modèle 3033 de la firme IBM, réalisant au moins la transformation et la remise en ordrc de traces converties associées à des paires sources-détecteurs dont les positions des coordonnées (SP,D) des stations sources-points-détecteurs sont connues, par l'utilisation d'une équation de transformation pour des points communs de réflexion associés (CRP) choisis dans le groupe comprenant: CRP = k D + ( 1-k)SP: pour des conversions d'ondes P en ondes Sv à la cible; CRP = ( 1-k)D + k SP: pour des conversions d'ondes Sv en ondes P à la cible o k (VP/V 5) l(Vp/Vs) + {l + (X/h)sin a}( 1 +fr 1/2 l f= l(VP/V) (X/h)2 cos 2 a; l(Vp/VS) + 1 l Vp et Vs sont les vitesses des ondes P et des ondes Sv, respectivement, dans le terrain de couverture; X est la distance du déport source-récepteur h est la profondeur du réflecteur cible; a est l'angle d'inclinaison du réflecteur cible; et SP et D sont les coordonnées de la source et du détecteur, respectivement, suivant la ligne de levé ou de prospection. A la suite d'une évaluation centrée sur la détermination de la constante "k" indiquée ci-dessus en fonction des rapports de vitesse dans la couverture, de l'angle d'inclinaison et de la profondeur (h) du réflecteur cible et de la distance de déport X source-récepteur, il apparaît aisément une transformation par laquelle sont produits ensuite des signaux d'une plus grande intelligi- bilité et d'une plus grande clarté pour une interprétation géophysique Il convient de noter que, bien que l'orga- nigramme 199 montre le procédé de l'invention d'une façon générale, un programmeur normalement expérimenté dans les techniques classiques de collecte et de traitement CMP, en plus des équations de transformation mentionnées ci-dessus, peut programmer la plupart des ordinateurs classiques en relativement peu de temps pour atteindre les buts et objectifs de l'invention Un langage typique de programmation pouvant être utilisé est le langage FORTRAN. Bien que l'instruction 200 soit la clé de la transformation des coordonnées de source (SP) et de récepteur (D) en coordonnées CRP par évaluation des équa- tions ci-dessus, on suppose que les-données ont été collectées comme indiqué à l'étape 201, précédem- ment, par des procédés à points milieux communs (CMP) Autre- ment dit, on considère comme un fait établi que les données collectées à l'étape 201 ont été produites séquentiellement par une source sismique classique placée en une série de positions de sources-points (SP), puis qu'elles ont été collectéesde façon redondante sous forme de traces de conversion par une série de détecteurs placés dans des positions connues de détecteurs (D) suivant la ligne de levé ou de prospection Etant donné que les traces converties sont associées chacune à une paire source- détecteur dont les coordonnées de position de la station source-pointdétecteur (SP,D) sont connues, un traitement ultérieur, effectué conformément à l'étape 200, provoque leur transformation aisée sous la forme de coordonnées réelles d'un point commun de réflexion Autrement dit, l'étape 200 établit les coordonnéos c OEP pour un groupement de traces converties afin de corriger dans le temps chaque trace convertie, la correction portant sur (i) toutes différences d'élévation qui ont pu exister entre ses coordonnées de la station source-point-détecteur au moment de la détermination de la trace (corrections statiques), et (iil des différences de courbure associées aux trajets de parcours de l'énergie (corrections dynamiques) Le- résultat est que les traces finales représentent une séquence de collecte imaginaire, comme si une source associée à une trace donnée était placée dans chaque CRP et activée, cette-opération étant suivie immédiatement d'un changement de position d'un détecteur au CRP et de la réception de phase de conversion de l'onde produite. De cette manière, les CRP associés à chaque groupement de traces converties (sous la forme de positions (SP,D) de chaque paire sourcedétecteur et d'une constante "k" évaluée) tiennent compte de trajets de parcours non symé- triques des rayons incidents et réfléchis, ainsi que de l'inclinaison et de la profondeur du réflecteur cible En- suite, les ondes converties font l'objet d'une sommation par une instruction 202, sous forme de coordonnées CRP, puis les groupements sont affichés conformément à une instruction 203-, sous forme d'une section sismique à déport zéro des traces converties. La figure 13 est également intéressante pour l'explication de l'instruction 200, figure sur laquelle est représenté un diagramme de sommation simplifié 204 qui sera décrit en association avec le système de collecte de la figure 7. Comme représenté, chaque point du diagramme 204 présente une coordonnée de source (SP) suivant l'axe 205 correspondant à la position de la source qui donne naissance à la trace correspondant à ce point Le même point présente également une coordonnée (D) suivant l'axe 206 correspondant à la position du récepteur dont le signal de sortie est également -associé à des traces de la même manière En d'autres termes, le diagramme 204 est une courbe des coordonnées SP et D de toutes les traces comprenant une séquence de collecte sismique CMP. La position de l'origine de la ligne sismique se trouve à l'intersection 207 des axes 205 et 206, c'est-à-dire à D = SP = 0 L'exemple montré sur la figure 13 est une ligne sismique constituée de six sismogrammes à 13 traces qui ont été enregistrés avec une géométrie de fin de déroulement, dans laquelle l'écartement incrémentiel des points tels que montrés sur le diagramme 204, à savoir - ASP et AD, est établi à une valeur constante. Pour un déport donné, tel que suivant l'axe 208, des lignes de groupement-communes classiques, géné- ralement indiquées en 209, coupent l'axe 208 perpendi- culairement Par conséquent, les coordonnées (c'est-à- dire "adresses") de traces groupables(alignées suivant de telles lignes 209 de groupement) sont aisément évaluées à l'aide de l'équation: CMP = k D + ( 1-k)SP o k est déterminé de façon à être égal à 0,5. Autrement dit, pour des réflecteurs plats classiques sans conversion de mode, les points de réflexion du réflecteur sont des projections verticales des points milieux entre les positions de paires sourcesrécepteurs associées, produisant la trace Par conséquent, un point commun de réflexion sur un réflecteur cible est associé verticalement à une paire de coordonnées connues (SP) et (D) Par exemple, dans le cas de coordonnées (SP) et (D) égales à 6,(en association avec l'axe de déport 208), les traces correspondantes pour un groupement approprié de- vraient être alignées suivant la ligne de groupements 209 a en traits mixtes, perpendiculairement à-l'axe 208. En outre, à l'intérieur du groupement de traces définies apparaissant sur la ligne 209 a, la distance de déport de la source au récepteur, associée à toute trace du groupe- ment, peut également être déterminée en raison du fait que les coordonnées de la trace de déport la plus longue sont une fonction directe des coordonnées de sa source- point (SP) et de la station de son récepteur (D) Dans l'exmple ci-dessus (portant sur la ligne de groupement 209 a), ceci apparaît évidemment à des coordonnées (SP) = 2 et (D) = 10. Dans le cas d'ondes de conversion arrivant sur des réflecteurs plats, les lignes de groupement ne sont plus orientées suivant les diagonales à 45 degrés indiquées en 209 sur la figure 13, mais elles sont modifiées suivant des jeux de lignes de groupement en relation avec l'évalua- tion de la constante "k" sur la base du type de conversion des ondes élastiques qui se produit et en fonction d'un rapport Vp/Vs choisi dans la couverture, au-dessus du réflecteur cible, c'est-à-dire soit (i) suivant des lignes 211 a, 211 b en traits pleins, associées à la conversion de l'onde P en onde Sv, soit (ii) suivant des lignes 212 a, 212 b, en traits pleins, associées à la conversion de l'onde Sv en onde P. Les pentes des lignes 211 a, 211 b et 212 a, 212 b pour un rapport Vp/Vs choisi, par exemple Vp/Vs = 24, sont identifiées par des évaluations "k" dans lesquelles cette dernière valeur s'avère être égale à 0,73; et 1 =k s'avère être égal à 0,24 Par conséquent, les-coordonnées pour des points communs de réflexion sur une cible plate sont trouvées par l'évaluation de: CRP = 0,73 D + 0,27 SP et CRP = 0,27 D-+ 0,73 SP pour des conversions de P à Sv et de Sv à P respectivement, afin de donner les jeux de lignes 211 a, 211 b et 212 a, 212 b de la figure 13. Il convient de noter que les jeux de lignes de groupement résultantes 211 a, 211 b et 212 a, 212 b, bien que ne coïncidant pas avec des lignes de groupement de points milieux classiques, peuvent néanmoins faire l'objet d'une sommation à l'aide d'un algorithme de sommation défini par les équations mentionnées ci-dessus, améliorées pour la recherche de coordonnées (SP en fonction de D) concernant chacune des lignes de groupement définies, par exemple suivant un trajet bidimensionnel "de passage". A cet égard, la tolérance bidimensionnelle suivante pour limiter la zone de recherche afin d'inclure uniquement les traces transformées dont les adresses les placent dans une "case" choisie de chaque trajet de passage défini, s'est avérée adaptée. Inclusion = données dans le tableau I ci-dessous. TABLEAU I SP = 7, D = 6 SP = 6, D = 6 SP = 5, D = 6 SP = 4, D = 7 SP = 3, D = 7 SP = 2, D = 8. Un processus, décrit ci-dessus, de "formation de cases" est également utile au traitement de données classiques en raison du fait que des irrégularités appa- raissent souvent sur le terrain par suite d'interférences de routes, de fleuves et d'obstacles de culture, par rapport aux positions souhaitées de tir et d'enregistre- ment Les premiers empêchent le développement de la séquence exacte de positions de tir et de détecteurs, montrée sur la figure 13, c'est-à-dire la séquence dans laquelle des lignes de groupement communes peuvent être établies dans les lignes 209 de la figure 13. Dans le cas de réflecteurs inclinés, plusieurs variables supplémentaires associées à l'évaluation de la constante "k" sont déterminées de la manière décrite précédemment, c'est-à-dire l'angle d'inclinaison du ré- flecteur et la profondeur du point de ce réflecteur, par évaluation des traces classiques collectées simultanément. Les figures 14 et 15 illustrent comment la présente invention peut être associée directement à des données sismiques CMP classiques pour tenir compte des réflecteurs à inclinaison. On peut à présent revenir à la description de la transformation des équations de coordonnées pour un réflecteur incliné Il ressort des figures 14 et 15 que des événements apparaissant dans les traces classiques enregistrées, en fonction des coordonnées de la source (SP) et du récepteur (D), sont d'abord évalués Puis l'inclinaison et la profondeur déterminées des réflecteurs cibles sont utilisées dans le traitement des traces converties intéressantes L'évaluation d'inclinaison et de profondeur d'événements intéressants fait l'objet d'une orientation-machine Il est important, pour ces détermina- tions, de maximiser la réponse du signal lorsque l'angle d'inclinaison de l'événement est modifié systématiquement. En référence à la figure 13, le groupement et la sommation des traces classiques provenant de cibles plates s'effectuent suivant les lignes diagonales 209. Par conséquent, l'inclinaison réelle de tout réflecteur d'une manière décrite précédemment en référence aux équa- tions théoriques formant la base de la présente invention, ne change la pente que des lignes correspondant à l'évé- nement particulier intéressant Il convient de noter en particulier que, étant donné que les projections normales des points de réflexion sont elles-mêmes codées et mémorisées, * conformément à des instructions de traitement, la profon- deur et l'inclinaison de ces événements peuvent de la même manière être aisément incorporées dans des codes mémorisés pour être utilisées ultérieurement dans le trai- tement des traces converties. A présent, d'une façon plus détaillée, comme indiqué sur la figure 14, les traces du terrain sont d'abord mises en mémoire Une instruction 300 assume alors la commande et trie les traces classiques en fonction des coordonnées des paires sources-récepteurs à partir desquelles les traces respectives sont déterminées Autrement dit, les traces sont mémorisées pour former la fonction de trace W(SP,D), mais les coordonnées d'adresses de points milieux et de déport de ces traces sont également enregistrées. Ensuite, une instruction 301 transforme les traces triées précédemment en termes de points de réflexion qui sont des projections verticales des points milieux de paires particulières sources-récepteurs associées aux traces intéressantes, sur une série de fenêtres de temps particulières w wl Autrement dit, chaque trace est transformée initialement comme une fonction d'une valeur de coordonnées de points milieux pour former la fonction trace W(y,t) pour une coordonnée de déport fixée (X) La série de fenêtres de traces identifiée avec le même déport fixé {X) est ensuite corrigée dynamiquement et statiquement par une instruction 302, après quoi les traces font l'objet d'une sommation initiale en fonction des points milieux communs (CMP) des paires sources- récepteurs respectives par l'intermédiaire d'une instruction 303. A la réponse négative à une étape de décision 304 en fonction de points d'échantillonnage individuels le long de la fenêtre de trace, une estimation d'incli- naison est fournie par une instruction 306 A l'aide d'une fonction vitesse, la profondeur est également estimée. Ensuite, à une instruction 308, les coordonnées des réflexions de la trace sont de nouveau estimées en fonction nées SP et D de la projection perpendiculaire du point de réflexion sur un horizon préalablement choisi; et (iii) de la profondeur h du point de réflexion Simultané- ment, on peut noter les opérations effectuées par l'ins- truction 309 qui est mise en action pour mémoriser les coordonnées précitées SP et D de la projection, l'angle d'inclinaison a et la profondeur h du réflecteur cible, pendant que la recherche pour la réponse maximale du signal se poursuit. Une itération suivant la boucle 310 se produit jusqu'à ce qu'une réponse affirmative soit donnée à l'étape de décision 304 basée sur le fait que la réponse au signal a été ou non maximisée Au cours de cette itération, l'instruction 309 place simplement en mémoire tampon les données provisoires d'inclinaison, de profondeur et de coordonnées de surface, jusqu'à ce que l'étape 304 libère la boucle 310 des opérations actives A ce moment, les dernières données mises en mémoire tampon sont libérées simultanément vers un enregistreur 307 sur bande pour aider à la poursuite du traitement des données converties, comme décrit ci-dessous Lorsque les dernières données sont enregistrées en 307, la mémoire tampon est effacée. Une instruction 311 prend ensuite en charge la commande L'instruction 311 détermine si la fenêtre de temps w wl est la dernière fenêtre pour l'évaluation de la trace Si une réponse négative est donnée à l'ins- truction 311, l'instruction 312 incrémente la fenêtre de temps d'une unité, et le processus parcourt la boucle 313 de façon à être répété Si une réponse affirmative est donnée à l'instruction 311, l'instruction 313 prend la commande du processus L'instruction 313 incrémente les opérations en faisant progresser le processus vers les coordonnées de groupement suivantes Si la coordonnée de déport précédente était (X), la coordonnée de déport fixée suivante dans le temps serait alors la coordonnée de déport X+ 1. On revient à la boucle 312 Les opérations effectuées par les instructions 301-312, comme décrit précédemment, sont déclenchées de façon séquentielle pour tenir compte des données d'inclinaison, de profondeur et de projection pour le nouveau groupement de traces sur les fenêtres w wl De cette manière, le processus peut être répété et peut revenir plusieurs fois sur lui- même pour générer une série de groupements à différentes coordonnées de déport, avec des données souhaitées d'in- clinaison, de profondeur et de projection. Dans le même temps que l'instruction 313 incrémente le processus pour générer de nouveaux groupe- ments à des coordonnées de déport augmentant progressi- vement, les instructions 314 et 315 font également l'objet d'opérations séquentielles pour (i) mémoriser les groupements en fonction du déport et (ii) produire ensuite une sommation de ces groupenents en tant que sections sismiques classiques Le processus s'achève avec une réponse affirmative à l'étape de décision 316. A présent, en référence à la figure 15, il convient de noter que le procédé de production de la section de déport zéro converti est similaire à celui montré sur la figure 14 Cependant, il existe quelques différences Par exemple, après que les traces converties présentes sur une bande 320 ont été introduites et triées par une instruction 321, l'interpréteur prérègle la source pour coder, par l'intermédiaire d'une instruction 322, la sélection de l'équation de transformation de coordonnées, lors d'opérations 323 ou 324 Comme indiqué précédemment, ces équations sont choisies dans le groupe comprenant CRP = k D + ( 1-k)SP: (pour des conversions d'onde P en onde Sv); et CRP = ( 1-k)D + k SP: (pour des conversions d'onde Sv en onde P). De façon plus détaillée, tandis que l'opération 321 trie les traces de conversion en fonction des coordon- nées des paires sources-récepteurs, c'est-à-dire des coordonnées SP et D, il convient de noter que des adresses de déport de ces traces sont également enregistrées Autre- ment dit, lorsque les traces sont triées pour former la fonction trace W(SP,D), des coordonnées de déport pour une série de fenêtres de temps particulières w wl font également partie intégrante de chaque adresse de trace. Les instructions 323, 324 transforment les traces triées précédemment entermes de coordonnées de points communs de réflexion (CRP) qui sont des projections de surfaces associées à des traces de paires sources- récepteurs particulières, sur la série de fenêtres de t Laps w *,*w et de coordonnées d'adresses de déport X XL. Après que les segments de traces identifiées avec les suivantes (c'est-àdire le même déport fixe X, la même fenêtre de temps W et les mêmes coordonnées CRP) ont fait l'objet d'une correction statique, ces segments de traces font l'objet d'une sommation à une instruction 325 Cette sommation tient évidemment compte des informations d'inclinaison du réflecteur, de profon- deur du réflecteur et de coordonnées de déport du réflecteur, ces données étant fournies par la bande 307 aux opérations 323 et 324 comme montré sur la figure 15. En ce qui concerne l'instruction 325, il convient de noter que les traces transformées sont additionnées suivant des lignes de groupement qui ont été modifiées par une évaluation du facteur "k", augmentée par un principe de recherche de "case" qui recherche dans deux dimensions, suivant le trajet choisi de groupements de passage à évaluation "k", comme indiqué précédemment De cette manière, des traces transformées, qui ont été déterminées par une évaluation du facteur "k" et dont les adresses placent ces traces dans la "case" examinée, sont additionnées Ces étapes de sommation ont été décrites précédemment. Une instruction 326 de décision prend ensuite en charge la commande pour décider soit de continuer le traitement des données par une bôucle d'itération 327, soit de permettre à une instruction 329 de prendre la commande Autrement dit, l'itération du processus par la boucle 327 se poursuit jusqu'à ce qu'une réponse affir- mative soit donnée à l'étape de décision 326, qui demande si les données associées à la dernière fenêtre wl ont été traitées ou non Autrement dit, l'instruction 326 détermine si les données traitées sont associées à la dernière fenêtre wl Si tel est le cas, une instruction 329 prend la commande; dans le cas contraire, c'est-à- dire en cas de réponse négative à l'instruction 326, une instruction 328 incrémente d'une unité la fenêtre de temps, et la boucle 327 est parcourue de façon que le processus soit répété Lorsque l'instruction 329 prend la commande, elle incrémente les opérations en faisant avancer le processus aux coordonnées de groupement suivantes Si la coordonnée de déport précédente est (X), la coordonnée de déport fixée, qui suit dans le temps, devrait donc être une coordonnée de déport X+ 1 La boucle 327 est de nouveau parcourue. Les opérations effectuées par les instructions 322-326, telles que décrites précédemment, sont donc effectuées de façon séquentielle pour produire des données converties de sommation, convenablement transformées. Autrement dit, le processus se répète plusieurs fois de façon à générer une série de groupements obtenus par sommation, à des coordonnées de déport différentes, en même temps que des données de coordonnées souhaitées d'inclinaison, de profondeur et de projection CRP. En même temps qu'une instruction 329 incrémente le processus pour produire de nouveaux groupements par sommation, à des coordonnées de déport augmentant progressi- vement, des instructions 330 et 321 sont également exécutées de façon séquentielle pour (i) mémoriser les groupements par sommation, en fonction du déport et (ii) générer une section de ces groupements en tant que section sis- mique convertie Le processus s'achève par une réponse affirmative à l'étape de décision 333. EXEMPLE On a entrepris un essai sur le terrain dans la Vallée de Sacramento, en Californie Les paramètres concernant les vibrateurs, les récepteurs et l'enregis- trement sont donnés ci-dessous dans le tableau I Io TABLEAU II Source d'ondes P: vibrateurs Mertz, modèle 9 A Diagramme:de balayage: 3 vibrations x 20 balayages sur 69 m B Espacement du vibrateur: 16,5 m C Mouvement: 1,8 m D Fréquences de balayage: 6-56 Hz, relevé E Temps de balayage: 16 S + 8 S d'écoute Configuration de l'enregistrement A Couverture: 2400 % B Dispositif: V P -288 m-1980 m C Intervalle de groupe: 36 m D Intervalles de géophone: 3 m E Rangée de géophones: alignement de camionnettes Instrumentation d'enregistrement A Gain: binaire B Intervalle d'échantillonnage: 4 ms C Longueur d'enregistrement: 24 s D Gain de préamplification: 2 ** 7 E Filtres: 1 Coupure inférieure: inutilisée 2 Coupure supérieure: 83,5 Hz 3 Eliminateur de bande: en service, 60 Hz F Nombre de canaux: 4 x ( 48 données + 4 Aux) Source d'ondes S: vibrateurs Mertz modèle 13/609 A Diagramme de balayage: 3 vibrations x 20 balayages sur 69 m B Ecartement des vibrateurs: 16,5 m C Mouvement: 1,8 m D Fréquences de balayage: 6-28 Hz, relevé E Temps de balayage: 16 S + 8 S d'écoute La figure 16 représente une section de déport zéro classique dérivée des données de terrain de cet essai. De cette section, on peut noter initialement que la production est constituée d'hydrocarbures gazeux provenant de grès Winters, comme indiqué globalement dans le cadre 403 Des interprétations indiquent en outre que le grès est recouvert de schistes delta Iques qui forment une fermeture impénétrable Les points lumineux 404, 405 (à l'intérieur du cadre 403) indiquent la position du gaz retenu. Des événements supplémentaires intéressants sont identifiés en 400,401 et 402, L'événement 400 est situé dans un complexe sédimentaire tertiaire L'événe- ment 401 est situé à la base de ce complexe L'événement 402 s'étend d'un schiste à faible vitesse juste au-dessus du grès Winters intéressant. Les figures 17 et 18 sont des sections par sommation de déport zéro des données converties; sur la figure 17, le son d'une onde P est utilisé; sur la figure 18, on utilise une source d'ondes Sv Les sections de conversion résultantes des ondes P-Sv et Sv-P, res- pectivementsont traitées conformément au procédé selon l'invention montré sur la figure 15. En résumé, il ressort de ces figures que le schiste à faible vitesse situé juste au-dessus du grès Winters, c'est-à-dire l'événement indiqué en 402 sur la figure 16, a été notablement étendu dans la direction latérale des sections, c'est-à-dire dans la région encadrée d'un trait plein 410 Il est encore plus remarquable de noter que les événements 400, 401 At 402 de la figure 18 pré- sentent une t ESQL Ution remarquable, bien que le niveau d'amplitude de la conversion des ondes Sv-P soit théori- quement très inférieur à celui de la conversion des ondes P-Sv montrées sur la figure 17. L'amélioration de la résolution des événements dans les sections données sur les figures 17 et 18 est une fonction directe du rapport Vp/Vs de la couverture. Par conséquent, une détermination correcte des vitesses des ondes P et Sv est nécessaire pour obtenir des sections précises à déport zéro réel. Même avec de telles données, des évaluations de la vitesse horizontale de sommation, nécessaires lors- que l'on utilise les équations de transformation de coordonnées selon l'invention exigent pne bonne estimation initiale de cette valeur A cet égard, l'utilisation de la moyenne géométrique de vitesses d'ondes connues P et Sv s'est avérée très utile, à savoir Vps = La concordance entre la moyenne géométrique des vitesses des ondes P et S et les vitesses de sommation finale pour les événements 401 et 402 des figures 16 à 18 est bonne Sur la figure 19, dans le cas de l'événement 401, seule une petite différence de quelques centaines de mètres par seconde apparaît De façon similaire, sur la figure 20, dans le cas de l'événement 402, de bons résultats analogues sont évidents Il convient également dc noter que les écarts semblent diminuer en passant de la figure 19 à la figure 20 Ceci est prévisible étant donné que le rapport du déport à la profondeur des réflecteurs concernés, à savoir l'événement 401 par rapport à l'évé- nement 402, diminue en fonction du déport. Il va de soi que de nombreuses modifications peuvent être apportées à l'appareil et au procédé décrits et représentés sans sortir du cadre de l'invention. Z 544869 REVENDICATIONS 1 Procédé pour améliorer la résolution de zon- es, sa mes lc'sctêe S aer des procédés classiques de collecte de points milieux communs (CMP), consistant à actionner séquentiellement une source sismique classique ( 28, 29) en une série d'emplacements de sources-points (SP) suivant une ligne de levé ( 13) et à collecter de façon redondante des réflexions à la fois des ondes classi- ques et des ondes qui en résultent par conversion, au moyen d'une série de détecteurs multicomposantes ( 31, 32, 33) placés en plusieurs stations (D) le long de la ligne de levé, chacune des traces résultantes classiques et de conversion étant associée à une paire sources-points/ station de détecteur de coordonnées connues (SP,D), le procédé, qui porte sur le traitement à la fois,des traces classiques et converties d'une manière systématique afin que des groupements classiques et convertis résultant de ces traces échantillonnent chacun un point de réflexion sur un réflecteur cible ( 23) placé dans le sous-sol et commun à chaque groupement, quelle que soit l'inclinaison et la profondeur dudit réflecteur cible, étant caractérisé en ce qu'il consiste ^ (i) à générer des enregistrements sismiques portant sur le terrain et comprenant au-moins des enregis- trements sismiques classiques et convertis séparés,par positionnement et utilisation d'un ensemble de sources et de détecteurs multicomposantes, tels que des coordonnées de sources-points/station de détecteur individuelles puis- sent être associées de façon redondante à un nombre choisi de traces classiques et converties desdits enregistrements, les traces classiques et converties étant les signaux simultanés de sortie des détecteurs, (ii) à établir, pour lesdites traces classiques et converties, une série de points séparés de réflexion commune (CRP), chaque point (CRP) étant associé aux traces de groupement qui lui sont communes, comme si une source associée à une trace donnée d'un groupement commun était placée à chaque point (CRP) et actionnée, suivie immédiatement par un changement de position d'un détecteur à chaque point (CRP) et de la réception d'ondes classiques ou conver- ties constituant ladite trace classique ou convertie; i-i chaque point ICRP) associé au groupement commun de traces classiques étant constitué pour tenir compte de différents trajets de parcours de l'onde dus à la profondeur et à l'inclinaison de chaque réflecteur cible; (iv) chaque point (CRP) associé audit groupement commun de traces converties subissant une transformation dynamique de coordonnées de la source-point et de la station du détecteur pour tenir compte des trajets de parcours non symétriques des rayons incidents et réfléchis, de l'in- clinaison ainsi que de la profondeur des réflecteurs cibles, l'équation de transformation de coordonnées étant choisie dans le groupe comprenant: CRP = k D + ( 1-k)SP: poursdes conversions d'ondes P en ondes Sv à ladite cible; CRP = ( 1-k)D + k SP: pour des conversions d'ondes Sv en ondes P à ladite cible; vv k = (Vp/Vs) k =Ps l(Vp/Vs) + { 1 + {X/h)sin a}(l+f)-l/2 l; o l{(Vp/vs)11 f= P (X/h) cos a; l(Vp/Vs) + 1 l Vp et Vs sont les vitesses de l'onde P et de l'onde Sv, respectivement, dans le terrain de couverture; X est la distance de déport source-récepteur h est la profondeur du réflecteur cible; a est l'angle d'inclinaison du réflecteur cible; et SP et D sont les coordonnées de la source et du détecteur, respectivement, suivant la ligne de levé 2 Procédé pour améliorer la résolution de données s Lsmiquas coilectées ar dcss proce_és c asscques de collecte à points milieux communs (CMP) consistant à actionner séquentiellement une source sismique classique ( 28, 29) en une série d'emplacements de sources-points (SP) suivant une ligne de levé ( 13) et à collecter de façon redondante des réflexions à la fois des ondes clas- siques et de leurs phases de conversion, au moyen d'une série de détecteurs multicomposantes ( 31, 32, 33) placés en plusieurs stations (D) de détecteurs le long de ladite ligne de levé, afin que des groupements classiques et convertis résultant de ces traces échantillonnent chacun systématiquement un point de réflexion d'un réflecteur cible ( 23) placé dans le sous-sol et commun à chaque groupement, quelle que soit l'inclinaison et la profondeur du réflecteur cible, le procédé étant caractérisé en ce qu'il consiste (i) à générer des enregistrements sismiques de terrain comprenant au moins des enregistrements sismiques classiques et convertis séparés, par position et utilisa- tion d'un ensemble de sources et de détecteurs multicompo- santes tel que les coordonnées individuelles de sources- points/stations de détecteurs (SP,D) puissent être associées de façon redondante à un nombre choisi de traces classiques et converties desdits enregistrements, ces traces converties et classiques constituant les signaux simultanés de sortie des détecteurs, (ii) à établir séparément, pour lesdites traces classiques et converties, une série de coordonnées de points communs de réflexion (CRP) associées chacune à un groupement de traces comme si une source associée à une trace donnée était placée à chaque point (CRP) et activée suivie immédiatement par un changement de position d'un détecteur en chaque point (CRP) et par la réception d'ondes classiques ou converties constituant ladite trace classi- que ou convertie (iii) chaque point (CRP) associé aux traces classiques étant établi de façon à tenir compte des trajets de parcours différents de l'onde, par suite de la profondeur et de l'inclinaison de chaque réflecteur cible; (iv) chaque point (CRP) associé aux traces converties subissant au moins une transformation dynamique des coordonnées des sources-points et des stations de détecteurs pour tenir compte des trajets de parcours non symétriques des rayons incidents et réfléchis, de l'incli- naison ainsi que de la profondeur des réflecteurs cibles, à l'aide d'une équation de transformation de coordonnées choisie parmi: CRP = k D + ( 1k)SP: pour des conversions d'ondes P en ondes Sv à la cible; CRP = ( 1-k) D + XSP: pour des conversions d'ondes Sv en ondes P à la cible; k (Vp/Vs) _ = l(Vp/Vs) + { 1 + (X/h)sin x}(l+f)-l/2 l; o f l(Vp/Vs) l (/h 2 Cs 2 f p (X/h)2 CO 52 a l(Vp/VS) + 1 l Vp et Vs sont les vitesses des ondes P et des ondes Sv, respectivement, dans le terrain de couverture; X est la distance de déport source-récepteur; h est la profondeur du réflecteur cible; a est l'angle d'inclinaison du réflecteur cible; et SP et D sont les coordonnées de la source et du détecteur, respectivement, suivant la ligne de levé 3 Procédé selon l'une des revendications 1 et 2, caractérisé en ce que les valeurs d'inclinaison et de profondeur utilisées par ladite équation choisie de transformation de coordonnées sont établies comme sous- produits de traitement de transformation et d'amplification des traces classiques. 4 Procédé selon la revendication 3, caracté- risé en ce que le traitement de transformation et d'ampli- fication des traces classiques consiste à établir des coordonnées de points milieux communs pour lesdites traces classiques sous la forme d'une fonction des coordonnées de la source-point et de la station de détecteur (SP,D) de chaque trace, à grouper lesdites traces suivant des signes établis de groupement des points milieux communs, à réaliser une sommation des traces groupées pour former des sections à déport fini, puis à estimer la profondeur et l'inclinaison à partir desdites sections d'événements et à estimer ensuite de nouveau les points (CRP) réels à partir des données d'inclinaison et de profondeur, à grouper de nouveau lesdites traces classiques suivant des lignes de groupement modifiées, notamment pour établir leurs points (CRP) réels, et à effectuer de nouveau une sommation des traces regroupées jusqu'à ce qu'une section de déport zéro réel soit formée, les estimations finales de profondeur et d'inclinaison des événements étant en- suite utilisées dans ladite équation de transformation de coordonnées de l'étape (iv). Procédé selon la revendication 4, caracté- risé en ce qu'il consiste également (v) à réaliser une sommation des traces conver- ties associées à chaque point (CRP) établi, en séquence, afin de former une section convertie zéro réelle présen- tant un étalement horizontal insignifiant des événements en haute intensité se produisant le long de cette section. 6 Procédé selon la revendication 4, caracté- risé en ce que l'étape consistant à grouper de nouveau et effectuer une nouvelle sommation des traces classiques jusqu'à ce qu'une section de déport zéro réelle soit formée est une opération d'itération qui se termine lorsque le signal produit en réponse audit événement intéressant est maximal. 7 Procédé selon l'une des revendications 1 et 2, caractérisé en ce que ladite étape (ii) consiste à corriger statiquement les traces classiques et converties pour normaliser les coordonnées des sources-points et des stations de détecteurs sur-un plan élévationnel horizon- tal commun, s'étendant le long de la liane de levé 8 Procédé selon l'une des revendications 1 et 7, caractérisé en ce que le réflecteur cible est sensi- blement parallèle au plan horizontal de référence, une équation de transformation de coordonnées pour les traces classiques étant de la forme: CR Pclassique = 0,5 SP + 0,5 D o Vp/Vs de la couverture est égal à l'unité, et SP, D sont les coordonnées de sources-points et de stations de détecteurs, respectivement. 9 Procédé selon la revendication 8, caracté- risé en ce que la constante (k) de l'équation choisie de transformation de coordonnées pour lesdites traces converties satisfait l'égalité suivante: k = (Vp/Vs) (Vp/Vs + 1) o Vp/Vs de la couverture peut notamment être déterminé expérimentalement ou par des données empiriques. Procédé selon l'une des revendications 2 et 3, caractérisé en ce que les valeurs d'inclinaison et de profondeur du réflecteur cible sont en relation mutuelle telle que la constante (k) de l'équation choisie de trans- formation dynamique des traces converties puisse prendre la forme approchée suivante: k = (YV/VS) l(Vp/Vs) + ( 1 + (X/h sin a)l 11 Procédé pour accroître la résolution d'évé- nements d'amplitude dans des enregistrements sismiques réalisés par des procédés de collecte (CMP), consistant à générer des ondes classiques à l'aide d'une source sismique classique ( 28, 29) en une série de positions de sources-points (s), et à collecter de -eçsri ec'c 4 a nte des phases converties des ondes classiques sous forme de traces converties en une série de détecteurs ( 31, 32, 33) positionnés dans des emplacements connus (d) le long d'une ligne de levé ( 13), le procédé, qui a trait au traitement des traces converties d'une manière efficace afin que des trajets de parcours non symétriques des rayons incidents et réfléchis desdites ondes classiques générées soient pris en compte avant la sommation des traces, étant caractérisé en ce qu'il consiste: (i) à activer séquentiellement ladite source classique à ladite série d'emplacements de sources-points (ii) à collecter de façon redondante au moins des phases converties des ondes classiques dans lesdites stations connues (d) de détecteurs par l'intermédiaire desdites séries de détecteurs afin de produire une série de traces converties associées chacune à une paire source- détecteur occupant des positions dont les coordonnées des sourcespoints/stations de détecteurs sont connues (s,d) (iii) à traiter lesdites traces converties afin de définir une série de coordonnées de points communs de réflexion (CRP) associées chacune à un groupement de traces converties, comme si une source associée à une trace donnée était placée dans chacun desdits points (CRP) et activée, suivie immédiatement par un changement de posi- tion d'un détecteur en chaque point (CRP) et la réception de phases converties de l'onde générée, constituant ladite trace, chaque point (CRP) subissant au moins une trans- formation dynamique des coordonnées des sources-points et de la station de détecteur afin de tenir compte des trajets de parcours non symétriques des rayons incidents et réfléchis, de l'inclinaison ainsi que de la profondeur du réflecteur cible, l'équation de transformation des coordonnées étant choisie parmi les suivantes: CRP = kd + ( 1-k)s: pour des conversions d'ondes P en ondes Sv à ladite cible; CRP = ( 1-k) + ks: pour des conversions d'ondes Sv en ondes P à ladite cible; o (vp,/vs) k= k = _V V 5 k = l(Vp/Vs) + { 1 + (X/h)sin a} (l+f)l-l/2; (Vp/Vs) + 1 f _(Vp/Vs) + 1+ (X/h)2 cos 2 a; (Vp/Vs) 1 Vp et Vs sont les vitesses des ondes P et des ondes Sv, respectivement, dans la couverture; X est la distance de déport source-récepteur; h est la profondeur du réflecteur cible; a est l'angle d'inclinaison du réflecteur cible; et S et d sont les coordonnées de la source et du détecteur, respectivement, le long de la ligne de levée. 12 Procédé selon la revendication 11, caracté- risé en ce que l'inclinaison et la profondeur du réflecteur cible sont en relation mutuelle telle que la constante "k" de ladite équation choisie de transformation dynamique présente la forme approchée suivante: k= (Vp/vs) t(Vp/Vs) + (l + X/h sin a)2 l 13 Procédé selon la revendication 11, caracté- risé en ce que le réflecteur cible est sensiblement parallèle au plan de référence horizontal, ladite constante "k" de l'équation choisie de transformation étant de la forme: k (VP/VS) U(Vp/V 5) + 1 l 14 Procédé selon la revendication 11, carac:ê- risé en ce qu'il consiste également - (a) à établir une série de trajets de passage de groupement commun imaginaires sur un diagramme de sommation classique qui identifie des traces par coordon- nées orthogonales sources-points/âtations de détecteurs, chaque trajet de passage ayant une pente sensiblement définie par ladite équation de transformation de coordon- nées; et (b) à grouper en séquence, le long de chaque trajet, lesdites traces corrigées statiquement, dont les coordonnées des sources-points et des stations de détecteurs les placent sur ce trajet afin d'établir une série de groupements de traces converties et corrigées. 15 Procédé selon la revendication 14, caracté- risé en ce qu'il consiste en outre à réaliser une sommation des traces de chacune desdites séries de groupements pour former une section sismique à sommation, déport zéro réel de traces converties. 16 Procédé pour accroître la résolution d'évé- nements d'amplitude dans des enregistrements sismiques produits par des procédés de collecte (CMP), consistant à produire des ondes classiques à l'aide d'une source sis- mique classique ( 28, 29) en une série d'emplacements (s) de sourcespoints, et à collecter de façon redondante des phases converties desdites ondes classiques sous forme de traces converties en une série de détecteurs ( 31, 32, 33) positionnés dans des emplacements connus (d) le long d'une ligne de levé ( 13), de manière que des trajets de parcours non symétriques des rayons incidents et réflé- chis des ondes classiques générées soient pris en compte avant la sommation des traces, le procédé étant caractérisé en ce qu'il consiste: (i) à activer séquentiellement ladite source classique dans ladite série d'emplacements des sources- points; (ii) à collecter de façon redondante au moins des phases converties desdites ondes classiques dans les- dites stations (d) connues des détecteurs par l'intermé- diaire de ladite série de détecteurs afin de produire une série de traces converties associées chacune à une paire source-détecteur dans des positions ayant des coordonnées connues (s,d) de sources-points et de stations de détecteurs; (iii) à établir, pour chacune des traces conver- ties, des coordonnées d'un point commun de réflexion (CRP) afin qu'un groupement de traces converties puisse être associé de la même manière à ce point, comme si une source indexée sur une trace convertie donnée était placée audit point (CPR) et activée, suivie immédiatement par un changement de position d'un détecteur audit point (CRP) et par la réception de phases converties de ladite onde générée comprenant la trace, chaque point (CRP) subissant au moins une transformation dynamique des coordonnées de la source-point et de la station de détecteur pour tenir compte des trajets de parcours non symétriques des rayons incidents et réfléchis, de l'inclinaison ainsi que de la profondeur du réflecteur cible, l'équation de la transformation de coordonnées étant choisie parmi les suivantes: CRP = kd + ( 1-k)s: pour des conversions d'ondes P en ondes Sv à ladite cible; CRP = ( 1-k) + ks: pour des conversions d'ondes Sv en ondes P à ladite cible; o (Vp/Vs) k = l(Vp/vs) + { 1 + (X/h)sin a} (l+f)l-1/2; (Vp/Vs) + 1 f p S + (X/h)2 cos 2; (V /Vs) 1 p' Vp et Vs sont les vitesses des ondes P et des ondes Sv, respectivement, dans la couverture X est la distance de déport source-récepteur h est la profondeur du réflecteur cible a est l'angle d'inclinaison du réflecteur cible; et s et d sont les coordonnées de la source et du détecteur, respectivement, le long de la ligne de levé