FR 2909448 A1 20080606 FR 0610581 A 20061201 10 La présente invention concerne le domaine des mesures pétrophysiques sur des échantillons de roches. Plus particulièrement, l'invention concerne le domaine des procédés permettant de caractériser la distribution tridimensionnelle de la perméabilité absolue d'un échantillon hétérogène issu d'une roche. État de la technique On connaît différentes approches pour caractériser une carte tridimensionnelle (3D) de perméabilité d'un échantillon d'une roche. Identification d'une carte 3D basée sur une identification "manuelle" 20 Cette approche consiste à utiliser un test de traceurs, c'est-à-dire un test miscible mais sans contraste de viscosité, pour obtenir une carte de perméabilité. Cette technique est décrite dans les documents suivants : - Dabbouk, C., Ali, L., Williams, G., I3eattie, G. :"Waterflood in vuggy layer in a Middle East reservoir ù Displacement physics understood", SPE 78530, Abu Dhabi 25 International Petroleum Exhibition and Conference, 13-16 October 2002. -Olivier, P., Cantegrel, L., Laveissière, J., Guillonneau, N. : "Multiphase flow behaviour in vugular carbonates using X-Ray", SCA 2004-13, Society of Core Analysts Symposium, Abu Dhabi, 5-9 October 2004. 1 15 2909448 2 Une ou plusieurs lois K((D) sont utilisées pour déduire d'une carte de porosité, une carte de perméabilité servant à la simulation de l'expérience de traceurs. Ces lois K((D) sont alors ajustées de façon manuelle jusqu'à obtenir un bon accord avec l'expérience. L'avantage de cette approche est sa simplicité. En revanche, un inconvénient majeur est qu'il existe un très 5 grand nombre de solutions à ce problème. Or, un calage manuel ne donne accès qu'à une seule solution. Identification d'une carte 3D basée sur une inversion "mathématique" Cette technique s'inscrit toujours dans le cadre des tests de traceurs. Selon cette méthode, il est possible de suivre, par exemple par tomographie, le déplacement d'un front à différents temps. A partir de cette information, la carte de perméabilité est solution d'un problème non linéaire avec des conditions aux limites connues. Cette méthode est décrite dans : - Zhan, L., Yortsos, Y.C. :"A direct method for the identification of the permeability 15 field of an anisotropie porous medium", SPE 62976, Annual Technical Conference and Exhibition, Dallas, 1-4 October 2000. Dans la mesure où les conditions aux limites (profil de pression ou perméabilité sur les frontières) et la carte de porosité sont connues, cette technique conduit à une solution unique, qui n'est valide que lorsque les contrastes de perméabilité sont faibles. La mise en oeuvre de 20 cette technique reste donc limitée en pratique, car les conditions aux limites ne sont pas toujours connues. Cette méthode ne s'applique pas aux milieux très hétérogènes. Profil ID basé sur un déplacement visqueux La mise en oeuvre d'un déplacement visqueux pour obtenir des informations sur le 25 profil ID de perméabilité a été proposée par Fincham et Gouth : - Fincham, A., Gouth, F. :"Improvements of coreflood design and interpretation using a new software", SCA 2000, Society of Core Analysts Symposium, Abu Dhabi, 9-12 October 2000. Ces auteurs proposent d'interpréter le signal de pression différentielle d'une injection 30 d'huile visqueuse à fort débit dans un massif initialement saturé en saumure pour déterminer le profil ID de perméabilité à l'huile à Swi. Grâce à cette technique les auteurs réussissent à 2909448 3 mieux caractériser l'hétérogénéité "locale" de leur échantillon et ainsi à affiner l'interprétation des expériences visant à la détermination des courbes de perméabilités relatives. Cependant, cette approche ne permet pas de remonter aux valeurs des perméabilités absolues que ce soit en 1D ou en 3D. L'objet de la présente invention concerne une méthode permettant de caractériser de façon rapide, précise et non destructive, la distribution tridimensionnelle de la perméabilité absolue d'un échantillon hétérogène quelconque. 10 L'invention concerne une méthode de détermination d'une distribution tridimensionnelle de la perméabilité absolue d'un échantillon hétérogène, à partir d'une carte de porosité tridimensionnelle dudit échantillon, dans laquelle on effectue les étapes suivantes : on réalise au moins une expérience de déplacement miscible visqueux au cours de laquelle on détermine une évolution de la pression différentielle DP(t)mes de part et d'autre 15 dudit échantillon ; - on détermine un profil 1D de perméabilité absolue le long dudit échantillon en utilisant ladite évolution de la pression différentielle ; - on construit une première carte 3D de perméabilité (K((D)) à partir de ladite carte 3D de porosité ; 20 - on estime une pression différentielle simulée DP(t)sim, en simulant ::iumériquement un test miscible visqueux à partir de ladite première carte de perméabilité et dudit profil 1D de perméabilité ; - on détermine ladite distribution tridimensionnelle de la perméabilité absolue dudit échantillon en modifiant au moins une fois ladite première carte de perméabilité de façon à 25 minimiser l'écart entre ladite pression différentielle simulée DP(t)sim et ladite pression différentielle mesurée au cours du temps DP(t),r,es. On peut déterminer ladite carte de porosité à partir de mesures statiques non destructives. On peut déterminer ladite carte de porosité à partir de mesures par rayons X en utilisant 30 des fluides de densités différentes saturant ledit échantillon. 5 2909448 4 On peut déterminer ladite carte de porosité à partir de mesures par résonance magnétique nucléaire. On peut déterminer ladite première carte de perméabilité (K(t)) à partir de cosimulations. 5 On peut déterminer ladite première carte de perméabilité (K((D)) à partir d'au moins une relation déterministe entre la porosité et la perméabilité. On peut modifier ladite première carte de perméabilité en appliquant une méthode de déformation graduelle. On peut déterminer en outre un profil de concentration ou des cartes 3D de 10 concentration à des temps différents, que l'on compare à des réponses simulées équivalentes Csim. 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 illustre le dispositif expérimental permettant de mettre en oeuvre la méthode selon l'invention. La figure 2 montre le principe du test miscible visqueux avec un échantillon hétérogène "composite". La figure 3 montre l'exemple d'un signal de pression différentielle DP attendu pour l'échantillon de la figure 2. La figure 4 illustre l'injection du fluide 2 (F2) dans l'échantillon initialement saturé par le fluide 1 (F1). La figure 5 montre une courbe de la pression différentielle (DP) mesurée en fonction du temps (T). La figure 6 illustre un profil de perméabilité (K) calculé en fonction du temps T et un profil de perméabilité théorique. 2909448 5 La figure 7 montre un profil de perméabilité calculé en fonction de la distance et un profil de perméabilité théorique. La figure 8 décrit les étapes utilisées pour déterminer une carte de perméabilité 3D en calant les pressions différentielles. 5 La figure 9 décrit les étapes utilisées pour déterminer une carte de perméabilité 3D en calant les pressions différentielles et les données de concentrations. Description détaillée de la méthode La méthode selon l'invention repose sur une approche intégrée combinant l'acquisition 10 de données expérimentales rapides et non destructives pour l'échantillon et leur interprétation à l'aide d'un outil numérique dédié pour accéder à la carte 3D de perméabilité absolue. Selon la méthode, deux types de mesures expérimentales sont nécessaires : une mesure statique de la carte 3D de porosité de l'échantillon: Cette carte peut être calculée par exemple à partir de mesures RX (mesures par rayons X) 15 obtenues avec des fluides de densités différentes saturant le milieu poreux. D'autres techniques comme la tomographie RMN permettent également d'obteni:r cette carte de porosité. La valeur locale de la porosité peut être alors calculée de manière déterministe en utilisant la différence d'atténuation entre les signaux (RMN et/ou X) obtenus sous les deux conditions de saturation, à savoir pour l'échantillon sec et l'échantillon saturé. Une telle 20 méthode est décrite par exemple dans le document suivant : Maloney, D. R., Wegener, D. C. "Significance absorption coefficients when determining in situ saturations by linear X-ray scans," International Symposium of the Society of Core Analysts, Socorro, USA, 2000. une mesure reposant sur au moins un déplacement miscible: Ce déplacement est réalisé en utilisant un fluide d'injection plus visqueux que le fluide en 25 place (rapport de viscosité par exemple supérieur à 50 cP). Il peut être avantageux de réaliser les tests en injectant le fluide visqueux dans les deux sens pour affiner la précision au niveau de l'interprétation des résultats en terme de perméabilité. Cette mesure peut être la différence des pressions mesurée de part et d'autre de l'échantillon. 2909448 6 La méthode comporte ensuite une étape d'interprétation des résultats. Celle-ci est réalisée en deux étapes distinctes. La première étape permet d'obtenir un. profil ID de perméabilité absolue le long de l'échantillon en utilisant l'évolution de la pression différentielle mesurée au cours des tests de déplacement miscible visqueux. Ce profil 1D de 5 perméabilité est utilisé dans la seconde étape d'interprétation en combinaison avec la carte 3D de porosité pour obtenir la carte 3D de perméabilité. Mise en oeuvre expérimentale 1-Détermination d'une carte tridimensionnelle de la porosité de l'échantillon 10 De nombreuses techniques permettant l'acquisition d'une carte 3D de porosité sont bien connues des spécialistes. On pourra se référer par exemple au document suivant : Maloney, D. R., Wegener, D. C. "Significance absorption coefficients when determining in situ saturations by linear X-ray scans," International Symposium of the Society of Core Analysts, Socorro, USA, 2000. 15 2- Déplacements miscibles visqueux Le dispositif nécessaire est schématisé sur la figure 1. Il est constitué principalement d'une cellule porte échantillon (CE) classique munie ou non de prises de pression locale, d'un système de pompage (SP), et d'un système d'acquisition automatique (SADP) de la pression différentielle (DP) relié à ordinateur (PC). Un des points critiques est l'utilisation d'un capteur 20 différentiel adapté à la gamme de perméabilité de l'échantillon testé. On peut positionner en sortie un conductimètre (CDT) pour suivre la production du fluide injecté en y ajoutant un traceur. Une burette de production (BP) peut être ajoutée en sortie du dispositif. L'évolution du fluide injecté au sein même de l'échantillon peut être obtenue par l'acquisition de mesures RX à l'aide d'un scanner (CTS) de type CT-scanner. Ce dernier dispositif peut en outre 25 permettre de faire l'acquisition de la carte 3D de porosité dans le même environnement. Néanmoins, l'acquisition de ces différents types de mesure (conductimètre en sortie et profils intermédiaires de concentration en fluide injecté) n'est pas obligatoire pour analyser les données finales. L'échantillon monté dans la cellule porte échantillon (CE) est saturé au départ avec un fluide de viscosité modérée (typiquement une saumure) que l'on appelle le fluide 1 dans la suite du document. On mesure alors sa perméabilité absolue moyenne à l'aide de la pompe 2909448 7 (SP) et du capteur différentiel (SADP) en utilisant la loi de Darcy. Cette valeur de perméabilité est importante, car c'est elle qui dimensionne le test miscible visqueux en fonction du matériel disponible. Si on dispose par exemple d'un capteur différentiel (SADP) de 10 bars pleine échelle, on calcule à partir de cette valeur de perméabilité, le débit 5 d'injection du fluide de viscosité élevée (appelé fluide 2 dans la suite) de manière à ce que la pression différentielle induite reste en deçà de la limite du capteur. Avant le démarrage de l'injection du fluide 2, on laisse volontairement un volume mort dans l'embout d'entrée (EE) saturé avec le fluide 1 en place dans l'échantillon. Ce volume est très utile pour mesurer de façon précise le niveau de référence de la pression différentielle du 10 fluide en place dans les conditions de l'expérience (débit). L'injection du fluide 2 se fait à débit imposé de préférence en utilisant une pompe sans pulsation pour éviter les artefacts au niveau du signal de pression. On utilise typiquement des débits d'injection de quelques centimètres cubes par heure. L'évolution de la pression différentielle au cours du test est mesurée de façon continue y compris pendant la période initiale de vidange du volume mort 15 qui fournit précisément le niveau de référence de la pression différentielle pour le fluide 1. On injecte en général environ deux volumes poreux de fluide 2 pour être sûr d'obtenir un régime stabilisé en terme de pression différentielle. L'évolution de cette pression est mesurée de manière continue au cours des tests. Une fois le test miscible visqueux terminé, le fluide 1 est remis en place dans 20 l'échantillon par injection avec la pompe. L'évolution de la pression ou (et) le: conductimètre permet(tent) alors d'évaluer le bon état de remise en place du fluide 1. Le déplacement étant ici fortement instable, on est beaucoup plus sensible aux hétérogénéités et il est nécessaire en général d'injecter un volume plus important pour retrouver l'état de base du fluide 1. On réalise alors un second test qui suit exactement les mêmes procédures que le 25 précédent, mais en procédant à l'injection dans le sens inverse. Interprétation des résultats Profil 1D de perméabilité L'identification du profil ID de perméabilité repose sur l'interprétation des tests 30 précédemment décrits. La viscosité plus élevée du fluide 2 par rapport au fluide 1 déplacé assure en théorie un rapport de mobilité favorable, ce qui se traduit par un déplacement de type piston dans l'échantillon. Le temps de déplacement du front depuis l'entrée vers la sortie 2909448 8 multiplié par le débit d'injection correspond alors au volume poreux de l'échantillon étudié. Ce calcul permet par ailleurs de convertir les temps d'expérience en positions du front de fluide 2 dans l'échantillon au cours de l'injection. La loi de Darcy implique que pour un échantillon parfaitement homogène, l'évolution 5 de la pression se fait de manière strictement linéaire avec une pente unique qui correspond à la valeur de la perméabilité. Pour un échantillon composite avec deux régions de perméabilité différentes, comme illustré sur la figure 2, on observe deux parties linéaires distinctes (figure 3). L'échantillon composite de la figure 2 est constitué d'une région R1 de perméabilité K1 inférieure à la perméabilité K2 de la seconde région R2. 10 La figure 3 illustre un exemple de signal de pression différentielle DP attendu pour l'échantillon composite de la figure 2. On note que la perméabilité est plus faible par rapport à la moyenne dans les 2 premiers tiers de l'échantillon. Dans le sens de l'injection du fluide miscible visqueux (IF), la courbe de pression différentielle DP est représentée en trait épais continu. La première pente est plus prononcée car c'est la région la moins perméable (R1) que 15 le front visqueux traverse en premier. Ensuite, on observe un infléchissement de la pente dès lors que le front visqueux arrive au niveau de la région 2 plus perméable. Dans cet exemple théorique, la cassure entre les deux pentes est atteinte pour un temps qui n'est pas symétrique par rapport aux deux plateaux car la région 1 est plus longue et les porosités ((D) sont égales. Le test inverse (SI) est représenté par une courbe discontinue. Ce test conduit à une signature 20 complètement symétrique en terme de pression différentielle. La première pente suivie est la plus faible tandis que la cassure survient plus tôt, les plateaux restant les mêmes. DP1 représente la pression différentielle lorsque l'échantillon est entièrement saturé par le fluide injecté (fluide 2), et DP2 représente la pression différentielle lorsque l'échantillon est entièrement saturé par le fluide en place (fluide 1). 25 La pression différentielle DP(t) est la différence de pression le long de l'échantillon, et, compte tenu de la nature du test, elle dépend du temps t. Cette pression différentielle est la somme de la pression différentielle dans la partie de l'échantillon saturée par le fluide 2 et de la pression différentielle dans la partie de l'échantillon saturée par le fluide 1 : 30 DP(t)= J Q~2 dx+ J Q~c dx o AK(x) x AK(x) (1 2909448 9 Xf est la position du front et dépend du temps. L est la longueur de l'échantillon. Q est le débit de fluide à travers la surface A, c'est-à-dire la surface de l'échantillon perpendiculaire à la direction d'écoulement. K est la perméabilité. Les viscosités respectives des fluides 1 et 2 sont notées i et 11.2. La figure 4 illustre l'injection du fluide 2 (F2) dans l'échantillon 5 initialement saturé par le fluide 1 (F1). Le fluide 2 est injecté (IF), il sature l'échantillon entre les coordonnées longitudinales de 0 àXf Le fluide 1 sature l'échantillon entre les coordonnées longitudinales de Xf à L (longueur de l'échantillon). Si on se réfère à la méthode de Buckley-Leverett (1942) développée pour les déplacements immiscibles, la position du front est donnée par : 10 Xf = t çbA où 0 est la porosité. La dérivée de Xf par rapport au temps est donc Q/(qA). On en déduit que la dérivée de DP par rapport au temps est : aDP Q2 at A20K(X f)(~z ùPi) L'équation 2 donne la relation entre la position du front et le temps. On peut alors 15 estimer la perméabilité en fonction du temps : 02 1 K(t)= -~u2 ù) aDP at ou de façon équivalente la perméabilité le long de l'échantillon. Cette première étape permet de déterminer le profil ID de perméabilité. Une valeur de perméabilité est représentative de la moyenne des perméabilités sur une tranche de l'échantillon, 20 perpendiculaire à la direction d'écoulement. Afin d'illustrer cette première étape du processus de caractérisation 1D de la perméabilité, un échantillon a été construit en mettant bout à bout quatre échantillons différents, quasi-homogènes et de perméabilités effectives connues. L'échantillon a été initialement saturé par une saumure (30 g/1 de NaCl). Puis, on y injecte de la glycérine à 25 raison de 15 cm3/h. La viscosité de la saumure est de 1,01 cP et celle de la glycérine de 60 cP. L'évolution de la pression différentielle (DP) au cours du temps (T) est représentée sur la figure 5. L'équation 4 permet d'en déduire le profil de perméabilité (K) en fonction du temps (T), illustrée sur la figure 6. Ce profil est comparé aux perméabilités effectives mesurées pour (2 (3 ( 4 2909448 10 les 4 bouts d'échantillons. Les morceaux d'échantillons avaient été sélectionnés pour leur caractère homogène. Le profil calculé reproduit raisonnablement bien le profil théorique. La figure 7 montre le même profil en fonction de la longueur (d) de l'échantillon. Sur la figure 6, les points représentent le profil de perméabilité calculé en fonction du temps, et la courbe 5 noire représente le profil de perméabilité théorique. De même, sur la figure 7, les points représentent le profil de perméabilité calculé en fonction de la distance, et la courbe noire représente le profil de perméabilité théorique. Carte 3D de perméabilité 10 Le principe de l'interprétation repose sur une combinaison de résultats dynamiques (profil 1D de perméabilité) avec des résultats statiques (carte 3D de porosité). En effet, en supposant qu'il existe une loi porosité/perméabilité (K(1i))) valable à l'échelle de la discrétisation retenue pour établir la carte 3D de porosité, il est alors possible de calculer une carte 3D de perméabilité a priori. En moyennant alors la perméabilité au niveau des 15 différentes tranches (calcul en parallèle), on obtient alors une information directement comparable avec le profil 1D de perméabilité mesuré directement par interprétation de tests dynamiques. L'enjeu est alors d'optimiser la loi K(q) de façon à obtenir la meilleure concordance possible entre les deux profils. La figure 8 illustre les différentes étapes suivies pour obtenir une carte 3D de 20 perméabilité cohérente avec les données disponibles, c'est à dire la carte 3D de porosité et la pression différentielle mesurée au cours du temps. 1) Dans un premier temps, on construit une carte 3D de perméabilité (K((D)) à partir de la carte 3D de porosité (carte 3D de porosité t déduite de données CT). On peut se placer dans un cadre stochastique et procéder par co-simulation, auquel cas la carte de perméabilité 25 dépend d'un germe et est corrélée à la porosité. On peut aussi employer une relation déterministe entre la porosité et la perméabilité. Par exemple, la loi K(çb) peut prendre la forme suivante : K = a (q û ó) où : a est un préfacteur, 30 la porosité résiduelle, i.e. une valeur de porosité pour laquelle l'écoulement ne peut plus se produire, 2909448 11 f un exposant traduisant l'évolution de la perméabilité en fonction de la porosité utile. Dans le cadre de l'étude d'un milieu poreux présentant une structure plus complexe (bimodale par exemple), il est possible d'introduire deux lois K1(4 ), K2((i)) mais aussi une valeur seuil de porosité 4c permettant de distinguer les mailles où appliquer chacune des lois. 5 2) La carte 3D de perméabilité étant connue, on simule numériquement le test miscible visqueux (SNT). Les résultats de la simulation, c'est-à-dire la pression différentielle simulée (DP(t)sim), sont alors comparés à la pression différentielle mesurée au cours du temps (DP(t)mes). On en déduit la valeur d'une fonction objectif (OF). La fonction objectif OF mesure l'écart entre les données réelles et les réponses correspondantes simulées. On peut 10 aussi envisager d'exprimer cette fonction différemment. Par exemple, elle pourrait mesurer l'écart entre le profil 1 D de perméabilité déduit des données de pression du profil 1 D de perméabilité déduit des pressions simulées. 3) La dernière étape consiste à perturber (PERT) les paramètres du problème afin de déterminer une carte 3D de perméabilité qui minimise la fonction objectif, c'est-à-dire qui 15 minimise l'écart entre la pression différentielle mesurée et le signal numérique correspondant. Dans le cas d'une co-simulation, on modifiera la carte de perméabilité en appliquant par exemple la méthode de déformation graduelle. Dans le cas d'une relation déterministe, on modifiera les paramètres intervenant dans la loi, par exemple a, 0) et fi. 20 Ce processus de calage peut être enrichi par des mesures de concentrations (C). En effet, lorsque l'expérience se fait sous scanner, il est possible d'obtenir un profile de concentration ou encore des cartes 3D de concentration à des temps différents (Cmes) à partir de mesures RX (Maloney et Wegener, 2000). Ces données pourront être intégrées à la 25 fonction objectif pour être comparées aux réponses simulées équivalentes Csim (figure 9). Les applications potentielles de cette méthode sont importantes en particulier en ce qui concerne la caractérisation des milieux poreux hétérogènes comme les milieux vacuolaires par exemple, ou encore les milieux poreux soumis à des écoulements réactifs (suite à l'injection 30 de CO2 ou à une opération de stimulation acide de puits par exemple). REVENDICATIONS 1) Méthode de détermination d'une distribution tridimensionnelle de :[a perméabilité absolue d'un échantillon hétérogène, à partir d'une carte de porosité tridimensionnelle dudit 5 échantillon, caractérisée en ce qu' elle comporte les étapes suivantes : - on réalise au moins une expérience de déplacement miscible visqueux au cours de laquelle on détermine une évolution de la pression différentielle DP(t)mes de part et d'autre dudit échantillon ; - on détermine un profil 1D de perméabilité absolue le long dudit échantillon en 10 utilisant ladite évolution de la pression différentielle ; - on construit une première carte 3D de perméabilité (K(cD)) à partir de ladite carte 3D de porosité ; - on estime une pression différentielle simulée DP(t)sim, en simulant numériquement un test miscible visqueux à partir de ladite première carte de perméabilité et dudit profil 1D de 15 perméabilité ; - on détermine ladite distribution tridimensionnelle de la perméabilité absolue dudit échantillon en modifiant au moins une fois ladite première carte de perméabilité de façon à minimiser l'écart entre ladite pression différentielle simulée DP(t)sim et ladite pression différentielle mesurée au cours du temps DP(t)mes• 20 2) Méthode de détermination de la distribution 3D de la perméabilité absolue d'un échantillon hétérogène selon la revendication 1, dans laquelle on détermine ladite carte de porosité à partir de mesures statiques non destructives. 3) Méthode de détermination de la distribution 3D de la perméabilité absolue d'un échantillon hétérogène selon la revendication 2, dans laquelle on détermine ladite carte de 25 porosité à partir de mesures par rayons X en utilisant des fluides de densités différentes saturant ledit échantillon. 4) Méthode de détermination de la distribution 3D de la perméabilité absolue d'un échantillon hétérogène selon la revendication 2, dans laquelle on détermine ladite carte de porosité à partir de mesures par résonance magnétique nucléaire. 2909448 13 5) Méthode de détermination de la distribution 3D de la perméabilité absolue d'un échantillon hétérogène selon la revendication 1, dans laquelle on détermine ladite première carte de perméabilité (Ke)) à partir de co-simulations. 6) Méthode de détermination de la distribution 3D de la perméabilité absolue d'un échantillon hétérogène selon la revendication 1, dans laquelle on détermine ladite première carte de perméabilité (K(1)) à partir d'au moins une relation déterministe entre la porosité et la perméabilité. 7) Méthode de détermination de la distribution 3D de la perméabilité absolue d'un échantillon hétérogène selon la revendication 5, dans laquelle on modifie ladite première carte 10 de perméabilité en appliquant une méthode de déformation graduelle. 8) Méthode de détermination de la distribution 3D de la perméabilité absolue d'un échantillon hétérogène selon la revendication 1, dans laquelle on détermine en outre un profil de concentration ou des cartes 3D de concentration à des temps différents, que l'on compare à des réponses simulées équivalentes Csim.