Brevet FR2920876A1 — publié le 20090313

La présente invention concerne le domaine pétrophysique et plus précisément la détermination de la saturation en saumure au sein d'un milieu poreux, tel qu'une fonnation géologique par exemple. La connaissance de la saturation en saumure permet de calculer les quantités d'hydrocarbures en place dans le cas d'un réservoir pétrolier, ou la quantité de CO2 stockée 15 dans le cas d'un stockage de CO2. Il est donc très important pour le spécialiste en charge de l'exploitation d'une formation souterraine, dans le cadre de production d'hydrocarbure ou de stockage de gaz acides par exemple, de connaître la saturation en saumure de cette formation. État de la technique 20 Il est connu de déterminer la saturation d'une formation géologique en réalisant des mesures électriques au sein de puits, en cours ou non de forage. Ces mesures au sein d'un puit sont appelées des diagraphies (� well logs �). Ces diagraphies électriques mesurent la résistivité électrique du milieu. Elles doivent être analysées de façon à déduire les valeurs de saturation au sein de la formation. 25 Une méthode classique pour réaliser une telle interprétation des diagraphies de résistivité consiste à utiliser les lois d'Archie, et notamment la relation suivante : RI = R = S, -n 0 avec : RI : indice de résistivité Sw saturation en saumure du milieu Rt : résistivité électrique du milieu Ro : résistivité électrique du milieu saturé à 100% par la saumure (Sw =1). n un nombre entier La résistivité électrique du milieu Rt est mesurée par diagraphies électriques, tandis que la résistivité électrique du milieu saturé à 100% par la saumure Ro, est déterminée à partir de mesures effectuées sur des carottes. Une étape importante de l'interprétation des diagraphies de résistivité (Rt) en terme de 10 saturation en saumure (SS,), consiste donc à définir le paramètre n. Pour ce faire, il est connu de mesurer d'une part l'indice de résistivité et la saturation en saumure sur des échantillons issus de la formation géologique. Ces mesures sont souvent réalisées dans les conditions de pression et de température rencontrées au sein de la formation. Ainsi, en mesurant plusieurs couples (RI, SH,), on peut estimer le paramètre n de la 15 relation d'Archie. La mesure de l'indice de résistivité d'un milieu poreux constitue l'étape expérimentale la plus difficile dans ce processus d'interprétation de la réponse électrique des milieux poreux. Une technique disponible pour mesurer l'indice de résistivité et la saturation en saumure 20 sur un grand nombre d'échantillons, consiste à placer ces échantillons sur une plaque poreuse semi-perméable à l'eau, l'ensemble étant disposé dans une enceinte pouvant être amenée en pression. Lorsque la pression P dans l'enceinte est augmentée progressivement, la saturation des échantillons diminue progressivement. A chaque palier de pression, on retire les échantillons de l'enceinte et on détermine la saturation par pesée. La résistivité est mesurée en 25 appliquant deux plaques métalliques connectées à un appareil de mesure aux deux faces de l'échantillon. Cependant, il est connu que la détermination de la saturation par cette technique est entachée d'erreurs dues à la perte de grain lors de la manipulation, et que la mesure de résistivité est entachée d'erreurs dues aux résistances de contact entre les plaques et le milieu 30 poreux. De plus, pour les échantillons peu perméables, une plaque poreuse de très faible perméabilité est nécessaire, nécessitant des durées d'expériences longues (de l'ordre de 2 mois au moins). La méthode selon l'invention Ainsi, l'objet de l'invention concerne une méthode alternative pour mesurer une saturation en un fluide conducteur et une résistivité d'un milieu poreux, initialement saturé par le fluide. La méthode comporte les étapes suivantes : a)- on extrait au moins un premier échantillon solide du milieu ; b)- on place le premier échantillon dans une cellule de centrifugation, à côté d'un second échantillon adapté à limiter un pied capillaire dans le premier échantillon, les échantillons étant disposés l'un à côté de l'autre de façon à optimiser un contact capillaire entre eux ; c)- on désature partiellement les deux échantillons en les soumettant à une centrifugation 15 selon une vitesse de rotation fixée et un temps de rotation suffisant pour permettre une désaturation du premier échantillons ; d)- on mesure la résistivité du premier échantillon, en le plaçant dans une cellule de mesure de résistivité, adaptée à injecter un courant par l'intermédiaire d'au moins un premier couple d'électrodes, et à mesurer un potentiel entre au moins un second 20 couple d'électrodes, ces électrodes étant disposées radialement autours du premier échantillon ; et e)- on détermine la saturation du premier échantillon en mesurant une aimantation nucléaire de l'échantillon à l'aide d'un appareil RMN. 25 Selon l'invention, le second échantillon peut être extrait directement du milieu poreux duquel est issu le premier échantillon. Le contact capillaire peut être optimisé en intercalant entre les deux échantillons une membrane mouillable au fluide. On peut utiliser la méthode pour estimer une relation permettant de déduire la saturation en fonction de la résistivité. Pour ce faire, on réitère les étapes b) à e), en modifiant la vitesse de rotation de façon à obtenir plusieurs couples de valeurs de résistivité et de saturation pour le premier échantillon. Puis on estime la relation à partir de ces couples. Enfin, on peut utiliser cette relation pour mesurer la saturation en un fluide conducteur d'une formation souterraine. La méthode comporte alors les étapes suivantes : - on extrait des échantillons de la formation que l'on sature avec le fluide ; - on détermine une relation permettant de déduire la saturation en fonction de la résistivité pour chacun des échantillons, à l'aide de la méthode selon la méthode décrite précédemment ; - on sélectionne des échantillons de façon à ce que chacun d'entre eux ait une relation différente des autres échantillons ; on mesure l'indice de résistivité et la saturation en fluide de chacun des échantillons sélectionnés, à l'aide d'une méthode adaptée à fonctionner selon des conditions de température et de pression propres à la formation souterraine ; on réalise au moins une diagraphie de résistivité eléctrique à travers la formation ; - on en déduit la saturation en fluide dans la formation en appliquant la relation. D'autres caractéristiques et avantages de la méthode selon l'invention, apparaîtront à la lecture de la description ci-après d'exemples non limitatifs de réalisations, en se référant aux figures annexées et décrites ci-après. Présentation succincte des figures - la figure 1 schématise les étapes de la méthode permettant de mesurer rapidement la saturation en un fluide conducteur et l'indice de résistivité d'un milieu poreux. la figure 2 schématise la mise en oeuvre de la désaturation par centrifugation, selon 25 l'invention. - la figure 3 illustre schématiquement un profil de saturation Si, dans les échantillons El et E2 en contact capillaire, dans le cadre d'une centrifugation dans un dispositif selon la figure 2.20 la figure 4 montre des profils de saturation S , dans les échantillons El et E2, pour différentes vitesses de rotation de centrifugation. les figures 5A et 5B montrent un dispositif de mesure radiale de la résistivité. - la figure 6 montre un exemple de résultat de mesure sur deux Grès de Fontainebleau de 5 grande perméabilité. - la figure 7 monte des mesures d'indice de résistivité et de saturation sur un calcaire à structure bimodale en drainage. la figure 8 monte des mesures d'indice de résistivité et de saturation sur un calcaire à structure bimodale en drainage. Description détaillée de la méthode La méthode selon l'invention permet de mesurer rapidement la saturation en un fluide conducteur et l'indice de résistivité d'un milieu poreux. Habituellement, ces mesures sont réalisées dans le but d'obtenir des mesures 15 extrêmement fiables, peu importe le temps nécessaires pour effectuer ces mesures. La méthode selon l'invention, au contraire, propose une méthode permettant d'obtenir rapidement ces mesures sans perte de précision. La méthode de mesure est décrite dans le cadre d'un milieu contenant les deux fluides suivant : saumure et air. Cependant, la méthode est également applicable au milieu contenant 20 de la saumure et de l'huile, comme un réservoir pétrolier par exemple, et même à tout milieu contenant au moins un fluide conducteur. Pour un milieu contenant de la saumure (fluide conducteur) et de l'air, et initialement saturé par la saumure, on applique la procédure suivante, décrite sur la figure 1 : - on extrait au moins un échantillon solide du milieu (en général une carotte) - EXT 25 - on désature partiellement l'échantillon par centrifugation - DES on mesure l'indice de résistivité par mesure radiale - RI on mesure la saturation à l'aide d'un appareil RMN - Sw 10 - Désaturation partielle de l'échantillon Le but de cette étape est d'obtenir une saturation moyenne Sw dans l'échantillon. La technique utilisée consiste à drainer l'échantillon par centrifugation, selon un temps fixé, et selon une vitesse de rotation co fixée, de façon à désaturer partiellement l'échantillon, jusqu'à obtenir une valeur particulière, mais inconnue, de saturation S. La figure 2 schématise la mise en oeuvre de la désaturation par centrifugation, selon l'invention. L'axe de rotation de la centrifugeuse (AR) est représenté à gauche. A titre d'exemple, les longueurs des échantillons sont identiques (0.03m), et leurs hauteurs (diamètres s'il s'agit de carottes cylindriques) également (he=0,04m). La distance R représente une distance par rapport à l'axe de rotation (AR). Deux valeurs particulières sont indiquées sur la figure 2 : la distance Ro indique la distance entre cette axe et l'échantillon placé au sein de la cellule de centrifugation, au plus proche de l'axe AR, c'est-à-dire l'échantillon 1 (E1) ; la distance R,nax représente la distance entre l'axe de rotation (AR) et l'extrémité la plus éloignée de l'échantillon 2 (E2). Selon un exemple, Rmax = 0,263m. Lorsqu'un échantillon saturé d'un fluide 1 est centrifugé dans un autre fluide 2 de densité différente, une différence de pression AP apparaît, produisant un changement de saturation de l'échantillon. La différence de pression générée est une fonction de la distance R par rapport au centre de rotation, de la vitesse de rotation w et des densités des fluides. Cette pression est donnée par AP(R)= 1 w2(Rx ûRZ)Pw û Pg) 2 où pH, et pg sont les densités des deux fluides (saumure, air dans ces expériences) et R,nax• Lorsque la pression générée (Eq 1) liée à la rotation est supérieure à la pression capillaire, la saumure s'écoule. La saturation en saumure diminue jusqu'à un équilibre, obtenu quand le débit est nul. A l'équilibre et pour une pression donnée, la saturation est donnée par la courbe de pression capillaire. En effectuant la centrifugation à différentes vitesses de rotation, on obtient ainsi différentes saturations. Il n'est cependant pas nécessaire d'atteindre l'équilibre capillaire. D'après l'équation (1), la pression maximale (la saturation la plus faible) est obtenue à la face d'entrée (Ro en drainage) du fluide de plus faible densité. Sur la face de sortie, la pression 6 (1) capillaire est nulle et ainsi, la saturation est théoriquement égale à un. La zone à saturation élevée constitue le pied capillaire. A cause de cette inhomogénéité de pression, le profil de saturation n'est pas uniforme le long de l'axe de rotation, comme le montre schématiquement la figure 3, qui illustre un profil de saturation SW dans les échantillons El et E2 en contact capillaire, en fonction du rayon R Pour éviter des problèmes d'interprétation des mesures, liés à la saturation non uniforme dans l'échantillon, on superpose par exemple deux échantillons voisins issus d'une même carotte. Puis, on mesure les valeurs moyennes de saturation et de résistivité de l'échantillon 1. Pour que l'échantillon 1 ait un profil de saturation uniforme, il est nécessaire d'avoir un bon contact capillaire avec l'échantillon 2. A cette fin, on peut utiliser une membrane mouillable à la saumure entre les deux échantillons. A cette fin également, il est préférable que l'échantillon 1 n'ait pas un diamètre supérieur à l'échantillon 2 (he, figure 2). On obtient typiquement des profils de saturation comme le montre la figure 4. Cette figure montre des profils de saturation S, dans les échantillons El et E2, en fonction d'une longueur normalisée ((R-Rmax) / (Rmax Ro)), pour différentes vitesses de rotation de centrifugation (de 180 tour/min (courbe du haut) à 4500 tour/min (courbe du bas)). On peut toutefois remplacer l'échantillon 2 par un autre échantillon issu d'un autre milieu, pourvu que sa porosité et sa perméabilité soient sensiblement identiques à celles de l'échantillon 1. Ainsi, selon l'invention, on désature l'échantillon 1 dans une cellule de centrifugation, à côté d'un échantillon 2 adapté à limiter le pied capillaire dans l'échantillon 1. Pour ce faire, les échantillons sont disposés l'un à côté de l'autre, l'échantillon 2 étant l'échantillon le plus éloigné du centre de rotation de la centrifugeuse. De plus, les échantillons doivent être disposés de façon à optimiser le contact capillaire entre eux. Ensuite, on désature partiellement ces deux échantillons en les soumettant à une centrifugation selon une vitesse de rotation fixée et un temps de rotation suffisant pour permettre une désaturation de l'échantillon 1 (et de l'échantillon par conséquent). 2-Mesure de l'indice de résistivité RI de l'échantillon Dispositif La mesure de résistance est réalisée dans une cellule de mesure de résistivité à électrodes radiale. Un tel dispositif est décrit dans les brevets EP 0 701 128 et FR 2 781 573. Il est composé d'électrodes positionnées radialement autour de l'échantillon. On peut par exemple utiliser un dispositif comportant six électrodes, tel que celui décrit en figure 5A. La figure 5B est une projection verticale de la cellule et des électrodes. La procédure est décrite en référence à ces figures. Deux couples d'électrodes C+ et C-injectent un courant dans l'échantillon (1) et un autre couple d'électrode V+ et V- mesure le potentiel électrique. Elles sont connectées aux entrées de potentiel et de courant d'un impédance mètre (6) via des fils électriques (3). Pour éviter le séchage de l'échantillon pendant la mesure, des mousses (8) saturées de saumure sont disposées en haut et en bas pour maintenir une humidité élevée de l'air au contact de l'échantillon. Une des difficultés pour utiliser ce type d'appareil, est la nécessité d'obtenir un bon contact électrique avec l'échantillon. De plus, il faut que ce contact soit reproductible de façon à ce que les mesures soient elles-mêmes reproductibles. Les électrodes sont moulées ou insérées dans une gaine en élastomère (4) permettant d'appliquer une pression de confinement. Il peut s'agir d'une gaine en Viton par exemple. Le fluide de confinement contenu dans le compartiment (2) est relié à une pompe (5) remplie d'une huile visqueuse servant à appliquer une pression de confinement fixée à 30 bar. Cette pression nous permet d'avoir un bon contact électrique entre les électrodes et l'échantillon, par suite les champs électriques dans l'échantillon restent quasi constants entre les différentes mesures sur un même échantillon. Ceci nous garantit une condition de mesure constante pour l'indice de résistivité. Un second problème couramment rencontré, est la perte de grain lors de la manipulation de l'échantillon, conduisant à des erreurs de mesures. Selon l'invention, pour retirer l'échantillon sans l'endommager, après la mesure électrique, on utilise avantageusement la gaine (4) et la pompe (5). En effet, la pompe peut retirer l'huile dans la cellule, donc produire un écartement des électrodes permettant de retirer sans glissement l'échantillon. Mesure On injecte un courant entre C+ et C-, on mesure un potentiel entre V+ et V-, et on déduit une résistance R=U/I. On mesure la partie réelle et la partie imaginaire de la résistance à plusieurs fréquences. La mesure est faite avec un impédance mètre qui détermine l'impédance complexe. On s'affranchit ainsi des effets capacitifs (partie imaginaire). On prend en général la partie réelle à 1kHz pour calculer l'indice de résistivité. La partie imaginaire sert à évaluer la fiabilité de la mesure. À 1 kHz, elle doit être environ 10 fois plus faible que la partie réelle. La mesure se stabilise pendant une dizaine de minutes après l'application de la pression de confinement. 3-Mesure de la saturation On retire l'échantillon du dispositif Puis, on mesure la saturation à l'aide d'un appareil RMN qui mesure l'aimantation nucléaire. Un tel appareil est bien connu des spécialistes. Un noyau atomique peut absorber les rayonnements électromagnétiques d'une fréquence spécifique en présence d'un champ magnétique. Ce phénomène est appelé la Résonance Magnétique Nucléaire (RMN). Dans la description de l'invention, on considère le noyau d'hydrogène est présent dans la saumure. Un proton peut passer à un niveau d'énergie plus haut s'il absorbe une quantité d'énergie 20 à la fréquence de Larmor co, cette dernière étant reliée au champ magnétique Bo dans lequel est placé l'échantillon: DE=hw et w=ïB0 (3) où h Est la constante Planck divisée par 2n. y est le rapport gyromagnétique du noyau considéré. Bo est un champ magnétique statique le long l'axe Z. Un autre champ magnétique 25 oscillant BI est appliqué dans le plan XY. Quand la fréquence de BI approche w, les protons irradiés passent au niveau d'énergie le plus haut. Le retour des protons vers le niveau d'énergie inférieur après la suppression du champ BI s'appelle la relaxation. Les temps caractéristiques nécessaires à un noyau pour revenir à l'équilibre sont le temps de relaxation spin-réseaux noté Tl, et le temps de relaxation spin-spin 30 T2. Nous mesurons T2 en général. 9 L'appareil RMN mesure l'aimantation nucléaire globale M(t) retournant à l'équilibre en fonction du temps à l'aide de la séquence CPMG. La relaxation spin-spin T2 est décomposée en une somme d'exponentielles telle que : 2n .r" M(t=2n~z)=1A;exp ù T n.=1...n (4) 2i Les principaux paramètres d'acquisition sont : i, demi-temps inter écho, par exemple 60 microsecondes, et n, nombres d'échos ou de points de mesure (il varie en fonction des échantillons). Dans la mesure du possible, on tente d'atteindre un rapport signal sur bruit supérieur à 100. En général la mesure de l'aimantation totale permet d'obtenir directement le volume poreux Vp en comparant le signal avec une quantité de saumure connue Verdon, V = Mi(Sw =100%)) x V e p Mi(saumure etalon) talon où Mi est l'aimantation M(t==0). Dans le système saumure-air, l'air ne produit pas de signal RMN, et on obtient directement les saturations Sw au cours de la centrifugation ou de l'imbibition en faisant le rapport des aimantations Mi selon: Mi(Sw)) Mi(Sw =100%)) La méthode selon l'invention consiste donc à mesurer, à l'aide d'un appareil RMN, les aimantations nucléaires, puis en utilisant la formule (6). Il faut donc connaître Mi(Sw =100%)) . Les échantillons étant saturés au départ de la méthode, il est nécessaire de réaliser une mesure de l'aimantation nucléaire avant l'étape de désaturation. Exemples de mesure (5) (6) 1) Grès de Fontainebleau: La figure 6 montre le résultat de mesures de résistivité et de saturations en saumure sur deux Grès de Fontainebleau de grande perméabilité (1000 mD). On observe une bonne cohérence entre les deux échantillons. La pente de la courbe RI = S,' est proche de 2 au dessus d'une saturation d'environ 0.20. La méthode proposée a permis de montrer qu'en dessous de 0.20, on observe un changement de pente, qui est interprété comme le début d'un régime de conduction dominée par les films. 2) Calcaire à structure bimodale, cas 1 Sur une structure complexe tel qu'un calcaire contenant deux tailles de pores, la courbe d'indice de résistivité (figure 7) montre un comportement particulier: on observe une déviation 10 vers le haut. Les profils de saturation montrent que la saturation est uniforme aux différentes étapes de saturation. Les distributions de temps de relaxation montrent deux modes correspondants chacun à une distribution de taille de pore. Lorsque le système poreux est partiellement désaturé, on voit que les grands pores sont vidés en premier, alors que les petits pores restent complètement saturés jusqu'à à une saturation moyenne de 0.38. 15 3) Calcaire à structure bimodale, cas 2 Dans le cas de ce calcaire, la courbe (figure 8) d'indice de résistivité ne présente pas d'anomalies marquantes. Il existe une hystérésis faible en imbibition. Les temps de relaxation montrent que les deux populations sont vidées en même temps dans ce cas. En imbibition, les gros pores restent vides. Modes de réalisation particuliers Selon un mode de réalisation, on complète les mesures effectuées sur les échantillons, par une mesure de distribution de temps de relaxation, à l'aide d'un appareil RMN. Ces mesures de distribution servent à déterminer les pores vidés pour une pression donnée. Cela 25 influence la réponse électrique. L'utilisation des techniques RIVIN complémente ainsi l'information obtenue et est utile pour l'interprétation de la réponse électrique. En effet, dans des structures complexes comprenant plusieurs populations de tailles de pore (par microporosité et macroporosité), la manière dont les différentes populations se drainent, influence la réponse électrique. 30 La distribution de temps de relaxation donne une information sur la distribution de taille de pore suivant la relation : 20 S p = T z û P2 Vp T2B où p2 est le coefficient de relaxivité de surface, et T2B le temps de relaxation propre de la saumure. Le temps de relaxation T2 d'un pore est lié à son rapport surface sur volume (Sp/Vp). La distribution de T2, calculée selon l'équation 4, traduit la distribution du rapport V/S dans le réseau poreux, appelé distribution de "tailles de pores". Ceci est vrai dans le cas d'un échange lent entre les pores. Ainsi, les temps de relaxation les plus longs correspondent aux pores de plus grandes dimensions. Dans le système saumure-air, le temps de relaxation T2 d'un échantillon partiellement saturé est lié à la surface mouillée par l'eau Sn, et à la saturation S, , de ce pore : ]l S. 1 û=P2 * +ù72 p Swp T2B Selon un autre mode de réalisation, on complète les mesures effectuées sur les carottes, par une mesure de profil de saturation. Ces mesures de profils de saturation permettent de savoir si la saturation est uniforme ou pas : si les échantillons 1 et 2 ne sont pas en contact 15 capillaires, il y aura un profil de saturation sur l'échantillon 1, semblable à 2. Il faut alors recommencer les mesures en améliorant le contact capillaire entre les échantillons. En utilisant les résultats d'une injection de mercure bien connue des spécialistes, on peut prévoir le profil de saturation en fonction de la vitesse de rotation. En supposant que l'échantillon est mouillable à l'eau, la relation de conversion est la suivante : YHg cos °Hg où y est la tension superficielle et 8 l'angle de contact, les indices désignant le système mercure - air (Hg) et le système air - saumure. A une position et vitesse de rotation donnée par l'équation (1), on calcule la pression obtenue et on cherche la saturation correspondante dans la courbe de pression capillaire air -saumure. Cela nous permet de trouver la saturation 25 locale à chaque position R, et d'obtenir ainsi les profils de saturation (figure 4). On voit que dans l'échantillon 1, la saturation est quasi-uniforme et que le pied capillaire est contenu dans l'échantillon 2. 12 (7) (8) 20 1p _ p Ysaumure-air cos °saumure--air saumure-air Hg Profil de saturation La mesure du profil de saturation est réalisée avec une séquence standard d'imagerie RMN, (une séquence d'écho de Hahn avec gradient pulsé). Lorsqu'un gradient de champ est appliqué, le champ magnétique varie linéairement le long la direction du gradient. La fréquence de Larmor des protons varie proportionnellement le long du gradient magnétique. w(y) = 7Bo (1 + Gy) (9) La fréquence locale code la position des protons suivant la direction du gradient. En reliant la densité de protons mesurée et les positions, on obtient le profil de densité protonique dp(z) le long la direction du gradient. Dans le système saumure-air, la mesure du profil de densité dp(z) permet d'obtenir le profil de saturation Sw(z): Sw(z) _ dp(z) dp(z,Sw =100%) Selon un autre mode de réalisation, on remplace le drainage par centrifugation par une imbibition spontanée. L'échantillon est alors posé dans un bécher dans lequel on ajoute une quantité contrôlée de saumure. La saumure se répand alors dans tout l'échantillon. En répétant cette opération, on peut obtenir les différents paliers de saturation en imbibition. La face d'entrée de la saumure en imbibition correspond à la face de sortie de la saumure en drainage. L'échantillon et la saumure sont à la pression atmosphérique. Les processus de mesure électrique et de mesure RMN sont les mêmes qu'en drainage. Applications de la technique de mesure selon l'invention Selon un mode de réalisation de l'invention, on peut utiliser la méthode de mesure rapide de la saturation et de l'indice de résistivité d'un milieu poreux, pour définir une relation permettant de déduire la saturation en fonction de la résistivité. Pour ce faire on réalise plusieurs mesures à partir du même échantillon issu du milieu. A la fin de la première série de mesures, l'échantillon est remis dans la centrifugeuse pour obtenir un autre palier de saturation. On modifie alors la vitesse de rotation o. Puis on réitère 13 (10) les étapes de mesures, de façon à obtenir différentes saturations, et pour chacune de ces saturations, on mesure l'indice de résistivité et la saturation. On obtient ainsi plusieurs couples de valeurs de résistivité (RI) et de saturation (Sw) pour l'échantillon étudié. A partir de ces couples, on estime une relation permettant de déduire la saturation en fonction de la résistivité. On peut par exemple utiliser la loi d'Archie, et déterminer, grâce aux couples, la valeur de l'indice n : avec : RI Selon un autre mode de réalisation, la méthode peut être appliquée à la détermination de la saturation en saumure (fluide conducteur) d'une formation souterraine, telle qu'un réservoir pétrolier, ou un réservoir de stockage de gaz, par l'interprétation de mesures électriques, telles que des diagraphies. En effet, cette nouvelle méthode de mesure permet de caractériser très rapidement un grand nombre d'échantillons à un coût faible, en raison de l'utilisation d'une centrifugeuse et d'un montage électrique simple. Bien qu'elle ne permette pas d'explorer des domaines de pression et température proches du réservoir, il est utile d'établir une courbe de référence dans des conditions simplifiées, pour évaluer, par la suite avec des techniques plus contraignantes, les effets de pression et de température. Ainsi, la méthode peut être utilisée dans le but de choisir les échantillons, avant une interprétation plus précise en conditions de pression et de température de la formation. La méthode comporte alors les étapes suivantes : On extrait des échantillons de la formation. Puis on les sature en saumure. Puis, en utilisant la méthode précédente, on détermine une relation permettant de déduire la saturation en fonction de la résistivité pour chacun des échantillons. On sélectionne alors des échantillons de façon à ce que chacun d'entre eux ait une relation différente des autres échantillons. Par exemple, si l'on choisit d'utiliser la loi d'Archie, on estime une valeur de l'indice �n� pour chaque échantillon. Tous les échantillons ayant une relation commune, indice de résistivité SW saturation en saumure du milieu n : un nombre entier c'est-à-dire avec une même valeur de �n �, sont regroupés. Puis, on en choisit un seul parmi cette famille d'échantillon. Ces échantillons sélectionnés sont représentatifs d'une corrélation particulière entre la saturation et la résistivité. Ensuite, on mesure l'indice de résistivité et la saturation en saumure de chacun de ces échantillons sélectionnés, à l'aide d'une méthode adaptée à fonctionner selon des conditions de température et de pression propres à la formation souterraine. On peut par exemple utiliser la méthode décrite dans Fleury M., "Advances in Resistivity Measurements using the FRIM Method at Reservoir Conditions. Applications to Carbonates". Proceedings of the International Symposium of the Society of Core Analysts, Septembre 2003, Pau, FRANCE. Enfin, on déduit la saturation de la formation souterraine, en appliquant la relation à au moins une diagraphie de résistivité eléctrique, mesurée à travers la formation. On applique les différentes relations en fonction de la profondeur de la diagraphie et de celle de l'échantillon correspondant, ou à travers un regroupement des différentes mesures appelé "Rock-Typing" non nécessairement reliées de manière simple à la profondeur. REVENDICATIONS 1. Méthode pour mesurer une saturation en un fluide conducteur et une résistivité d'un milieu poreux, initialement saturé par ledit fluide, caractérisée en ce que : a)- on extrait au moins un premier échantillon solide dudit milieu ; b)- on place ledit premier échantillon dans une cellule de centrifugation, à côté d'un second échantillon adapté à limiter un pied capillaire dans ledit premier échantillon, les échantillons étant disposés l'un à côté de l'autre de façon à optimiser un contact capillaire entre eux ; e)- on désature partiellement les deux échantillons en les soumettant à une centrifugation selon une vitesse de rotation fixée et un temps de rotation suffisant pour permettre une désaturation dudit premier échantillon ; d)- on mesure la résistivité du premier échantillon, en le plaçant dans une cellule de mesure de résistivité, adaptée à injecter un courant par l'intermédiaire d'au moins un premier couple d'électrodes, et à mesurer un potentiel entre au moins un second couple d'électrodes, lesdites électrodes étant disposées radialement autours dudit premier échantillon ; et e)- on détermine la saturation du premier échantillon en mesurant une aimantation nucléaire de l'échantillon à l'aide d'un appareil RMN. 2. Méthode selon la revendication 1, dans laquelle ledit second échantillon est extrait dudit milieu poreux. 3. Méthode selon l'une des revendications précédentes, dans laquelle on optimise le contact capillaire en intercalant entre les deux échantillons une membrane mouillable audit fluide. 4. Méthode selon l'une des revendications précédentes, dans laquelle on réitère les étapes b) à e), en modifiant la vitesse de rotation de façon à obtenir plusieurs couples de valeurs de résistivité et de saturation pour ledit premier échantillon, puis on estime une relation permettant de déduire la saturation en fonction de la résistivité. 5. Méthode pour mesurer la saturation en un fluide conducteur d'une formation souterraine, dans laquelle : - on extrait des échantillons de ladite formation que l'on sature avec ledit fluide ; - on détermine une relation permettant de déduire la saturation en fonction de la résistivité pour chacun desdits échantillons, à l'aide de la méthode selon la revendication 4 ; on sélectionne des échantillons de façon à ce que chacun d'entre eux ait une relation différente des autres échantillons ; on mesure l'indice de résistivité et la saturation en fluide de chacun desdits échantillons sélectionnés, à l'aide d'une méthode adaptée à fonctionner selon des conditions de température et de pression propres à ladite formation souterraine ; on réalise au moins une diagraphie de résistivité eléctrique à travers ladite formation ; on en déduit la saturation en fluide dans la formation en appliquant ladite relation.