FR 2844355 A1 20040312 FR 0211283 A 20020911 Désignation du domaine technique La présente invention concerne une méthode de mesure de la mouillabilité d'échantillons de roches par la mesure du temps de relaxation par résonance magnétique nucléaire à faible champ. La méthode trouve des applications notamment pour l'analyse de roches prélevées dans une formation souterraine recelant ou susceptible de receler des hydrocarbures La connaissance de différents paramètres et notamment de la mouillabilité des roches, est aussi utile notamment quand on doit procéder à une récupération assistée d'une formation, par une injection d'un fluide sous pression, et que l'on doit déterminer par des 10 tests préalables le fluide (liquide ou gaz) qui convient le mieux pour déplacer les effluents. L'invention trouve aussi des applications en génie civil pour faire de l'hydrologie de terrains pour évaluer leur degré de pollution par exemple, ou encore dans le bâtiment pour tester des matériaux de construction afin notamment de décider de traitements hydrofuges par exemple. Présentation de l'art antérieur On connaît des méthodes pour déterminer la mouillabilité des roches vis à vis de l'eau et de l'huile qui peuvent y être contenues, comportant la réalisation de cycles de drainage de la roche c'est-à-dire de déplacement des fluides visant à diminuer la saturation en eau, suivi d'une imbibition, en désignant par ce terme un déplacement des fluides 20 permettant d'augmenter la saturation en eau (Sw) de la roche. La pression capillaire Pc en un point se définit comme la différence à l'équilibre entre la pression P(huile) de l'huile et celle P(eau) de l'eau. Ce paramètre n'a de sens que si les deux fluides sont en phase continue dans le milieu poreux. Pour un milieu mouillable à l'eau, seules les valeurs positives ont un sens. Par contre, lorsque le milieu a une mouillabilité mixte, les fluides 25 peuvent rester en phase continue aussi bien pour les pressions capillaires (Pc) positives que négatives. Pour une application de ce type, un cycle complet de mesure de la pression capillaire comporte généralement (Fig. 1): 2 2844355 a) un drainage primaire positif d'un échantillon saturé initialement en eau à 100% (courbe 1); b) une imbibition positive (courbe 2); c) une imbibition négative (courbe 3); d) un drainage négatif (courbe 4); et e) un drainage secondaire positif (courbe 5). l existe différents types de dispositif capables de dresser les courbes de la Fig. 1. Dans un dispositif dit à "plaque poreuse", décrit notamment par le brevet US 4 506 542, l'échantillon de roche poreuse contenant deux fluides en phases continues, est 10 placé dans une cellule allongée terminée à ses deux extrémités opposées par des barrières capillaires perméables à un premier fluide. Ce premier fluide est injecté sous pression au travers de la première membrane et on mesure la différence de pression entre la pression d'injection et la pression du fluide évacué à l'autre extrémité. Les pressions des deux fluides ainsi que la pression capillaire Pc sont constantes le long de l'échantillon, et la 15 saturation est supposée uniforme. Il est également connu de procéder par centrifugation à vitesse progressive au moyen de dispositifs de centrifugation tels que ceux décrits par exemple par les brevets FR 2 772 477 (US 6 185 985) ou FR 2 763 690 ou encore dans les brevets EP 603 040 (US 5 463 894) ou FR 2 798 734) du demandeur. L'échantillon initialement saturé avec un premier liquide (de la saumure par exemple) est placé (Fig.10) dans un récipient ou godet allongé A contenant un deuxième fluide de densité différente de l'huile par exemple) . Le godet A est fixé au bout d'un bras tournant B, et on lui applique une force centrifuge de façon à étudier les déplacements des fluides dans l'échantillon au cours d'au moins deux phases distinctes. Durant une première 25 phase de drainage, on soumet alors l'ensemble à une force centrifuge dirigée suivant la longueur du récipient de manière à exercer sur lui une force d'expulsion qui tend à faire sortir une partie du premier fluide. Le champ de pression créé par la centrifugation s'exprime en fonction de la densité r, du rayon R et de la vitesse angulaire w, par la relation: 1/2 w2**-. p (Rmax2- R2), pour chaque fluide. On impose que la pression des 30 deux fluides à la sortie de l'échantillon soit la même et qu'elle s'annule à la sortie. Dans le 3 2844355 même temps, le deuxième fluide pénètre à l'intérieur de l'échantillon. Les deux fluides se déplacent à l'intérieur de l'échantillon jusqu'à une position d'équilibre o la force due à la pression capillaire dans les pores, compense la force centrifuge exercée. Une sonde de mesure, est placée dans le godet, sur le côté de l'échantillon. La sonde peut être de type 5 capacitif par exemple et comporter une tige métallique, isolée par une gaine en céramique. On mesure la capacité entre la tige et le fluide conducteur (la saumure), qui est proportionnelle à la hauteur immergée. Avec ce moyen de mesure, la précision de mesure est de 1,5% du volume poreux. La sonde détecte la position de l'interface entre les deux liquides dans le godet et transmet les signaux de mesure à un automate E de commande et 10 d'acquisition de signaux de mesure regroupant des moyens hydrauliques de circulation des liquides et un appareil d'acquisition Dans la phase de ré-imbibition, on diminue la vitesse de façon à étudier la réintégration du fluide initial dans l'échantillon. Les saturations locales mesurées avec ce type de dispositif, sont calculées par un programme d'inversion à partir de la quantité totale 15 d'eau expulsée hors de l'échantillon. Suivant une autre méthode dite "dynamique", on place un échantillon dans une cellule allongée terminée à ses deux extrémités par des membranes perméables à l'eau. A une première extrémité, on injecte directement de l'huile sous pression dans l'enceinte. De l'eau y est également injectée mais cette injection est effectuée au travers de la membrane 20 et à une pression inférieure. A l'extrémité opposée, l'huile est évacuée directement alors que l'eau sort au travers de la membrane terminale. En réglant les débits d'injection d'huile et d'eau, on s'arrange pour que la pression capillaire soit la même à l'entrée de l'enceinte qu'à sa sortie, ce qui entraîne une saturation uniforme pouvant être déduite du bilan des fluides. La pression capillaire est obtenue par exemple en mesurant la différence entre les 25 pressions de l'huile et de l'eau à la sortie de l'enceinte. Une telle méthode est décrite notamment par Brown H.W. in "Capillary Pressure Investigations" Petroleum Transaction AIME, vol. 192, 1951. Des exemples de mise en oeuvre sont décrits par exemple dans les brevets EP 729 022 (US 5 698 772) ou EP-A- 974 839 (US 6 229 312) du demandeur. On connaît également une méthode dite semi dynamique dans laquelle on confine 30 un échantillon de roche imbibé avec un premier fluide dans une cellule close, on injecte un autre fluide sous pression à une première extrémité de l'enceinte et on balaie l'extrémité opposée par un fluide à basse pression mis en circulation par des moyens de pompage qui 4 2844355 entraîne au dehors le fluide drainé. Le dispositif comporte des moyens de mesure de la pression et la saturation de l'échantillon différents de sa longueur, la quantité de fluide évacué et la résistivité électrique de l'échantillon. Cette méthode est mise en oeuvre par exemple dans le brevet FR 2 708 742 (US 5 679 885) du demandeur. Ces courbes de drainage et d'imbibition ayant été établies, il est connu de calculer les indices de mouillabilité WI à partir des aires Ad et Ai délimitées par les courbes de pression capillaire positives et négatives, comme le montre la Fig.1, par la relation WI(USBM) = Lg Ad Ai La méthode selon l'invention La méthode de mesure de la mouillabilité d'un échantillon de roche poreuse en présence d'eau et d'huile, selon l'invention comporte la détermination de la surface poreuse mouillée *par l'en eau et de la surface poreuse mouillée par l'huile quand l'échantillon est saturé en eau et en huile, et le calcul de l'indice de mouillabilité par combinaison des valeurs obtenues des dites surfaces. La détermination de la surface poreuse mouillée en eau et de la surface poreuse mouillée par l'huile quand l'échantillon est saturé en eau et en huile, est obtenue par exemple par le biais de mesures de temps de relaxation de l'échantillon placé dans un appareil à résonance magnétique nucléaire. On détermine par exemple l'indice de mouillabilité par la relation: IMN = SM -SM ou la relation IN = îocr SMw SMNS + SM, Mv C10 SMo dans lesquelles SM. est la surface poreuse mouillée par l'eau et SM, est la surface poreuse mouillée par l'huile lorsque le milieu poreux est saturé en eau et huile. Suivant un mode de mise en oeuvre, on détermine l'indice de mouillabilité par les opérations suivantes: a) on mesure les temps caractéristiques de relaxation de l'échantillon saturé en eau; 2844355 b) on mesure les temps caractéristiques de relaxation de l'huile dans l'échantillon en présence d'eau, dans une zone proche de la saturation (Swir) de l'échantillon; c) on mesure les temps caractéristiques de relaxation de de l'eau dans l'échantillon en présence d'huile, dans une zone proche de la saturation résiduelle (Sor); d) on mesure les temps de relaxation de l'échantillon dans un état o son point de saturation à 100% en huile est atteint; et e) on combine les mesures des différentes temps de relaxation pour obtenir le dit indice. Suivant un mode préféré de mise en oeuvre, on détermine les temps de relaxation des étapes a) à c) après application à l'échantillon d'une centrifugation. Suivant un mode préféré de mise en oeuvre, on détermine les temps de relaxation de l'étape d) après déplacement forcé des fluides dans l'échantillon placé dans une cellule de confinement. De préférence, on choisit une huile dont le temps de relaxation intrinsèque (TB) est le plus grand possible et le plus proche de l'eau, du dodécane par exemple. Les temps caractéristiques de relaxation sont par exemple ceux qui correspondent soit au maxima des courbes de saturation soit à des valeurs moyennes des dites courbes. On vérifie à l'expérience que la mesure de IRMN obtenue par la méthode est aussi sensible mais nécessite bien moins de temps à établir et qu'elle est applicable à un grand nombre d'échantillons. Présentation succincte des figures Les caractéristiques et avantages de la méthode et du dispositif selon l'invention, apparaîtront plus clairement à la lecture de la description ci-après d'un exemple non limitatif de réalisation, en se référant aux dessins annexés o: - la figure 1 montre des exemples de courbes de pression capillaire nécessaires à la 25 détermination de l'indice de mouillabilité classique USBM (United States Bureau of Mines); 6 2844355 - les figure 2a et 2b montrent respectivement la distribution des temps de relaxation T2 pour une roche mouillable à l'eau saturée en eau et en huile, et une représentation de la distribution de phase (matrice en hachuré, eau en gris clair et huile en gris plus foncé); - les figure 3a et 3b montrent respectivement la distribution des temps de relaxation T2 5 pour une roche de mouillabilité intermédiaire saturée en eau et en huile (même partie centrale que pour la figure 2), et une représentation de la distribution de phase (matrice en hachuré, eau en gris clair et huile en gris plus foncé); - la figure 4 montre sur des courbes de pression capillaire, les valeurs de saturation o les mesures de relaxation RMN sont faites - la figure 5 montre un exemple de distribution des temps de relaxation T2 aux valeurs de saturation indiquées à la figure 4 (de bas en haut: 100% eau, Swi, Sor et 100% dodecane); et - la figure 6 montre la relation entre l'indice IRMN et l'indice IUSBM standard; - la figure 7 montre schématiquement une cellule de confinement avec un circuit 15 permanent de chauffage par circulation de fluide; et - la figure 8 montre schématiquement l'agencement d'un dispositif de centrifugation utilisé pour la mise en oeuvre de la méthode. Description détaillée de la méthode On rappelle que la technique d'analyse RNIN consiste essentiellement on le 20 rappelle, à appliquer à un objet à tester un premier champ magnétique statique Bo de polarisation destiné à aligner les noyaux des protons d'hydrogène initialement orientés au hasard, suivant la direction du champ puis à un deuxième champ magnétique impulsionnel oscillant à la fréquence de Larmor perpendiculaire au premier créé par des bobinages excités par un signal de commande pour réaliser une expérience de résonance magnétique 25 nucléaire. Quand ce champ impulsionnel cesse, le retour des noyaux vers leur état initial ou relaxation, engendre des signaux électromagnétiques (échos) qui sont détectés et analysés. A partir des caractéristiques d'amplitude de ces signaux, on détermine la présence de paramètres physiques de l'objet. 7 2844355 On propose de définir un nouvel indice de mouillabilité construit en combinant des valeurs de la surface poreuse SM, mouillée par l'eau et de la surface poreuse SM, mouillée par l'huile lorsque le milieu poreux est saturé en eau et huile. L'indice peut être calculé par exemple par la relation suivante: RMN SMW + SMo (la) ou encore par la relation: IRPMN =logo SM (lb) Les quantitées SM, et SMo sont obtenues par exemple par la mesure du temps de relaxation dominant dans une expérience de résonance magnétique nucléaire à faible 10 champ. En effet, dans une telle expérience, la relation fondamentale liant le temps de relaxation longitudinal T1 ou transversal T2 à la surface S et au volume de ce pore est la suivante: 1 1 S ( T,2 TIB,2B V ou T1B,2B est le temps de relaxation du fluide hors du milieu poreux. Fondamentalement, cette relation provient du fait que les molécules dans le pore subissent des mouvements effets de diffusion et interagissent avec la paroi poreuse pendant le temps de mesure (le temps de relaxation typique est de l'ordre de 100 ms). Les interactions de surface sont rassemblées dans le coefficient Pî,2 appelé relaxivité de surface. Nous avons négligé par 20 ailleurs un terme de gradient qui est important lorsque la mesure est faite à haut champ magnétique (>0. lT). La relation 2 est strictement valable pour un pore saturé d'un seul fluide. En général, les milieux poreux naturels présentent une distribution de tailles de pore et donc du rapport S/V. On observe ainsi en général une distribution des temps de relaxation mas ceci ne modifie pas la méthode décrite ici. Lorsque deux fluides sont présents dans un pore au sein du milieu poreux, le même mécanisme physique de diffusion est valable, mais la distribution des deux fluides vis à vis 8 2844355 de la surface poreuse a une importance cruciale. Par exemple, lorsque le milieu est mouillable à l'eau, l'eau se trouve en surface et interagit avec elle, tandis que l'huile se trouve à l'intérieur du pore et n'interagit pas avec la surface. Lorsque l'on mesure la distribution des temps de relaxation pour un tel système, on observe le résultat de la figure 5 2a. L'eau est marquée par des temps de relaxation plus courts qu'à la saturation 100% en eau, car le volume d'eau V a diminué (voir équation 1). L'huile est caractérisée par un temps de relaxation qui est le même que si la mesure était faite hors milieu poreux (TIB,2B) car il n'y a pas d'interaction avec la surface poreuse. fl est par ailleurs connu que les interactions à l'interface eau/huile ne donnent pas lieu à une relaxation différente de celle 10 de l'huile hors milieu poreux, à condition que l'huile considérée soit légère (par une huile raffinée tel le dodecane). Schématiquement, la distribution de phase est représentée à la figure 2b. Lorsque le même milieu poreux a une mouillabilité différente, non franche à l'eau, la distribution des temps de relaxation est modifiée (figure 3a). Essentiellement, le temps caractéristique de l'huile sera plus faible du fait d'une interaction de l'huile à la 15 surface poreuse car l'eau n'est plus en premier voisin de la surface. On représente la distribution possible des deux fluides à la figure 3b. Un tel effet a été obtenu en utilisant les procédures standards permettant de reproduire les propriétés de surface de roches provenant d'un réservoir pétrolier: nettoyage puis vieillissement dans l'huile réservoir à la température du réservoir à la saturation irréductible en eau. Ainsi, on généralise la relation (2) à un système diphasique eau-huile. Pour l'eau, nous avons: i i Sm T = p1 V (3) Tl,2w TBw,2Bw W et pour l'huile: 1 1 + SM, (4) T1T,20 lo,2Bo 10,20 VO On peut donc déterminer les quantités SM, et SM, par la mesure des temps de relaxation Tl ou T2 en milieux poreux, des temps de relaxation TIB ou T2B des fluides hors milieux poreux, des volumes liquides VO ou V, présents dans le milieux poreux. Les constantes d'interaction p peuvent être déterminées séparement, mais nous verrons que seul le rapport de ces constantes est nécessaire. 9 2844355 En introduisant les relations 3 et 4 dans la relation de base 1, nous obtenons: Sw( - T) - CpSo( - T) r _ w bw o b IRMN T Tb (5) sw( 1- -) + CpSo(- - -) T T T T I w lbw so lbo o Sw=Vw/Vp est la saturation en eau, So=Vo/Vp est la saturation en huile Vp le volume poreux Cp= P- est le rapport des relaxivités pour les surface mouillées par l'eau (w) et par l'huile (o). Pour alléger la notation, les indices 1 et 2 ont été supprimés. Il reste donc à déterminer les saturations les plus appropriées pour mesurer les quantités SM, et SM0. Pour des raisons liées au calcul des distributions de temps de 10 relaxation, on choisira la saturation irréductible en eau (Swi, figure 4) pour déterminer la surface mouillée par l'huile en présence d'eau (SMw), et la saturation résiduelle en huile (Sor, figure 4) pour déterminer la surface mouillée par l'eau en présence d'huile. On choisit également une huile dont le temps de relaxation intrinsèque (TB) est le plus grand possible et le plus proche de l'eau. En effet, si TB est trop faible, la méthode sera limitée à 15 des milieux poreux dont le rapport S/V est grand (des tailles de pore petites), ou à des milieux dont la relaxivité de surface est grande. Le dodecane par exemple est une huile raffinée qui est adaptée à la mesure proposée car son temps de relaxation intrinsèque (TBO) est de 1 s, proche du temps de relaxation de l'eau (TBW, environ 2.7 s). Il faut en général exclure les huiles pétrolières brutes car elles ont des temps de relaxation trop faibles et 20 possèdent également une distribution de temps de relaxation qui gêne considérablement l'analyse. Les détails pratiques du calcul sont détaillés ci-dessous. A partir de la distribution des temps de relaxation aux différentes saturations (voir l'exemple de la figure 5), on peut facilement déterminer le temps de relaxation correspondant au pic dominant (Tw à Sor, To à Swi), et ce sont ces valeurs que l'on utilise 25 dans la relation 5. Il est connu que la détermination du temps de relaxation dominant est robuste et dépend peu du procédé de calcul des distributions. On note que le pic dominant correspond à celui de l'huile à Swi, et à l'eau à Sor. La présence du deuxième fluide à ces 2844355 différentes saturations est importante du point de vue physique mais perturbe peu le calcul. On peut montrer que quelque soit le temps de relaxation du deuxième fluide présent en faible quantité (l'eau à Swi, l'huile à Sor), le temps de relaxation du fluide dominant fluctue peu. Cet aspect est important pour la robustesse du calcul. La figure 5 présente la 5 distribution des temps de relaxation transversal T2. On peut également utiliser le temps de relaxation longitudinal Ti mais cette mesure est moins favorable en général car T1>T2 mais T1B=T2B pour un fluide donné. Pour déterminer le rapport de relaxivité de surface Cp, on utilise les temps de relaxation dominants Twloo et T.Ioo respectivement aux deux saturations Sw=100% et 10 So=100%. La formule utilisée est la suivante: C p= 1/TW100 -1TBW pO 1/TT100 -1/TBO (6) La figure 5 illustre les deux mesures effectuées. Ainsi, au total, on doit donc mesurer les distributions de temps de relaxation à quatre saturations: 100% saumure, 100% huile et les saturations extrêmes Swi et Sor. Pour cela on pourra utiliser la technique 15 de la centrifugation pour atteindre, en partant de Sw=100%, Swi puis Sor. Un palier de rotation unique à vitesse maximale permet de les atteindre rapidement. Pour la mesure à So=100%, on utilise une succession de déplacements miscibles en injectant un alcool puis l'huile raffinée. Mise en oeuvre Pour la mise en oeuvre de la méthode, on utilise par exemple un appareil de mesure RMN qui comporte classiquement (Fig.7) des aimants 1 avec un entrefer dans lequel sont placés des bobinages 2. Les bobinages sont connectés à un boîtier électronique pour engendrer les signaux excitateurs (signaux dans la gamme des radiofréquences par exemple) créant un champ électrique oscillant et acquérir la réponse de l'échantillon S aux 25 signaux excitateurs. Un dispositif de ce type est mis en oeuvre par exemple dans la demande de brevet FR-A-2 823 308. il 2844355 Exemple pratique de mise en oeuvre L'indice de mouillabilité est obtenu en réalisant par exemple la succession d'étapes suivantes avec un échantillon de réservoir dont la mouillabilité est intermédiaire: 1 On sature l'échantillon à 100% avec de la saumure de gisement et on le place dans 5 l'appareil de mesure RMN (celui de la Fig.7-8 par exemple) pour mesurer les distributions des temps de relaxation TI, T2 et en déduire le paramètre TwlOO de l'équation 6; 2 On le place ensuite dans le dispositif de centrifugation (tel que celui schématisé à la Fig.10) et on le centrifuge à la vitesse maximale dans un godet rempli d'huile jusqu'à atteindre la saturation irréductible en eau Swir; après quoi on mesure les distributions des 10 temps de relaxation T1, T2 dans l'appareil de RMN, de façon à en déduire To de l'équation ; 3 Placé de nouveau dans le dispositif de centrifugation, on centrifuge l'échantillon à vitesse maximale jusqu'à atteindre le point de saturation résiduelle en huile Sor; puis on le transfère dans l'appareil RMN de nouveau pour mesurer les temps de relaxation TI, T2 et en 15 déduire T. 4 Ceci fait, on place l'échantillon dans une cellule de confinement telle que celles décrites dans les brevets précités et on injecte un solvant miscible à l'eau et à l'huile puis de l'huile (par exemple du dodécane) jusqu'à atteindre la saturation 100% en huile. Ce point étant atteint, on mesure de nouveau les temps de relaxation TI, T2 dans l'appareil RMN 20 d'o l'on déduit la valeur de Toloo.nécessaire à la détermination de Cp. On connaît ainsi tous les paramètres pour calculer le coefficient de mouillabilité IRMN. Comparaison avec l 'indice IuSBM La comparaison entre l'indice de mouillabilité lUSBM et le nouvel indice IRMN 25 montre une bonne corrélation entre ces deux quantités (figure 6) pour des roches réservoir. On peut donc avantageusement remplacer la mesure de IUSBM par une mesure de IRMN aussi sensible mais nécessitant bien moins de temps à établir et applicable à un grand nombre d'échantillons. 12 2844355 On a -décrit un exemple de mise en oeuvre o l'on utilise une méthode de relaxométrie de type RMIN pour mesurer la surface et le volume des pores. Cette méthode n'a cependant rien de limitatif. On peut utiliser toute autre méthode d'analyse et notamment la cryomicroscopie. 13 2844355 REVENDICATIONS 1) Méthode de mesure de la mouillabilité d'un échantillon de roche poreuse en présence d'eau et d'huile, caractérisée en ce qu'elle comporte la détermination de la surface poreuse mouillée en eau et de la surface poreuse mouillée par l'huile quand l'échantillon est 5 saturé en eau et en huile, et le calcul de l'indice de mouillabilité par combinaison des valeurs obtenues des dites surfaces. 2 Méthode selon la revendication 1, caractérisée en ce que la détermination de la surface poreuse mouillée en eau et de la surface poreuse mouillée par l'huile quand l'échantillon est saturé en eau et en huile, est obtenue par le biais de mesures de temps de 10 relaxation (TI, T2) de l'échantillon placé dans un appareil à résonance magnétique nucléaire 3) Méthode selon la revendication 1 ou 2, caractérisée en ce que l'indice de mouillabilité est obtenue par la relation: = SM" -SM" RMVSMW + SM0 o SM, est la surface poreuse mouillée par l'eau et SM, est la surface poreuse mouillée 15 par l'huile lorsque le milieu poreux est saturé en eau et huile. 4) Méthode selon la revendication 1 ou 2, caractérisée en ce que l'indice de mouillabilité est obtenu par la relation: RMN = 1og1O SM0 o SM, est la surface poreuse mouillée par l'eau et SM, est la surface poreuse mouillée 20 par l'huile lorsque le milieu poreux est saturé en eau et huile. ) Méthode selon l'une des revendications précédentes, caractérisée en ce que l'on détermine l'indice de mouillabilité par les opérations suivantes: a) on mesure les temps caractéristiques de relaxation de l'échantillon saturé en eau; b) on mesure les temps caractéristiques de relaxation de l'huile dans l'échantillon en 25 présence d'eau, dans une zone proche de la saturation (Swir) de l'échantillon; - 14 2844355 c) on mesure les temps caractéristiques de relaxation de de l'eau dans l'échantillon en présence d'huile, dans une zone proche de la saturation résiduelle (Sor) d) on mesure les temps de relaxation de l'échantillon dans un état o son point de saturation à 100% en huile est atteint; et e) on combine les mesures des différentes temps de relaxation pour obtenir le dit indice. 6) Méthode selon la revendication 5, caractérisée en ce que l'on mesure les temps de relaxation des étapes a) à c) après application à l'échantillon d'une centrifugation. 7) Méthode selon la revendication 5, caractérisée en ce que l'on mesure les temps de relaxation de l'étapes d) après déplacement forcé des fluides dans l'échantillon placé 10 dans une cellule de confinement. 8) Méthode selon l'une des revendications précédentes, caractérisée en ce que l'on choisit une huile dont le temps de relaxation intrinsèque (TB) est le plus grand possible et le plus proche de l'eau. 9) Méthode selon la revendication 5, caractérisée en ce que les temps 15 caractéristiques de relaxation sont ceux qui correspondent soit au maxima des courbes de saturation soit à des valeurs moyennes des dites courbes.