Brevet FR2602059A1 — publié le 19880129

q La présente invention concerne des procédés de diagraphie électrique de fond et des appareils pour mettre en oeuvre ces procédés, et plus particulièrement un procédé pour traiter systématiquement des infor5 mations de tension et de courant de diagraphie préliminaires de manière à aider à l'amélioration de l'interprétation de la résistivité d'une formation terrestre autour d'un trou de sondage d'o ont été recueillies les données préliminaires, comme par l'intermédiaire O10 d'une électrode unique de courant émettant de façon continue et d'une multiplicité d'électrodes potentielles d'un ensemble mobile de diagraphie. Selon un des aspects de la présente invention, le traitement des données de terrain implique l'utilisation 15 de matrices d'impédance réciproques précises de façon surprenante indexées sur des jeux d'incréments prédéterminés mais qui se recouvrent d'analyses de profondeur de manière à produire une série de réponses de synthèse qui peuvent être associées à un nombre étonnamment 20 grand d'ensembles simulés par ordinateur. A partir des réponses de..thèse, on peut faire facilement une comparaison avec les réponser.- produite, par des ensembles similaires dans des conditions connues de trou de sondage et de formation et, en conséquence, 25 on peut déterminer la résistivité vraie (Rt) de la formation ainsi que la résistivité de la zone envahie (Rxo). En accord avec un autre aspect de la présente invention, une telle comparaison implique la manipulation 30 des matrices d'impédance réciproque de manière à indiquer de façon plus précise des modifications latérales et verticales dans les résistivités de la formation et de la zone envahie de manière similaire au cas de plusieurs ensembles modernes de diagraphie focalisés mais 35 pour lesquels une base connue de données a été développée précédemment à des fins de comparaison. D'une connaissance de la distribution des tensions dans la formation terrestre que traverse un trou de sondage, résultant de l'écoulement de courant 5 imposé dans la formation, on peut déterminer la saturation en hydrocarbures de la formation. Les matrices rocheuses généralement ne conduisent pas l'électricité. Cependant, si la formation est oreuse et contient des fluides, on peut faire pénétrer du courant à travers la formation 10 et la distribution des tensions peut être mesurée le long du sondage. L'impédance de la formation est liée à sa,:capacité à empêcher l'écoulement de courant à travers cette formation et elle est mesurée en ohms. La résistivité de la formation est liée également à 15 la capacité de cette dernière d'empêcher l'écoulement du courant mais elle est mesurée non pas en ohms mais en fonction d'ohms par mètre carré par mètre ou d'ohms par mètre. En d'autres termes, la résistivité d'une formationest l'impédance (en ohms) d'un mètre par un mètre par un mètre cube de la formation lorsque le courant s'écoule entre: faces opposées du cube. Les résistivités ombent dans la gamme de 0,2 à 1000 ohms-mètres dans les formations terrestres les plus perméables connues de la Demanderesse. Comme la formation devant subir la diagraphie est traversée par un sondage contenant un fluide ayant une résistivité autre que celle de la formation adjacente, la résistivité apparente obtenue (Ra) peut différer de la résistivité réelle (Rt) de la formation. En d'au30 tres termes, la présence du sondage rempli d'un fluide ayant une résistivité Rm différente de celle de la formation, le fait que le filtrat de fluide de perçage envahit la formation à un degré limité et balaye l'eau de la formation et quelques hydrocarbures de façon 35 à avoir une résistivité Rxo différente à nouveau de celle de la formation, et le fait que les électrodes de mesure peuvent traverser des formations voisines, perturbent tous les résultats finals. Certains procédés de diagraphie électrique 5 remédient à de telles perturbations à cause de nouvelles conditions de sondage. Par exemple, des diagraphies classiques de résistivité (diagraphies non focalisées) fournies par des outils de surveillance électrique classiques (ES) fournissent des bonnes estimations 10 de résistivité vraie seulement dans des lits homogènes épais ayant des porosités supérieures à 15%. Pour des conditions de lit plus mince, de tels outils peuvent fournir des résultats fiables si l'invasion du filtrat est peu profonde, si la résistivité vraie est faible 15 jusqu'à modérée et si la résistivité de la zone rincée est égale ou inférieure à la vraie résistivité à mesurer. Des diagraphies additionnelles plus avancées ont été développées pour se concentrer sur l'augmentation despropriétés de focalisation des outils électriques afin de pallier les perturbations m ntionnées ci-dessus. Par exemple, des familles d'outi le résistivité ont été développées dans ce d;ier quart de siècle; ils utilisent des courants de focalisation pour commander les trajets suivis par le courant de mesure: parmi 25 de tels outils on peut citer les outils de diagraphie de focalisation comportant l'outil focalisé sphériquement. De tels outils utilisent des courants fournis à des électrodes spéciales de l'outil et les valeurs lues sont moins affectées par les conditions de sondage 30 et la présence de lits adjacents. Cependant, à un degré sensible, les deux types de diagraphie n'ont pas été assez souples dans les conditions variables de sondage rencontrées dans les champs de production actuels, soit sur terre soit en mer. Par exemple, les diagraphies ES classiques ont une structure trop large pour fournir un procédé permettant à l'utilisateur de déterminer la réponse de focalisation d'outils électriques indépendants du système des électrodes. Au contraire, des diagraphies électriques focalisées sont formulées trops strictement pour fournir de tels résultats indépendants. En d'autres termes, des mesures insuffisantes en résultent, fournissant des résultats de caractéristiques de focalisation allant au-delà de la configuration d'origine. En plus, des 10 facteurs de calibrage pour des outils focalisés profondément et peu profondément apparaissent être choisis de manière que leurs réponses soient égales à la résistivité vraie de la formation dans des formations non envahies ayant des contrastes de résistivité formation/ 15 boue dans une gamme de 10/1 à 100/1 normalisée sur un trou de sondage de 20 cm. De là, pour permettre à l'utilisateur d'avoir la faculté de tester différentes réponses de focalisation, indépendamment du système des électrodes, on devait dévelooner un procédé de diagraphie entièrement différent. Une telle proposition est itenue dans le brevet des Etats-Unis d'Amérique N 3 076 '38 intitulé "Electrical Logging" de R.B. Stelzer dans lequel un outil à ensembles multiples d'électrodes est utilisé 25 pour fournir des mesures de tension et de courant qui peuvent être disposées dans un format de matrice dans un ordinateur numérique en fonction de la profondeur le long du trou de sondage. Dans ce brevet, la genèse du format de 30 la matrice est décrite sous la forme d'un ensemble de 2 x 3 divisé en six sous-matrices dont l'une est un ensemble carré dont les entrées de façon surprenante apparaissent être indépendantes de tout système ultérieur d'électrodes à synthétiser. La sous-matrice carrée 35 mentionnée ci-dessus a des rangées qui peuvent être remplies de valeurs de données de champ brutes par exemple pour identifier la tension à une position commune de profondeur et des colonnes de valeurs qui peuvent identifier la réponse en tension au niveau d'une série d'électrodes de tension (y compris l'électrode d'émission de courant) en fonction d'un emplacement commun d'électrodes de courant. On pense que cette proposition est la première à reconnaître que des données de diagraphie d'électro10 des (c'est-à-dire courant, résistance et tension) pouvaient être combinées dans un tel format de matrice. Dans des opérations sur le terrain, une électrode de courant montée au fond est excitée de façon continue à mesure que la sonde est déplacée dans le 15 trou de sondage. Des mesures absolues de tension à chacune de la série d'électrodes de détection de tête (y compris l'électrode d'émission de courant) sont détectées et enregistrées par rapport à une électrode de référence éloignée de tension de tête. Une élect ode 20 de courant de retour est également montée sur la 1b ide de l'outil en un emplacement approprié par ra, - aux autres électrodes et le courant de l'électrode émettrice est également mesuré et enregistré. En divisant les tensions absolues mesurées par le courant correspon25 dant mesuré en accord avec la loi d'Ohm en format de matrice, on peut établir une matrice de résistance R entre des valeurs arbitraires de synthèse de tension et de courant. (Dorénavant, les quantités de matrice seront soulignées). En principe, une telle matrice 30 de résistance est appropriée pour synthétiser sensiblement les réponses d'outils de diagraphie électrique classiques par la manipulation des éléments de la matrice. De telles opérations impliquent spécifiquement une sousmatrice telle qu'expliquée avec plus de détail ci35 dessous et ceci est très important lors de la réalisation du schéma en accord avec la proposition à cause de sa propriété de base de permettre aux courants de synthèse d'être déterminés de façon unique à partir des tensions correspondantes, ou vice versa. Il est également important pour la mise en oeuvre pratique de la proposition de reconnaître qu'il sera généralement nécessaire de résoudre des systèmes d'équations impliquant la sous-matrice mentionnée ci-dessus ou ce qui est équivalent de calculer de 10 façon précise l'inverse mathématique de la sous-matrice pour simuler les réponses d'outils modernes de focalisation,:En d'autres termes, la solution de la réciproque d'une telle sous-matrice sera généralement nécessaire pour la synthèse d'outils modernes focalisés et spéciale15 ment pour la synthèse de nouvelles combinaisons d'électrodes inconnues jusqu'ici nécessitant des relations de tension et de courant arbitraires. Ainsi, la proposition, ci-dessus est applicable seulement aux situations pour lesquelles il est possible de produire l'inverse 20 de la sous-matrice avec une précision suffisante. iais l'expérience a indiqué que dans beaucoup d'applicat-ions sur le terrain, de tels résultats ne sont pas possibles. Ce problème est lié aux contraintes numériques imposées par le procédé de mesure qui provoque finalement une 25 précision finie limitée des mesures de tension et il est apparu régulièrement dans les situations sur le terrain pour lesquelles le contraste de la résistivité de la formation par rapport à celle de la boue est plus grand que 100 à 1 (à savoir, dans des situations 30 o on utilise des fluides de forage salés; o la formation non envahie a une porosité faible; et o il y a une saturation d'hydrocarbures de modérée à élevée). On pense que lorsque la proposition ne fournit pas des résultats précis, ceci provient du fait que dans 35 de tels constrates élevés, le potentiel tend à changer très lentement d'électrode à électrode. Ainsi, il a été impossible de préserver la précision nécessaire pour résoudre de façon précise la variation graduelle appliquée. En conséquence, lorsqu'on manipule subséquem5 ment des valeurs de potentiel de matrice, par exemple dans le cas o des calculs de virgules flottantes spécifient des différences dans le potentiel entre des électrodes voisines, le procédé de cette proposition s'écroule. Plus récemment, une seconde 10 proposition a été avancée dans le brevet des EtatsUnis d'Amérique N 4 087 741 intitulé "Downhole Geoelectric Remote Sensing Method", de I.R. Mufti, dans lequel on décrit un outil à ensembles d'électrodes multiples pour la détection d'anomalies résistives latérales 15 éloignées du trou de sondage. De façon typique, de telles anomalies sont des d6mes de sel. Ce système utilise le principe de superposition pour réaliser la synthèse de divers outils à quatre électrodes à la manière d'outils de diagraphie électrique espacés o20 très longs (ULSEL). Voir R.J. RUNGE ET COLLABORATEURS, "Ultra-long Spaced Electric Log (ULSEL)", THE LOG ANALYST, Vol. 10, N 5, septembre-octobre 1969. De façon particulière, dans cette composition, un ensemble d'électrodes de courant monté au centre 25 (à savoir une électrode de courant ayant des électrodes de détection de tension disposées symétriquement audessus et en dessous de l'électrode de courant) est disposé sur une bride d'une longueur très importante. L'électrode de courant est excitée de façon continue 30 à une fréquence basse à mesure que la bride est déplacée à travers le trou de sondage. Des différences de tension entre des électrodes de détection adjacentes au-dessus et au-dessous de l'électrode de courant sont mesurées et enregistrées. Le but exclusif de l'outil est de 35 synthétiser divers outils à grande distance et à quatre électrodes pour la détection d'anomalies latérales. Comme les électrodes de détection de la tension ne sont pas espacées uniformément, et comme les quantités liées à la résistance du point de commande, (c'est5 à-dire l'impédance du point de commande au niveau des électrodes d'émission de courant) ne sont pas mesurées, la proposition n'a pas pour résultat le type de formulation de matrice que propose soit la première proposition soit la présente invention. En d'autres termes, tandis 10 que la seconde proposition permettra des calculs de potentiels au niveau des électrodes données en présence de c.ertains courants arbitraires au niveau d'autres électrodes, il ne permettra pas les calculs inverses, c'est-à-dire le calcul du courant à une position donnée 15 d'électrodes pour des potentiels donnés à d'autres positions d'électrodes par l'intermédiaire d'une matrice d'impédance mesurée. On ne Deut donc pas l'utiliser non plus en principe ou en pratique pour synthétiser d'autres types d'outils de diagraphie ayant généralement de 20 l'intérêt. Dans une demande de brevet parallèle, la Demanderesse décrit en détail comment utiliser la matrice d'impédance améliorée telle que fournie par une électrode de courant unique centrale 25 activée en continu en conjonction avec une multiplicité d'électrodes potentielles d'un ensemble mobile. Ces étapes ont compris une sous-étape permettant de déterminer les valeurs de résistivité synthétique qui ont été formulées en utilisant les tensions synthétiques 30 produites à partir des valeurs de tension et de courant mesurées préliminaires. Comme les valeurs de potentiel synthétique sont basées en partie sur des diffences potentielles entre des paires adjacentes d'ensembles d'électrodes, elles sont étonnamment précises pour tout type de conditions de trous de sondage y compris celles o des contrastes de résistivité de la formation et de la boue de forage dépassent 100:1. Cependant, il y a une autre nécessité à corréler de façon précise les différences dans les réponses des ensembles autres que par la production d'une série de valeurs de résistivité synthétique pour prendre en compte d'une manière systématique la variation dans les réponses de tels ensembles en fonction par 10 exemple de l'invasion du filtrat de la boue de forage et des contrastes de résistivité différente dans le trou de forage et autour de celui-ci qui traverse la formation contrôlée. En accord avec la présente invention, on 15 développe un procédé pour synthétiser les caractéristiques de réponse vraies d'une combinaison de différents outils de diagraphie électrique centrés dans le trou dans des conditions variées difficiles de trous de sondage, le procédé consistant à (1) déterminer les 20 valeurs d'impédance d'une formation de terrain traversée par un trou de sondage rempli d'une boue de forage d'une résistivité (Rm), et (2) manipuler sélectivement les valeurs d'impédance en tant qu'entrées d'impédance d'une série de matrices réciproques de manière à synthé25 tiser les fonctionnements de différents outils centrés dans le trou sur un incrément de profondeur associé avec une précision surprenante. Chacune des matrices d'impédance réciproque est associée à un ensemble de matrices réciproques indexé sur un parmi une série 30 d'incréments finis de balayage de profondeur qui se recouvrent, ces incréments appartenant à la formation et étant mesurés le long du trou de sondage. Chaque incrément de balayage dépend de la longueur d'ensemble L de l'ensemble des électrodes de manière à définir 35 des marqueurs de profondeur profonds et peu profonds, et sont indexés centralement sur la profondeur du trou de sondage d'un ensemble d'électrodes centrales du système au moment du recueil des données. De façon détaillée, pour fournir une indica5 tion vraie de la résistivité de la formation (Rt) même lorsque cette dernière est espacée du trou de sondage par une zone envahie de.résistivité (Rxo) d'une longueur latérale inconnue en raison de l'invasion du filtrat de la boue de forage, les étapes suivantes sont effec10 tuées en séquence: (i) un assemblage de M ensembles d'électrodes centrées sur le trou ayant un espacement "a" incrémentiel égal entre les électrodes est tout d'abord calibré pour obtenir des jeux de facteurs d'étalonnage normali15 sés sur des configurations connues d'amorçage de tension et une réponse de courant associée dans une zone de réponse de résistivité connue, les ensembles d'électrodes ayant un indice connu de numérotation d'ordre interne, lesdits jeux de facteurs d 'étalonnage pouvant être chacun adressés en fonction de conditions de trous de forage y compris les conditions difficiles de trous de forage ainsi que- par un type d'assemblage synthétique focalisé par ordinateur; (ii) ensuite, un assemblage sur le terrain 25 d'ensembles d'électrodes centrées sur le trou d'espacement incrémentiel égal entre électrodes "a" est placé dans le trou de forage, les ensembles d'électrodes ayant un indice connu de. numérotation d'ordre interne allant de 1 à 2M-1, un ensemble d'électrodes centrales comprenant une électrode de courant et une électrode de potentiel ayant un indice de numérotation M et des électrodes d'extrémité ayant des indices 1 et 2M-1, tous les ensembles d'électrodes sauf l'ensemble central comprenant seulement une électrode de mesure 35 de potentiel, et la profondeur absolue d'au moins un 1d ensemble d'électrodes étant connue de façon continue par rapport à un niveau de données de profondeur prédéterminé mesuré depuis la surface; (iii) ensuite, l'assemblage est déplacé de façon continue le long du trou de forage et le courant de l'électrode de courant de l'ensemble d'électrodes centrales d'une valeur connue est injecté dans la boue de forage de résistivité (Rm) et la zone de résistivité Rxo envahie et de là dans la formation ayant une résisti10 vité Rt; (iv) pendant l'injection du courant, le potentiel absolu à chaque électrode de potentiel ainsi que les différences de potentiel entre des paires adjacentes d'électrodes de potentiel sont mesurés d'une manière 15 rapide, les mesures de potentiel pouvant être indexées sur des stations communes de diagraphie d'une série de stations de diagraphie espacées, également situées le long du trou de forage et ayant une distance incrémentielle d'espacement de "a", "a" étant la distance entre 20 les ensembles d'électrodes; (v) ensuite des valeurs d'impédance sont calculées à partir des valeurs de potentiel absolues et de différence de potentiel absolu mesurées et leur courant d'injection associé, chaque valeur étant indexée 25 sur les nombres d'indexation interne connus d'électrodes de courant actif et de potentiel utilisés dans les mesures; (vi) ensuite, les valeurs d'impédance sont réindexées pour former des entrées d'impédance d'une 30 série d'ensembles de matrices modifiées se recouvrant AZ, chaque ensemble AZ étant associé à un segment prédéterminé de ladite formation égale en dimension verticale à M stations de diagraphie, et comprenant M x M entrées d'impédance, M étant le plus grand nombre de l'indice 35 de numérotation de l'ensemble des électrodes centrales et le rapport entre le nombre des entrées d'impédance de différence et des entrées absolues étant d'environ (M-1) 1:1; (vii) ensuite, chaque ensemble de matrices 5 modifiées AZ est inversé pour former son ensemble de matrices réciproques AZ 1 en accord avec les techniques classiques d'inversion des matrices; (viii) des paramètres de réponse focalisés par ordinateur sont ensuite produits en utilisant l'en10 semble de matrices réciproques AZ 1 de l'étape (vii) en conjonction avec les mêmes configurations de réponse de courant et d'amorçage de tension de l'étape (i) (ix) finalement, des jeux de facteurs d'étalonnage de l'étape (i) sont recherchés jusqu'à 15 ce que le produit d'un jeu particulier de facteurs détalonnage et des paramètres de réponse de l'étape (viii) pour tous les types d'assemblage d'outils synthétiques soit essentiellement une constante, de sorte que la condition difficile de trou de forage en est ?0 déduite même en présence de contrastes élevés entre la résistivité vraie et la résistivité de la boue et nonobstant le fait que des jeux synthétiques de configurationsde potentiel ont été utilisés comme initiateurs des paramètres de réponses focalisées par ordinateur 25 et produits subséquemment. Sur les dessins: - la figure 1 est une élévation latérale partielle d'une formation de terrain traversée par un trou de forage contenant un assemblage de diagraphie 30 électrique de la présente invention et illustrant un procédé pour traiter des estimations de la distribution des impédances de la formation du terrain autour du trou de forage y compris une forme de la présente invention dans laquelle les valeurs de courant et de tension préliminaires pour les estimations sont obtenues par l'intermédiaire d'un assemblage de diagraphie qui se déplace de façon continue le long du trou de forage lorsqu'ont lieu des mesures de courant et de potentiel 5 à une série de stations de diagraphie indexées sur des positions de profondeur désignant des profondeurs fixes en unités d'espacement d'électrodes; - la figure 2 est un détail de l'ensemble d'électrodes centrales M de la ficure 1 illustrant des 10 détails de structure des électrodes de courant et de potentiel comprenant cet ensemble. - la figure 3 est une section prise le long de la ligne 3-3 de la figure 1 illustrant comment des mesures détectées fournies par l'assemnblage peuvent être utilisées pour en déduire 15 les caractéristiques de résistivité de la formation en fonction de distances décalées même si la formation a subi l'invasion du filtrat de la boue de forage; - la figure 4 est une représentation schématiaue du procédé de la présente invention pour indiquer 20 si tématiquement la distribution des impédances de 1 mation de terrain voisine en utilisant un assemblage,odifié d'ensembles de neuf électrodes o les calractéristiques d'indexation de désaffichage final sont montrées en détail; - les figures 5 et 6 sont des représentations schématiques d'éléments de circuit de fond et de tête pour mettre en oeuvre le procédé de la présente invention; - la figure 7 est un réseau artifiel équiva30 lent simulant une formation de terrain pour illustrer la correspondance unique entre les admittances de ce réseau et la matrice d'impédance réciproque AZ- 1 de la présente invention; - la figure 8 est une courbe de facteur 35 pseudo-géométrique pour différents assemblages focalisés par ordinateur en fonction du diamètre d'invasion (Di) utilisant une configuration normalisée de réponse de courant et d'amorçage de tension et une série de conditions connues de trou de forage; - les figures 9 à 27 représentent des courbes des facteurs d'étalonnage par rapport aux contrastes de résistivité illustrant comment des jeux de tels facteurs sont uniques pour une condition choisie de trou de forage, de manière à permettre la détection 10 en accord avec le procédé de la présente invention. Dans la description qui va suivre, les mêmes,. éléments sont désignés dans le mémoire et sur les dessins par les mêmes numéros de référence. Les dessins ne sont pas à l'échelle et certaines caractéris15 tiques de l'invention peuvent être exagérées en échelle ou en fonctionnement pour obtenir plus de clarté ou de concision. En référence à la figure 1, on voit un trou de forage 8 pénétrant dans une formation souterraine 20 9. On su, ose que le trou de forage 8 contient une bo-e de f. 10. Le but de la boue de forage 10 est d'entraîner les débris hors de la tête de perçage pendant le forage du trou 8 et également d' empêcher l'affaissement du trou de forage pendant le perçage et la pro25 duction de la formation 9. Un câble conducteur 12 d'une structure classique est suspendu à l'intérieur du trou de forage 8; il s'agit de bras multiples d'acier souple entrelacés avec plusieurs conducteurs électriques. Le câble 12 est suspendu dans le trou de forage 8 par 30 l'intermédiaire de la poulie 13 attachée au derrick 14 et de là à l'unité de levage 16. Le but de cette unité de levage 16 est d'enrouler ou de dérouler le câble 12 dans le trou de forage 8. A la surface de la terre 15, des signaux se 35 trouvant sur les conducteurs électriques du câble 12 traversent un ensubie de baaues de coulissement (non illustrées) se trouvant sur l'unité de levage 16 en direction d'un circuit 17 de traitement es de commande à l'intérieur de la fourgonnette 18 par exemple par 5 l'intermédiaire du conducteur 19. Au fond, de tels signaux prennent leur origine dans les caractéristiques opérationnelles de l'assemblage de diagraphie électrique 21 de la présente invention ou sont une fonction de celles-ci. ASSEMBLAGE DE DIAGRAPHIE ELECTRIQUE 21 Au fond, le câble 12 est attaché à l'assemblage de diagraphie 21 par l'intermédiaire d'un bouchon fileté 22. Au-dessus du bouchon 22, il y a un centreur de tête 23. En dessous de l'assemblage 21, il y a un centreur de fond 24. Les centreurs 23 et 24 ont la même structure et chacun comprend une série de bras 25 sollicités par ressorts qui sont rappelés radialement vers l'extérieur pour entrer en contact avec la paroi latér le 5 du trou de forage 8. L'arrange20 ment des centre Is 23 et 24 ressemble à celui des centralisateur -.. ilisés dans -divers types d'outils de diagraphie et d'inspection et- a pour but 'e localiser l'assemblage de diagraphie 21 en coïncidence avec l'axe de symétrie A-A du trou de forage 8. L'assemblage de diagraphie 21 a une structure cylindrique ayant un support défini par un mandrin 26. Le mandrin 26 a une surface externe 26a. Au niveau de cette surface, il y a une série d'ensembles d'électrodes ayant un indice de numérotation interne E1, E2,.... 30 EM...E2M_1 (M=2N+I1 étant un nombre entier impair). Pour que les ensembles voisins soient isolés électriquement les uns des autres, on insère une matière isolante entre chaque ensemble électrique E1, E2...E2M_ et la surface 26a du mandrin 26 comme expliqué ci-dessous. 35 Horizontalement, l'assemblage 21 est divisé en trois sections: (i) une section de tête 27 a l'intérieur de laquelle se trouvent une électrode de courant de retour 28 et une électrode de potentiel de référence 29, le butde telles électrodes 28 et 29 étant de compléter la boucle de courant et de normaliser les mesures de potentiel, respectivement, comme il est habituel dans la technique de diagraphie électrique. Si on le désire, les brins métalliques du câble 12 peuvent également être utilisés comme électrodes de référence ainsi 10 que comme électrode de courant de retour, comme il est habituel dans cette technique; (ii) la section centr.ale 30 qui comporte des filets 31 à son extrémité supérieure pour entrer en contact avec la section de tête 27. L'intérieur de la section centrale 30 comporte 15 une série d'éléments de servo-commande et de mesure devant être décrits en détail ci-après, et sous la commande du circuitmètre de l'unité 17 de traitement et de commande à la surface 15; et (iii) la section de fond 33 sur laquelle repose la série des ensembles 20 d'électrodes E1, E2,...E2M, mentionnés précédemment. La figure 2 mo la construction de l'ensemble d'électrodes centrales EM avec plus de détail. Comme on le voit, un tel ensemble comprend un élément d'électrodes de courant 34 et une électrode 25 de détection de potentiel 35, les deux étant annulaires de manière à s'adapter aux deux extrémités de l'élément isolant 42. Entre les éléments 34 et 35, se trouve un élément isolant 41qui fait partie de l'élément 42 et sert à garder les électrodes 34 et 35 hors contact 30 physique l'une avec l'autre. Les éléments 34, 35 et 42 sont situés autour de la surface externe 26a du mandrin 26 concentrique par rapport à l'axe de symétrie B-B. Il est à noter que lorsque les éléments 35 et 34 sont attachés de la manière montrée à la figure 1 de manière à fonctionner comme électrodes de potentiel et de courant, respectivement, ils ont des largeurs dans le sens de la verticale qui sont faibles en comparaison de l'espacement axial "a" entre les ensembles 5 d'électrodes et en plus ils sont logés si étroitement les uns par rapport aux autres (c'est-à-dire que l'épaisseur de la lèvre 41 de la figure 2 est faible en comparaison de l'espacement axial "a" vers l'ensemble d'électroces espacé au voisinage de cet ensemble) que pour 10 la résolution axiale leur position axiale est équivalente à une seule profondeur fixe, à la figure 1. En ce qui concerne les ensembles restants d'électrodes E 1E2e... EM_1 et EM+.,E2M_, l'élément d'électrode de courant 34 n'est pas directement utilisé 15 dans la présente invention. Seule l'électrode de détection de potentiel 35 doit être présente. Ainsi, l'élément de courant 34 peut être considéré comme retiré des ensembles restants E1,E2,...,EM_1 et EM+1' E2MEn ce qui concerne l'o:dre de dénombrement 20 interne des ensembles d'électrod:s, il est à noter que l'ensemble peu profond --st l'ensa.._. E1 et l'ensemble le plus profond est E2M-1 et que la profondeur de chaque ensemble d'électrodes E1,E2....E2M_1 peut être déterminée sur la base de profondeur = dk + (p-1)a, o p est 1,2...2M-1; dk est la profondeur absolue de l'ensemble d'électrodes le moins profond de l'assemblage; et 2M-1 est l'indice de l'ensemble d'électrodes le plus profond de l'assemblage. 30 OPERATIONS En bref, selon la présente invention, le but de l'assemblage de diagraphie 21 est de fournir une distribution d'impédance de la formation 9 de telle manière que ces valeurs préservent une relation de 35 un à un entre chaque combinaison subséquente de tensions et de courants qui pourraient être utilisées -éventuellement pour synthétiser les caractéristiques de réponse à partir de tout nombre d'outils électriques différents dans une grande variété de conditions difficiles et 6 de différence dans les trous de forage et les formations de terrains. Pour atteindre ce but, l'ensemble de diagraphie 21 fonctionne de la manière montrée à la figure 1 de sorte qu'à mesure que des valeurs de courant et de potentiel sont recueillies systématiquement en fonc10 tion de la profondeur, le circuit de traitement et de commande 17 à la surface peut manipuler sélectivement de telles valeurs pour former les entrées d'une série d'ensembles de matrices, chaque ensemble de matrice étant utile pour faire la synthèse des opérations des 15 différences d'outils de diagraphie. Le but en est de fournir une indication vraie des caractéristiques de la formation 9 même si celle-ci peut avoir été envahie par du filtrat de boue de forage dans le trou de forage ou autour de celui-ci à un degré inconnu, c la manière 20 de la figure 3, ou bien si des jeux synt hétiques ou non de potentiel ou des jeux de valeur de. ant sont utilisés subséquemment comme initi ateurs u assemblage synthétique d'outils focalisés par ordinateur. En d'autres termes, comme le montre la 25 figure 3, il est à noter que le trou de forage 8 devant subir la diagraphie par l'assemblage de diagraphie 21 peut contenir de la boue de forage 10 d'une résistivité Rm, que la boue de forage 10 peut constituer un gâteau de boue 6 au niveau de la paroi latérale 30 5 du trou de forage 8; et que le filtrat de boue peut envahir la formation 9 à une distance D du centre du trou de forage 8, ce qui crée différents niveaux de résistivité dans le trou de forage 8 et autour de celuici, à savoir une zone de gâteau de boue de résistivité 35 Rmc, une zone rincée 7 de résistivité Rxo et une zone 1 9 de transition 11 de résistivité qui peut varier entre la résistivité de la zone rincée et celle de la formation 9 à savoir entre Rxo et Rt et qui dans la discussion est appelée Rapp. Pour indiquer les modifications de 5 résistivité latérale vis-à-vis de la boue de forage , la zone rincée 7, et la formation 9, il faut non seulement que les valeurs de diagraphie de courant et de potentiel telles que fournies par l'assemblage de diagraphie 21 (en fonction de la profondeur) soient systématiquement recueillies et indexées mais encore que les valeurs d'impédance calculées à partir de ces valeurs soient indexées de façon précise de sorte qu'ensuite, à savoir après que la diagraphie a été achevée, elles peuvent être manipulées sélectivement pour fournir 15 différents degrés de résolution latérale que les initiateurs de réponse pour de telles manipulations soient ou non des jeux artificiels de valeur de courant ou de potentiel. De cette manière, les modifications de résistivité latérale mentionnées ci-dessus peuvent 20 être déterminées aisément. Pendant le prélèvement des données en ac avec la igure 1, l'assemblage de diagraphie 21 se déplace dans le trou d'assemblage 8 à une variation de vitesse constante. Les courbes de courant de l'ensem25 ble central EM sont activées en continu. Des mesures ont lieu à des instants pendant lesquels les ensembles d'électrodes,E1,E2..E2M_1 se sont déplacés de façon précise d'une distance "a" égale au facteur d'espacement de l'ensemble des électrodes comme mentionné ci-dessus 30 et ceci a lieu assez vite par rapport à la vitesse de l'assemblage 21 dans le trou de forage 8 pour qu'on puisse les appeler "instantanées" par rapport à une série de stations communes de diagraphie équidistantes d,d2...dv le long du trou de forage 8. Les mesures de potentiel commencent de façon typique avec l'ensemble E1 et sont faites dans les séquences ordonnées à travers les ensembles restants E2,E3, etc.,et se terminent avec l'ensemble E2M_. A mesure que l'assemblage de diagraphie 21 roule vers le haut ou vers le bas du trou de forage, une station de diagraphie, par exemple à l'endroit o l'ensemble E1 descend de la station de diagraphie d1 à la station d2, le processus de prélèvement est répété. Il en résulte une série de valeurs de courant et de potentiel qui sont recueillies systématiquement en fonction 10 de la profondeur pour être manipulées par la suite comme indiqué ci-dessous. A noter cependant que les valeurs de courant et de potentiel à manipuler n'ont une intégrité de formation de terrain que si elles sont liées au même local fixe de prélèvement. En d'au15 tres termes, les valeurs doivent être indexées à un intervalle particulier de balayage de profondeur fixe (égal à M stations de diagraphie comme expliqué cidessous) ce qui empêche l'entrelacement de valeurs semblables de différents intervalles de balayage de 20 profondeur, d'une manière expliquée également avec plus de détail ci-après. La figure 4 illustre comment font lieu le prélèvement et l'indexage systématiques pendant les opérations. A des fins de description, on suppose que 25 le nombre des ensembles d'électrodes comprenant l'assemblage de diagraphie 21 a été grandement diminué, à savoir par exemple calculé à partir du grand assemblage de la figure 1 jusqu'à un assemblage à 9 électrodes comprenant les ensembles d'électrodes E1,E2...E9. On in30 jecte continuellement du courant au moyen de l'électrode de courant (non montré) de l'ensemble d'électrodes centrales E5. Ainsi, on suppose que l'électrode de courant au niveau de l'ensemble central 5 de l'assemblage à 9 électrodes est activée et que des potentiels absolus 35 et de différence sont mesurés au niveau des quatre ensembles au-dessus de l'électrode de courant, y compris l'ensemble des électrodes de courant (c'est-à-dire sur les ensembles ayant J 'ordre de numérotation 1 2,....,5, en commençant par l'ensemble le moins profond par rapport à la surface 15). Ensuite, on suppose que pendant la seconde moitié du cycle de prélèvement, la même électrode de courant (qui occupe la cinquième position de l'électrode, selon son nombre interne) reste activée et que des potentiels absolus et de différence sont également 10 mesurés sur les quatre ensembles les plus profonds ayant l'ordre de numération de 6,7,...,9. r En associant les valeurs mesurées de potentiel et de courant dans 5x5 matrices, on voit à l'évidence que des quantités de potentiel (à la fois absolu et 15 différentiel) recueillies à des électrodes de potentiel qui ne se trouvent pas en dessous de l'électrode de courant fourniront des entrées de matrice sur sa diagonale et au-dessus de celle-ci, tandis que celles recueillies en dessous de l'électrode de courant fourniront 20 des entrées en dessous de leur diagonale. Mais comme l'assemblage se déplace en continu, il sera nécessaire de former chaque ensemble de matrice à partir de plusieurs cycles différents de prélèvement. La figure 4 montre comment le recueil systè25 matique et l'indexage ont lieu pendant de telles opérations dans lesquelles quatre cycles séparés de prélèvement, à savoir cyles 1,2,...5, pour des positions de diagraphie A,B,C,D et E sont sont décrits en détail. A la figure 4, l'ordonnée de la courbe indique des unités 30 de profondeur et l'abscisse indique des unités de temps incrémentiel 1,2...5. L'espacement entre les ensembles E1, E2,..., E9 est égal au facteur d'espacement "a" et il en est de même de la distance entre des stations de diagraphie adjacentes d,d2,d2...d14. Bien que l'as35 semble 21 se déplace continuellement le long du trou de forage 8, chaque emplacement A,B,...,E marque un moment pendant lequel le prélèvement des valeurs de potentiel de phase et de courant a lieu. Il est à noter à ce propos que pendant le recueil des données en accord 5 avec la figure 4, l'assemblage est roulé continuellement vers le bas. Le mouvement de l'assemblage 21 a lieu à cause du déroulement du câble 12 par l'intermédiaire de l'unité de levage 16 à la surface 15. Les valeurs recueillies sont transmises en tête par l'intermédiaire 10 du câble 12 et de là depuis l'unité de levage 16 jusqu'au circuit de traitement et de commande 17. Etant donné la grande masse de données, l'indexation des valeurs de diagraphie est plutôt importante et elle dépend des positions de profondeur absolue ainsi que 15 relative de l'électrode de courant émettrice ainsi que de celle des électrodes de mesure de potentiel comprenant les ensembles d'électrodes E1, E2,...E9. Par exemple, pour des mesures prélevées lorsque l'assemblage 21 est à la position A à la figure 20 4, l'électrode de courant de l'ensemble d'électrodes est à un marqueur de profondeur dk+4a qui coincide avec la station de diagraphie (d5). Pour l'assemblage 21, chaque cycle de mesure 1,2,...,5 nécessite le prélèvement des valeurs analogiques suivantes: (1) huit valeurs 25 de différence de potentiel, (2) neuf valeurs de potentiel absolu, (3) une valeur d'intensité de courant et (4) deux paires de valeurs de commande liées à l'indication de la distorsion de phase, c'està-dire l'indication de la distorsion par l'intermédiaire d'une différence 30 de temps entre le courant à l'électrode de courant de l'ensemble E5 et le potentiel aux deux électrodes de potentiel le plus éloignées. Ces valeurs sont transmises en tête par l'intermédiaire du câble 12 et de là à la surface 15 de l'unité de levage 16 au circuit 35 de traitement et de commande 17 pour le stockage et la manipulation en accord avec le procédé de la présente invention. Afin de s'assurer que les adresses des valeurs recueillies de courant et de potentiel sont 5 complètes, les indices suivants sont enregistrés relativement aux valeurs recueillies de courant et de potentiel à savoir: (i) par des marqueurs de profondeur dk,dk+a, dk+ 12a o le facteur "a" et l'espacement incrémentiel entre des ensembles d'électrodes et dk et la profon10 deur absolue de l'ensemble d'électrodes E1 au début du prélèvement des données, à savoir lorsque l'assemblage est placé à la position A; (ii) par des stations de diagraphie à électrodes dénombrées consécutivement (d1,d2,d3,...d13) associées à toute l'opération de 15 diagraphie tout comme lorsque la position relative de chaque station présente de l'intérêt; (iii) par le numéro de la station de profondeur de balayage (Sd1, Sd2, etc) associé à la profondeur de l'ensemble d'électrodes centrales correspondant à des ensembles particu20 liers de matrices intéressants dont sept sont montrés à la figure 4, à savoir dk+6a, dk + 7a,...,dk + 12a. Ces valeurs peuvent être indexées dans plusieurs formats différents en fonction du prélèvement des données, et dont des exemples typiques sont les affichages 46,47, 25 48, 49 et 50, et ensuite elles sont réindexées en format d'ensembles de matrices comme montré à l'affichage 51. Il est également à noter que les affichages 46,47... ont une autre étiquette d'annotation, à savoir le fait que les valeurs illustrées formant chacun des 30 affichages doivent être encore indexées pour indiquer la profondeur de l'électrode de courant de l'ensemble d'électrodes E5 pendant chacun des cycles de prélèvement 1,2,3...5 qui donnent naissance aux affichages 46,47...50. Un tel système d'annotation peut égalenent être transféré 35 dans l'affichage 51 d'ensembles de matrices réindexées des valeurs d'impédance associées à ces mesures, comme expliqué ci-dessous. En d'autres termes, on suppose que des profondeurs absolues des. stations de diagraphie numéro5 tées sont connus; si bien que lorsque l'assemblage 21 est situé à la position A, les ensembles d'électrodes E1,E2... E9 sont associés à l'indice de numéroation interne 1,2,...,9 d'ordre consécutif; de là, lorsque l'électrode de courant de l'ensemble d'électrodes E5 est 10 à la profondeur dk+4A et que des mesures sont prises sur les ensembles associés d'électrodes, les valeurs absolues et différentielles de potentiel et l'intensité du courant seraient indiquées par les quantités suivantes: - Vl,5(dk-4a) V2,5(dk+4a) V9,5(dk+4a) AV2,5(dk+4a) AV3,5(dk+4a AV9,5(dk+4a) J5(dk+4a) Il est à noter en ce qui concerne les indices du potentiel absolu que le premier symbole d'indice concerne le numéro d'indice interne de l'ensemble d'électrodes auquel a lieu la mesure de potentiel et que 35 le second symbole d'indice identifie le numéro d'indice interne de i'!ecnrcE ce curant suLs:san une exc -_ -i on, tandis que le terme entre parenthèses concerne la profondeur absolue a partlr,ar exemple, de 1a surface 15 ju s' ' 1la position de l'électrode de courant. En ce qui concerne l'étiquette 5 d'adresse mentionnée en dernier, on voit que la s-ation de diagraphie de l'électrode de courant, à savoir la station de diagraphie (d.),doit également être utilisée comme élément de substitution car la profondeur absolue peut être calculée ensuite. Il est à noter que les différences de potentiel sont mesurées entre les paires d'ensembles d'électrodes, c'est-à-dire entre les ensembles d'électrodes 1 et 2; 2 et 3; 3 et 4; etc. Ces valeurs sont également indexées de la même manière que ci-dessus, c'est15 à-dire en accord avec ce qui suit: AVi (dk+4a) o i = 2,3,...,9. Il est à noter que ci-dessus, le premier indice concerne la position du plus profond des ensembles d'électrodes d'une paire et on suppose que la valeur 20 de normalisation pour former la valeur de potentiel de différence concerne l'ensemble d'électrodes ordonné descendant. En d'autres termes, la valeur AV2,5 (dk+ 4a) indique que la différence de potentiel est mesurée 25 entre les électrodes de potentiel des ensembles E1 et E2 numérotés intérieurement 1 et 2 respectivement, et que l'électrode de courant est placée à l'ensemble 5 numéroté selon l'ordre interne, tandis que la valeur AV4 5(dk+4a) indique que la différence de potentiel est mesurée entre les ensembles E3 et E4 de numéro interne d'ensembles 3 et 4, l'émetteur de courant étant associé avec l'ensemble 5 de numéro interne. Il est à noter que la représentation des valeurs mentionnées ci-dessus 35 telles qu'indiquées comprend des entrées de colonne 46a et 46b de l'affichage 46. L'intensité de courant est montrée comme étant l'entrée de la colonne 46c tandis que les mesures de temps Tl(dk+4a),T9(k+4a) associées à l'indication de la distorsion de phase éventuel5 le sont indiquées comme étant des entrées de la colonne 46d. La prochaine étape du procédé selon la présente invention est de répéter les mesures décrites ci-dessus aux positions B,C, D et E de la figure 4, 10 à savoir l'électrode de courant étant aux emplacements de profondeur dk+45a,dk+6a, dk+7a et dk+8a avec la paire de valeurs de commande dans des coordonnées appropriées de temps de manière à indiquer la présence (ou l'absence) de distorsion de phase, de la manière montrée ci-dessus. 15 Ces valeurs occupent des entrées des colonnes 47a,47b... 47d de l'affichage 47; des colonnes 48a,48b...48d de l'affichage 48, des colonnes 49a,49b...49d de l'affichage 49; et des colonnes 50a,50b...50d de l'affichage 50. 0 Le tableau I ci-dessous résume les mesures sous forme de tableau pour plus de clarté. T.ABLEAU I C46a C_6b C46c 1- -----__ Vl,5(d5) V2,5(d5) V3,5(ds) V4,5(d5) V5,5(d5) V6,1(d5) V7,1(d5) V8,1(d5) V9,1(d5) C47a V1,5(d6) V2, 5(d6) V3 (d6) V6') V4. 5i'6) V5,5(d6) V6,5(d6) V7,5(d6) V8,5(d6) V9,5(d6) C48a AV2,5(d5) AV3,5(d5) AV4,5 (d5) AV5,5(d5) AV6,5-(d5) AV7,5(d5) AV8, 5 (ds) AV9,5(d5) C47b AV2,5 (d6) AV3,5(d6) LV4, 5 (c6) AV5, 5(d6) AV6,5 (d6) AV7,5(d6) AV8, 5(d6) AV9,5 (d6) C48b J5(d5) Ti(d>) 3:( 5) = AFFICHAGE 46 - C47c J5(d6) Tl(d6) To(d6) = AFFICHAGE 47 C48c C48d Vl,5(d7) V2,5(d7) V3,5(d7) V4,5(d7) V5,5(d7) V6,5(d7) V7,5(d7) V,5(d7) V9,5(d7) AV2,5(d7) AV3,5(d7) AV4,5(d7) AV5,5(d7) AV6,5(d7) AV7, 5(d7) LVQ AVs,5(d7) AV9,5(d7) J5(d7) Tl(d7) T9g(d7) : AFFICHAGE 48 C49a Vl,5 (d8) 5 V2,5(d8) V3,5(d8) V4,5(d8) V5,5 (d8) V6,5(d8) 10 V7,5(d8) V8,5(d8) V9,5(d8) C49b AV2,5(d8) V3 2, 5 ( d 8) AV3,5(d8) AV4,5 (d8) AV5, 5(d8) AV6,5(d8) AV7,5(d8) AV8, 5(d 8) AV9,5(d8) C50b C49c C49d J5(d8) Tl(d8) T9(d8) = AFFICHAGE 49 C50a C50c C50d V1,5(d9) V2,5(d9) V3, 5 (d9) 20 V4,5(d9) V5 5 (dg9) - V6, 5(d9) V7,5(d9) V8,5(dg) 25 V9,5(d9) AV2,5(d9) AV3,5(d9) AV4, 5(d9) AV5 5 (d9) z, 5 ( dg9) G"7,5(d9) AV7, 5(d9) AV8,5 (d9) AV9,5(d9) J5(d9) Tl (d9) T9g(d9) = AFFICHAGE 50 Légende: d5 = dk+4a; d6= dk+5a; d7= dk+6a; d8 = dk+7a; et d9 = dk+8a A partir des valeurs mesurées ci-dessus de potentiel et d'intensité de courant, on peut déterminer le rapport des valeurs mesurées associées avec le même jeu de variables électriques des affichages 35..50 en utilisant les indices et les équations 46,47,48... 50 en utilisant les indices et les équations suivants à savoir pour l'affichage 46: Zi,5(dK+ 4a) = Vi,5(dk+4a)/J5 (dk+4a), i=1,2,...,9 LZi,5(dk+4a) = AVi,5(dk+4a)/J5 (dk+4a), i=2,3,...,9 Pour l'exemple de la figure 4, de telles entrées ont été mises sous forme de tableau sur le tableau II ci-dessous: TABLEAU II 21,5(d5) Z2,5(d5) Z3,5(d5) Z4,5(d5) Z5,5(d5) Zl, 3(d6) Z2,5(d6) Z3,5(d6) Z4,5(d6) Z5,5(d6) Z6, 5 (d5) Z7,5(d5) Z8,5(d5) 29,5(d5) Z6,5(d6) Z7,5(d6) Z8,5(d6) Z9,5(d6) Z6,5(d7) Z7,5(d7) Z8,5(d7) Z9,5(d7) Z6,5(d8) Z7,5(d8) Z8,5(d8) Z9,5(d8) AZ2,5(d5) LAZ3,5(d5) AZ4,5(d5) AZ5,5(d5) AZ6,5(d5) LZ2, 5(d6) AZ3,5(d6) LZ4 5(d6) AZ5,5(d6) LZ6,5(d6 AZ2,5(d7) AZ3,5(d7) AZ4,5(d7) AZ5,5(d7) AZ6,5(d7) AZ2,5(d8) AZ3,5(d8) AZ4,5(d8) AZ5,5(d8 AZ6,5(d8) &Z7, 5(d5) AZ8,5(d5) DEPUIS LZ9,5(d5) L'AFFICHAGE 46 LZ7,5(d6) AZ8,5(d6) DEPUIS DEPUIS ( L'AFFICHAGE 47 AZ7,5(d7) AZ8,5(d7) DEPUIS Zg5(d7) L'AFFICHAGE 48 AZ7,5(d8) AZ8,5(d8) DEPUIS AZ9,5(d8) L'AFFICHAGE 49 30 Zl,5(d7) Z2,5(d7) Z3,5(d7) Z4,5(d7) Z5, 5(d7) Z1 5(d8) Z2,5(d8) Z3,5(d8) 5 Z4,5(d8) Z5,5(d8) Z1,5(d9) 22,5(d9) Z3, 5(d9) 24,5(d9) Z5,5(d9) Z6,5(d9) Z7 5(d9) Z8,5(d9) Z9 5(d9) AZ2,5(d9) AZ3, 5(d9) AZ4,5(d9) AZ5,5(d9) AZ6,5(d9) AZ7,5(d9) AZ8,5(d9) AZg,5(d9) LEGENDE: d5= dk+4a; d6 = dk+5a; d7 = dk+6a; d8 = dk+7a; d9 = dk+8a Ainsi, en général, pour une matrice fixe ayant M x M entrées (l'ensemble d'électrodes ayant le numéro interne le plus grand ayant l'indice 2M1, M étant égal à 2M+1 = 5 dans l'exemple ci-dessus) 20 l'ensemble autour du marqueu de profondeur d (M+N-1)a k sera construit à partir. sous-jeu des quantités suivantes: Vi,M(dk+(M+r-2)a), i=1,2....,2M-1 r=l,2,...,M AVi, Ml(dk+(M+r-2)a), i=2,3....,2M-l r=l,2,....,M JM(dk+(M+r-2)a), r=1,2,...., M 30 Dans les quantités ci-dessus, l'indice "r" représente un indice d'affichage numérique identifiant des enregis35 trements particuliers appliqués par la production de l'ensemble pour chaque activation de courant par affichage. Par exemple, pour les affichages 46,47...50 pour former l'ensemble de matrice 51 de la figure 4, on voit que jusqu'à cinq affichages, il y a des enregistre5 ments insuffisants pour produire un ensemble. Ensuite, les résultats ci-dessus peuvent être réindexés en format d'ensemble de matrices pour produire l'affichage 51 tel que mentionné précédemment. A noter à cet égard que les entrées de matrice mention10 nées dans l'affichage 51 préservent la relation de un à un des valeurs de courant et de potentiel recueillies-avec l'assemblage de diagraphie 21 aux positions de diagraphie différentes de la figure 4. Ces entrées sont résumées sous forme de tableau au tableau III 15 et ont été annotdes au tableau II à des fins de discussion. En comparant les entrées des tableaux II et III, il est à noter que la profondeur de balayage (Sd1) de l'ensemble de matrices illustré coincide avec le marqueur de profondeur (dk + 6a)qui se trouve deux marqueurs 20 de profondeur en dessous de l'endroit o l'ensemble central a été placé initialement loisque le prélèvement a eu lieu (c'est-à-dire au cycle 1), tandis que la profondeur de balayage prochaine dans le temps (Sd2) est à la profondeur de dk + 7a qui se trouve être une station 25 de diagraphie en dessous de Sd. TABLEAU II- C1 C2 C3 C4 C5 Z5, 5(d5) AZ6,5(d5) AZ(Sdl)=AZ7, 5(d5) AZ8,5(d5) AZ9,5(d5) Z4,5(d6) A5,5(d6) AZ6,5(d6) AZ7,5(d6) 7zS,5(d6) Z3,5() ^Z- 5(d-) LZ4,5(d. 7) aZ5,5(d-7) AZ6,5(d7) AZ7,5( 7) Z2,5(d8) AZ3,5(d8) AZ4,5(d8) AZ5, 5(d8) AZ6, 5(d8) Zl 5(d9) R LZ2 5(d9) R LZ3,5(d9) R3 AZ4,5(d9) R4 AZ5,5(d9) R5 LEGENDE: d5 = dk + 4a; d6 d8 = dk + 7a; d9 = dk + 5a; d7 = dk + 6a; = dk + 8a; Sd1 = d7 = dk + 6a. Il est à noter au tableau III que la première rangée (R1) peut également être obtenue par réciprocité. Ainsi, comme la relation suivante se vérifie, 20 c'est-à-dire pour M=5: Zi,5 (d2M-i) = Z2M-i,5(ds), i = 1,2,.. les matrices des tableaux III et IV sont équivalentes. TABLEAU IV Cl C2 C3 C4 C5 Z5,5(d5) Z4,5(d6) Z3,5(d7) Z2,5(d8) Z1,5(d9) AZ6,5(d5) LZ5, 5(d6) AZ4,5(d7) AZ3,5(d8) A2,25(d9) _Z(Sdl)=AZ7,5(d5) AZ6,5(d6) LZ5,5(d7) AZ4,5(d8) AZ3,5(d9) AZ8,5(d5) AZ7,5(d6) AZ6,5(d7) AZ5,5(d8) AZ4,5(d9) AZ9, 5(d5) AZ8,5(d6) AZ7,5(d7) AZ6,5(d8) AZs5,5(d9) Rl R2 R3 R4 R5 LEGENDE:d5 = dk + 4a; d6 = dk + 5a; d7 = dk + 6a; d8 = dk + 7a; d9 = dk + 8a; Sd1 = d7 = dk + 6a. Les matrices des -tableaux III et IV ont des colonnes indexées C1,...,.C5 (en allant de gauche à droite) et des rangées indexées Ri,....R5 (en commiençant du haut en bas) ce qui est la notation normale 5 des séquences d'indexation des matrices. Suivant cette convention, pour le cas général o les dimensions de la matrice sont M x M au lieu de 5x5, et si on désigne les éléments de la matrice modifiée AZ(Sd1) par les quantités AZ(dk + (M+N-1)a)) j, les éléments sont donnés 10 par les formules suivantes (o dk+(M+N-1)a = Sd1, M+N-1 = 3N): (les éléments de la première rangée) kZ(dk+3Na))l, = ZM+1_ j M(dk+(M+j-2) a), J=1,...,M; (éléments de la seconde à Mième rangée) (LZ'dk+3Na))i,j = AZM+i-j M(dk+(M+j-2)-a, i=2,....M j=1,...,M Il est à noter que les éléments ci-dessus de la première rangée peuvent être obtenus au moyen de la réciprocité comme cela a été fait en passa' t 20 du tableau III au tableau IV ci-dessus pour le c. s o M=5. Ainsi, les entrées recueillies comprenant l'affichage 51 ont les caractéristiques communes suivantes: (i) elles sont liées à des valeurs préliminaires 25 mesurées à l'intérieur de la formation au voisinage du trou de forage; (ii) elles sont indexées sur une zone particulière à l'intérieur du terrain, le symbole d'indice se référant à la station d'électrodes de l'électrode centrale de l'assemblage, et les deux valeurs 30 de profondeur étant les limites de la zone intéressante, et (iii) elles sont également indexées sur un numéro particulier de stations de profondeur de balayage associé à la profondeur de l'électrode centrale. On peut également voir dans les formules 35 ci-dessus que les éléments sur la diagonale principale et au-dessus d'elle (c'est-à-dire les éléments en diagonale et en supradiagonale) résultent de mesures de potentiel au niveau des électrodes qui ne sont pas situées en dessous de l'électrode EM qui délivre le cou5 rant. En outre, on peut voir que les éléments qui sont situés en dessous de la diagonale principale résultent de mesures de potentiel au niveau d'électrodes qui ne sont pas situées au-dessus de l'électrode EM qui délivre le courant. Si l'on se réfère au tableau III en tant que modèle, on peut voir que dans chacune des colonnes, la profondeur de l'électrode de courant est constante et s'accroît de colonne à colonne de gauche à droite. En plus, en partant de la colonne 1, on voit que dans 15 toute rangée donnée, l'indice le plus à gauche (qui est lié à l'ensemble d'électrodes à numéro interne impliqué dans la mesure correspondante de potentiel) décroît de gauche à droite en commençant par 5 (en général indice M) dans la première rangée et en commen20 çant par 9 (en général indice 2M-1) dans la dernière rangée. Il est à noter également que les entrées en diagonale du tableau III en parallèle sur la diagonale principale ont des numéros d'indexation communs de 25 courant et de potentiel mais dans un ordre croissant avec la profondeur. En conséquence de l'indexation des valeurs mesurées d'impédance, telle que montrée au tableau III, le spécialiste qui interprète peut aisément synthétiser les opérations des différents 30 types d'assemblage de diagraphie de manière à fournir une multiplicité de confiarations différentes de pénétration du courant et de là il peutdéterminer plus aisément les modifications de résistivité en fonction de la distance latérale dans la formation sous surveillance. Il est à noter en outre que tandis que l'affichage 51 ne représente qu'une matrice 5 x 5, lors des opérations réelles, il y a beaucoup plus d'éléments par ensemble de matrice. Par exemple, si l'ensemble d'électrodes numéroté en dernier de l'assemblage de 5 la diagraphie est désigné par le nombre entier 2M1 (M étant égal à 5 comme dans l'exemple ci-dessus) et que l'électrode centrale est numérotée M, il y a donc M x M entrées par ensemble de matrice. Cependant, il est à noter qu'une rangée, habituellement la rangée 10 R1, sera occupée par des valeurs d'impédance absolues et que les rangées restantes, à savoir R2,R2...RM seront composées de valeurs de différence de potentiel déterminées en prenant la différence entre les mesures de paires voisines d'électrodes de potentiel de l'assembla15 ge de diagraphie. A cause de l'existence et de la possibilité de mise en oeuvre du principe de superposition linéaire, les matrices montrées aux tableaux III et IV ont la possibilité de lier systématiquement tout jeu arbitraire de courants émetteurs de l'assemblage 20 de diagraphie avec une distribution de potentiel correspondante de tout nombre choisi d'électrodes de potentiel de l'assemblage (ou vice versa) en accord avec AV (dk + 3Na) = AL2 (dk + 3Na) * J (dk + 3Na) o (-) indique une matrice ou une quantité vecteur; la 30 profondeur dk + 3Na correspond à celle du numéro de profondeur de balayage associé à une matrice particulière; et AV et J sont des vecteurs de colonne Ml donnés par AV(dk+3-a) = V1 (--2a) LV2 (dk+(2N,+i)a) I.I I.I 1 VM(dk+4Na) et J(dk+3Na) = Jl(dk+2NLa) I _ J2(dk+a) 1 I. I | I. I. I JM(dk+4Na) l MATERIELS DE FOND ET DE.TETE Pour permettre une commande précise de l'assemblage de diagraphie, la présente invention selon la figure 1 envisage d'utiliser un circuit de commande de surface pour surveiller les opérations de fond, 25 c'est-àdire d'utiliser un dispositif de commande à l'intérieur du circuit 17 de commande et de traitement de la fourgonnette 18 à la surface 15 pour synchroniser les opérations d'une servo-commande à l'intérieur de la section centrale 30 se trouvant dans l'assemblage 30 de diagraphie 21 de la figure 1. Les figures 5 et 6 illustrent de telles opérations en détail, le circuit de tête étant représenté à la figure 5 et le circuit de fond à la figure 6. Comme le montre la figure 5, le circuit 35 17 comprend une horloge maître 60 permettant de commander le circuit de réception et de transmission 61, un circuit logique d'entrée maître 63, et un sortie d'entrée sortie d'ordinateur numérique 64 par l'intermédiaire du circuit logique de synchronisation 65. Des données de diagraphie ayant un format et un caractère devant être décrits ci-.après montent à partir de l'assemblage de diagraphie 21 à travers le câble 12 sur la poulie 13 du derrick 14. Ensuite, les données sont transmises par l'unité de levage 16 par l'intermédiaire du conducteur 19 vers 10 le circuit de réception et de transmission 61 et de là par le circuit logique d'entrée maitre 63 vers l'ordinateur numérique 64. Dans cet ordinateur, les données peuvent être affichées en temps réel sur l'unité d'affichage 66 (par exemple pour vérifier la distorsion de 15 phase) et ensuite après des calculs d'impédance la matrice finale est enregistrée dans l'enregistreur 67. Comme les données de diagraphie sont initialement dans un format qui est incompatible avec le fonctionnement de l'ordinateur, elles doivenrt tout d'abord être dérmulti20 plexées dans le circuit logique d'entrée maître 64 (c'est-à-dire changées d'un format sériel à un format parallèle) et ensuite elles sont indexées avec des marqueurs de profondeur appropriés en provenance du codeur de profondeur 62 attaché à l'unité de levage 25 16. Pour fournir des longueurs appropriées de mots et de blocs aux données compatibles avec le traitement à l'intérieur de l'ordinateur 64, le circuit logique 63 d'entrée maître est commandé soigneusement en utilisant le circuit logique de synchronisation 65 en conjonc30 tion avec l'horloge maître 62. En ce qui concerne le codeur de profondeur 62, il est à noter que pendant les opérations, le codeur de profondeur fournit la profondeuer absolue d'un emplacement de référence de l'assemblage 21 par rapport 35 à la surface 15 (de préférence sur la base de la profon- deur de l'ensemble d'électrodes le moins profond de l'assemblage de diagraphie à chaque position de diagraphie). La profondeur associée aux mesures prenant son origine sur les ensembles d'électrodes restants de 5 l'assemblage 21 de la figure 21 est déterminée à partir de l'espacement connu "a" entre des ensembles d'électrodes voisins E1,E2,..,E2M-l' En fonctionnement, l'horloge maître 60 produit une série d'impulsions de synchronisation (horlo10 ge) qui sont appliquées au circuit logique de synchronisation 65 et de là aux circuits restants pour commander le moment de fonctionnement de tels circuits. Chaque opération nécessite habituellement un certain nombre d'impulsions d'horloge et, par conséquent, l'unité de 15 temps pour pouvoir achever l'une des diverses opérations est un multiple exact des impulsions d'horloge. Parexemple, la lecture du circuit logique d'entrée maitre 63 est réalisée pendant un intervalle spécifique de temps qui est un multiple exact des impulsions horloge 20 de l'horloge maître 60. Lorsqu'un sous-jeu de circuit est bloqué, un nouveau sous-jeu est habilité par les impulsions de temps de manière à pouvoir exécuter des nouvelles opérations. Il est à noter à ce propos que l'unité de commande de format 68 est susceptible de modifications manuelles pendant le transfert des données vers l'ordinateur 64. De cette manière, le format des données peut être modifié pour satisfaire à différentes exigences d'interprétation qui ont lieu au niveau de l'unité 66 d'affichage en temps réel et sur l'enregis30 treur de données 67 de la manière expliquée précédemment. Comme indiqué ci-dessus, la figure 6 montre en détail le circuit de fond. Tel que la figure le montre, une horloge 80 en conjonction avec le circuit logique 81 de synchro35 nisation et de commande est utilisée pour commander les opérations des ensembles d'électrodes E1,E2...,E2M_1 de l'assemblage de diagraphie indiqué de façon générale par le numéro 21 à la figure, en conjonction avec le circuit de synchronisation de tête du circuit 17 de commande de traitement et en réponse à celui-ci. Pour plus de détail, on suppose que l'horloge 80 tourne librement et un signal d'initialisation 79 du circuit de tête 17 apparaît dans le circuit logique 81 et amorce les opérations. On suppose également qu'au moyen du 10 signal de début 79, l'initialisation a eu lieu à l'intérieur des autres éléments de circuit de l'ensemble de circuit de fond. En d'autres termes, en utilisant un signal d'initialisation en provenance du circuit logique 81 de commande et de synchronisation, les élé15 ments suivants: dispositif de format/mémoire tampon d'émission 85; convertisseur 86 analogue/numérique logique de commande de gain 87; et multiplexeur 89 sont initialisés correctement pour débuter l'opération. Simultanément, le compteur 90 associé au 20 générateur de courant 91 est initialisé et les opérations commencent pour permettre l'injection de courant par l'intermédiaire de l'ensemble d'électrode EM. En d'autres termes, après que le générateur de courant 91 a été activé pour injecter du courant 25 dans la formation voisine par l'intermédiaire de l'électrode de courant EcM, l'intensité de courant, ainsi que les potentiels particuliers absolus et de différence sur les ensembles d'électrodes de potentiel E1,E21..., EM1 sont mesurés et passés aux circuits d'acquisition M-i des données. Plus particulièrement, les potentiels absolus et les potentiels de différence sont tous passés au multiplexeur 89, au convertisseur analogique/numérique 86 et au formateur/mémoire tampon d'émission 85 par l'intermédiaire des circuits d'échantillonnage/maintien 35 88a, 88b,...,88m et 84b,84c,...84m en utilisant des amplificateurs binaires particuliers de gain d'une série de tels amplificateurs indiqués généralement par les numéros 97 et 98. Il est à noter que de tels comparateurs 83b et 83c sont reliés en permanence aux 5 sorties des amplificateurs 97a et 97m respectivement. De cette manière, des mesures de phase entre l'électrode de courant d'excitation et les électrodes de potentiel de l'ensemble éloigné peuvent être obtenues. Intensité de courant A la figure 6, on voit que l'intensité du courant est mesurée par l'intermédiaire d'un amplificateur de gain binaire 100 dont le gain est également commandé par la logique de commande de gain 87 et dont le signal de sortie est appliqué à un multiplexeur 15 89 par l'intermédiaire du circuit 88 de l'échantillonnage et de maintien. Pour une gamme utile de courants de 500 micro-ampères à 10 ampères, la résistance 101 se trouvant dans le trajet de courant ne devrait pas dépasser 0,1 ohm, ce qui donne une entrée de tension vers l'amplificateur 100 qui se trouve dans une gamme de 500 microvolts à 1 volt. Par conséquent, son gain programmable est formé par des étapes binaires allant de 5 à 100 000 et nécessite au moins un code de gain de 15 bits. Mesures de phase Pour mesurer la phase des potentiels sur certains ensembles d'électrodes, les périodes intermédiaires entre le passage par zéro du signal en fonction de la référence de phase, c'est-à-dire le début du 30 courant sinuso3dal,sont mesurées. Le contenu du compteur sert de référence de phase et est chargé dans les registres de phase 92b et 92c à l'instant précis o les comparateurs 83b et 83c détectent un passage par zéro des signaux de potentiel correspondants. Dans 35 la présente forme de réalisation, les phases des poten- tiels sur les deux ensembles d'électrodes éloignés c'est-à-dire les phases à E1 et E2M-1 par rapport au courant qui est injecté par l'intermédiaire de l'ensemble EM sont recherchées. Toute modification dans les compta5 ges indiquées dans un registre de phase donné est directement proportionnelle à la distorsion de phase, ce qui fournit une indication directe de la fiabilité des mesures associées. Mesures des potentiels absolus et différentiels Pour fournir une intensité de courant sur l'électrode EcM de l'ensemble d'électrodes EM, la séquence suivante doit avoir lieu: tout d'abord le compteur 90 est remis à zéro par l'intermédiaire d'un signal approprié provenant de la logique de synchronisation 15 et de commande 81. Les impulsions d'horloge à l'entrée du compteur 90 fournissent des incréments à son contenu jusqu'à ce que soit produit un cycle complet de courant. Le convertisseur numérique/analogique 102 de contrôle de la sinusoïde convertit ensuite le contenu du compteur 20 90 pour produire une série de valeurs discrètes de co' dont-les amplitudes individuelles varient sinusocalement en fonction du temps. Après amplification par l'intermédiaire d'unamplificateur 103, le courant variant sinuso dalement est discriminé par le commutateur 25 de courant 95 vers l'électrode de courant EcM de l'ensemble d'électrode EM et de là, dans la formation voisine de la manière décrite précédemment. A cet égard, on suppose que les ensembles d'électrodes qui constituent l'assemblage général de diagraphie qui est décrit actuel30 lement ont des impédances internes essentiellement infinies de sorte qu'elles ne prélèvent aucun courant appréciable du milieu environnant et ils sont comme décrit précédemment des petits anneaux de peu de volume de sorte que leur présence n'a pas tendance à altérer 35 de façon importante le champ de potentiel au voisinage de la surface externe de l'assemblage. Additionnellement, le courant en provenance de l'électrode de courant doit bien sûr effectuer son retour pour fermer le circuit de mesure et ceci est accompli au moyen de l'électrode de retour 28 située à l'écart, comme lemontre la figure 1. L'électrode de retour 28 apparaît située à l'infini pour tous buts pratiques éventuels. Ensuite, on fait des mesures des potentiels absolus de tous les ensembles (par rapport à l'électrode 10 de référence 29 de la figure 1) ainsi que tous les potentiels différentiels sur toutes les paires voisines d'électrodes de potentiel Epl,EP2,...,EP2M-1 des ensembles d'électrodes E1,E2'...,E2M-1 comprenant l'assemblage 21. En d'autres termes, les potentiels absolus aux 15 électrodes de potentiel Epl,-Ep2,...nEp2M_1 sont indiqués par l'intermédiaire des amplificateurs 97a,97b,....,97m tandis que les différences de potentiel de paires adjacentes d'électrodes de potentiel Ep1, Ep2,....nEp2M-1 sont- mesuré-s au moyen d'amplificateurs différentiels 20 à gain bina- e 98b,98c,...,98m. iite, le processus de prélèvement est répète en utilisant le circuit de discrimination 95 pour réactiver l'ensemble d'électrodes EM par du courant. Le potentiel absolu est mesuré sur toutes les électrodes 25 Epî,EP2.... EP2M 1 au moyen d'amplificateurs 97a,97b,..., 97m tandis que des différences de potentiel des paires voisines d'électrodes Ep1l,Ep2,..., Ep2M_1 sont mesurées au moyen d'amplificateurs 98b,98c,...,98m d'une manière semblable, comme indiqué ci-dessus. Tout le processus 30 de prélèvement est répété chaque fois que l'assemblage de diagraphie s'est déplacé d'une station de diagraphie soit vers le haut, soit vers le bas. Il est à noter que le procédé de prélèvement décrit ci-dessus peut être répété entre les stations 35 de diagraphie si on le désire. Le fait d'acquérir con- tinuellement des données de -cette manière permet divers schémas de renforcement des signaux, ce qui donne des estimations améliorées d'impédance par rapport aux bruit environnants. Il est à noter également que lorsque le cycle de prélèvement est répété, les mesures sont à noter à la surface 15 et ensuite traitées pour fournir des entrées d'impédance d'une série d'ensembles finis de matrices qui se chevauchent associés à une série 10 d'intervalles de profondeur. En d'autres termes chaque ensemble d'entrées d'impédance est indexé sur un intervalle." sélectionné de profondeur central dépendant de la longueur active de l'assemblage des électrodes (entre électrodes E1 et E2M_) comme décrit précédemment. Il reste maintenant à indiquer brièvement comment les données peuvent être mesurées et ensuite transmises en tête sur la base d'enregistrements de données divisés en mots et blocs compatibles avec l'ordinateur 64. Une brève description du format de prélèvement 20 et de transmission st dans l'ordre des choses; elle a lieu ci-dessous. Formats de prélèvement et de transmission Un enregistrement de données consiste en informations de profondeur fournies par le codeur de 25 profondeur 62 de la figure 5, suivies par des données recueillies par l'assemblage 21 de la figure 1, comme esquissé précédemment. De façon plus spécifique, l'outil se trouvant à une profondeur présélectionnée, les informations 30 de profondeur constituent l'entrée de l'ordinateur numérique 64 par l'intermédiaire de la logique d'entrée mâitre 63 et un signal de départ vers l'outil est transmis par l'intermédiaire du circuit 61 de réception et d'émission. Comme décrit précédemment, ce signal 35 de départ amorce le processus de prélèvement des données de sorte que (i) des données de potentiel absolu, (ii) des données de potentiel différentiel, (iii) des données d'intensité de courant et (iv) des données de phase sont transmises séquentiellement par l'intermédiaire 5 de la mémoire tampon 85 de production de format et de transmission en une séquence prédéterminée au circuit de traitement et de commande 17 comme indiqué à la figure 6. Chaque mot de données consiste en une sortie 10 de 16 bits en provenance du convertisseur analogique/ numérique plus un maximum de 16 bits pour le code de gain,correspondant. Les mesures sur les électrodes sont indexées en gain. Les gains doivent être réglés dans des étapes binaires de 10 à 3500 pour des mesures 15 de potentiel absolu (nécessitant un code de gain de 9 bits); de 10 à 200 000 pour des mesures de différence de potentiel (nécessitant un code de gain de 15 bits) -et de 5 à 100 000 pour des mesures de courant (nécessitant un code de gain de 15 bits). La mesure dephase 20 ne nécessite pas un code de ú-in. La quantité de - _ es contenues dans un enregistrement de données tel qu'esquissé ci-dessus est déterminée par le nombre des ensembles d'électrodes actives sur l'outil. Par exemple, en supposant que 25 145 ensembles d'électrodes actives sont utilisés (c'està-dire M = 73, N = 36, de sorte que 2M-1 = 4N + 1 = 145), chaque fois que le processus de prélèvement décrit ci-dessus est répété, ceci donnera 145 potentiels absolus, 144 potentiels différentiels, une intensité de courant, 30 deux mesures de phase et en correspondance, 290 codes de gain. Il en résultera: (292 + 290)*(16) = 9312 bits d'informations par enregistrement de données, ou cycle de prélèvement. Simultanément, avec le traitement en temps réel des données acquises, il est possible de tranférer les informations à un dispositif de mémorisation tel qu'une bande magnétique à laquelle on peut accéder 5 ensuite pour le traitement. Le format particulier qu'on choisit dans ce but est quelque peu arbitraire mais devrait être éventuellement dicté par des raisons pratiques relevant des caractéristiques du dispositif d'enregistrement et de la nature du traitement subséquent. A des fins d'identification complètes, les enregistrements de données sont précédés par des informations de tête en plus des informations habituelles décrivant le temps, le site, etc. et devraient également contenir des données décrivant des paramètres tels 15 que l'espacement entre électrodes "a", le diamètre du mandrin, le nombre d'électrodes actives utilisées lors du prélèvement des données, l'emplacement de l'électrode de référence, etc. Pour plus de clarté, les points suivants 20 sont répétés: (1) Même si i'ensemble d'.:trodes EM bifurque pour former des premiers et seconds petits anneaux, ceux-ci occupent essentiellement les mêmes coordonnées de profondeur. En d'autres termes, les distances axiales 25 entre les petits anneaux et les dimensions axiales de ces derniers sont assez faibles en comparaison de la distance entre des ensembles voisins d'électrodes pour que,en pratique, l'ensemble d'électrodes EM ait les mêmes caractéristiques que les ensembles restants qui 30 sont constitués seulement par des électrodes de potentiel. (2) Pour manipuler les quantités en tantque fonctions des différentes variables d'une manière rapide, les ensembles d'électrodes E1,E2,..., E2M_1 sont 35 numérotés en utilisant le schéma d'indexation défini précédemment en commençant en haut du mandrin et en terminant au point le plus bas du mandrin dans le trou de forage. L'ensemble d'électrodes centrales est numéroté M et l'ensemble d'électrodes placées le plus bas est 5 numéroté 2M-1 (M = 2N+1). En conséquence, une fois que N est établi, la numérotation des ensembles d'électrodes en résulte directement. Par conséquent, en supposant que N=36, il -en résulte que M est 73 et que les ensembles d'électrodes au-dessus de l'ensemble central 10 seraient E1,E2,. * ,E72, et ceux en dessous seraient..DTD: E74'E75.. E145À (3) Les quantités de potentiel mesurées sont également indexées, selon les systèmes décrits plus haut. En d'autres termes, le courant d'électrodes 15 étant émis depuis l'ensemble d'électrodes EM à la profondeur dk, les potentiels absolu et différentiel des ensembles E1E2... E2_1 sont indexés sur la base de l'ordre de numérotation interne et l'ensemble actif comme décrit précédemment, dk apparaissant Un tant qu'argument pour désigner la profondeur d.: station occupée par l'électrode d'émission de courant On peut déterminer facilement à partir des valeurs obtenues ci-dessus de courant et de potentiel, leurs rapports d'impédance,et on peut les indexer en fonction 25 de la position de l'électrode de courant et de la position de l'électrode de potentiel, comme décrit précédemment. -En continuant à déplacer l'assemblage vers le haut ou vers le bas le long du trou de forage à 30 mesure que l'assemblage est forcé de répéter le processus de prélèvement à chaque station de diagraphie, on voit qu'il est possible de développer une série de mesures d'impédance indexées sur la position de l'électrode de courant et la position d'électrode de potentiel. 35 Ensuite, ces valeurs peuvent être réindexées pour former 4 7 un ensemble de matrice d'entrées d'impédance associées avec la formation voisine du trou de forage qui traverse la formation en étude, et normalisées sur la profondeur de la zone de traversée de la matrice; sur le nombre 5 des stations de balayage; ainsi que sur la profondeur de l'ensemble d'électrodes centrales. Il est à noter que par suite du principe de la superposition linéaire, les ensembles de matrices d'impédance fournis par le procédé de la présente inven10 tion ont la capacité de lier un jeu arbitraire de courants en provenance d'un assemblage d'électrodes M à la- distribution correspondante de potentiel absolu et de potentiel différentiel apparaissant sur les électrodes M, par rapport à l'électrode de référence éloi15 gnée, telle que mentionnée précédemment. Il est à noter finalement que les ensembles d'impédance de la présente invention peuvent comprendre systématiquement un très grand nombre d'entrées, à savoir que si la dernière électrode est indexée 2'20 1, il y a M x M entrées dans chaque ensemble, av-c une rangée constituée par des valeurs d'impédance dér. vées dC.o mesures de potentiel absolu sur les ensembles d'électrodes de l'assemblage. Par conséquent, le plus grand nombre des entrées sont constituées par des valeurs 25 d'impédance calculées à partir de différences de potentiel entre les paires adjacentes d'ensembles d'électrodes, à savoir égales à M-1 rangées d'entrées. Ces entrées définissent une matrice modifiée d'impédance de la formation, indexée comme décrit précédemment. 30 EMPLOI FINAL DES DONNEES Pour pouvoir déterminer de façon précise la résistivité de la formation en fonction de la distance latérale depuis le trou de forage, on suppose que les entrées d'impédance améliorées de chaque ensemble 35 de matrices ont été déterminées en accord avec les étapes ci-dessus. On suppose également que la matrice réciproque comprenant chaque ensemble a été produite. Alors, sur la base du principe de superposition, tout jeu arbitraire de différence de tension entre des ensem5 bles d'électrodes voisins par rapport à une réponse de courant correspondante peut être lié par une telle matrice déterminée d'impédance réciproque. C'est-àdire que pour une série de différences de tension de commande entre un nombre sélectionné d'ensembles d'élec10 trodes simulant différents assemblages d'outils focalisés, on peut exprimer des réponses correspondantes de courant en tant que systèmes d'équations linéaires dont les termes comprennent des vecteurs de colonne AV, J et la matrice réciproque AZ 1 à savoir: 15 J = AZ-1 AV o AVT = (V1(V2-V1)(V3-V2)...(VM-VM-1)) jT =(J1 J2 ' ' JM) o AV et JT désignent les transposées des vecteurs 20 AV et J, respectivement. Tandis que l'art antérieur a suggéré (i) comment utiliser les réciproques de la relation cidessus, c'est-à-dire comment lier une série de tensions synthétiques calculées en produisant un jeu de vecteurs 25 de courant avec la matrice de résistance de la formation en étude (voir Stelzer, ouvrage cité) et (ii) comment construire un réseau équivalent simulant une formation de terrain appliquant la matrice de résistance réciproque pour aider au calibrage des outils électriques (voir 30 THE LOG ANALYST, mai-juin 1979, "A Theory of Equivalent Artificial Networks Simulating The Subsurface Formations, and Their- Application In Electric Well Logging", Zoltan Barlai), personne à la connaissance de la Demanderesse n'a envisagé un procédé pour lier systématiquement une caractéristique 35 déterminable de la- formation autre que des modèles de tension synthétiques, en général ou en particulier, pour une utilisation lors de la détermination d'une caractéristique de la formation (ou des séries de caractéristiques) en fonction de l'invasion du filtrat de boue et de contrastes de résistivités différentes dans le trou de forage et aux environs. Par exemple, dans Stelzer, ouvrage cité, la résistivité apparente d'un assemblage synthétique sélectionné est déterminée par la combinaison linéaire 10 des valeurs normalisées de tensions synthétiques divisées par une quantité produite avec une combinaison linéaire des réciproques des distances linéaires effectives absolues parmi les électrodes. De tels résultats n'ont eu lieu seulement qu'après que les tensions synthétiques 15 ont été calculées en se basant sur la production de différents jeux de modèle de courant avec une matrice de résistance précédemment déterminée pour la formation en étude. De plus, en supposant que les rayons effectifs des électrodes de l'assemblage de diagraphie sont négli20 geables, on peut dire que la résistivité est une fonction directe des entrées de tensions synthétiques, Mais l'expérience a indiqué que souvent d'autres facteurs ont un effet important sur la réponse de l'assemblage spécialement lorsqu'une invasion importante du filtrat 25 de boue est ae grands contrastes de resistivité différents sont présents dans le trou de forage et autour de celui-ci. La figure 7 est un réseau équivalent artificiel de bornes qui simule de façon originale une forma30 tion de terrain sur une longueur équivalente égale à M stations de diagraphie le long d'un trou de forage traversant la formation en étude et illustre la correspondance unique des admittances de ce réseau et des entrées de chaque ensemble de matrice d'impédance réci35 proque telles que résultant de l'assemblage d'électrodes et du circuit de tête de la figure 1, même là o l'invasion du filtrat de boue dans la formation est importante. En d'autres termes, la figure 7 montre que chaque ensemble de matrices d'impédance réciproque représente un réseau équivalent de bornes qui simule de façon précise la formation en étude même dans des conditions difficiles de diagraphie. Comme illustré, le réseau comprend un assemblage d'adrittances localisées 140 et 144 - reliant ensemble 10 les bornes P1,P2...P2N+ ainsi que respectivement à la borne de masse 145. Les bornes de ce réseau correspondent à M stations de diagraphie qui définissent l'intervalle de balayage intéressant, P1,correspondant à la station la plus élevée et P2N+ 1 correspondant à la 15 station la plus basse. Le nombre total d'admittancesdistinctes dans le réseau est de (N+1)(2N+1). Il en résulte que s'il y a M stations d'électrodes dans la zone intéressante, il y a M (=2N+1) bornes qui interconnectent les admittances localisées à la borne de masse 20 145. Si des courants arbitraires d'intensités J1,J2....JM sont appliqués aux bornes P1,P2,...P2N1 et que des tensions correspondantes de noeuds sont définies comme étant V1,V2,...,VM, l'application de la loi de Kirchoff à chaque borne fournit le jeu suivant d'équations: J1 = Yi,lVl + Y1,2(V1-V2) +...+ Y1,M(Vi-V.M1) J2= Y2,1(V2-V1) +Y2,2V2 + + Y2,M(V2-V.). ....DTD: . ...DTD: J =Yl,l(VN-Vl) + YM,2(VZl -V2) + '+ r,MV o la définition Y. = Y. est impliquée. Ili ji] Ce jeu d'équations peut être arrangé d'une manière évidente pour fournir le système suivant: J1 = Ylv1 - Y l,2V2 - ' - YlmVI J2 = -Y2,lVl + Y2V2 -.... - Y2,MVI JM = -,1lV - YM1,2V2 -..... + YMVI4 10 o les termes en diagonale de la matrice sont obtenus en additionnant toutes les admittances de liaison à partir de la borne correspondant à la rangée (ou colonne) donnée, c'est-à-dire Yp = Yp,1 + Yp,2 + + YpM p=1,2,...,M Pour la simplicité et la clarté de la notation qui va suivre, les coefficients yi et Yij seront remplacés en outre par des X selon le schéma suivant, c'està-dire Xij = -Y i non égal j X. = Y. 1,1 1 de sorte que le système ci-dessus d'éauations peut être réécrit de façon appropriée comme suit: J1 = Xlv1 + X1,2V2 +... + X1 MVM J2 = X2,lVl + X2, 2V2 + ' ' À + X2,MVM J30M= XM lVl + XM,2V2 + + XM,MVM Ce jeu d'équations concerne les tensions absolues; cependant, on désire les placer sous une forme impliquant des tensions de différence entre des bornes adjacentes. 35 Dans ce but, et à la suite de la notation utilisée plus haut, il est possible d'exprimer ces tensions d.e différence AV selon les relations AVj = V = Vj, j=2,3,...,M Vp = V = AV2 + AV3 +... + AVp, p=2,3,.... M Ceci a son tour permet au système ci-dessus d'équations concernant les tensions et courants absolus d'être réécrits comme suit: = U1,jVl + Ul, 2AV2 + '.. + U1,MAVM J2 = U2,1V1 + U2,2AV2 + * * ' + U2,MAVM JM = Ut,1 V1 + UM,2AV2 +. . + UM,MAVM o UPq = +X +. + Xpq + M pq=1,2,..., M UXp, q p, q +l ', ''' Le système précédent d'équations donne le résultat 20 désiré, c'est-à-dire qu'en notation de matrice, on a la forme suivante: J = U AV o AVT = (V1 (V2 - V1) (V3 - V2).. (V - VMi)) JT = (J1 J2 J M) et les éléments de la matrice U sont des combinaisons simples linéaires uniques des admittances localisées du circuit des paramètres Yi. (par l'intermédiaire des X). Par conséquent, en supposant que l'inverse de la matrice d'impédance modifiée AZ fournie par l'assemblage de diagraphie de la figure 1 existe, on peut alors déterminer à partir de la relation U = AZ1 de façon uniique les admittances ii pour mettre le circuit équivalent dans une correspondance un à un avec la réciproque (l'inverse) de la matrice d'impédance modifiée. Ceci signifie qu'une discussion des diverses propriétés de l'assemblage de diagraphie peut être effectuée en utilisant le réseau équivalent. C'està-dire qu'en accord avec la présente invention, la correspondance parfaite et unique de la matrice d'admit10 tance du réseau de la figure 7 et la matrice d'impédance réciproque fournie par l'appareil de la figure 1 Deut être utilisée de la manière suivante pour fournir un procédé systématique de détermination d'une série de caractéristiques de réponse de manière à prendre 15 en compte les contrastes des différentes résistivités et l'invasion du filtrat dans le trou de sondage et autour de celui-ci. En d'autres termes, la présente invention décrit un procédé pour faire la synthèse des caractéristiques de réponse vraies d'une combinaison 20 de différents outils de diagraphie électrique centrés sur le trou dans des conditions variées et difficiles de trous de sondage, le tout premièrement en déterminant les valeurs d'impédance d'une formation de terrain traversé par un trou de sondage rempli d'une boue de 25 forage d'une résistivité Rm et (2) en manipulant sélectivement les valeurs résultantes d'admittance en tant qu'entrées d'une matrice réciproque de manière à faire la synthèse des opérations des différents outils centrés sur le trou sur un incrément associé de profondeur 30 avec une précision surprenante, chacune des matrices d'impédance réciproque étant associée à un ensemble de matrices indexé sur un incrément parmi une série d'incréments finis chevauchant le balayage de profondeur de la formation mesurés le long du trou de sondage. 35 Chaque incrément de balayage dépend de la longueur L de l'assemblage des électrodes pour définir des marqueurs de profondeur peu profonds et profonds et se trouvent indexés au centre sur la profondeur dans. le trou de sondage d'un ensemble d'électrodes centrales de l'assemblage au moment du prélèvement des données. Afin de fournir une indication vraie de la résistivité Rt de la formation même lorsque la formation est espacée du trou de sondage par une zone envahie d'une résistivité Rxo d'une étendue latérale inconnue 10 à cause de l'invasion du filtrat de la boue de forage, le procédé de la présente invention spécifie les étapes suivantes: (i) tout d'abord, un assemblage de M ensembles d'électrodes centrés sur le trou ayant un espace15 ment "a" entre électrodes incrémentiel égal est étalonné de manière à obtenir des jeux de facteurs d'étalonnaae normalisés sur des modèles connus d'amorçage de tension dans une zone de réponse de résistivité connue, les ensembles d'électrodes ayant un indice connu de numérota2 tion d'ordre interne et chacun comprenant une électrode de courant et une électrode de potentiel, les jeux de facteurs de calibration étant chacun adressables en fonction des conditions du trou de sondage y compris les conditions difficiles de trous de sondage ainsi 25 qu'au moyen d'un type d'assemblages synthétiques focalisés par ordinateur; (ii) ensuite, un assemblage sur le terrain est placé dans le trou de sondage, cet assemblage ayant un ensemble d'électrodes centrales comprenant une élec30 trode de courant et une électrode de potentiel, et une série d'ensembles d'électrodes additionnels dont chacun comprend seulement une électrode de mesure de potentiel.-La profondeur absolue d'au moins un ensemble d'électrodes est connue de façon continue par rapport 35 à un niveau de données de profondeur Prédéterminée mesurée depuis la surface terrestre; (iii) ensuite encore, l'assemblage est déplacé de façon continue le long du trou de sondage et un courant de valeur connue est injecté depuis 5 l'électrode centrale de courant par la boue de forage d'une résistivité Rm et la zone envahie de résistivité Rxo et de là dans la formation de résistivé Rt; (iv) pendant l'injection du courant, le potentiel absolu de chaque électrode de potentiel ainsi 10 que les différences de potentiel entre des paires adjacentes d'électrodes de potentiel sont mesurés d'une manière rapide, les mesures de potentiel pouvant être indexées sur des stations communes de diagraphie d'une série de stations de diacraDhie éauidistantes localisées le long du trou de forage et ayant une distance d'espacement incrémentiel de "a", cette distance étant la distance entre des ensembles voisins d'électrodes; (v) ensuite, des valeurs d'impédance en provenance des potentiels mesurés absolus et différen20 tiels et de leur courant d'injection associé sont calculées -- ae valeur étant indexée sur les éléments d'indexation interne connus d'électrodes actives de courant et de potentiel utilisées dans les mesures; (vi) ensuite encore, les valeurs d'impédance 25 sont réindexées dans des entrées d'impédance d'une série d'ensembles de matrices modifiés qui se chevauchent AZ, chaque ensemble AZ étant associé à un segment prédéterminé de ladite formation, égal en étendue verticale à M stations de diagraphie, et comprenant M x M entrées 30 d'impédance, M étant l'indice de numérotation interne de l'ensemble d'électrodes de courant dans ledit assemblage et le rapport du nombre des entrées d'impédance de différence par rapport aux entrées absolues étant d'environ M-1:1; (vii) après cela, chaque ensemble de matrices AZ est inversé pour former un ensemble modifié réciproque de matrice LZ- 1 en accord avec des techniques classiques d'inversion de matrice (viii) ensuite, des paramètres de réponse focalisées par ordinateur sont produits en utilisant la matrice réciproque AZ- 1 de l'étape (vii) en conjonction avec les mêmes modèles d'amorce de tension de l'étape (i); (ix) finalement, des jeux de facteurs d'étalonnage de l'étape (i) sont recherchés jusqu'à ce que les produits d'un jeu particulier de facteurs d'étalonnage et des paramètres de réponse de l'étape (viii) pour tous les assemblages d'outils synthétiques 15 soientessentiellement une constante afin que la condition de sondage difficile soit déduite même en présence de contrastes marqués entre la résistivité réelle et celle de la boue et indifféremment du fait que des jeux synthétiques de modèles de potentiel aient été utilisés comme 20 initiateurs d-s paramètres de réponse, focalisés par ordinateur, ite produits. Comme décrit brièvement dans la présente invention, de la manière exposée ci-dessus, les étapes nécessitant des commentaires additionnels sont expliqués 25 avec plus de détail ci-dessous. ETAPE (i) Dans cette étape, il est à noter que pour calibrer un assemblage de M ensembles d'électrodes centrées sur le trou sur une série de différentes condi30 tions de trous de sondage et de réponsesde l'assemblage, le nombre des ensembles d'électrodes de l'assemblage doit correspondre à celui utilisé lors du prélèvement des données de diagraphie sur le terrain. Dans cet exemple, M est supposé être égal à 73 (2M-1=145), le facteur d'espacement "a" est de 7,5 cm, le diamètre du mandrin est de 9,5 cm et le diamètre du trou de sondage est de 20,3 cm. Il en résulte que la longueur totale de l'assemblage actif L est de 2Ma = 18,24 m. Additionnellement, différents contrastes de résistivité (et différentes distances d'invasion de filtrat) doivent être supposés et des jeux de facteurs d'étalonnage doivent être produits, ces facteurs étant indexés sur différentes 10 conditions de trous de sondage, ainsi que sur des modèles normalisés d'amorçage de tension et une réponse de courant qui y est associée. Spécifiquement, une série de contrastes entre la résistivité de la formation et celle de la 15 boue est tout d'abord supposée dans un trou de 20,3 cm avec des conditions connues d'invasion et des réponses caractéristiques connues d'assemblages d'outils synthétiques calculés en appliquant des techniques appropriées de résolution au prob'.ème bien connu de la valeur limite 20 décrivant les modèles de courant et de tension dans des milieux conduct. A cet égard, voi- V. N. Dakhnov, (1962), "Geophysical Well Logging"!', traduction par G. V. Keller, Magazine Colorado School of Mines, volume 57, N 2, chapitre 3. Ensuite, des jeux de facteurs d'étalonnage sont déterminés et ainsi le produit des réponses calculées et du jeu de-facteurs d'étalonnage est égal à la résistivité caractéristique désirée qui peut généralement être prise comme étant Rt. On suppose à cet égard que les facteurs d'étalonnage sont désignés par le symbole kq, et alors les résultats ci-dessus sont définis par une série de relations de la forme Raq = k Rq, q = 1,2,...,N+1 = K, 35 dans laquelle Ra représente maintenant la réponse q de résistivité apparente de la formation en question et R est la valeur de résistance brute calculée. Les q facteurs d'étalonnage peuvent alors être déterminés en posant: Ra = Rt, q = 1,2,....,N+1, par exemple, pour donner la réponse désirée de résistivité pour la formation particulière considérée. En accord avec la présente invention, les assemblages d'outils synthétiques désirés peuvent être décrits en supposant 10 que le réseau de la figure 7 simule la formation en étude et que l'assemblage donné d'outils à synthétiser est positionné et fonctionne de manière qu'initialement l'ensemble d'électrodes central soit choisi comme l'électrode de mesure du courant, c'est-à-dire sur l'ensemble 15 central (N + 1) et l'intensité du courant à l'intérieur de la formation voisine en question est déterminée lorsqu'il y des jeux choisis de potentiels aux bornes P1I P2,', 'PM du circuit de la figure 7. Les réponses dépendant de la profondeur seront donc ainsi caractéri20 sées initialement par le rapport VN+ /JN+, et le courant JN+1 sera déterminé par -elation suivante: JN+=l Y+i,i(VN+1 -V)+...+YN+i,N+lVN+l+.+'..+YN+,M(V N+lV-) On note ainsi que la détermination cidessus du courant JN+1 implique seulement l'admittance reliant la (N+1)ième électrodes à la borne de masse 145 de la figure 7, ainsi que les admittances entre 30 le (N+1)ième ensembles d'électrodes et les ensembles restants. De là, on peut isoler diverses admittances ou combinaisons en imposant des distributions appropriées de tension sur les ensembles d'électrodes de l'assemblage. Par exemple, on suppose qu'une série de 35 vecteurs de potentiel linéairement indépendants S1, S2, SNSN+ sont imposes sr ies er.s.'.es ci'électrodes de l'assemblage de diagraphie ayant la distribution montrée au tableau V. TABLEAU V 0 0 1 1 0 1 o si =. I.0 I.1 = 10 1* o | 0 I * 1 À '1 I. I l o 20 1 0 I I. I. _9= I l 1 1 --N+i-- 1 1 0. la la a0 0 1 1- il 0 0. I b= I. 1-... 1 1I 11 1 -N+l--| l |--N 1 1 1 1 I. I À.I lai l 0 1 I I l sl=l +1-- - 0i 1h * I il il il il 1h * I * I h1 ol 1 l 1 1 1 SN+1= À * I 11 1 1 1 1 1 1 1 1 1 1 1 1 1 1 LEGENDE: S1 contient 1 valeur unitaire de potentiel S2 contient 3 valeurs de potentiel unitaires SN contient 2N-1 valeurs de potentiel unitaires SN+ 1 contient 2N+1 valeurs de potentiel unitaires Comme on le voit, chaque vecteur est symétrique autour de l'ensemble N+1 et présente des amplitudes unitaires individuelles soit de ZERO soit de UN. Le nombre des UNS dans toute colonne M x 1 de tout vecteur S est en accord avec 2p-1 o p=1,2...N+1. 2p-1 oû p=1,2...N+1. Tandis que le vecteur S1 a un seul potentiel d'amplitude UN (à l'ensemble N+1) pour une distri bution en tiroir de 1 x 1 autour de l'ensemble N+1, le vecteur S2 a une longueur de distribution de trois définissant une distribution en tiroi r de 3 x 1 autour du même ensemble, tandis que le vecteur S3 a une distribution en tiroir de x 1 au même endroit. Par ailleurs, pour les vecteurs SN et SN+1 10 o N=36, si M est 73, les distributions en tiroir sont 71 x 1 et 73 x 1, respectivement, comme le mpntre le tableau V. A mesure que chacun des vecteurs de potentiels caractéristiques du tableau V est appliqué, un 15 à la fois, en partant du vecteur S N+1, les valeurs correspondantes de courant peuvent être déterminées séquentiellement pour l'assemblage d'électrodes synthétisées. A des fins d'annotaticnen accord avec la présente invention. la composante de courant résultante est appleée 20 le courant de mesure pour le vecteur de potentiael correspondant de la distribution désignée, à sav - appelé J(SN+1) o SN+1 désigne le vecteur de potentiei à partir duquel on détermine le courant. De là, pour calculer le courant de mesure à l'électrode N+1 en partant du vec25 teur SN+ 1, c'est-à-dire pour des tensions identiques UN sur tous les ensembles d'électrodes, la (N+1)ième rangées de la matrice d'admittance produite et du vecteur de potentiel SN+ l fournit: J(S)Y J(-N+1) YN+1i,N+1 o J(SN+1) indique que le courant de mesure est en réponse au vecteur de potentiel SN+I' La quantité YN+1,N+1 peut être considérée comme l'auto-admittance de l'ensemble (N+1I) d'électrodes dans la formation particulière en question. De même, pour la distribution de potentiel décrite par SN, la (N+i)ième rangée de la matrice d'admittance produite et du vecteur SN fournit J(SN) y= N+1,1 + YB+1,N+1 + YN+I,M Et en outre, pour le vecteur S N-, la (N+1) ième rangée de la matrice d'admittance produite et du vecteur de potentiel SN-1 est égale à J(SN-1) YN+1,1+YN+1,2+YN+1 +YN+1+1,M-1+YN+ MM Le modèle peut se discerner aisément, c'est-à-dire que 10 chaque courant de mesure est simplement la somme de l'auto-admittance de l'électrode de N+1 plus toutes les admittances qui la relient aux autres électrodes ayant un potentiel de ZERO (sur la base des ZERO dans l'ordre de distribution du vecteur de potentiel). Il 15 en résulte facilement que l'expression générale pour le courant de mesure est donnée par J(SS) = YN+ 1,N+'+1 q=N+1l, et J(s) = N-(q-l) N J) SUl4 y + SUÂI Y + Y, q=l,2,..,N i=l N+l,i i=q N+l,i+N+l N+l,N+l Comme la distribution du vecteur de potentiel caractéristique est telle que toutes les composantes 25 qui ne sont pas nulles ont une valeur UN, la réciproque du courant de mesure J(Sq) peut être considérée comme équivalent à la résistance de la formation de terrain pour laquelle la synthétisation a lieu, à savoir de telle façon que Rq = l/J(Qq) et Ri une série de Rq, des jeux de paramètres de réponse peuvent être produits pour lesquels des jeux de facteurs d' étalonn age kq peuvent à leur tour être calculés. Il en résulte que le spécialiste interprétant ces données 15 bénéficie d'un procédé systématique par lequel la résistivité de la formation et la profondeur de l'invasion de la boue de forage dans le trou de sondage et autour de celui-ci sur la base de telles caractéristiques peuvent être aisément déterminées. La figure 8 illustre graphiquement la relation entre les numéros d'indexation croissants de la réponse synthétisée des assemblages d'outils et la profondeur de la réponse depuis l'intérieur de la formation de terrain environnante. Comme illustré, le facteur pseudo-géométrique G pour les divers assemblages synthétiques est étudié en fonction du diamètre d'invasion (Di) d'un filtrat de boue selon les modèles d'amorçage de tension indiqués au tableau V. Un simple profil d'étape sépare la zone 30 envahie Rxo de la zone non envahie Rt comme il est habituel dans les calculs de ce type. A cet égard, le facteur pseudo-géométrique G est défini de la manière habituelle par la relation G (Ra -Rt)/(Rxo - Rt) q q o les termes Rt et Rxo sont tels que définis précédemment. Le terme Ra dénote la résistivité apparente q des assemblages synthétisés après que les facteurs d'étalonnage ont été choisis. Dans cet exemple, dans 5 le but de construire les courbes de la figure 8, des multiplicateurs constants k ont été déterminés, de q manière que la résistivité apparente de chacque assemblage puisse être rendue égale à Rt dans la situation pour laquelle le rapport Rt/Rm est 100/1, avec un trou 10 de forage de 20, 3 cm et aucune invasion. Les courbes de la figure 8 montrent clairement la tendance des réponses produites par ordinateur à être sélectives en profondeur, avec les courbes à accroissement plus lent 150a,150b,...,150d caractérisant 15 les réponses les plus profondes et les courbes 150p,150q, z caractérisant les réponses peu profondes. Il est à noter que les courbes 150a et 150s se rapprochent des réponses des diagraphies respectives de Laterolog profonds (LLd) et de Laterolog peu profonds (LLS) couvertes 20 par des marques présentées par, et constituant, le produit d'outils focalisés mis en service par Schumlumber, Inc., Houston, Texas actuellement utilisés en commun dans la technique de diagraphie électrique. De telles courbes 150a et 150s avec les courbes additionnelles 25 de la figure 8 démontrent la supériorité de l'étendue de couverture et de la résolution du procédé de la présente invention. Les figures 9 à 27 sont des études des facteurs d' étalonnage en fonction d'une série de contrastes de différentes résistivités illustrant comment des jeux de tels facteurs sont uniques par rapport à une condition choisie de trou de forage. A cet égard, les courbes de la figure 9 couvrent un cas très simple qui sert à illustrer le 35 principe. Dans la figure, une série de courbes 170 est montrée; il s'agit du facteur d'-uaicnnaze logarithmique (kq)par rapport au logarithme du rapport Rt/Rm. La formation modèle sur laquelle sont basés les calculs consiste simplement en un trou de sondage 5 de 20, 3 cm dans un milieu par ailleurs homogène. C'està-dire qu'il n'y a pas d'invasion. En outre, la longueur de l'assemblage est de 9 m et consiste en M=73 ensembles d'électrodes. Les courbes 170 indiquent avec plus de 10 détail qu'il y a des jeux de facteurs d'étalonnage tels que les jeux 170a, 170b,... qui peuvent être associés de façon unique avec des paramètres intéressants de la formation. Pour ce cas illustratif très simple, il n'y a bien sur seulement qu'un paramètre en question 15 à savoir Rt qui peut être déterminé de façon unique si la résistivité de la boue Rm 'st connue. A noter également que l'ordonnée de droite est désignée avec des nombres entiers 1 à 37 représentant les outils de lecture la moins profonde et de lecture la plus 20 profonde de l'assemblage décrit précédemment qui ont été simulés synthétiquement en utilisant les modèles de réponse et d'initiation mentionnés au tableau V. Les facteurs d'étalonnage kq ont été calculés sur la base de la supposition suivante; si les réponses brutes des outils, c'est-à-dire les k sont donnés dans q un contraste de résistivité connue, leurs produits avec les k (qui correspondent au contraste de résistivité connue) fourniront des résistivités apparentes qui sont exactement égales au Rt correspondant dans un 30 trou de sondage de 20,3 cm. Il est clair alors que si les données brutes sont disponibles pour une condition inconnue, c'est-à-dire si les R sont disponibles à q partir d'une formation de terrain qui est considérée comme étant non envahie,on n'abesoin dans cette situation 35 que de chercher les courbes de la figure. 9 pour détermi- ner le jeu de facteurs de calibration qui fait en sorte que les résistivités apparentes résultantes soient toutes égales à une valeur constante d'o l'on déduit la condition de trou de forage, c'est-à-dire Rt. A cet égard, on voit qu'au-delà de contrastes de 1/1 (entrée du contraste de diagraphie égale à 0,0), les réponses des outils focalisés par ordinateur à lecture plus profonde, par exemple des indices d'outils de 10 à 37, ont des facteurs d'étalonnage qui ne dépendent pas fortement du contraste. Par ailleurs, les réponses des outils focalisés par ordinateur à lecture peu profonde indiquent la présence de facteurs d ' étalonnagequi dépendent fortement du contraste. Les figures 10 à 27 indiquent en outre 15 les aspects d'interprétation de la présente invention, par exemple dans des conditions de trou de forage impliquant des degrés croissants d'invasion de filtrat de boue. Si on désire tenir compte de la simple 20 invasion dans le modèle de la formation sur lequel les facteurs d'étalonnage doivent être basés, le nombre des paramètres est alors accru, c'est-à-dire que les facteurs d'étalonnage dépendront désormais non seulement de Rt mais également de Rxo, Di et Rm, 25 pour un diamètre fixe de trou de forage. Dans une telle situation, les courbes des figures 10 à 27 sont représentatives, car elles sont des courbes de facteurs d'étalonnage (logarithmique) par rapport au logarithme du rapport Rt/Rxo pour divers diamètres d'invasion 30 (profil en degrés). Il est à noter que le rapport de Rxo à Rm est fixé à 10/1-. De telles courbes représentent alors une base de données qui peut être utilisée de la manière indiquée ci-dessus pour en dériver les caractéristiques d'invasion de la formation. C'est-à-dire 35 que la condition du trou de forage est déduite en recher- chant les facteurs de ca-libraticn pour déterminer le jeu qui, lorsqu'il est multiplié par les réponses brutes correspondantes,donne des résistivités apparentes qui sont toutes égales à la même valeur de constante, les 5 facteurs d'étalonnage ayant été déterminés de manière à fournir la valeur de Rt dans une telle circonstance. Le procédé suppose donc de façon plus générale la production d'une base de données dans laquelle des jeux de facteurs d' étalonnage associés à différen10 tes conditions de trou de sondage (ressemblant au jeux a,170b,. de la figure 9 ou plus généralement des jeux du type contenu dans les figures 10 à 27) ont été calculés pour un assemblage particulier de diagraphie d'une certaine configuration par lequel des entrées 15 de matrice d'impédance doivent être obtenues sur le..DTD: terrain. A cet égard, l'assemblage d'électrodes (2M1) a été indiqué précédemment comme étant la configuration normale d'outils pour les opérations sur le terrain. On -suppose également que de tels jeux de facteurs d'étalonnage peuvent être emmagasinés dans une série de tableaux de facteurs indexés par des indices des outils 1 à 37 associés aux modèles uni taires en ti ro ir de lxl, 3x1, 5x1, lxl...71x1 et 73x1 respectivement, dans l'ordinateur numérique du processeur 25 contrôleur 17 de la figure 1. De tels tableaux sont indexes de manière que les entrées dans les tableaux soient annotées par configuration d'outils et modèle de réponse de courant et d'amorçage de tension, c'està-dire que la configuration des outils a été fixée 30 en accord avec la distribution des vecteurs de potentiel autour de l'assemblage de la manière mentionnée précédemment. On voit donc que de telles entrées de tableaux, à savoir les facteurs d'étalonnage doivent prendre en compte une variété de conditions différentes de 35 trou de forage et de formation,par exemple des conditions variables d'invasion ainsi que des conditions de noninvasion comme mentionné précédemment. ETAPES (ii) à (vi) Après qu'on a produit la série des tableaux 5 de facteurs d' étalonnage dans lesquels les jeux de facteurs d'étalo nnage sont adressés et annotés comme indiqué ci-dessus en accord avec le procédé de la présente invention, l'assemblage de diagraphie de la présente invention est situé dans le trou de forage 8 de la 10 figure 1. La prochaine diagraphie de la formation de la manière décrite précédemment a lieu, ce qui produit une série d'ensembles modifiés de matrices AZ, chacun étant associé à un segment prédéterminé de la formation de terrain en étude. Comme indiqué précédemment, chaque 15 ensemble de matrice comprend M x M entrées, M étant l'indice de numérotation interne de l'ensemble d'électrodes de courant de l'assemblage. Dans l'exemple d'un assemblage (2M-i) de 145 électrodes, M est 73. ETAPES (vii) à (ix) 2C Après détermination des ensembles modifiés de matrice, chacun est inversé pour former un ensemble -1 modifié réciproque de matrice AZ 1 en utilisant les techniques classiques d'inversion de matrice. Par exemple, l'expérience a indiqué que les techniques classiques 25 bien connues de résolution des équations telles que l'élimination Gaussienne sont suffisantes pour résoudre de façon satisfaisante le problème de l'inversion. A cet égard, voir par exemple "Linear Algebra And Its Applications", Gilbert Strang, Academic Press, 1976. 30 Ensuite, on crée une série de réponses produites par ordinateur en utilisant les matrices -1 obtenues modifiées réciproques AZ 1 normalisées sur des modèles particuliers de réponse de courant et de tension ressemblant à ceux utilisés en produisant les tableaux 35 de facteurs de l'étape (i) comme décrit précédemment. En d'autres termes, en utilisant les données fournies de l'assemblage sur le terrain, on produit une série de valeurs de réponse R qui sontassociées à un segment q de la formation égal à M stations de diagraphie le 5 long du trou de forage lorsque la diagraphie de la formation a eu lieu. Finalement, les valeurs R obtenues q sont comparées aux tableaux de facteurs contenant les jeux de facteurs d' étalonnage. La recherche est terminée (ce qui indique les paramètres intéressants de 10 la formation) lorsque certains critères de meilleure sélection sont satisfaits. Par exemple pour un jeu donné, de valeursRq, la série des tableaux contenant les facteurs d ' étalonnage sont recherchés séquentiellement jusa'à ce que la production d'un jeu porédéterminé de 15 facteurs d' étalonnage à l'intérieur d'un tableau particulier de facteurs fournisse le résultat désiré, c'est-à-dire lorsque les produits des valeurs R et le q jeu sélectionné de facteurs d'étalonnage sont constants et égaux pour toutes les 37 réponses d'outils focalisées 20 par ordinateur. Si la sortie de la recherche a comme résultat un format de courbe, les valeurs Ra étant conq sidérées comme les ordonnées et le numéro d'indice des outils 1 à 37 comme l'abscisse, le meilleur jeu de données en accord avec le procédé de la présente 25 invention a lieu lorsqu'un jeu donné de facteurs d'étalonnage multiplié par rapport aux valeurs Rq (=Raq) s'approche d'une ligne droite horizontale. Il est à noter cependant que les tableaux de facteurs d'étalonnage doivent être augmentés par l'addition d'autres 30 conditions de la formation et/ou du trou de forage. A cet égard, on peut annoter des tableaux additionnels par d'autres variables de trous de forage telles qcue montrées aux figures 10 à 27. La figure 10 par exemple illustre qu'un 35 tableau de facteurs peut être construit facilement pour une formation de terrain ayant des lits épais avec une résistivité de boue de 0,1 ohm par mètre (Rm) et une résistivité de la zone envahie de 1 ohm par mètre (Rxo). Les facteurs d'étalornaze sont considérés comme l'ordonnée et le contraste de résistivité de la formation et de la zone envahie est considéré comme l'abscisse pour un diamètre d'invasion de 20,3 cm (c'està-dire aucune invasion). Les figures 11 à 27 illustrent comment 10 on peut construire même plus de tableaux de facteurs pour le même trou de forage mais pour différentes conditions de la formation par rapport à ceux de la figure , des variations de diamètre d'invasion étant cependant prises en considération, c'est-à-dire dans le cas o 15 ces dernières augmentent progressivement par exemple de 25,4 cm à la figure 11 à 63 cm à la figure 15, et finalement à 300 cm à la figure 27. Il est à noter que pour une invasion allant de modérée à profonde, les ourbes indiquent une forte dépendance des facteurs 20 d 'é alonnage illustrés par rapport à l'amplitud. du,ontraste de résitivité de la formation et de 1zone envahie, c'est-àdire Rt/Rxo. Egalemient, comme la distance latérale s'accroit dans les figures, les réponses des outils focalisés par ordinateur à lecture 25 plus profonde jouent un rôle progressivement plus important et ajoutent largement au nombre des facteurs d'étalonnage disponibles dans l'étape de production et de là augmentent de façon importante la fiabilité des résultats finals. L'essentiel de toute analyse 30 est le fait que les jeux particuliers de facteurs d'étalonnage pour différents assemblages des outils peuvent être liés en combinaison de façon unique à des paramètres de la formation qui sont intéressants et ainsipeuvent être utilisés pour fournir une indication precise surprenante de tels parametres. Bien sûr, les facteurs d'étalonnace sont anciens dans la technique de diagraphie électrique. Par exemple, pour les outils classiques focalisés qui 5 fournissent des courbes de réponse tels que 150a et s de la figure 8, o les facteurs d'étalonnace sont choisis de manière que leurs réponses soient égales approximativement à la résistivité vraie de la formation (Rt) dans des formations non envahies ayant des contras10 tes de résistivités de la formation et de la boue dans la gamme de 10/1 à 100/1, normalisés sur un trou de forage de 20,3 cm, il est à noter que de tels outils sont câblés en dur. De là, une fois que leurs facteurs d' étaionnage sont déterminés, ils restent fixes à 15 moins qu'une modification de la conception nécessite une modification, comme lorsqu'un changement a lieu par exemple dans l'espacement ou la dimension des électrodes de tels outils. Mais comme les caractéristiques de fonctionnement en accord avec la présente invention 20 ne dépendent pas de la conception mais sont facilement calculées e:.cilisant les étapes du procédé de la présente invention comme mentionné ci-dessus, ils peuvent être plus facilement associés à des conditions différentes de trous de forage et de là offrent une solution 25 pour déterminer les paramètres intéressants de ia formation. Bien que des formes de réalisation préférées de la présente invention aient été décrites en détail, il est à noter que cette invention n'est pas limitée 30 à ces formes de réalisation car beaucoup de modifications apparaîtront facilement aux spécialistes de cette technique et ainsi l'invention doit bénéficier de l'interprétation la plus large possible selon les termes des revendications suivantes. REVENDICATIONS 1. Procédé pour faire la synthèse des caractéristiques de réponse vraies d'une combinaison de différents outils de diagraphie électrique centrés 5 sur un trou dans des conditions difficiles variées de trous de sondage, consistant (1) à déterminer des valeurs d'impédance d'une formation de terrain traversée par un trou de forage rempli d'une boue de forage d'une résistivité (Rm) et (2) à manipuler sélectivement 10 les valeurs d'impédance en tant qu'entrées d'impédance d'une matrice de manière à faire la synthèse des opérations des différents outils centrés sur le trou sur un incrément de profondeur associé, le tout avec une précision surprenante, chacune des matrices d'impédance 15 étant associée à un ensemble de matrices indexées sur un incrément parmi une série d'incréments finis chevauchants de balayage de profondeur de ladite formation mesurés le long du trou de forage, chaque incrément de balayage dépendant de la longueur L de l'assemblage 20 d'électrodes de man -re à définir des marqueurs de profondeur peu profonds et profonds en étant indexés centralement sur la profondeur du trou de forage d'un ensemble d'électrodes centrales de l'assemblage au moment du prélèvement des données, de manière à fournir 25 une indication vraie de la résistivité (Rt) de la formation même si celle-ci est espacée du trou de forage par une zone envahie de résistivité (Rxo) d'une étendue latérale inconnue en raison de- l'invasion du filtrat de la boue de forage, le procédé étant caractérisé 30 en ce qu'il comprend les étapes consistant: (1) à calibrer un assemblage de M ensembles d'électrodes centrés sur le trou ayant un espacement incrémentiel d'électrodes égal "a" pour obtenir des jeux de facteurs d'étaionn age normalisés sur des 35 modèles connus de tension dans une zone de résistivité connue de répornse, ies ensembles d'électrodes ayant un indice de numerotation connu d'ordre interne, les jeux de facteurs d'6talonnaae étant chacun adressables en fonction des conditions du trou de forage y compris les conditions difficiles de trous de forage ainsi que par un assemblage particulier du type synthétique focalisé par ordinateur; (ii) à positionner un assemblage de champ dans le trou de forage, cet assemblage comprenant un ensemble d'électrodes centrales comportant une électrode' 10 de courant et une électrode de potentiel et une série d'ensembles d'électrodes additionnelles dont chacun comprend seulement une électrode de mesure de potentiel, la profondeur absolue d'au moins un ensemble d'électrodes étant connue en continu par rapport à un niveau de 15 données de profondeur prédéterminé mesuré depuis la surface du sol; (iii) à déplacer continuellement l'assemblage le long du trou de forage et injecter en continu du courant d'une-valeur connue depuis l'électrode de courant 20 de l'ensemble central, à..:,-rs la boue de forage de résistivité (Rm) et la zone envahie de résistivité (Rxo) et de là dans la formation de résistivité Rt; (iv) pendant l'injection du courant, à mesurer le potentiel absolu de chaque électrode de 25 potentiel ainsi que les différences de potentiel entre des paires adjacentes d'électrodes de potentiel d'une manière rapide, les mesures de potentiel pouvant être indexées des stations communes de diagraphie d'une série de stations de diagraphie équidistantes situées 30 le long du trou de forage et ayant une distance incrémentielle d'espacement de "a", la valeur "a" étant la distance entre des ensembles voisins d'électrodes. (v) à calculer la valeur d'impédance à partir des valeurs mesurées de potentiel absoiu et différentiel et leurs courants associés d'ir.jection, chaque valeur étant indexée sur lesdits numéros d'indexation connus internes d'électrodes actives de courant et de potentiel utilisées dans les mesures; (vi) à *réindexer les valeurs d'impédance en entréesd'impédance d'une série d'ensembles chevauchant de matrice, chaque ensemble étant associé à un segment prédéterminé de ladite formation égal en étendue vertica10 le, sur M stations de diagraphie et comprenant M x M entrées d'impédance, la valeur 2M-2 étant le plus grand nombre de l'indice de numérotation des ensembles d'électrodes comprenant ledit assemblage et le rapport du nombre des entrées d'impédance de différence aux 15 entrées absolues étant d'environ M-1:1; (vii) à inverserchaque ensemble de matrice modifiée; (viii)à produire des paramètres de réponse focalisées par ordinateur en utilisant des modèles 20 de tension de l'étape (i); (ix)à rechercher des jeux dt facteurs a'étalonnage de l'étape (i) de sorte que la condition difficile de trou de forage en est déduite même en présence de constrastes élevés entre la résistivité 25 vraie et la résistivité de la boue et nonobstant le fait que des jeux synthétiques de modèles de potentiel ont été utilisés comme initiateurs des paramètres de réponse focalisés par ordinateur et produits subséquemment. 2. Procédé selon la revendication 1, caractérisé en ce que l'assemblage de M ensembles d'électrodes comprend des ensembles El, E2..EM et le numéro d'indexation interne commence à l'ensemble E1 et se termine à l'ensemble EM, la vale'ur E1 étant l'ensemble le moins 35 profond de l'assemblage et EM étant l'ensemble le plus profond par raptrz a t a s-fa ^_r-, -t!es ensembles E1 et EM comportant des électrodes de courant. 3. Procédé selon la revendication 2, caractérisé en ce que la longueur L de l'assemblage comprenant 5 des ensembles d'électrodes E1, E2...E2MI est la distance verticale séparant l'ensemble d'électrodes le moins profond E1 et l'ensemble le plus profond E2M_1 le long du trou de forage. 4. Procédé selon la revendication 2, caracté10 risé en ce que la longueur L de l'ensemble définit également une zone de formation d'une extension verticale L/2 qui est liée par la profondeur aux entrées calculées d'impédance de chaque matrice M x M de l'étape (vi). 5. Procédé selon la revendication 1, caractéz 15 risé en ce que les entrées de chaque matrice- M x M comprennent M colonnes et M rangées d'entréesd'impédance, les colonnes pouvant être constituées d'entrées d'impédance absolue Z et d'entrées d'impédance de différence modifi ée AZ et les entrées d'impédance de d'fférence 20 modifiée A Z de chaque colonne étant chacune associées à une profondeur commune à l'intérieur de la -j nation de terrain, désignée par un argument de profondeur entre parenthèses, associé à la profondeur instantanée de l'électrode de courant centrale à mesure que les 25 données sont prélevées, chaque entrée dans la première rangée étant une impédance absolue dont l'indexation interne est basée sur le même modèle que les rangées restantes, ou sinon, est une valeur d'impédance absolue associée à une profondeur fixe commune de l'ensemble 30 d'électrodes de courant central, désigné également par un argument de profondeur entre parenthèses et résultant du principe de la réciprocité. 6. Procédé selon la revendication 5, caractérisé en ce que chacune des colonnes de chaque matrice 35 M x M est associée à une électrode de potentiel dont l'indice de numerotation interne varie dans un ordre donné de rangée à rangée le long de chaque colonne. 7. Procédé selon la revendication 5, caractérisé en ce que chacune des rangées de chaque matrice M x M est associée à une électrode de potentiel dont l'indice de numérotation interne change selon un ordre connu de colonne à colonne le long de chaque rangée. 8. Procédé selon la revendication 5, caractérisé en ce que M a la valeur de 5. 9. Procédé selon la revendication 8, caractérisé en ce que les entrées de la matrice 5 x 5 sont constituées d'entrées d'impédance absolue Z et d'entrées d'impédance de différence modifiée WZ de sorte que l'indice le plus à gauche de numérotation 15 interne va cycliquement de 1 à 9, en décroissant de gauche à droite le long de chaque rangée et en augmentant de haut en bas le long de chaque colonne, en commençant et en finissant à 5 dans les entrées à gauche en haut et à droite en bas, l'indice de numérotation 20 interne le plus à droite étant fixé à 5 à toutes lee ent ées. 10. Procédé selon la revendication 9, caractérisé en ce que les entrées d'impédance de la matrice 5 x 5 sont chacune identifiées par un premier sous25 indice associé à l'électrode de potentiel de numérotation interne d'o la valeur de potentiel absolu ou de différence par laquelle l'entrée d'impédance a été calculée était associée à l'origine, suivi par une virgule et ensuite un second sous-indice associé à l'électrode 30 de courant central d'ordre interne d'o est parti le courant de commande associé au calcul d'impédance. 11. Procédé selon la revendication 10, caractérisé en ce que les sousindices des entrées d'impédance absolue Z et les entrées d'impédance de 35 différence AZ de la matrice 5 x 5 pour des ensembles de numérotation interne,... sont selon 'ordre suivant en forme matricielle, la première forme comprenant des entrées d'impédance absolue et de difference modifiée avec un modèle uniforme d'indexation interne et la seconde forme utilisant le principe de la réciprocité permettant à des impédances absolues dérivées d'une profondeur fixe commune d'électrodes de courant central d'être utilisées en formant la premaLère rangée: (première forme) Z5,5(d5) Z4,5(d6) Z3,5(d7) AZ6,5(d5) AZ5,5(d6) AZ41,5(d7) AZ(Sdl)=AZ7, 5(d5) AZ6,5(d6) AZ5,5(d7) AZ8,5(d5) AZ7,5(d6) AZ6,5(d7) A29,5(d5) AZ8, 5(d6) AZ7,5(d7) Z2,5(d8) AZ3, 5(d8) AZ4, 5(d8) aZ5,5(d8) aZ6,5(d8) Z1, 5(d9) AZ2,5(d9) LZ3,5(d9) AZ4 5(d9) AZ5,5(d9) (seconde forme) Z5,5(d5) AZ6,5(d5) -AZ(Sdl)=Az7,5(d5) AZ8,5(d5) AZ9,5(d5) 26,5(d5) Z7,5(d5) AZ5, 5(d6) AZ4,5(d7) AZ6,5(d6) AZ5,5(d7) AZ7,5(d6) AZ6,5(d7) AZ8,5(d6) AZ7, 5(d7) Z8,5(d5) AZ3; (d8) AZ4,5(d8) AZ5,5(d8) AZ6,5(d8) Z9,5(d5) LZ2,5(d9) AZ3,5(d9) AZ4,5(d9) AZ5,5(d9) o d5, d6, d7, d8 et d9 sont les profondeurs de l'électrode de courant central activée à mesure que les données. sont prélevées. 12. Procédé selon la revendication 5, caractérisé en ce que les entrées de la matrice M x M sont 35 constituées par des entrées d'impédance absolue Z et des entrées d'-!-édanrce diffrentieile modifiée AZ de telle manière que l'indice de numérotation interne le plus à gauche s'étende cycliquement de 1 à 2M-1, en décroissant de gauche à droite le long de chaque ran5 gée et en augmentant de haut en bas le long de chaque colonne, en débutant et en terminant à M dans les entrées en haut à gauche et en bas à droite, l'indice de numnérotation interne le plus à droite étant fixé à M pour toutes les entrées. 13. Procédé selon la revendication 12, caractérisé en ce que les entrées d'impédance Z et WZ sont chacune identifiées par un premier sous-indice associé à l'électrode de potentiel à numérotation interne d'o la valeur de potentiel absolue ou différentielle 15 par laquelle l'entrée d'impédance a été calculée était associée à l'origine, le tout suivi d'une virgule, et ensuite un second sous-indice associé à l'électrode de courant central à numérotation interne o à pris son origine le courant de cornande pour le calcul d'impé20 dance. 14. Procédé selon la revendication 1, caractérisé en ce que l'étape (i) est caractérisée en outre par la production de jeu de facteurs d'étalonnage pour différentes combinaisons d'assemblagesen supposant 25 une configuration d'assemblage comprenant M ensembles d'électrodes également espacés commandés par une série de modèles indépendants de tension en fonction des constrastes différents de résistivité entre la boue et la formation. 15. Procédé selon la revendication 14, caractérisé par l'étape additionnelle consistant à emmagasiner les facteurs d' étalonnage produits en fonction des contrastes de résistivité entre la boue et la formation et des combinaisons différentes d'assem35 blages. 16. Procédé selon la revendication!5, caractérisé par l'étape additionnelle consistant à emmagasiner les facteurs d 'étalonnace produits en fonction des contrastes de résistivité entre la boue et la zone envahie, des constrastes de résistivité entre la formation et la zone envahie, des diamètres de la zone envahie et des diamètres du trou de forage et desdites combinaisons différentes d'assemblages. 17. Procédé selon la revendication 14, 10 caractérisé en ce que l'étape (ix) comporte la multiplication de chaque jeu de facteurs d'étalonnage de l'étape (i) sur une base de un à un par des paramètres de réponse individuels de l'étape (vii) normalisés sur des modèles semblables d'amorçage de tension de 15 l'étape (i) pour former une série de valeursnormalisées et la sélection d'une série particulière qui contient des valeurs d'amplitude essentiellement constantes en tant que meilleur résultat. 18. Procédé selon la revendication 1 compre20 nant une étape consistant à fournir les impédances d'entrée pour une série d'ensembles de matrices, chaque ensemble de matrices étant indexé sur un incrément parmi une série d'incréments finis chevauchant le balayage de profondeur de la formation de terrain, mesurés 25 le long du trou de forage, chaque incrément de balayage dépendant de la longueur d'assemblage L de l'assemblage des électrodes de manière à définir des marqueurs de profondeur peu profonds et profonds et étant indexés centralement sur la profondeur dans le trou de forage 30 d'une électrode unique émettrice de courant d'un ensemble d'électrodes centrales de l'assemblage au moment du prélèvement des données, de manière à fournir une indication vraie de la résistivité de la formation (Rt) même lorsque la formation est séparée du trou de forage 35 par une zone envahie de résistivité (Rxo) d'étendue latérale inconnue à cause de invasion du filtrat de boue de forage, cette étape étant caractérisée par les sous-étapes consistant à: (i) attacher rigidement un assemblage d'en5 sembles d'électrodes centrés sur le trou d'espacement d'électrodes incrémentielles éga 1 "a" à un câble de diagraphie, les ensembles d'électrodes ayant un indice connu de numérotation d'ordre interne, un ensemble d'électrodes centrales comprenant une électrode de 10 courant et une électrode de potentiel ayant un indice de numéro M et un ensemble d'électrodes terminal ayant un indice de numéro de 2M-1 lorsque M=2N+1, tous les ensembles d'électrodes sauf l'ensemble central comprenant une électrode de mesure de potentiel seulement; 15 (ii) positionner l'assemblage d'ensembles d'électrodes dans le trou de forage, la profondeur absolue d'au moins un ensemble d'électrodes étant connue en continu par rapport à un niveau de données de profondeur prédéterminé mesuré depuis la surface du sol; (iii) déplacer en continu l'assemblage le long du trou de forage et injecter en continu du courant depuis l'électrode de courant de l'ensemble d'électrodes central d'une valeur connue à travers la boue de forage de résistivité (Rm) et la zone envahie 25 de résistivité (Rxo) et de là dans la formation de résistivité Rt; (iv) mesurer le potentiel absolu à chaque électrode de potentiel des ensembles d'électrodes ainsi que les différences de potentiel entre des paires voisi30 nes d'électrodes de potentiel à une série de stations de diagraphie équidistantes situées le long du trou de forage, et ayant une distance incrémentielle d'espacement de "a", cette valeur "a" étant la distance entre les ensembles d'électrodes. (v) calculer les valeurs d'impédance à partir des valeurs mesurées de potentiel absolu et différentiel et leur courant d'injection associé, chaque valeur étant indexée sur les numéros connus d'indexation interne des électrodes actives de courant et de potentiel utilisés dans les mesures; (vi) réindexer les valeurs; d'impédance en entrées d'impédance d'une série d'ensembles chevauchants de matrices, chaque ensemble étant associé à un segment prédéterminé de ladite formation égal en étendue verticale, sur M stations de diagraphie et comprenant M x M entrées d'impédance, 2M-1 étant le numéro le plus grand de l'indice de numérotation de l'assemblage des ensembles d'électrodes et le rapport entre le nombre des entrées d'impédance de différence 1-5 et celui des entrées absolues étant M-1:1, de sorte que chaque matrice M x M puisse être utilisée de façon surprenante pour déduire la résistivité vraie de ladite formation même en présence de contrastes élevés entre la résistivité vraie et la résistivité de la boue et 20 nonobstant le fait de savoir si oui ou non des jeux synthétiques de potentiel ou des jeux de valeurs de courant sont utilisés subséquemment comme initiateurs de réponses produites subséquemment d'assemblagesd'outils synthétiques focalisés par ordinateur. 19. Procédé selon la revendication 18, caractérisé en ce que l'assemblage des M ensembles d'électrodes comprend des ensembles E1,E2,...E2M_1 et que l'ordre d'indexation interne commence à l'ensemble E1, comporte un ensemble central EM et se termine à l'en30 semble E2M1, la valeur E1 étant l'ensemble le moins profond de l'assemblage et E2M étant l'ensemble le plus profond par rapport à la surfacedu sol. 20. Procédé selon la revendication 19, caractérisé en ce que la longueur L de l'assemblage 35 comprenant les ensembles d'électrodes El, E2. EM *. E2M! est ïa distance verticale séparant l'ensemble d'électrodes le moins profond E. et l'ensemble le plus profond E2M-1 le long du trou de forage. 21. Procédé selon la revendication 20, caractérisé en ce que la longueur L de l'ensemble définit également une zone de formation, d'étendue verticale L/2, qui est liée par la profondeur aux entrées calculées d'impédance de chaque matrice M x M de l'étape (vi). 22. Procédé selon la revendication 18, 10 caractérisé en ce que les entrées de chaque matrice M x M comprennent M colonnes et M rangées d'entrées d'impédance, les colonnes étant constituées par des entrées d'impédance absolue Z et des entrées d'impédance différentielle modifiées ALZ et ces dernières étant 15 chacune associéesdans chaque colonne avec une profondeur commune à l'intérieur de la formation de terrain, désignée par un argument de profondeur entre parenthèses, associé avec la profondeur instantanée de l'électrode de courant central à mesure que les données sont préle20 vées, et chaque entrée dans la première rangée pouvant consister soit en une valeur d'impédance différentielle modifiée basée sur le même modèle que les rangées restantes soit en une valeur d'impédance absolue associée avec une profondeur commune d'ensemble d'électrodes 25 de courant central désignée également par un argument de profondeur entre parenthèses et résultant du principe de la réciprocité. 23. Procédé selon la revendication 22, caractérisé en ce que chacune des colonnes de chaque 30 matrice M x M est associée avec des électrodes de potentiel dont l'indice de numérotation interne change dans un ordre connu de rangée à rangée le long de chaque colonne. 24. Procédé selon la revendication 22, 35 caractérisé en ce que chacune des rangées de chaque matrice M x M est associée à des électrodes de potentiel dont l'indice de numérotation interne change séquentiellement dans un ordre connu de colonne à colonne le long de chaque rangée. 25. Procédé selon la revendication 22, caractérisé en ce que M est égal à 5. 26. Procédé selon la revendication 25, caractérisé en ce que les entrées de la matrice 5 x 5 sont constituées par des entrées d'impédance absolue 10 Z et des entrées d'impédance différentielle modifiée AZ de telle manière que l'indice d'alimentation interne le plus à gauche s'étende cycliquement de 1 à 9, en décroissant de gauche à droite le long de chaque rangée et en augmentant de haut en bas le long de chaque colonne, 15 en débutant et en se terminant à 5 dans les entrées en haut à gauche et en bas à droite, l'indice de numérotation interne le plus à droite étant fixé à 5 pour toutes les entrées. 27. Procédé selon la revendication 26, 20 caractérisé en ce que les entrées d'impédance de la matrice 5 x 5 sont chacune identifiées par un premier sous-indice associé à l'électrode de potentiel de numérotation interne d'o la valeur de potentiel absolu ou différentiel par laquelle l'entrée d'impédance a 25 été calculée était associée à l'origine, suivi par une virgule et ensuite un second sous-indice associé à l'électrode de courant central de numérotation interne d'o est parti le courant de commande associé au calcul d'impédance. 28. Procédé selon la revendication 27, caractérisé en ce que les sous-indices des entrées d'impédance absolue Z et les entrées d'impédance modifiée àZ de la matrice 5 x 5 pour des ensembles à numérotation interne E1,E2,...,E9 ont l'ordre suivant en forme matri35 cielle, la première forme comprenant des entrées d'impé- dance absolue et différentielle modifiées avec un modèle uniforme d'indexage interne, et la seconde forme utilisant le principe de la réciprocité permettantà des impédances absolues dérivées d'une profondeur fixe commune d'électrodes de courant central d'être utilisées pour former la première rangée, c'est-à-dire: (première forme) Z5, 5(d5) AZ6,5(d5) aZ_(Sdl)=gZ7,5(d5) Z28,5(d5) AZ9,5(d5) Z4,5(d6) AZ5,5(d6) AZ6,5(d6) LZ7, 5(d6) AZ2,5(d6) Z3,5(d7) AZ4,5(d7) AZ5,5(d7) AZ6,5(d7) AZ7,5(d7) Z2,5(d8) LZ3, 5(d8) AZ4,5(d8) LZ5,5(d8) AZ6,5(d8) Z2,5(d9) AZ2,5(d9) AZ3,5(d9) LZ4, 5(d9) LZ5,5(d9) (deuxième forme) Z5,5(d5) aZ6,5(d5) AZ(Sd1)= Z7,5 (d5) Z8, 5 (d5) AZ9,5(d5) Z6,5(d5) AZ5,5(d6) AZ6,5(d6) AZ7,5(d6) AZ8,5(d6) Z7,5(d5) àZ4,5(d7) AZ5,5(d7) AZ6,5(d7) AZ7,5(d7) 28,5(d5) LZ3,5(d8) LZ4,5(d8) A25, 5(d8) 2Z6,5(d8) Z-,5(d5) LZ",5(d9) AZ3,5(d9) AZ4,5(d9) AZ5,5(d9 o: d5, d d7, d8 et d9 sont les profondeurs de l'électrode de courant central activée à mesure que les données 30 sont prélevées. 29. Procédé selon la revendication 22, caractérisé en ce que les colonnes de la matrice M x M sont constituées par des entrées d'impédance absolue Z et des entrées d'impédance différentielle modifiée 35 AZ de telle manière que l'indice de numérotation interne le plus à gauche s'étende cycliquement de i à 2M-1 en décroissant de gauche à droite le long de chaque rangée et en augmentant de haut en bas le long de chaque colonne, en débutant et en se terminant à M dans les entrées de gauche en haut et de droite en bas, l'indice de numérotation interne le plus à droite étant fixé à M pour toutes les entrées. 30. Procédé selon la revendication 29, caractérisé en ce que les entrées d'impédance Z et AZ sont chacune identifiées par un premier sous-indice associé à l'électrode de potentiel de numérotation interne d'o la valeur de potentiel absolue'ou différentielle par laquelle l'entrée d'impédance a été calculée était associée à l'origine, suivie par une virgule 15 et ensuite un second sous-indice associé à l'électrode de courant central de numérotation interne d'o est parti le courant de commande pour le calcul des impédances.