FR2900239A1 - Methode pour determiner la viscosite de fluides par resonances magnetique nucleaire - Google Patents

Abstract

- Méthode pour déterminer la viscosité d'un fluide à partir de mesures par résonance magnétique nucléaire.- On définit une courbe de corrélation entre viscosité et moyenne logarithmique de distribution de temps de relaxation transversale. Pour se faire, on prend en compte une distribution type de temps de relaxation transversale T2, en fonction du fluide étudié. Cette distribution est complète. On en déduit une valeur de viscosité. On modifie ensuite cette distribution en fonction de caractéristiques de l'appareil de mesure, pour déterminer une distribution équivalente à une distribution qui serait obtenue par mesure RMN réalisée avec cet appareil. On calcule la moyenne logarithmique de cette distribution équivalente. Puis l'on simule plusieurs baisses de la viscosité en translatant la distribution type, pour obtenir plusieurs couples de valeurs viscosité / moyenne logarithmique, dont on déduit la courbe de corrélation. Cette courbe permet alors de déterminer la viscosité du fluide à partir de mesures par résonance magnétique nucléaire.- Application à l'exploration pétrolière.

Description

10 La présente invention concerne une méthode pour déterminer la viscosité

d'un fluide à partir de mesures RMN, quelque soit la viscosité du fluide. En particulier, la méthode permet de prendre en compte l'effet d'atténuation des composantes courtes dans une distribution de temps de relaxation transversale T2. En cela, la méthodologie proposée tient compte : 15 - des caractéristiques, telles que le temps inter écho de l'appareil RMN mesurant des temps de relaxation transversale ; - de la forme de la distribution des temps de relaxation transversale T2 du fluide dont on souhaite déterminer la viscosité. 20 L'utilisation des mesures RMN (Résonance Magnétique Nucléaire) dans l'industrie pétrolière est largement répandue depuis de nombreuses années. Les appareils de mesure à résonance magnétique nucléaire (RMN), comportent classiquement des aimants et des bobinages placés dans l'entrefer des aimants. L'échantillon, sur lequel on souhaite effectuer des mesures, est disposé à l'intérieur des 25 bobinages. On applique à cet échantillon, soumis au champ magnétique des aimants, un champ électrique oscillant (généralement dans la gamme des radiofréquences (RF)) en connectant le bobinage à un générateur de signal à fréquence variable, et l'on enregistre la réponse de l'échantillon aux signaux excitateurs. Pour stabiliser le champ magnétique créé par les aimants, on stabilise leur température en établissant une circulation forcée d'air et en interposant entre les aimants et le bobinage excitateur, un circuit de refroidissement où circule un fluide. Dans le cadre d'opérations de mesure, il est utile de porter l'échantillon au préalable à une température définie. C'est le cas par exemple quand il s'agit d'étalonner en laboratoire des mesures faites par des sondes de diagraphie par RMN à différentes profondeurs de puits d'exploration de zones souterraines. Présentation de l'art antérieur La mesure de relaxation RMN consiste à exciter les spins des protons du fluide contenu dans l'échantillon étudié en leur imposant une aimantation, et à observer leur retour à l'équilibre. Deux mécanismes permettent ce retour à l'équilibre : - Un échange d'énergie entre les spins et leur environnement, cela s'appelle la relaxation longitudinale, et est caractérisé par un temps Tl. - Une perte de cohérence entre les différents spins, correspondant à une relaxation transversale, et caractérisée par un temps T2. Ces mesures RMN permettent d'obtenir une courbe de variation de l'aimantation en fonction du temps. De ces courbes on extrait les constantes de temps TI et T2. La mesure du temps TI peut être effectuée par exemple par la séquence dite d'inversion récupération (IR). La mesure du temps T2 peut être effectuée par exemple par l'utilisation de la séquence dite CPMG. Cette séquence est illustrée sur la figure 1 et elle est décrite par exemple dans le document suivant : H. Carr, E. Purcell, (1954) Physical Review 94: 630. Dans le cas ou la relaxation est décrite par un seul temps (décroissance mono exponentielle), la théorie proposée par Abragham (Abragham, A. (1961). Principles of Nuclear Magnetism . Clarendon Press Oxford.) permet de relier ces temps aux temps de corrélation des mouvements : 1 - J rC + 4rC +w rc 1+4ufôr2 ~ 1 =C(3 +5 rC + rc T2 2 21+w-ârc 1+4urôr~ ~

C=3/~o\Zh274 10 ~42r1 b6 où r~ est le temps de corrélation des interactions dipolaires entre deux protons de rapport gyromagnétique y et distants de b. h est la constante de Planck divisée par 27r, uo est la perméabilité magnétique du vide et co0 est la fréquence de Larmor (fréquence de résonance). Selon cette théorie, on relie la relaxation RMN de deux protons sur une sphère de rayon a, à la viscosité rl du milieu dans lequel évolue cette sphère par la relation de Stokes Einstein : où k est la constante de Boltzmann et TK la température en Kelvin. Le temps de 10 corrélation rC est donc proportionnel à la viscosité. La figure 2 illustre l'évolution des temps de relaxation TI et T2 en fonction de Tc selon la théorie proposée par Abragham : pour de faible valeur de rC, les temps de relaxation TI et T2 sont égaux. Ensuite, TI présente un minimum T Imrn. La courbe de Tl remonte ensuite. Par contre T2 est inversement proportionnel à rc, sur une très grande 15 gamme de temps de corrélation, il n'atteint un plateau que pour les temps de corrélation très longs, pour lesquels il apparaît un effet solide (ES). Ce type d'expérience de mesures RMN est devenu important en pétrophysique pour caractériser, par exemple, la taille des pores dans les roches : il a été montré que les spins d'un échantillon saturé en eau relaxent avec une distribution de temps T2, c'est- 20 à-dire que tous les protons ne perdent pas leur aimantation en même temps. Une telle distribution est représentée figure 3. Elle indique le nombre de protons (NP) qui 6Dra, kTK Drot 87rr7a 3 (1)

relaxent avec un temps de relaxation transversal T2 donné. Cette distribution peut être interprétée comme la distribution des tailles de pores dans la roche considérée. De plus, si l'on mesure les corrélations entre TI et T2, on peut savoir si le liquide est de l'eau ou un hydrocarbure.

Il est également connu d'utiliser la mesure de T2 pour déterminer la viscosité des huiles brutes. D'ailleurs, la détermination de la viscosité des huiles brutes a été la première application de la relaxation RMN dans l'industrie pétrolière, après que des travaux aient démontré une corrélation globale entre la viscosité des huiles et le temps de relaxation mesuré. Dans le document suivant : Morris, C., (1994), Hydrocarbon saturation and viscosity estimation from NMR logging in the belridge diatomite , SPWLA 35th Annual Logging Symposium, Morris et al. ont mesuré les distributions de temps de relaxation transversal (T2) de plusieurs huiles. Ils ont remarqué que les huiles raffinées présentent des distributions de T2 étroites, alors que celles des pétroles bruts s'étalent sur plusieurs décades. Ils remarquent également que plus l'huile est visqueuse, plus sa distribution est large et plus le poids des composantes courtes est important. L'existence de ces distributions a amené au choix d'un paramètre représentatif : la moyenne logarithmique de la distribution T2,ni. Dans le document suivant : Straley, C., Rosini, D., Vinegar, H., Tutunjian, P., and 20 Morris, C. (1994). Core analysis by low field NMR . In SCA Proceedings. Straley a obtenu une corrélation empirique en faisant des mesures sur des huiles mortes: 29877 17 Il propose un ajustement avec une pente 0.9: 1200 T2ml = 0.9 où T2m1 est en milliseconde, TK est la température en Kelvin et 77 la viscosité en centipoises. T _ 1200TK = 4.03 TK 2ml 25 Cependant, d'autres auteurs ont montré que la mesure de laboratoire peut être influencée par la présence d'oxygène. On montre que le préfacteur de la corrélation est modifié, la pente restant inchangée. En effet, l'oxygène moléculaire est un centre paramagnétique qui accélère la relaxation RMN. Les expériences réalisées sur des échantillons désoxygénés conduisent à des temps de relaxation T2 plus longs, confirmant l'effet accélérateur de l'oxygène sur la relaxation RMN. Cette problématique n'intervient que pour les mesures de laboratoire, sur des échantillons non préservés. En effet, l'huile in-situ ne contient pas d'oxygène. Problème de détection des composantes courtes : Une étude de LaTorraca et al. sur des huiles lourdes montre un effet important. Cette étude est présentée dans le document suivant : LaTorraca, G., Dunn, K., Webber, P., and Carlson, R. (1998). Low field NMR determinations of the properties of heavy oils and water-in-oil emulsions. Magnetic Resonance Imaging, 16 :659ù662.

LaTorraca remarque que lorsqu'on augmente le temps inter écho (TE) dans la séquence de mesure CPMG (figure 1), la valeur moyenne T2i,,l du temps de relaxation détecté augmente, et l'amplitude totale du signal détecté diminue. En effet, lorsque le temps inter écho augmente, les composantes courtes de relaxation ne sont plus détectée correctement et sont sous estimées. Il propose une corrélation empirique pour tenir compte de cet effet : T2ntl ù (TE + 0.5) 298 où TE et T2,,,I sont en milliseconde. D'autres auteurs (Bryan et al.) ont proposé la corrélation suivante: 1.15 17 RHI4.SSI 2 ml où R-II est un indice exprimant la perte de signal liée aux composantes courtes. On retrouve une telle corrélation dans le document suivant par exemple : 2200 + 470TE2 TK x ~-

Bryan, J., Kantzas, A., and Bellehumeur, C. (2002). Viscosity predictions for crude oils and crude oil emulsions using low field NMR. Society of Petroleum engineers.

Problème des distributions: Hirasaki et al. ont compilé de nombreuses données de la littérature dans le document suivant : Hirasaki, Lo, and Zhang (2003). NMR properties of petroleum reservoir fluids. Magnetic Resonance Imaging, pages 269-277.

Ils mettent en évidence la non validité des corrélations T2-viscosité pour les huiles très visqueuses. Ils observent également le comportement du temps de relaxation longitudinal Tl. Ce temps présente un plateau inattendu pour les huiles visqueuses. Ces auteurs ouvrent la voie en proposant d'utiliser un effet de distribution de temps de relaxation pour expliquer ce comportement. En effet, dans le cas de Tl, des travaux ont montré que la prise en compte d'une distribution modifie fortement la forme de la courbe TI f(rr) (Figure 2). Cependant, l'effet d'une distribution n'a pas été envisagé pour T2. Finalement, si de nombreuses corrélations sont disponibles pour faire le lien entre les temps de relaxation mesurables par RMN et la viscosité de fluides, la littérature met en évidence deux problèmes majeurs : - Problème de détection des composantes courtes des temps de relaxation transversale T2 : Ce problème est directement lié au type d'appareil utilisé pour effectuer la mesure RMN, et plus particulièrement au temps inter écho (TE) spécifique de cet appareil. Lorsqu'on augmente le temps inter écho dans la séquence de mesure CPMG (figure 1), les premières valeurs de temps de relaxation transversale T2, c'est-à-dire les temps courts, ne sont pas mesurés. En effet, la première valeur mesurée correspond au temps inter écho comme l'illustre la figure 1. Les mesures de temps de relaxation transversale T2 sont réalisées au temps TE, puis 2.TE, puis 3.TE, ... Ces mesures sont illustrées par des carrés noirs sur la figure 1. PDCC indique l'intervalle de temps pour lequel il y a un problème de détection des composantes courtes. Du fait de cette perte de signal, la valeur de la moyenne logarithmique (T2in1) de la distribution des temps T2 augmente, et l'amplitude totale du signal détecté diminue. Il est donc important, surtout pour des fluides visqueux tels que les huiles lourdes, de ne pas sous estimer ces composantes courtes, pour ne pas surestimer la moyenne logarithmique (T2,,I) de la distribution. - Problème lié à la forme de la distribution des temps de relaxation transversale T2 Des études ont montré que plus l'huile est visqueuse, plus sa distribution est large et plus le poids des composantes courtes est important. Ceci est dû à la grande diversité de composés chimiques présents dans les pétroles bruts visqueux. De plus les huiles très visqueuses comportent des fractions de composés lourds (asphaltènes) très importantes, ce qui se traduit sur la distribution de T2 par une fraction de composantes courtes importante. Plus la distribution de temps de relaxation est large plus sa moyenne logarithmique (T2,,I) est élevée. La linéarité classiquement observée entre cette moyenne et la viscosité ne l'est plus, comme l'illustre la figure 7. Il est donc impératif de prendre en compte cet effet de distribution pour déterminer de façon précise et fiable la viscosité en fonction des mesures RMN.

La méthode selon l'invention permet de déterminer la viscosité d'un fluide, tel qu'une huile brute, à partir de mesures RMN, quelque soit la viscosité du fluide, tout en prenant en compte l'effet d'atténuation des composantes courtes dans la distribution des temps de relaxation transversale T2. En cela, la méthodologie proposée tient compte : - des caractéristiques, telles que le temps inter écho, de l'appareil RMN réalisant la mesure pour estimer les temps de relaxation transversale T2 ; - de la forme de la distribution des temps de relaxation transversale T2 d'un échantillon étudié.

La méthode selon l'invention L'invention concerne une méthode pour déterminer la viscosité d'un fluide à partir de mesures par résonance magnétique nucléaire, dites mesures RMN, réalisées par un appareil de mesures RMN, caractérisée en ce qu'elle comporte les étapes suivantes : • on définit une courbe de corrélation entre viscosité et moyenne logarithmique de distribution de temps de relaxation transversale en réalisant les étapes suivantes : a) on prend en compte une distribution type de temps de relaxation transversale T2, en fonction dudit fluide, et comportant des composantes courtes desdits temps de relaxation ; b) on prend en compte des caractéristiques dudit appareil parmi au moins les caractéristique suivantes : un bruit aléatoire, un temps inter écho, un temps de relaxation transversale minimal à partir duquel ledit appareil ne peut plus effectuer de mesures ; c) on modifie ladite distribution type en fonction desdites caractéristiques dudit appareil, pour déterminer une distribution équivalente à une distribution qui serait obtenue par mesure RMN réalisée avec ledit appareil ; d) on détermine à partir de ladite distribution type une valeur de viscosité, et on détermine une moyenne logarithmique de ladite distribution équivalente ; e) on translate ladite distribution type, puis l'on réitère les étapes c) et d) pour obtenir un nouveau couple de valeurs viscosité /moyenne logarithmique ; f) on obtient plusieurs couples permettant d'inférer ladite courbe de corrélation entre viscosité et moyenne logarithmique, en translatant plusieurs fois ladite distribution type selon e) ; • on réalise sur ledit fluide une mesure RMN à partir dudit appareil pour déterminer une distribution de temps de relaxation transversale, dont on détermine une moyenne logarithmique caractéristique dudit fluide ; et • on détermine ladite viscosité dudit fluide à l'aide de ladite courbe de corrélation et de ladite moyenne logarithmique caractéristique. 8

