FR 2572189 A1 19860425 FR 8416149 A 19841022 SCHLUMBERGER Aff. 20 Procédé pour la détection de fractures du sous-sol à l'aide d'un 'outil acoustique. L'invention concerne la détection de fractures que présente le sous-sol, à partir d'un puits foré dans celui-ci, La détection de fractures est-c intéressante à de nombreux ti- tres. En prospection pétrolière, elle est particulièrement importante: les zones de fracture ont de fortes chances de constituer des zones de production d'hydrocarbures, car le système de fractures procure un trajet de haute perméabilité, et constitue par conséquent un mécanisme de drainage efficace, 1 à travers une roche réservoir qui, pour sa part, est dense et possède une faible perméabilité. Beaucoup d'outils actuellement connus sont susceptibles, après interprétation convenable des mesures, de donner des indications sur la présence de fractures. Aucun d'entre eux ne le fait de manière véritablement satisfaisante. Parmi les outils disponibles, une attention particulière a été portée aux outils à émission/réception acoustique (Sonic Tool)o Des études ont été menées dans le but d'évaluer la possibilité de repérer en pratique des fractures avec un tel outil. A cet égard, il est fait référence aux publica- tions suivantes R.oLo Morris, DoRo Grine and ToEo Arkfeld, 1964. "The use of compressional and shear acoustic amplitudes for the location of fractures": J. Pet. Tech., 16, 623-635. E.A. Koerperich, 1978, "Investigation of Acoustic Boundary Waves and Interference Patterns as Techniques for Detecting Fractures": J. Pet. Tech., 30, 1199-1207. F.L. Paillet, 1980, "Acoustic Propagation in the Vicinity of Fractures which intersect a Fluid-filled Borehole": Proceedings SPWLA Symposium, 1980, paper DD. Les procédés connus jusqu'à présent font largement appel à l'interprétation par un spécialiste d'une représentation con- venable des mesures effectuées par l'outil sonique. L'homme de l'art connaît les inconvénients d'un tel procédé. C'est pourquoi il est vivement souhaitable d'accéder à une détection de fractures qui puisse être automatisée, fiable, et autant que possible dépourvue d'ambigulté. La présente invention a essentiellement pour but de fournir des moyens qui permettent d'atteindre cet objectif. L'un des buts de la présente invention est d'obvier aux inconvénients des techniques antérieures de détection de fractures par atténuation d'amplitude des ondes sonores. Un autre but de l'invention est d'obtenir une identification de fractures avec une résolution de profondeur qui soit bien meilleure que l'espacement entre émetteur et récepteur de sons. Un troisième but de l'invention est de permettre l'évalua- tion objective de configurations à trajets multiples qui permettent l'identification de fractures. Un autre but encore de l'invention est de fournir aux spé- cialistes des informations et/ou une représentation amélio- rées des mesures afin de bénéficier de la plus grande préci- sion qu'apporte l'invention. Le procédé proposé est du type comprenant les étapes sui- vantes: a) descendre dans le puits, de manière repérée, un outil émetteur/récepteur acoustique possédant des transducteurs séparés pour l'émission et la réception avec au moins un espacement entre émetteur et récepteur; b) au cours de la'descente, produire répétitivement I'émis- sion acoustique de l'outil, et enregistrer à chaque fois la- forme d'onde d.échos WFn reçue par l'outil,en relation avec le ou les espacements, ces échos résultant du couplage acous- tique de l'outil au sous-sol à travers le fluide emplissant le puits. De façon générale, l'invention est remarquable en ce que la détection des fractures est effectuée par traitement de portions de formes d'onde WFn qui sont représentatives de l'un au moins des deux effets locaux corrélatifs des frac- tures que sont les pointes d'amplitudes dues à une élévation locale du couplage acoustique, et les configurations croisées (en anglais "CrissCross patterns") dues à une conversion de mode. Selon un premier aspect de l'invention, le procédé comporte en outre les étapes suivantes c) extraire de chaque forme d'onde WFn, dans une fenêtre choisie de celle-ci, contenant les échos en ondes de com- pression, une grandeur d'amplitude de compression An, ce qui fournit au moins une courbe d'amplitude de compression en fonction de la profondeur An(z), associée à un demi- espacement sn entre émetteur et récepteur; d) isoler de la courbe d'amplitude de compression An(z), les pointes d'amplitude de largeur voisine d'une valeur pré- déterminée, ce qui fournit au moins une courbe de pointes d'amplitude Pn (z); et e) déterminer la présence de fractures à partir des courbes Pn(z+s) d'une part et Pn(z-s) d'autre part. Ceci peut être appelé "analyse par les pointes d'amplitude". De préférence,: - l'étape a) consiste à descendre dans le puits un outil apte à opérer conjointement selon un nombre m d'espacements, Sn désignant le demi-espacement (avec n = 1 à m); - l'étape b) comprend l'enregistrement sélectif des formes d'ondes FWn reçues à chaque émission et pour chaque espacement - les étapes c) et d) sont exécutées sélectivement pour cha- que espacement, fournissant un jeu de courbes d'amplitude de compression An(z), puis un jeu de courbes de pointes d'amplitude Pn(z) respectivement associées à chacun des dèmi- espacements sn; et - l'étape e) comprend: 25.el) la détermination de deux courbes de coefficients de corrélation croisée T(z) et R(z), en fonction de la profon- deur, selon les relations: n = m I-T (z)] = Prod Pn(z+sn)] et n = 1 n = m 1-R (z)]m = Prod Pn(z-sn)] n = 1 e2) la détermination de la présence de fractures à partir des crêtes de la courbe T(z) d'une part, de la courbe R(z) d'au- tre part. Très avantageusement, l'étape d) comprend dl) la détermination de courbes lissées Mn(z) par filtrage en moyenne mobile de chacune des courbes d'amplitude An(z);et d2) la détermination de la courbe de pointes Pn(z) par sous- traction stricte de Mn(z) dans An(z) avec Pn(z)=An(z)-Mn(zo) si An(z) > Mn(z) et Pn(z) = 0 dans le cas contraire. Dans un mode de réalisation particulier de l'étape dl), le filtrage en moyenne mobile est réalisé à l'aide d'un filtre de m6dianie opérant sur 5 points0la répétition des émissions acoustiques se faisant tots les 15 cm environ. En pratiqcue, la valeur prédéterminée pour la largeur des pointes dVamplitude correspond à environ une longueur d'on- de de l'onde acoustique de compression dans le sous-sol, en particulier 0, 3 mètre pour une onde acoustique de 16 kHzo Par ailleurs, il est souhaitable que m soit au moins égal à 3o Selon d'autres particularités de l'inventions à l'étape c), la grandeur Anl est définie par le maximum de l'amplitude crête à crête de la forme d'onde FWn, sur une partie au moins de sa fenêtre contenant les échos en ondes de com= pressiono La grandeur d'amplitude de compression est de nature loga- rithmicqueo Pour une observation par un spécialiste, l'étape e) comprend l'établissement d'une représentation conjointe des courbes T(z), R{z) et de préférence Pn(z), en fonction de la pro- fondeuro Selon un autre aspect de l'invention, le traitement relatif aux configurations croisées comprendladétermination de l'amplitude ou de l'énergie propre à celles-ci, relativement à l'amplitude ou l'énergie des échos directs. - Très avantageusement, les portions de formes d'ondes prises en compte sont celles qui sont dues aux ondes de compression, et se terminent au premier écho du type onde de cisaillement ou ondes pseudo-Rayleigh. De manière plus précise, le traitement relatif aux configura- tions croisées comprend la:sommation d'amplitudes tirées des formes d'ondes, et relatives à des positions attendues pour les deux bras d'une configuration croisée, à l'intérieur d'une fenêtre qui recouvre l'espacement émetteur/récepteur, afin de déterminer l'énergie cumulée EX des échos sur les deux bras de la configuration croisée, la présence d'une fracture étant déterminée d'après l'écart de cette énergie cumulée par rapport à zéro. L'énergie cumulée EX est déterminée de préférence comme la moyenne géométrique des énergies E+, E.- relevées sur les deux bras. On note maintenant par S(t,p) le carré de la somme des ampli- tudes selon une ligne qui intersecte le temps t avec une pente p. Les énergies sur les bras ascendants E+ et sur les bras descendants E- de la configuration croisée peuvent être déterminées par les relations: E+ = somme S(t,p) /N+ (p >O) o N+ est le produit du nombre de niveaux ou échantillons dans la fenêtre d'accumulation par le nombre de lignes d'accumulation 'ascendantes' considérées et g -;f -z; :0:.: 2572189 : 97-:- E- Somme S(t,p) IN- - (p o N est le produit du nombre de- niveaux ou échantillons dans la fenêtred'accumulation par le nombre de lignes d'accumulation 'descendantes'= considérées. - - i 'L - -:f --3 - 7 C-. n Dans ces conditions, la-gamme des valeurs-de t utilisées coïncide avec la fenêtre des-échos-de compre ssion, tandis que la gamme- des valeurs de p coincide avec la pente des lignes joignant les échos:de compression aux échos de cisaillement.. Dans un mode de réalisat-ionparticulier -de -'inve-ntioni le traitement relatif aux configurati-ns croisées-comprend les étapes suivantes: - - c') définir une--fenêtre de sommation couvrant- 'espacement émetteur/récepteur et s'étendan-t- entre les- échos de com-: pression et la montée du premier écho pseudo-Rayleigh ou de cisaillement;: - --:: : : - :: - - ?.0 ' d') normaliser l'amplitude de chaque forme d'onde par rap- port à. l'écho direct le plus -intense qu'"elle contient; : e') éliminer les échos directs- par-filtragede.médiane dans la direction profondeur;;.::: : - - f f) déterminer les -amplitudes somméescomme -défini plus: haut pour un jeu choisi de L-ign-es de sommation possédant : des pentes et intersections temporelles correspondant -à une configuration croisée -dans ladite fenêtre,:- - - g') déterminer les grandeurs E+, E- et EX comme défini - plus haut, et h')- répéter les étpfcf- = =: lfen -ê: sommation. h') répéter lees étapes c'-):Jg') pour Tles différentes poSi- tions admissibles de la fenêtre de sommation.-.-::- 25721-89 - X 0 X f f8 Le procédé peut alors se compléter-d'une visualisation; de l'intensité- des configurations croisées, ou mieux d'une - visualisation conjointe des informations tirées- de l'éléva- __ tion locale du couplage acoustique et des configurations- croisées. D'autres caractéristiques- et avantages de l'invention appa- raltront=-àl'examen.de la description-détaillée ci-après,: ainsi que des dessins-annexés, sur lesquels:( - - la figure 1 est-une illustration schématique d'un outil dit LSS, propre à Ia mise-en oeuvre de l'invention; - {. - - les figures 2A et 2B sont des-chronogrammes illustrant- respectivement une-forme d'-onde d'émission FE et une forme - d'onde d'écho WFn; - - la figure 3 est un diagramme illustrant les différentes étapes de la mise en oeuvre du procédéad'analyse parf les pointes d'amplitude, dans unmode de réalisation particulier de celui-ci;::: - la figure 4 illustre, sous forme d'une représentation plane. à effet tridimensionnel,-un jeu de formes d'ondes reçues WFn,- : - ::.: - les figures 4a et 4b illustrent des formes d'ondes-reçues, ainsi que le traitement que celles-ci subissent par la mise en oeuvre du procédé d'analyse des pointes d'amplitude; -: - la figure 5 est un dessintrès schématique -illustrant certains phénomènes qui se déroulent au niveau d'une fracture; - les figures 5a et 5b sont des représentations-de signaux captés qui mettent en évidence certains phénomènes util-es- selon l'invention,en relation avec les configurations croisées; - la figure 6 est un schéma illustrant le procédé de traite-. ment des configurations croisées selon l'invention, dans un exemple particulier; - la figure 7 montre certains résultats obtenus par le procédé de la figure 6; et - les figures 8a et 8b illustrent des résultats obtenus par mise en oeuvre conjointe des procédés qui forment les deux aspects de la présente invention. Les dessins annexés sont porteurs d'informations qu'ils peu- vent seuls délivrer complètement Ils sont donc incorporés à la présente description, non seulement pour faire compren- dre l'invention, mais aussi pour contribuer-le cas échéant à la définir. La figure 1 illustre un outil émetteur/récepteur acoustique du type dit LSS (pour "Long Spacing Sonic Sonde", c'est-à- dire sonde acoustique à grand espacement). Axialementcette sonde comporte deux émetteurs acoustiques TI et T2, en principe omnidirectionnels latéralement, es- pacés de deux récepteurs acoustiques également omnidirec- tionnels R1 et R2. Au-dessous des émetteurs est prévu un dispositif de déclenchement de l'émission et d'enregistrement des échos, noté globalement DEo Cet outil permet de définir plusieurs espacements, à savoir: - un premier espacement 2os1 entre'll'émetteur T l'et le récepteur Ri; - un autre espacement 20s2 entre l'émetteur Tl et le récep- teur R2, de même qu'entre l'émetteur T2 et le récepteur R1;et - un troisième espacement 20s3 entre l'émetteur T2 et le récepteur R2. Les demi-espacements sl, s2 et s3 ont par exemple pour valeurs 1,22 mètre (4 pieds), 1,52 mètre (5 pieds) et 1,83 mètre (6 pieds). La figure 2A illustre une forme d'onde d'émission qui sera appliquée à l'un et l'autre des émetteurs T1 et T2. La forme d'onde est appliquée successivement aux deux émetteurs, à des instants suffisamment distants pour éviter toute interférence au niveau des récepteurs R1 et R2, mais aussi suffisamment rapprochés pour que l'on puisse considé- rer que les échos obtenus sont relatifs à la même profon- deur, compte tenu de la vitesse de déplacement de l'outil. Par ailleurs, la cadence de ce déclenchement séquentiel des deux émetteurs est choisie telle que l'on obtienne des infor- mations sensiblement tous les 15 centimètres, dans l'échelle des profondeurs. Les émetteurs T1 et T2 produisent des ondes acoustiques de compression, également dénommées "longitudinales", par oppo-. sition aux ondes de cisaillement ou "transversales". En ef- fet, seules les ondes longitudinales sont de nature à tra- verser le fluide qui remplit un puits sans atténuation excessive. De même, les récepteurs R1 et R2 vont recevoir à travers le fluide du puits des échos sous forme d'ondes longitudinales. Mais, par l'interaction des parois du puits avecle sous-sol, ces échos pourront correspondre à différents autres types d'ondes qui se seront propagées dans le sous-sol. Ainsi, la figure 2B montre qu'après l'impulsion d'émission, la forme d'onde WFn d'un signal reçu comprend tout d'abord une portion P1 qui est relative à des échos de compression, dûs notamment à l'aller-retour des ondes de compression à l'intérieur du fluide contenu dans le puits, après réflexion par les parois de celui-ci. 2 2572189 i: : ': 7: - A:. a Il intervient ensuite, dans une section p2-de ia forme d'on- de WFn des échos de cisaillement, dits également de Rayleigh. Une partie significative du signal émis sous-forme d'onde de compression va entrer dans lee sous-sol,-et -se propager dans celui-ci sous la forme d-'ondes de -cisaillement, qui -peu- vent pénétrer l' intérieur de la matière solide,- et égale- ment sous la forme d'ondes dîites de Rayleigh,; qui,-pour leur part,se propagent au voisinage d'un- interface solide/liquide (ou solide/gaz). Ces ondes se distinguent les unes des autres par des vitesses de propagation différentes, -d'oi une durée différente de leurs alternances. Il existe également des ondes dites pseudo-Rayleigh, qui apparaissent du fait que -l'énergie acoustique tend à. s 'écarter -de- 1 interface:(donnant lieu aux ondes de Rayleigh) dès lors que des ondes de com- pression s'établissent dans la boue de forage.- La portion p3 de la forme d'onde WFn fait ensuite apparaitre-des.échos dits de boue. L'homme de l'art sait que, à l'inverse des ondes de-compres- sion et de cisaillement, quifse propagent dans un milieu infiniles ondes de Rayleigh et pseudo-.Rayleigh ne,se pro- pagent que dans un milieu fini, c'est-à-dire qu'il s'agit d'ondes guidées. - t Un troisième type d'ondes guidées apparait dans un forage, et se manifeste-dans la portion P4 de la forme-d'onde WFn. Il s'agit des échos ou ondes de Stoneley, qui s'établissent dans la boue par interaction-entre celle-ci et les formations du sous-sol. Les fréquences les plus basses des ondes de Stoneley sont également appelées ondes tubulaires. - 0: Il a déjà été proposé d'utiliser l'atténuation qui se'mani- feste dans les formes d-'ondes enregistrées par un outil - acoustique, dans le but de repérer les fractures présentes dans le sous-sol.- - - Mais on comprend aisément, compte tenu notamment du grand 2 S72189 nombre de types différents d'ondes qui interviennent, que--;: l'interprétation des mesures est très délicate. La Demanderesse a observé{notaamment que l'atténuation du- signal reçu peut être aussi bien produite par des fractures que par d'autres inhomogénéités dans):la formation. De plus, la résolution en profondeur que l'on pouvait obtenir par- l'étude de l'atténuation est relativement modeste. On décrira maintenant en référence.-à la figure 3 un procédé- qui permet d'améliorer considérablemènt la détection de fractures, en se fondant-au contraire sur -'1élévation locale du couplage acoustique.entre le fluide du forage et les -for- mations-du sous-sol que t'on a-pu observer au voisinage de-- fractures. Pour la mise en oeuvre du procédé,-ondescend dans le puits, de manière repérée, un outil émetteur/récepteur acoustique tel que celui illustré sur la figure 1, ou en variante un autre outil, par exemple celui connu sous le nom de Digital Sonic Tool. Il-est considéré comme souhaitable que cet outil permette d'accéder à au moins trois espacements entre émet- teur et récepteur, pour chaque-niveau- de profondeur o 1l'on va rechercher une forme d'onde, et que ces espacements- soient supérieurs -au mètre. Au cours -de la descente, on-produit répétitivement l émission acoustique de l'outil, soit pour l'un de ces transductèurs, soit pour chacun- de -ces-transducteurs, en succession, comme - défini plus haut. On enregistre à chaque fois la forme d'on- de d'écho WFn reçue par les récepteurs de l'outil, en rela- tion avec les espacements. Il y a donc en principe une for- me d'onde pour chaque espacement et pour chaque profondeur, d'o l'écriture WFn(z).. - Comme décrit en relation aux figures-:2A et 2B,: les échos résultent du couplage acoustique de-l'outil au sous-sol à travers le fluide ou boue emplissant le puits. Compte tenu du fait que les données reçies sont échantil- lonnées en profondeur, on notera maintenant les formes d'on- des par WFn(zi). L'écart entre deux profondeurs d'échantil- lonnage, zi+1 - zi = Az est de l'ordre de 15 centimètres et l'indice n qui représente les espacementsrestvariable de 1 à m avec m au moins égal, à 3. Ces informations sont en principe transmises en surface, o les formes d'ondes sont alors acquises, comme le rappelle l'étape 10 de la figure 3. Intervient ensuite une étape 12, qui consisté, pour chacune des formes d'ondes WFn(zi), à déterminer l'amplitude crête à crête maximale que présente cette forme d'onde dans sa fenêtrede compression. La fenêtre de compression est aisément repérable comme on l'a vu en référence à la figure 2Bo On peut prendre l'amplitude maximale crête à crête soit sur tous les cycles compris dans cette fenêtre de compression, soit, plus simplement, déterminer l'amplitude maximale à partir de celle de la seconde ou troisième alternance pré- sente dans-la fenêtre de compression (P1 figure 2B)-, ou en- core la somme de celles-ci, De préférence, l'information d'amplitude ainsi obtenue, qui est rapportée à un espacement et à une profondeur, est transformée en son logarithme, par exemple décimal. Il en résulte une suite d'échantillons logarithmiques An(zi) ce qui termine l'étape 12 de la figure 3. Une telle suite d'échantillons s'analyse comme un tableau de grandeurs numériques. En l'espèce, on a m suites d'échan- tillons numériques associées respectivement à des profondeurs croissantes. Pour simplifier l'exposé, on appelle ici courbe chaque suite d'échantillons définie en relation avec la profondeur. La suitedu procédé consiste alors à isoler de chaque courbe d'amplitude de compression An(zi) , les pointes d'amplitude qui ont une largeur voisine d'une valeur prédéterminée. Il semble actuellement souhaitable de définir cette largeur prédéterminée comme correspondant à environ une longueur d'onde de l'onde acoustique de compression dans le sous-sol, ce qui donne en particulier 0,3 mètre (1 pied) pour une onde acoustique de 16 kHz. Ceci fournit des courbes de pointes d'amplitude Pn(zi) associées aux différents espacements. La présence de fractu- res est alors déterminée à partir des courbes Pn(z+sn) d'une part et Pn(z-sn) d'autre part, s désignant le demi-espacement. Ceci sera maintenant décrit de manière plus précise en se référant à nouveau à la figure 3. Après l'étape 12, l'étape 14 consiste à lisser par filtrage en moyenne mobile, chacune des courbes d'amplitude An(zi), ce qui fournit autant de courbes Mn(zi).-Ce lissage s'effec- tue en tenant compte des pointes d'amplitude--que l'on veut conserver dans le signal. Typiquement, pour un échantillonnage en profondeur tous les 15 centimètres,. il est apparu convenable d'effectuer le lissage de l'étape 14 à l'aide d'un filtre de médiane opérant sur cinq points de profondeur, c'est-à-dire deux points de part et d'autre du point concerné, en plus de celui-ci. Des fil- tres de médianes sont décrits dans l'article de J.B. Bednar, 1983, "Applications of Median Filtering to Ceconvolution, Pulse Estimation, and Statistical Editing of-Seismic Data": Geophysics, 48, pp. 1598 à 1610. L'étape 16 consiste alors à déterminer la différence algébri- que, provisoirement notée Pn(zi) entre An(zi) et Mn(zi). * Ensuite, un test 18 détermine si cette différence est posi- I 5_ '= P ' 3 --0 - 25 - 7- 2189 - 15 tive ou. négative ou nulle. Si elle est positive-, on conser- ve à Pn(zi) la définition que l:tuidonne l'tape 16. Si la différence est négative ou nulle- on affecteà-Pn(zi) la- valeur 0, à l'étape 20. Ceci peut être considéré comme la soustraction stricte, c'est-à-dire arithmétique, de Mn(z) X(i)'0::0X - - i-: - à An(z.) :- - f Ainsi a.pu être déterminée la:courbe de-pointes-d'amplitu- de- Pnzi). - Comme on dispose de:préférence de-plusieurs espacements, l'étape 22 va non seulement calculer tout d'abord les courbes Pn(zi-sn) et Pn(zi+sn), mais aussi utiliser ensuite ces cour- bes pour effectuer une comparaison entre les-résultats rela- tifs aux différents espacements. De préférence-,- cette compa- raison s'effectue-sous la forme du calcul de deux-coefficients de corrélation, fonction de la profondeur, et notés respecti- vement T(zi) et R(ii). Ces coefficients sont donnés par les relations suivantes: : : - ; - = n- m- - - ET(z)m -= Prod IPn(z%+s-)] -= n_1 - = n= m- - Ep,(z)3m =Prod ÀPn(z-Sn)j:':. -: - - n= i Comme on le verra mieux plus 1oin sur un-exemple, la Demande- resse a observé-qgue les crêtes dans les courbes T(zi) et R(zi) indiquent les emplacements-de couplage acostique augmentés,-c'est-à--dire de-fractures. La résolutIon-vertica- le de ce procédé est de l'ordre de la largeur- d'une crête, -c'est-à-dire typiquement: 0,3 mètre (1 pied). On note- que, par opposition à la technique antérieure, Le procédé -de l'invention aboutit à une-résolutionen-profondeur qui est bien meilleure que l'espacement entre émetteur et récepteur. t6- Le procédé peut alors se terminer par une étape 24, qui con- siste à visualiser conjointement lese deux.grandeurs T (zi _et R(zi), de préférence avec les cóurbes de pointes-d'amplitude Pn(zi) valables pour les m espacements que-permet l'outil. - _ =:: - _; Un exemple particulier- sera maintenant discuté èn référence. aux figures 4, 4a et 4b. - La figure 4 donne, en représentation tridimensionnélle, des formes d'ondes recueillies par un outil du genre-précité.,0 L'échelle- horizonàale est la profondeur,- l'échelle verticale: est -l'amplitude de l'onde, tandis que-des formes d'ondes =- relatives à plusieurs temps différe:nts, échantillonnées -de: 600 à 800 microsecondes-, sont illustrées-- conjointement. On- voit-immédiatement la difficulté-d' interprétation de telles courbes pour l'homme de:l'art, même expérimenté.- -: Sur la figure -4a, la-partie basse illustre l'amplitude des-: quatre formes d'ondes-obtenues-avec l'outil LS$S,pour les es- pacements de 8 et 10, 10-et 12 pieds en échelle- logarithmique. La partie intermédiaire-de laifigure 4a illustre les courbes- de crêtes:n(i.. - Pn=:z:) -: --:-- La partie-haute illustre-les courbes T(zi) et R(zi) en: regard l'une de l'autre.::::: :; : : La section illustrée sur la figure -4 correspond en fait à: - la région-qui se trouve autour de la profondeur-2655-mètres: sur la figure 2A. Il-.apparalt que l'amplitude des quatre for- mes d'ondes manifeste le comportement décrit plus-haut, rà- savoir que:: -: --l'amplitude montre deux pointes positives à des positions:: o l'émetteur ou le -récepteur se trouve adjacent I la-frac:: ture (ces points -sont: marqués T et. R respectivement sur la, figure 4); - le niveau général d'amplitude se trouve diminué entre ces deux pointes. Effectivement, pour chaque espacement, on obtient deux poin- tes à 2655,3 mètres, plus ou moins le demi-espacement, et les courbes T et R montrent deux pointes importantes (3 déci- bels) à cette profondeur. Ces pointes manifestent donc une fracture horizontale isolée. Mais, dans la section comprise entre 2680 et 2687 mètres sur 1.a figure 4a, il y a une forte atténuation d'amplitude, qui s'étend surune zone plus large que l'espacement émetteur/ récepteur. Trois pointes de récepteur (deux intenses, l'une plus faible) sont observées. Ceci suggère une zone fracturée qui contient au moins trois emplacements de couplage acousti- que augmenté. La section qui va des profondeurs 2688 à 2693 suggère la pré- sence d'une fracture verticale. En effet, et c'est là un autre aspect avantageux de l'invention, le procédé décrit per- met aussi la détection de fractures verticales, que la tech- nique antérieure considérait comme pratiquement inopérantes quant à l'atténuation sur les échos reçus. Bien que les phénnomènes ne soient pas simples à modéliser, la Demanderesse pense que, très schématiquement,des pointes intenses sont en principe indicatrices de fractures horizon- tales, tandis qu'une série de petites pointes semble indiquer une fracture verticale, dès lors qu'elles ne sont pratique- ment pas accompagnées d'une atténuation des échos reçus. On se réfère maintenant -à la figure 4b, qui-vient se ratta- cher à la figure 4a du côté des profondeurs faiIles. Il;.ap- paraît, dans la section comprise entre 2570 et 2610 mètres, une zone o l'amplitude change de manière tout à fait erra- tique. Toutefois, l'analyse de crête réduit ce magma d'infor- mations en un petit nombre de pointes intenses, qui sont indi- catrices d'une collection de zones de fracture isolées. D'autres aspects de ce procédé seront discutés plus loin. Auparavant, il est fait référence à la figure 5, laquelle explicite les conversions de modes de propagation qui peuvent se produire au-niveau d'une fracture FS du sous-sol. Il s'agit là uniquement d'illustrer les effets d'une frac- ture, étant observé qu'une conversion de mode se produit chaque fois que l'onde rencontre une frontière qui sépare des zones de fort contraste acoustique. Toutefois, il est estimé que les fractures vont produire des conversions de modes particulièrement intenses, tout spécialement entre les modes relatifs à la formation et ceux relatifs au forage. Les rai- sons en sont les mêmes que celles qui président au couplage. acoustique augmenté que l'on a décrit plus haut. Dans la partie gauche de la figure 5, un front d'onde de - compression FOC rencontre la fracture FS, située entre l'émetteur T et le récepteur R. Une partie de l'énergie incidente FOC va demeurer sous la forme d'une onde de com- pression,pour arriver au récepteur R. Une autre partie va se trouver convertie en onde de cisaillement OS,ou encore en une onde OR de l'un des modes guidés par le forage, à savoir Rayleigh, pseudo-Rayleigh ou Stoneley. Ces ondes sont détectées par le récepteur R. Les ondes de cisaillement et les ondes guidées vont arriver après les ondes de compression OC qui sont demeurées telles de l'émetteur au récepteur, mais avant les modes correspondants de cisaillement ou guidé qui auront voyagé directement de- l'émetteur au récepteur. Lorsque l'outil se déplace vers le haut, le temps d'arrivée ?57 2189 1-9 de ces échos 'à conversion de mode" va augmenter -à partir-: du temps de pure compression, pour-aller jusqu'-au temps de pur mode cisaillementou guidé,- au fur et: mesure que -la fraction du trajet des ondes qui'aura.,fait l'objet d'une conversion s'accroît de 0 à 1.: - - - -:7, X .-i- i:: i Les échos qui- résultent d'ủne: conversion- d'un mode rapide à un mode lent (en termes de vitesse de -propagation)-cons- - !tituent donc les bras sud-ouest à0nord-est-des X que forme-. la configuration croisée= Les bras nord-ouest à sud-est desX-sont explîu:s à l'aide d'un processus analogue illustré'dans 1a partie-de d- roite de la figure 5.; Ici, 1 "énergie est- convertie depuis un-mode lent à un mode -rapide- au niveau de la fracture. Plus- préi- sément,un front d'ondes de-cisaillement,-noté-FOS, qui - arrive au niveau de la fracture FS, va se--trouver pour- par- tie transformé en ondes de compression-:OC, donc plus rapi- des. Une partie de l'énergie peut demeurer sous la forme-. d'une onde de cisaillement OS plus lente. On observe donc bien une conversion d'un mode -lent à un mode rapide.,.- Si maintenant on-note par v et -u-les viteszses-respectives de propagation des-mode-s rapide et lent, et que celles-ci. sont constantes à I' intérieur de l'espacement émetteur/ - récepteur, alors la-pente p (la dérivée du temps par rapport à la profondeur)- des braes de la=configuration-en X possé- dera l'amplitude- 1/u - 1/v, qui est indépendante de l'espacement. -:: : Les figures 5a et 5b -font apparaitre des manifestations typiques de configurations croisées. La Demanderesse a pu observer que leurs traits-principaux sontles suivants:i 2 57721889 z r e à l: s - la configuration s'étend sur une profondeur égale à l'es- pacement émetteur/récepteur; :, 0 i i- '- i: = 'i - i i:: = -i -:chaque configuration est constituée de p-lusieurs croix-: ou X qui sontsuperposés aux échos directs;: -;; - dans la direction-de profondeur,les bras des X recouvrent- - toujours l'espacement-émetteur/récepteur;: - - -10 - dans-la direction temps, les-bras:desmX commencent soitz- dans la, région des échos-decompression, soit dans cel-le des échos de cisa'illement, et se terminent dans -la région des échos de cisaillement ou pseudo-Rayleigh, ou -bien des:échos -de Stoneley, respectivement. - -Bien entendu, tdes configurations croisées-, comme illustrées sur les figures 5a et 5b,sont accessibles à l'observation visuelle directe. Mais cela est apparu insuffisant pour dif- férentes raisons que 1'homme de-l'art connaît bien:-.. : : - l'influence subjective liée notamment à la qualité, de la-: visualisation à-densité variable, ainsi que du choix des - niveaux de gris; i - - := - 7:=i:: i - z : - -difficulté pour l-'observateur, meme entraîné, de faire - abstraction du côté relatif des configurations croisées, qui doivent être rapportées au niveau du signal de-fond. 0Or, en _ fait, les configurations croisées se&présentent souvent lors-- que les échos-directs sont faibles. Il en résulte que les i configurations croisées sont souvent sous-estimées, et même qu'elles passent quelquefois inaperçues.: - Ainsi, un second aspect de l'invention prévoit que l'on effectue sur les signaux relatifs aux configurations croisées, un traitement-qui comprend la détermination de l'amplitude et/ de l'-énergie propre à ces configurations croisées, relative- ment à l'amplitude ou l'énergie des échos-directs. - Il est préférable, mais pas absolument nécessaire, que les portions des formes d'ondes WFn qui sont prises en compte soient celles qui sont dues aux ondes de compression, et se terminent au premier écho du type onde de cisaillement ou onde pseudo-Rayleigh. Il est maintenant fait référence à la figure 6. La première étape du procédé relatif aux configurations croi- sées,étape notée 30, consiste à définir une fenêtre de somma- tion qui recouvre l'espacement émetteur/récepteur, et se trouve comme on vient de le dire définie entre les échos de compression, et la montée du premier écho de type pseudo- Rayleigh ou cisaillement (voir figure 2B). En pratique, pour chaque profondeur, on inclut dans la fenê- tre de sommation les parties de formes d'ondes qui se situent de chaque côté des limites, jusqu'à concurrence de la moitié de l'espacement émetteur/récepteur. Par exemple, pour un espacement de 3 mètres (10 pieds) , on utilisera 21 formes% d'ondes, échantillonnées à 15 centimètres, pour effectuer les différents calculs à une profondeur donnée, jusqu'à arri- ver à la grandeur EX définie plus loin. A noter que l'infor- mation ainsi obtenue n'est pas équivalente à l'amplitude An(z) mentionnée plus haut. Ensuite, l'étape 32 consiste à normaliser l'amplitude de chaque trace sur son maximum, qui est arbitrairement fixé à 1. Ceci assure que l'on donnera sensiblement le même poids à toutes les traces dans la sommation, et en même temps que toutes les amplitudes dans une trace donnée seront rappor- tées à l'écho direct le plus intense. - L'étape 34 consiste alors à appliquer un filtrage en profon- deur pour éliminer les échos directs. Ce filtrage est de préférence un filtrage de médiane, qui est moins critique que le précédent,et peut porter sur cinq ou sept points en profondeur. On passe ensuite à l'étape 36, qui consiste à déterminer les interceptions temporelles T et pente p d'une configuration croisée qui se trouverait dans la fenêtre. A partir de là, on peut construire des sommes de la forme - S(T,P) =oi Vil t+p-zi)o2 o Vi(t) est l'amplitude de la forme d'onde FWn en fonction du temps à la profondeur Zi, Tr est l'interception (ou inter- section) de la ligne de sommation avec la trace ou forme d'onde. considérée dans le temps, tandis que p est la pente de cette même Iigne de sommation. L'homme de l'art comprendra que cette sommation d'amplitudes correspond aussi à une sommation d'énergies. En pratique, on utilise un jeu limité de lignes de sommations qui pré- sentent les pentes et intersections temporelles appropriées, tout en s'étendant sur la fenêtre de sommation définie à l'étape 20. Chaque somme correspond donc à un bras qui est soit ascen- dant, soit descendant, dans l'un des X que comporte la con- figuration çroisée, s'il en existe une. L'étape 38 va alors consister à faire la somme des énergies E+ relatives aux bras montants, c'est-à-dire ceux dont la pente p est positive, ainsi que la somme E-, relative aux bras descendants, dont la pente est négative. Les énergies sont dans chaque cas rapportées au nombre d'échantillons considérés, qui sont notés N+ et N- respectivement. Chacune de ces grandeurs N est définie comme le produit du nombre de niveaux situés dans la fenêtre de sommation, par le nombre de lignes de sommation utilisées. 23. - '-;> On peut alors déterminer l'énergie totale EX:contenue-du fait de la conversion de mode dans-la configuration croisée, -en faisant la moyenne géométrique-des -grandeurs E+. et E-. -5 En présence d'une configuration croisée*: l'énergie EX ne sera pas nulle. Son intensité représente-en;fait-l aproportion -de l'énergie acoustique totale qui aura fait l'objet d'une con- version -de mode, d'o d'une part lapossibilité de détecter une fracture, d'autre part celle- d'"en apprécier éventuelle- ment l'importance.-- -: : = - -t-:, 4- --.e Bien entendu, une étape 40:va tester s"il 'existe une autre fenêtre de sommation- possible, auquel cas les étapes 32 -à -38 sont répétées avec cette ou ce-s autres-sfenëtres-de sommation. Au terme du procédé, une étape 42 --va permettre.!une visuali- sation de l'énergie EX déterminée- à chaque fois,en fonc.tion _ - _ du paramètre restant, à savoir la profondeur., La figure 7 fait apparaître l'-énergie-S sommée en--oblique,. pour une position de fenêtre qui e t à'-peu près centrée sur la configuration -illustrée à la profondeur 2655-mètres, pour - un espacement de 3,60 mètres (12 pieds)-. _-Cette figure-couvre une gamme étendue de valeurs de la pente et-de-l'interception temporelle.::' ] --:- : Il apparaît en fait que S-est concent-ré- sur deux bandes= étroites, qui présententl1-es-/caractéristiques suivantes-: - pente de plus ou moins 5-0 microsecondes par;pied-, qui est la pente de la ligne -joignant l'écho de compression au pre- m-ier écho du type pseudo-Rayleigh/cisaillemen-t:- - - - - valeur de l'intersection temporell e qui se situe -pratique- ment entre.900 et:1-000 microsecondes, ce qui correspond- au temps d'arrivée des échos de compression.--: -- - -. 242 -: :- --- Des mesures -;ont été faites avec-une fenêtre déplacée de:;- centimètres (6 pouces) vers- le -haut. Le résultat (non-: représenté) a été que= les-lobes situés-s-plus ou moins-Q50 -: microsecondes par pied sont décalés horizontalement vers la = gauche et la droite respectivement de 25 microsecondes, - Lorsque la fenêtre est-déplacée vers -le bas de 15 centimè- ' tres, ils sont à nouveau décalés hor0izontalement-de -25 micro- secondes, mais dans-la direction opposée. -_ -- Ces dernières observations:montrent bien quë-les lobesobser-: vés dans la grandeur S de la figure -7 corre-spondent effecti-; - vement aux bras d'un X dans la confïiguratien croisée. - - 0 z -f - =:;-é a nt- en pratiq X0 f- 0; -- es00-- - - Cecirindique égalementr en pratique, dans quelle mesure la -grandeur EX s'écartera de- sa valeur maximum lorsque la-fenê-: tre -de sommation n-' est-pas en alignement-exact avec: la:::- configuration croisée. - --- - Les figures 8a et 8b qui complètent- la précédente aux pro'- fondeurs plus faibles, montrent-de.manière combinée les- ré-- sultats des procédés qui'forment:les deux aspects de - 1 finvèntion. Tout à fait en bas sont iIlustrées- es-gr-andeurs T (z) et - R(z) obtenues -comme p-récédemment:décrit, -; -- :: Immédiatement- au-dessus:est- illustrée- -i évolution deè::la- -: - grandeur d'énergie-de- configuration croisée EX obtenue:- - pour l'espacement de 3 mètres -(10 pieds).i: -- Au-dessus sont illustrées d'autres mesures sur lesquelles;: 3 on reviendra plus loin. - - -: -: : : : Les figures 8a et 8b:correspondent par ailleurs, respecti- [: vement, -aux figures S5a et 5b. En ce qui concerne -les:confi- gurations-croisées, il y- a un -bon accord-entre -l-es figures - et 8, là o les -configurations croisées étaient d"ores e-t'- déjà clairement visibles dans la représentation à densité variable. Aux endroits o les échos directs sont-faibles,- et o les configurations croisées-éventuelles sont diffici- les à discerner, on observe malgré tout une réponse très nette au niveau de la grandeur EX. Cela constitue un avanta- ge important du procédé selon le second aspect de l'invention. Par ailleurs, les résultats des procédés qui forment les deux aspects de'l'invention sont en bon accord l'un avec l'autre, puisque tout pic dans la grandeur EX coincide d'une manière générale avec un pic individuel ou une rafale de-piés dans les grandeurs T et R. On observe que la grandeur EX,toutefois,n'a pas.une résolu- tion verticale (en profondeur) aussi bonne que les pics des grandeurs T et S. La résolution de la grandeur EX est alors estimée de l'ordre de 60 centimètres (2 pieds) pouri{ininter- valle d'échantillonnage de 15 centimètres (6 pouce). Toute- fois, on peut obtenir une meilleure résolution en prenant un intervalle d'échantillonnage en profondeur plus fin. Du bas vers le haut, les figures 8a et-8b comportent égale-- ment quatre autres pistes d'enregistrement, qui représentent respectivement: - pour la piste T3 rayonnement gamma naturel - total corrigé (c'est-à-dire uranium supprimé) - pour la piste T4 résistivité de la formation en mode respectivement profond (pointillés), superficiel (tiretés) et ponctuel (continus). - pour la piste T5 - - densité RHOB de la formation, par outil à rayonnement gamma, - porosité NPHI de la formation,par outil à rayonnement gamma. - pour la piste T6 - mesure du diamètre du puits, par outil à rayonnement gamma (pointillés longs). - densité relative DRHO mesurée à 2 espacements différents (pointillés courts).: L'homme de l'art comprendra que les procédés selon la présente invention sont tout à fait complémentaires, et qu'ils apportent des informations que les moyens de la tech- nique antérieure ne savaient pas fournir, mais aussi qu'ils peuvent néanmoins; comme cela est courant dans l'étude du sous-sol, se compléter avec-eux. - -- f- - 02572 189 - 27 - Revendications. -:!! i::::: 1. Procédé pour -la- détection de fràc.ture dans -le sous-sol à partir d'un puitscompxenant les étapes suivantés a)- descendre dan-s le puits, de fmanière repérée,.un-outil émetteur/récepteur acoustique possédantdes-transducteurs séparés (Tl, T2,R1, R2) pour 1'émission et-la réception avec au moins un espacement!entre émetteur-et récepteur; b) au cours de-la descente, produire répétitivement (DE) -- :l' émission--acoustique de l'outil, -et en-registrer- à= chaque -foies la forme d'onde d-'chos i-Fn,:reçue- parl'outil, en rela- tion-avec le ou-l-,esesepacements,. ces-!échos r-ésuItant du- cou- plage acoustique.de- l'-outil -au sous-sol à travers le flluide -emplissant le puits; --- X -0f z= f:Lt X-f - Ef V:, V X S: -f-- te 0f e c caractérisé e-n ce que la détection des fract-ures est-eúec- tuée par traitement de portions-ds:formes d'ondes-WEn qui sont représentatives de-l'un au moins-des deux effets locaûx corrélatifs des fractures que-sont- les -pointes d-amplitude- dues à une élévation- locale du coup-lage acoustique,-et les configurati on-scroises dues àune- conversion de mode. - 2. Procédé selon la revendication 1-, caractérisé -en Ce qu- iI comporte en- outre les:tapes- suivantes: i -: _ -:; =-i:.i!- *:: - _: c) extraire-(12) de chaque forme d'onde WFn.,.dans une- fenêtre choisie de cel-le-ci, contenant-Ies échos en- ondes de comm- 30.pression, une grandeur d'-amplitude-- de compression A (z):- - n ce qui fournit au moins une courbe d'-amplitude de compression en fonction de la profondeurfAfl(z), a-ssociée à un- demi- espacement s n entre émetteur -et rcepteur;- - -; d) isoler de la courbe d'amplitude -de compression An(z)les pointes-d'amplitude de- largeur voisine d'une valeur prédé- terminée (14-20)- et ce qui fournit au moins-une courbe de pointes d'amplitude Pn(z); et; :: /: e) déterminer la présence-de fractures (22) _-àpartir des--: courbes Pn (Z+s) d'une-part- et Pn(z-s)- d'autre part. - 3. Procédé selon la revendication 2,-caractérisé en-ce que - rl'étape. a) consiste à descendre -dans le' puits un;outil apte à opérer conjointement selon un nombre m d'espacements, sn désignant le demi-espacement- (avec n: mi)-; -r:1._ MI - l'étape b) comprend i 'enregistrement sélectif-des formes-:- d'ondes WFn reçues à chaque émission et- pour chaque espace-: n -.ment; = - -:-::=: - les étapes c) et d) sont exécutées -sélectivement pour cha- que espacement, fournissant un -jeu de courbes d'amplitudes - de compression An(zj,puis un jeu-de courbes de pointes d'am- plitude Pn(z), respectivement associées chacun des demi- espacements.sn-S et: - - - -= - : - l'étape e) comprend (22):;;'; :--, - a f =:E--, X E - f S= -C ': - el) la détermination de deux courbes-de-coefficients:de -:. corrélation croisée T(-z) et R(Z)-, enfonction-de la -profon' - deur, selon les relations: - n =m [-T(z)]t =Prod 1 Pn-(z-is] n= - = n- n =; :::- X S:S 0 -: -2 _ -: ' % :n =m [: --r:.: _: ' - -: _: - [R(z)]-m = Prod 'Pn(z-sn) - n = 1, -: - [: ':- : nOi : -f -: e2) la détermination de la présence de fractures à partir - des cretes de la courbe T(z) d'une part,de la courbe R({z) d'autre part. 4. Procédé selon la revendication 3, caractérisé en ce que l'étape d) comprend: dl) la détermination (-14) de.,courbes lissées Mn(z) par fil- trage en moyenne mobile de chacune des courbes d'amplitude An(z); et d2) La détermination de la courbe de pointes Pn(z) par sous- traction stricte de Mn(z) dans An(z), avec Pn(z) = An(z) - Mn(z) si An(z) > Mn(z) et Pn(z).= 0 dans le cas contraire (18, 20). 5. Procédé selon la revendication 4, caractérisé en ce qu'à l'étape dl) le filtrage en moyenne mobile (14) est réalisé à l'aide d'un filtre de médiane opérant sur cinq points, la répétition des mesures acoustiques se faisant tous les X cm environ. 6. Procédé selon l'une des revendications 2 à 5, caractérisé en ce que la valeur prédéterminée pour la- largeur des pointes d'amplitude correspond à environ une longueur d'onde de l'onde acoustique de compression dans le sous-sol,en parti- culier 0,3 mètre pour une onde acoustique de 16 kHz. 7. Procédé selon l'une des revendications 2 à 6, caractérisé en ce que m est au moins égal à 3. 8. Procédé selon l'une des revendications 2 à 7, caractérisé en ce qu'à l'étape c) la grandeur An est définie par le maximum de l'amplitude crête à crête de la forme d'onde WFn sur une partie au moins de sa fenêtre contenant'les échos en ondes de compression (12). 9. Procédé;selon l'une des revendications 2 à 8,caractérisé en ce que la grandeur d'amplitude de compression est de nature logarithmique (12). 10. Procédé selon l'une des revendications 2 à 9,prise en combinaison avec la revendication 3, caractérisé en ce que l'étape f) comprend l'établissement d'une représentation conjointe (24) des courbes T(z), R(z) et de préférence Pn(z), en fonction de la profondeur. 1l. Procédé selon la revendication 1, caractérisé en ce que le traitement relatif aux configurations croisées comprend la détermination de l'amplitude ou de l'énergie propre à celles-ci, relativement à l'amplitude ou l'énergie des échos directs (30-34). 12. Procédé selon la revendication 11, caractérisé en ce que les portions de formes d'ondes prises en compte sont celles qui sont dues aux ondes de compression, et se termi- nent au premier écho du type onde de cisaillement ou ondes pseudo-Rayleigh (30). 13. Procédé selon l'une des revendications 11 et 12, carac- térisé en ce que le traitement relatif aux configurations croisées comprend la sommation d'amplitudes, tirées des for- mes d'ondes, et relatives à des positions attendues pour les deux bras d'une configuration croisée, à l'intérieur d'une fenêtre qui recouvre l'espacement émetteur/récepteur (30-36), afin de déterminer l'énergie cumulée EX des échos sur les deux bras de la configuration croisée, la présence d'une fracture étant déterminée d'après l'écart de cette énergie cumulée par rapport à zéro. 14. Procédé selon la revendication 13, caractérisé en ce que l'énergie cumulée EX est déterminée comme la moyenne géomé- trique des énergies E+,E- relevées sur les deux bras (38). 2 2572189 - - '- = /i:: -c::C _- - 31 1 15. Procd-dé selon l'une des revendications. 1-3 et 14, carac- térisé en ce que, s('r,p) notant le carré-- de la somme des am- plitudes-selon une ligne qui intersecte lie temps--Tr avec-une pente p. on détermine les énergies sur les bras ascendants E+ et-sur les bras desceidants:E- de la configuration croi- sée par les relations (36,-38) - E+= Somme IS(',p)] /N+ - - - L p. >-O- f= 0 X C 0 X f D- - -- o N est le produit du nombre de niveaux ou échantillons- dans -la fenêtre diaccumulation--par le nombre de-lignes d'ac- cumulation 'ascendantes' -:considérées et.-- :. - - f s =: = E- So: -: _ -: o N est le produit du- nombre de niveaux-ou, chantillons-1. dans la fenêtre d'accumulation par-:le nombre de lignes d'accumulation '. "descendantes' considérées. X 16. Procédé selon la revendication 15, caractérisé en ce-que la gamme des valeurs de -T utilisées coincide avec la fenêtre des échos de compression, tandis que la gamme des-valeurs de p coincide avec la pente:des lignes joignant les échos de- compression. -- -: = -- 17. Procédé selon l'une- des revendications 1-5- -et 16, carac- térisé en ce que le traitement relatif aux confi-gurations { croisées comprend -les étapes su-ivantes -: - - _: -.: =:: _ : c') définir une fenêtre-de sommation (30-)'couvrant l'espace- ment émetteur/récepteur- et s'étendant entre -les échos de compression et la montée du premier-écho pseudo-Rayleigh ou de cisaillement; -:: -: - d') normaliser l'amplitude/de chaque forme d'onde par rapport à l'écho direct le plus intense qu'elle contient (32); -2572189 - e') éliminer.leséchos directs par filtrage de médiane dans la direction profondeur (34) - - : X X D-t =---: =; - -.:.: f')- déterminer les amplitudes sommées comme défini plus haut- pour un jeu choisi de lignes de sommation possédant des pentes et intersections temporelles correspondant à une - configuration croisée dans ladite fenêtre (36); - g') déterminer les grandeurs E+, E- et EX comme-défini plus;- haut (38);e -t - h') répéter les étapes c') à g') pour les differentes -; positions admissibles de la fenêtre de sommation.:- 18. Procédé selon la -revendication 17, caractérisé en ce-; qu' il comprena en outre une visualisation de l' intensité des; configurations-croisées (42). - - - - 19. Procédé selon l'une des revendications précédentes, - caractérisé par une visualisation conjointe des informations tirées de l'élévation locale du couplage acoustique-et des - configurations croisées.-:- :'-