Selon la méthode, la distribution type peut être issue d'une mesure RMN réalisée sur un échantillon de fluide, à une température telle que la plus petite valeur de temps de relaxation T2 mesurée est supérieure audit temps inter écho d'un facteur au moins égal à cinq. Cette distribution type peut également être issue d'une modélisation analytique dans laquelle on définit une largeur, une position et une asymétrie définissant ladite distribution. On peut déterminer la distribution équivalente en réalisant les étapes suivantes : - on tronque ladite distribution type en mettant à zéro des valeurs de ladite distribution type pour des temps de relaxation transversale T2 inférieurs audit temps de relaxation transversale T2 minimal ; - on reconstruit, à partir de ladite distribution tronquée, un signal RMN qui serait obtenu dans des conditions réelles, en prenant en compte ledit bruit et ledit temps inter écho ; - on détermine ladite distribution équivalente en réalisant une décomposition exponentielle dudit signal RMN reconstruit. Enfin, selon la méthode, le fluide peut être une huile brute. Présentation sommaire des figures Les caractéristiques et avantages de la méthode selon l'invention, apparaîtront plus clairement à la lecture de la description ci-après d'exemples non limitatifs de mise en oeuvre, en se référant à la figure annexée et décrite ci-après. - La figure 1 illustre la séquence dite CPMG pour mesurer les temps de relaxation transversale T2 ; - La figure 2 illustre l'évolution des temps de relaxation Ti et T2 en fonction de Tc selon la théorie proposée par Abragham, avec illustration de l'effet solide (ES) ; - La figure 3 montre une distribution des temps de relaxation transversale T2, pour laquelle il n'y a aucune perte des composantes courtes (DT2) ; - La figure 4 illustre l'effet solide sur la distribution synthétique de la Figure 3 (DT2`') ; -La figure 5 représente un exemple de signal RMN reconstruit (M(t)) à partir de la distribution DT21r de la figure 4 ; - La figure 6 représente un exemple de décomposition du signal illustré sur la figure 5 (DT2 P) ; -La figure 7 illustre le résultat de la méthode selon l'invention avec la courbe déduite des couples (T2in1, i) ainsi qu'une comparaison avec des données expérimentales disponibles dans la littérature ; - La figure 8 illustre les étapes de la méthode permettant de déterminer une courbe de corrélation entre la viscosité 17 et la moyenne logarithmique T2,n1• Description détaillée de la méthode La méthode permet de déterminer une courbe de corrélation entre une mesure représentative de la distribution des temps de relaxation transversale T2 et la viscosité d'un fluide. Plus précisément, à partir d'un temps inter écho donné (TE), et un fluide 1.5 visqueux donné (caractérisé par sa distribution des temps de relaxation T2), la méthode détermine un ensemble de couples (T2,,,lexp,r1), c'est-à-dire des couples comportant une valeur de moyenne logarithmique d'une distribution T2 (T2m1e,,p ), et une valeur de viscosité (77) correspondante. La méthode comporte sept étapes, illustrées sur le schéma de la figure 8 et 20 détaillées ci-après. 1) Obtention d'une distribution complète des temps de relaxation (DT2). La première étape consiste à définir une distribution type des temps de relaxation transversale T2, non altérée par des problèmes expérimentaux, c'est-à-dire une 25 distribution type pour laquelle il n'y a aucune perte des composantes courtes. Cette distribution des temps de relaxation transversale T2 non altérée est notée DT2. On qualifie cette distribution de type , car elle est spécifique de la famille de produits (fluides) que l'on souhaite étudier. Pour construire une telle distribution DT2, deux solutions peuvent être envisagés :10 - Lorsqu'un échantillon du fluide (huile lourde par exemple) est disponible, on réalise une mesure RMN de laboratoire à une température telle que tous les temps de relaxation transversale T2 soient bien déterminés (pas de perte des composantes courtes). Par exemple, la plus faible valeur de T2 mesurée devra être telle que T2 = Q.TE avec (3 supérieur où égal à 5. Classiquement (3 est pris égal à 5. Si ce n'est pas le cas, si des valeurs T2 ne peuvent être mesurées, on chauffe l'échantillon. Le fait de chauffer l'échantillon diminue la viscosité du fluide et décale donc la distribution des temps de relaxation vers des valeurs élevées de T2. On chauffe l'échantillon jusqu'à ce que la mesure RMN permette de mesurer une valeur de T2 telle que T2 = - Dans le cas où la détermination expérimentale de la distribution n'est pas possible, on peut utiliser un modèle analytique reproduisant les formes de distribution des temps de relaxation transversale T2 rencontrées habituellement. Par exemple, pour des huiles pétrolières, on génère une distribution DT2 ayant les caractéristiques principales d'une distribution expérimentale, c'est-à-dire : 2 à 3 décades de large et un épaulement aux temps courts. Pour cela, une distribution est définie comme suit pour un échantillonnage de T2 régulier sur une échelle logarithmique sur 4 décades : - -2 log( T2 ) +B T2m1 26-2 La largeur de la distribution générée peut être ajustée par le paramètre 6, sa position par le paramètre T2m1 et son asymétrie par le paramètre pa. Un exemple d'une telle distribution synthétique est donné figure 3, où u = 0.2 et pa=113. Le calcul de la moyenne logarithmique T2,n1 donne : 0,32. 36 8 B = ' ) Calcul de la viscosité (iq) La distribution DT2, étant générée de façon à être totalement résolue, c'est-à-dire avec toutes les composantes courtes de T2, on peut en tirer la moyenne logarithmique T2m1 exempte d'artéfacts de mesure et on peut utiliser les corrélations existantes pour 5 déterminer la viscosité correspondant à cette distribution DT2, notée ri. Elle est ainsi définie : avec Q et w des constantes réelles, T2.1, la moyenne logarithmique de la distribution DT2. 10 Les valeurs de Q et w peuvent être obtenues dans la littérature. Par exemple 1200 Straley propose T2m1 = 0,9 . Ces constantes peuvent également être déterminée par une 77 série de mesures RMN et viscosité de laboratoire ayant pour but de caler le modèle ci dessus. 3) Modification de la distribution DT2 pour prendre en compte l'effet solide 15 L'effet solide (ES) est illustré sur la figure 2. Cet effet signifie que les valeurs de temps de relation transversale T2 ne peuvent être inférieures à une certaine valeur, notée T2/imite, fonction uniquement du fluide étudié. Cette valeur limite est donnée par Abragham, par exemple dans le document suivant : Abragham, A. (1961). Principles 20 of Nuclear Magnetism . Clarendon Press Oxford. Cette valeur T2limtte correspond à l'apparition d'un signal RNVIN type solide qui ne peut être mesuré par le spectromètre de l'appareil de mesure RIvIN. T211mtte est donnée par la formule: 12 T2lim ite ( 1 \ 2 ù w0 \, Tlmin / wo est la fréquence de Larmor. La valeur de Tlmin est obtenue par minimisation de la fonction T, = f (r) (équation 1), donnée par Abragham dans le document cité ci-dessus. Ainsi en dessous de ce temps limite, la distribution est nulle. Un exemple d'une telle distribution est donné par la figure 4, qui illustre un exemple de l'effet solide à partir de la distribution synthétique de la Figure 3. On a choisi b=2,5.10-10m (unique paramètre qui apparaît dans l'équation 1), ainsi T2limtté 30 s. Cette nouvelle distribution (DT2 avec prise en compte de T2limite) est qualifiée de tronquée . Elle est notée : DT2'' 4) Reconstruction d'un signal RMN (M(t)). A partir de cette nouvelle distribution (DT2t'), on recalcule le signal RMN qui serait obtenu dans les conditions réelles. Pour ce faire, on génère un signal de la même qualité qu'un signal réel, en prenant en compte le bruit et le temps inter écho (TE), 15 caractéristiques de l'appareil de mesure RMN utilisé. On échantillonne le temps comme dans une séquence CPMG (voir figure 1) : t = n*TE. Le signal RMN recomposé est alors défini par: M(t) = D2 2" (TZ;) exp ût + br Tz; où t est le temps et br est une fonction aléatoire qui permet d'ajouter un bruit au 20 signal ainsi généré. Un exemple de signal RMN reconstruit à partir de la distribution DT2t' de la figure 4 est donné figure 5 : le rapport signal sur bruit est de 100 et le temps inter écho choisi est TE=400 s. Une séquence de type CPMG a été utilisée, mais la méthode selon l'invention reste valable pour toute autre séquence permettant la mesure de ces temps de relaxation. 5) Décomposition multi exponentielle du signal (DTf') Le signal RMN ainsi généré, M(t), est ensuite décomposé comme s'il s'agissait d'un signal expérimental. On recherche une distribution que l'on appelle 1.0 25 expérimentale , et que l'on note DT2exp. En fait il s'agit d'une distribution pseudo expérimentale : elle correspond à la distribution que l'on obtiendrait à partir de mesures réalisées par l'appareil RMN, caractérisé par son bruit, son temps inter écho et son temps de relaxation transversale minimal. Il s'agit donc d'une distribution équivalente . Selon un mode de réalisation, cette décomposition consiste à ajuster le signal recomposé par l'expression: M(t) _ E DT2exp (T2i)exp(ù nT/T I avec i= l,...,80 l où n est le nombre de points en temps espacés du temps inter écho TE lors de la reconstruction d'un signal RMN (étape 4), T21 est le ième échantillon de 80 valeurs prédéfinies espacées sur une échelle logarithmique. Le calcul de la distribution DT2exp aux valeurs T21 est bien connu pour être un problème mathématique mal posé. Afin d'éviter les oscillations de la solution, on utilise une méthode de régularisation, dans laquelle la fonction x2 suivante est minimisée: n 80 x2 = E M(ti ) ù DT2e" (T2; )exp( ;=i où a est le paramètre de régularisation choisi de manière à obtenir un lissage optimum, c'est la valeur la plus grande qui produit le plus faible écart entre le modèle et le signal, et les tj représentent un jeu d'intervalles de temps espacés sur une échelle logarithmique. On obtient une distribution que l'on appelle expérimentale , et que l'on note DT2exp. Un exemple de décomposition du signal illustré sur la figure 5 est donné figure 6. La moyenne logarithmique T2mtexp de cette distribution est extraite. Sur la figure 6, elle est égale à 0,76. D'autres méthodes de décomposition multi exponentielle du signal RMN peuvent être envisagées selon l'invention. \ z 80 + a 2 E (DT2exp (T2, ))2 i=125 6) Translation de la distribution complète DT2 (TRANS). Le couple de valeurs ri (obtenu à l'étape 2) et T2mrexp (obtenue à l'étape 5) sont stockées dans un graphique dont l'axe des abscisses représente la viscosité (1) et l'axe des ordonnées la moyenne logarithmique d'une distribution T2 (TZmi ). Afin d'obtenir un autre couple de valeurs, on translate la distribution DT2. Pour ce faire on modifie les abscisses de cette distribution sans modifier les ordonnées. A chaque ordonnée on associe non plus l'abscisse T2i, mais l'abscisse translatée TZ;a"S, avec : TZ,""s =, eXP(log(T2r) ù P)) avec p=0.1 par exemple

Ceci permet de translater la distribution sur l'axe des T2 sur une échelle logarithmique sans modifier la forme de la distribution. On obtient une nouvelle distribution DT2 à laquelle on applique les étapes 2,3,4,5 et 6. 7) Détermination d'une courbe de corrélation r~ / Ti 2,n(CORS On réitère la translation décrite à l'étape 6 précédente, ainsi que la détermination d'un couple (T2,nlexp, ri) autant de fois que nécessaire pour obtenir une bonne

approximation d'une loi , permettant de relier T2,nt à ri : r)=Y(T2",r) Le fait de translater la distribution permet d'obtenir un jeu de données où l'on a pour chaque position de la distribution, la viscosité i (calculée à l'étape 2) et la valeur T2,nrexp (calculée à l'étape 5). Cette itération est arrêtée lorsque la distribution DT2`r est complètement tronquée par l'effet solide (DT2r (Tu) = 0, V T2i ). Selon l'exemple de la figure 7, on peut par exemple estimer cette loi de corrélation `Y , à partir des différents couples (T2,nrexp, i1) générés, par la formule suivante : 1+ CT2 "Z 11=W(T2mr)=KJ 1+C où K=5300, C=1, nl =1/0.9 et n2=1/0.5 dans l'exemple de la figure 7. Cette figure 7 illustre le résultat de la méthode selon l'invention : la courbe 2 représente la fonction `l' déduite des couples (T2,n/, ri) ainsi obtenus. Sur cette même figure on compare ces résultats aux données expérimentales disponibles dans la littérature ou aux autres corrélations disponibles : les ronds représentent les données expérimentales issues de Hirasaki et al., et la courbe 1 illustre la corrélation proposée par Straley et al.. La courbe 'Y ainsi obtenue peut être utilisée pour prédire la viscosité du fluide étudié. On rappelle que l'obtention de cette courbe dépend : - du produit (le fluide), par l'intermédiaire de la distribution de temps de relaxation T2 choisie ; 10 - de l'appareil de mesure par l'intermédiaire du niveau de bruit des mesures et du temps inter écho (TE) lors de la génération de signaux composites ; La courbe ainsi créée peut être utilisée selon différentes façons : En laboratoire 15 Pour effectuer de nombreuses mesures sur des échantillons pétroliers par exemple, la méthode peut être appliquée de la façon suivante : - On construit une courbe `P (q = (T,.!)) à partir d'un échantillon, puis l'on réalise des mesures RMN sur tous les échantillons.On déduit de ces mesures et de la courbe la viscosité de chacun des fluides contenus dans chacun des échantillons. 20 Cette méthode sera notamment utilisée lorsque aucune courbe n'est disponible pour le type d'appareil de mesures RMN utilisé, ou si les échantillons présentent a priori des distributions de temps relaxation T2 bien particulières. - On peut également utiliser directement une courbe déjà établie pour le type d'appareil de mesures RMN utilisé. 25 In situ, au sein d'un réservoir pétrolier par exemple Pour déterminer la viscosité des fluides au sein d'un réservoir pétrolier, on réalise des mesures RMN par diagraphies.5 .7 - On peut alors construire une courbe `Y (r~ = Y(T2ni,)) à partir d'un échantillon prélevé, ou à partir de formules analytiques déduites des connaissances générales du réservoir, puis l'on réalise des mesures RMN par diagraphie. On déduit de ces mesures et de la courbe `Y la viscosité de chacun des fluides contenus dans le réservoir. Cette méthode sera notamment utilisée lorsque aucune courbe n'est disponible pour le type d'appareil de mesures RMN utilisé, ou si les fluides présentent a priori des distributions de temps relaxation T2 bien particulières.

- On peut également utiliser directement une courbe déjà établie pour le type d'appareil de mesures RMN utilisé. 1.0

Claims (5)

REVENDICATIONS

1) Méthode pour déterminer la viscosité d'un fluide à partir de mesures par résonance magnétique nucléaire, dites mesures RMN, réalisées par un appareil de mesures RMN, caractérisée en ce qu'elle comporte les étapes suivantes : • on définit une courbe de corrélation entre viscosité et moyenne logarithmique de distribution de temps de relaxation transversale en réalisant les étapes suivantes : a) on prend en compte une distribution type de temps de relaxation transversale T2, en fonction dudit fluide, et comportant des composantes courtes desdits temps de relaxation ; b) on prend en compte des caractéristiques dudit appareil parmi au moins les caractéristique suivantes : un bruit aléatoire, un temps inter écho, un temps de relaxation transversale minimal à partir duquel ledit appareil ne peut plus effectuer de mesures ; c) on modifie ladite distribution type en fonction desdites caractéristiques dudit appareil, pour déterminer une distribution équivalente à une distribution qui serait obtenue par mesure RMN réalisée avec ledit appareil ; d) on détermine à partir de ladite distribution type une valeur de viscosité, et on détermine une moyenne logarithmique de ladite distribution équivalente ; e) on translate ladite distribution type, puis l'on réitère les étapes c) et d) pour 20 obtenir un nouveau couple de valeurs viscosité / moyenne logarithmique ; f) on obtient plusieurs couples permettant d'inférer ladite courbe de corrélation entre viscosité et moyenne logarithmique, en translatant plusieurs fois ladite distribution type selon e) ; • on réalise sur ledit fluide une mesure RMN à partir dudit appareil pour 25 déterminer une distribution de temps de relaxation transversale, dont on détermine une moyenne logarithmique caractéristique dudit fluide ; et • on détermine ladite viscosité dudit fluide à l'aide de ladite courbe de corrélation et de ladite moyenne logarithmique caractéristique.

2) Méthode selon la revendication 1, dans la quelle ladite distribution type est issue d'une mesure RMN réalisée sur un échantillon de fluide, à une température telle que la plus petite valeur de temps de relaxation T2 mesurée est supérieure audit temps inter écho d'un facteur au moins égal à cinq.

3) Méthode selon l'une des revendications précédentes, dans la quelle ladite distribution type est issue d'une modélisation analytique dans laquelle on définit une largeur, une position et une asymétrie définissant ladite distribution.

4) Méthode selon l'une des revendications précédentes, dans laquelle on détermine ladite distribution équivalente en réalisant les étapes suivantes : - on tronque ladite distribution type en mettant à zéro des valeurs de ladite distribution type pour des temps de relaxation transversale T2 inférieurs audit temps de relaxation transversale T2 minimal ; - on reconstruit, à partir de ladite distribution tronquée, un signal RMN qui serait obtenu dans des conditions réelles, en prenant en compte ledit bruit et ledit temps inter écho ; - on détermine ladite distribution équivalente en réalisant une décomposition exponentielle dudit signal RMN reconstruit.

5) Méthode selon l'une des revendications précédentes, dans la quelle ledit fluide est une huile brute.