CA2605741A1 - Methode et dispositif pour evaluer des parametres d'ecoulement et des parametres electriques d'un milieu poreux - Google Patents
Methode et dispositif pour evaluer des parametres d'ecoulement et des parametres electriques d'un milieu poreux Download PDFInfo
- Publication number
- CA2605741A1 CA2605741A1 CA002605741A CA2605741A CA2605741A1 CA 2605741 A1 CA2605741 A1 CA 2605741A1 CA 002605741 A CA002605741 A CA 002605741A CA 2605741 A CA2605741 A CA 2605741A CA 2605741 A1 CA2605741 A1 CA 2605741A1
- Authority
- CA
- Canada
- Prior art keywords
- sample
- fluid
- volume
- production
- pressure
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Abandoned
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
- G01N33/24—Earth materials
- G01N33/241—Earth materials for hydrocarbon content
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume, or surface-area of porous materials
- G01N15/08—Investigating permeability, pore-volume, or surface area of porous materials
- G01N15/082—Investigating permeability by forcing a fluid through a sample
- G01N15/0826—Investigating permeability by forcing a fluid through a sample and measuring fluid flow rate, i.e. permeation rate or pressure change
Abstract
Méthode et dispositif pour évaluer dans le même temps et avec le même appareillage, la résistivité électrique et des paramètres d'écoulement d'un milieu poreux, à partir de mesures expérimentales issues d'au moins une expérience de déplacement de type multi-débits. On mesure de façon continue en fonction du temps la pression différentielle entre deux extrémités dudit échantillon, la résistivité électrique dudit échantillon, et pour chaque palier de stabilisation on mesure également au moins deux valeurs de volume de production du fluide déplacé. On en déduit ensuite une courbe continue de volume de production du fluide déplacé en fonction du temps à partir de laquelle on estime des paramètres d'écoulements à l'aide d'un simulateur numérique d'écoulement ; Application à l'exploitation de gisements pétrolier s par exemple.
Description
MÉTHODE ET DISPOSITIF POUR ÉVALUER DES PARAMETRES
D'ÉCOULEMENT ET DES PARAMETRES ELECTRIQUES D'UN MILIEU
POREUX
La présente invention concerne une méthode et un dispositif pour évaluer des paramètres d'écoulement et des paramètres électriques d'un milieu poreux, avec le même appareillage, la même expérience et à partir d'un mêine échantillon de roclle prélevé dans ce milieu. Il peut s'agir par exemple de roches d'une zone d'un gisement souterrain de fluides (gisement hydrocarbures, aquifère, etc.).
Plus particulièrement, l'invention concerne la mesure de la perméabilité
relative et de paramètres électriques (résistivité) d'un milieu poreux en soumettant un échantillon de ce milieu successivement à une ou plusieurs phases de drainage ou d'imbibition.
Le procédé et le dispositif peuvent trouver des applications dans de nombreux domaines, notamment dans l'étude de gisements pétrolifères. L'invention concerne alors la mesure des courbes de perméabilité relative qui contrôlent le comportement dynamique des réservoirs pétroliers et les mesures électriques qui servent à
la calibration des outils de diagraphie pour évaluer la quantité d'hydrocarbure dans le réservoir. L'évaluation de la richesse et de la productivité d'un gisement souterrain est obtenue par exemple par des simulations numériques à partir de valeurs, mesurées ou déterminées, de certains paramètres des roches tels que la pression capillaire et leur perméabilité relative pour les fluides qui y sont présents. Ces paramètres sont
D'ÉCOULEMENT ET DES PARAMETRES ELECTRIQUES D'UN MILIEU
POREUX
La présente invention concerne une méthode et un dispositif pour évaluer des paramètres d'écoulement et des paramètres électriques d'un milieu poreux, avec le même appareillage, la même expérience et à partir d'un mêine échantillon de roclle prélevé dans ce milieu. Il peut s'agir par exemple de roches d'une zone d'un gisement souterrain de fluides (gisement hydrocarbures, aquifère, etc.).
Plus particulièrement, l'invention concerne la mesure de la perméabilité
relative et de paramètres électriques (résistivité) d'un milieu poreux en soumettant un échantillon de ce milieu successivement à une ou plusieurs phases de drainage ou d'imbibition.
Le procédé et le dispositif peuvent trouver des applications dans de nombreux domaines, notamment dans l'étude de gisements pétrolifères. L'invention concerne alors la mesure des courbes de perméabilité relative qui contrôlent le comportement dynamique des réservoirs pétroliers et les mesures électriques qui servent à
la calibration des outils de diagraphie pour évaluer la quantité d'hydrocarbure dans le réservoir. L'évaluation de la richesse et de la productivité d'un gisement souterrain est obtenue par exemple par des simulations numériques à partir de valeurs, mesurées ou déterminées, de certains paramètres des roches tels que la pression capillaire et leur perméabilité relative pour les fluides qui y sont présents. Ces paramètres sont
2 généralement mesurés en laboratoire à partir d'échantillons de roches prélevées in situ par carottage.
État de la technique La mise en production d'un champ pétrolifère nécessite donc une étude précise des conditions de récupération des hydrocarbures. Il s'agit de déterminer d'une part, la quantité d'hydrocarbure en place, et d'autre part, de déterminer les paramètres d'écoulement polyphasiques dans la roche qui constitue le réservoir, et notamment lors de phase d'injection d'eau ou d'huile.
Les roches qui constituent un réseivoir peuvent être de nature différente en terme de propriétés pétrophysiques et il est courant de faire des regroupements, appelés "rock types", pour travailler sur des échantillons appartenant à la même famille ("rock type").
Pour chacune des familles, il est ainsi nécessaire de réaliser l'ensemble des mesures pétrophysiques pour paramétrer le simulateur, calculer les quantités en place... car ces différents paramètres varient d'une famille à l'autre. Pour un "rock type"
donné, il est ainsi courant de choisir précautionneusement plusieurs échantillons, l'un étant utilisé
pour la détermination des courbes de pression capillaire (Pc), l'autre pour les courbes de perméabilités relatives (krs), un autre pour les courbes de résistivité, etc.
La détermination des paramètres d'écoulement inultiphasique tels que les pernéabilités relatives et la pression capillaire constitue aussi un enjeu majeur dans le cadre de l'exploitation pétrolière que ce soit dans la phase d'évaluation du gisement ou de sa mise en production. Les perméabilités relatives (krs) sont utilisées dans les simulateurs numériques pour décrire les écoulements polyphasiques dans la roche, en particulier quand on injecte de l'eau ou du gaz pour maintenir la pression et récupérer plus d'huile. Ce sont ces courbes qui, typiquement, permettent le contrôle du profll de production dans les gisements. Des expériences de laboratoire, menées en conditions de pressions et températures de gisement, sont couramment conduites pour déterminer ces paramètres de manière représentative. En particulier, les perméabilités relatives sont classiquement obtenues par l'intermédiaire d'expériences de déplacement telles que décrites dans le document FR 2 708 742 (US 5 679 885). Ces expériences de déplacements multi-débits consistent à injecter un fluide (appelé fluide injecté) dans un
État de la technique La mise en production d'un champ pétrolifère nécessite donc une étude précise des conditions de récupération des hydrocarbures. Il s'agit de déterminer d'une part, la quantité d'hydrocarbure en place, et d'autre part, de déterminer les paramètres d'écoulement polyphasiques dans la roche qui constitue le réservoir, et notamment lors de phase d'injection d'eau ou d'huile.
Les roches qui constituent un réseivoir peuvent être de nature différente en terme de propriétés pétrophysiques et il est courant de faire des regroupements, appelés "rock types", pour travailler sur des échantillons appartenant à la même famille ("rock type").
Pour chacune des familles, il est ainsi nécessaire de réaliser l'ensemble des mesures pétrophysiques pour paramétrer le simulateur, calculer les quantités en place... car ces différents paramètres varient d'une famille à l'autre. Pour un "rock type"
donné, il est ainsi courant de choisir précautionneusement plusieurs échantillons, l'un étant utilisé
pour la détermination des courbes de pression capillaire (Pc), l'autre pour les courbes de perméabilités relatives (krs), un autre pour les courbes de résistivité, etc.
La détermination des paramètres d'écoulement inultiphasique tels que les pernéabilités relatives et la pression capillaire constitue aussi un enjeu majeur dans le cadre de l'exploitation pétrolière que ce soit dans la phase d'évaluation du gisement ou de sa mise en production. Les perméabilités relatives (krs) sont utilisées dans les simulateurs numériques pour décrire les écoulements polyphasiques dans la roche, en particulier quand on injecte de l'eau ou du gaz pour maintenir la pression et récupérer plus d'huile. Ce sont ces courbes qui, typiquement, permettent le contrôle du profll de production dans les gisements. Des expériences de laboratoire, menées en conditions de pressions et températures de gisement, sont couramment conduites pour déterminer ces paramètres de manière représentative. En particulier, les perméabilités relatives sont classiquement obtenues par l'intermédiaire d'expériences de déplacement telles que décrites dans le document FR 2 708 742 (US 5 679 885). Ces expériences de déplacements multi-débits consistent à injecter un fluide (appelé fluide injecté) dans un
3 échantillon de roche initialement saturé avec majoritaireinent un autre fluide (appelé
fluide déplacé). Le fluide injecté draine alors une partie du fluide saturant hors de l'échantillon. La quantité de fluide déplacé ainsi produite par l'échantillon est mesurée sous forme d'un volume. Ce volume est appelé volume de production du fluide déplacé.
Au fur et à mesure de l'injection, la quantité de fluide déplacé produite augmente, c'est-à-dire que le volume de production du fluide déplacé augmente. Puis, on atteint une stabilisation de ce volume, c'est-à-dire, que le volume n'augmente plus. Ceci cotrespond à un arrêt de la production du fluide déplacé hors de l'échantillon. Le débit d'injection est alors augmenté pour drainer une partie du fluide saturant encore en place dans l'échantillon, jusqu'à obtenir un nouveau palier de stabilisation. Le débit d'injection est donc augmenté plusieurs fois (après chaque palier de stabilisation), et pour chacun de ces paliers, l'on mesure l'évolution des paramètres suivants, qui constituent les données expérimentales : la pression différentielle, pour un des fluides, de part et d'autre de l'échantillon, notée AP' (t) , le volume de production du fluide déplacé, notée V(ts), et la saturation locale pour l'un des deux fluides (injecté ou déplacé) et pour chaque palier d'injection, notés S(t, ), S(t2 ), ... Le fluide injecté peut, par exemple, être de l'eau et le fluide déplacé de l'huile. L'obtention des perméabilités relatives passe alors par une étape d'interprétation des données expérimentales de manière analytique, mais aussi numérique à l'aide d'un simulateur d'écoulement, de manière à prendre en compte tous les phénomènes physiques (pression capillaire, gravité) qui influencent les données expérimentales obtenues. Une telle interprétation est décrite par exeinple dans le document FR 2 708 742 (US 5 679 885).
On connaît également, par ce document, un dispositif permettant de réaliser ce type d'expériences de déplacement multi-débits.
En ce qui concerne la détermination des paramètres électriques, l'approche la plus courante consiste également à réaliser des expériences de déplacement à débit ou à
pression fixées. La résistivité d'un échantillons solide peut être mesurée au moyen d'électrodes qui sont placées en contact avec sa surface, à des emplacements choisis, entre lesquels on fait passer un courant électrique. La mesure de la différence de potentiel apparaissant entre les emplacements des électrodes, donne directement la
fluide déplacé). Le fluide injecté draine alors une partie du fluide saturant hors de l'échantillon. La quantité de fluide déplacé ainsi produite par l'échantillon est mesurée sous forme d'un volume. Ce volume est appelé volume de production du fluide déplacé.
Au fur et à mesure de l'injection, la quantité de fluide déplacé produite augmente, c'est-à-dire que le volume de production du fluide déplacé augmente. Puis, on atteint une stabilisation de ce volume, c'est-à-dire, que le volume n'augmente plus. Ceci cotrespond à un arrêt de la production du fluide déplacé hors de l'échantillon. Le débit d'injection est alors augmenté pour drainer une partie du fluide saturant encore en place dans l'échantillon, jusqu'à obtenir un nouveau palier de stabilisation. Le débit d'injection est donc augmenté plusieurs fois (après chaque palier de stabilisation), et pour chacun de ces paliers, l'on mesure l'évolution des paramètres suivants, qui constituent les données expérimentales : la pression différentielle, pour un des fluides, de part et d'autre de l'échantillon, notée AP' (t) , le volume de production du fluide déplacé, notée V(ts), et la saturation locale pour l'un des deux fluides (injecté ou déplacé) et pour chaque palier d'injection, notés S(t, ), S(t2 ), ... Le fluide injecté peut, par exemple, être de l'eau et le fluide déplacé de l'huile. L'obtention des perméabilités relatives passe alors par une étape d'interprétation des données expérimentales de manière analytique, mais aussi numérique à l'aide d'un simulateur d'écoulement, de manière à prendre en compte tous les phénomènes physiques (pression capillaire, gravité) qui influencent les données expérimentales obtenues. Une telle interprétation est décrite par exeinple dans le document FR 2 708 742 (US 5 679 885).
On connaît également, par ce document, un dispositif permettant de réaliser ce type d'expériences de déplacement multi-débits.
En ce qui concerne la détermination des paramètres électriques, l'approche la plus courante consiste également à réaliser des expériences de déplacement à débit ou à
pression fixées. La résistivité d'un échantillons solide peut être mesurée au moyen d'électrodes qui sont placées en contact avec sa surface, à des emplacements choisis, entre lesquels on fait passer un courant électrique. La mesure de la différence de potentiel apparaissant entre les emplacements des électrodes, donne directement la
4 mesure de la résistivité. Pour que les mesures soient représentatives, il iinporte que le contact des électrodes avec la surface soit le meilleur possible. Un procédé
connu consiste par exemple à placer l'échantillon à tester dans une gaine souple élastique. Les électrodes sont placées entre l'échantillon et la gaine et reliées au travers elle par des conducteurs électriques à un système de mesure de la conductivité électrique.
La gaine est placée dans une enceinte de conflnement. Un fluide sous pression est admis dans l'enceinte ce qui à pour effet de plaquer la gaine et par conséquent les électrodes contre l'échantillon. Un tel procédé utilisant un fluide sous pression pour plaquer une gaine contre un échantillon, est utilisé par exemple dans les outils de pétrophysique décrits dans les documents FR 2 708 742 (US 5 679 885) et FR 2 724 460 (US 5 610 524) ou encore dans le document US 5 105 154.
Un dispositif permettant de réaliser de telles expériences de déplacement à
débit fixé est décrit par exemple dans le document FR 2 708 742 (US 5 679 885), ou encore à
pression imposée par centrifugation dans le document FR 2 758 881 (GB 2 322 942).
On connaît également par le document FR 2.778 449 un dispositif permettant de réaliser de telles mesures de résistivité en conditions de réservoir tout en mesurant et en contrôlant les niveaux de pression capillaire dans la roche pendant les phases de drainage. Cela permet de tracer l'évolution de l'indice de résistivité en fonction de la saturation en saumure de l'échantillon. Une céramique placée en sortie d'échantillon permet de descendre à des niveaux de saturation extrêmement bas ce qui est particulièrement utile pour couvrir la plage de saturation d'un nombre important de réservoirs. Les données de résistivité sont alors utilisées directement par les diagraphistes pour estimer les saturations en huile initiales et donc les quantités d'1luile en place. Par ailleurs, ce dispositif est aussi utilisé pour obtenir le même type de résultats pendant les phases d'imbibition (spontanée et forcée). La saumure est alors injectée dans l'échantillon à débit imposé à travers la céramique. Plusieurs phases à
débits d'injection croissants sont en général réalisées de manière à diminuer la saturation en huile dans l'échantillon jusqu'à sa valeur résiduelle. La quantité
d'hydrocarbure est alors déterminée à partir de ces mesures de résistivité électrique. Le principe repose sur le fait que la résistivité d'un échantillon saturé avec de la saumure ou de l'huile dépend de la saturation en eau car l'liuile n'est pas un fluide conducteur.
L'évaluation des parainètres d'écoulements (krs) et des parainètres électriques requiert donc un appareillage différent, et donc des écllantillons différents pour chaque appareil. Or, l'évaluation de la richesse et de la productivité d'un gisement souterrain nécessite une mise en cohérence entre les mesures réalisées sur ces différents
connu consiste par exemple à placer l'échantillon à tester dans une gaine souple élastique. Les électrodes sont placées entre l'échantillon et la gaine et reliées au travers elle par des conducteurs électriques à un système de mesure de la conductivité électrique.
La gaine est placée dans une enceinte de conflnement. Un fluide sous pression est admis dans l'enceinte ce qui à pour effet de plaquer la gaine et par conséquent les électrodes contre l'échantillon. Un tel procédé utilisant un fluide sous pression pour plaquer une gaine contre un échantillon, est utilisé par exemple dans les outils de pétrophysique décrits dans les documents FR 2 708 742 (US 5 679 885) et FR 2 724 460 (US 5 610 524) ou encore dans le document US 5 105 154.
Un dispositif permettant de réaliser de telles expériences de déplacement à
débit fixé est décrit par exemple dans le document FR 2 708 742 (US 5 679 885), ou encore à
pression imposée par centrifugation dans le document FR 2 758 881 (GB 2 322 942).
On connaît également par le document FR 2.778 449 un dispositif permettant de réaliser de telles mesures de résistivité en conditions de réservoir tout en mesurant et en contrôlant les niveaux de pression capillaire dans la roche pendant les phases de drainage. Cela permet de tracer l'évolution de l'indice de résistivité en fonction de la saturation en saumure de l'échantillon. Une céramique placée en sortie d'échantillon permet de descendre à des niveaux de saturation extrêmement bas ce qui est particulièrement utile pour couvrir la plage de saturation d'un nombre important de réservoirs. Les données de résistivité sont alors utilisées directement par les diagraphistes pour estimer les saturations en huile initiales et donc les quantités d'1luile en place. Par ailleurs, ce dispositif est aussi utilisé pour obtenir le même type de résultats pendant les phases d'imbibition (spontanée et forcée). La saumure est alors injectée dans l'échantillon à débit imposé à travers la céramique. Plusieurs phases à
débits d'injection croissants sont en général réalisées de manière à diminuer la saturation en huile dans l'échantillon jusqu'à sa valeur résiduelle. La quantité
d'hydrocarbure est alors déterminée à partir de ces mesures de résistivité électrique. Le principe repose sur le fait que la résistivité d'un échantillon saturé avec de la saumure ou de l'huile dépend de la saturation en eau car l'liuile n'est pas un fluide conducteur.
L'évaluation des parainètres d'écoulements (krs) et des parainètres électriques requiert donc un appareillage différent, et donc des écllantillons différents pour chaque appareil. Or, l'évaluation de la richesse et de la productivité d'un gisement souterrain nécessite une mise en cohérence entre les mesures réalisées sur ces différents
5 écliantillons.
La méthode et le dispositif selon l'invention permettent l'évaluation des paramètres d'écoulement polyphasique avec une prise en compte complète et précise des effets capillaires ainsi que des paramètres électriques d'un milieu poreux, avec un seul type d'appareillage, une même expérience et à partir d'un même échantillon de roche prélevé dans ce milieu.
La méthode et le dispositif selon l'invention L'invention concerne une méthode pour évaluer conjointement la résistivité
électrique et des paramètres d'écoulement d'un inilieu poreux, à partir d'au moins une =
expérience de déplacement de type multi-débits. Ce type d'expérience consiste à injecter un premier fluide, appelé fluide injecté, dans un échantillon du milieu initialement saturé avec un second fluide, appelé fluide déplacé, dans laquelle le débit d'injection du fluide injecté est augmenté au moins une fois après un palier de stabilisation du volume de la production du fluide déplacé. La méthode selon l'invention comporte les étapes suivantes :
- on mesure de façon continue en fonction du temps la pression différentielle entre deux extrémités dudit échantillon ( OP' (t) ) et la résistivité
électrique dudit échantillon (R(t)) - on mesure le volume de production du fluide déplacé à chaque palier de stabilisation ( V (t, ) , V (t2 ) ) ;
- on estime une relation entre la résistivité électrique de l'échantillon et le volume de la production du fluide déplacé à partir desdites mesures de résistivités (R(t)) et desdites valeurs de volume de production du fluide déplacé ( V(ti ), V(t2 ));
La méthode et le dispositif selon l'invention permettent l'évaluation des paramètres d'écoulement polyphasique avec une prise en compte complète et précise des effets capillaires ainsi que des paramètres électriques d'un milieu poreux, avec un seul type d'appareillage, une même expérience et à partir d'un même échantillon de roche prélevé dans ce milieu.
La méthode et le dispositif selon l'invention L'invention concerne une méthode pour évaluer conjointement la résistivité
électrique et des paramètres d'écoulement d'un inilieu poreux, à partir d'au moins une =
expérience de déplacement de type multi-débits. Ce type d'expérience consiste à injecter un premier fluide, appelé fluide injecté, dans un échantillon du milieu initialement saturé avec un second fluide, appelé fluide déplacé, dans laquelle le débit d'injection du fluide injecté est augmenté au moins une fois après un palier de stabilisation du volume de la production du fluide déplacé. La méthode selon l'invention comporte les étapes suivantes :
- on mesure de façon continue en fonction du temps la pression différentielle entre deux extrémités dudit échantillon ( OP' (t) ) et la résistivité
électrique dudit échantillon (R(t)) - on mesure le volume de production du fluide déplacé à chaque palier de stabilisation ( V (t, ) , V (t2 ) ) ;
- on estime une relation entre la résistivité électrique de l'échantillon et le volume de la production du fluide déplacé à partir desdites mesures de résistivités (R(t)) et desdites valeurs de volume de production du fluide déplacé ( V(ti ), V(t2 ));
6 - on établit une courbe continue de volume de production du fluide déplacé en fonction du temps (V(t)) à partir de ladite relation ;
- on estime des paramètres d'écoulements à partir de ladite mesure continue de pression différentielle ( AP' (t) ) et de ladite courbe continue (V(t)) à
l'aide d'un simulateur numérique d'écoulement.
La relation entre la résistivité électrique et le volume de la production du fluide déplacé peut être estimée en outre à partir d'une courbe de la résistivité en fonction de la saturation de l'échantillon en fin de palier.
Selon l'invention, pour mesurer la résistivité électrique, il est possible de disposer radialement des électrodes de mesure de la résistivité électrique autour de l'échantillon et d'utiliser un échantillon de faible longueur, typiquement de l'ordre de 3 cm.
L'invention concerne également un dispositif pour réaliser sur un échantillon solide poreux (S) des expériences de déplacements multi-débits. Ce dispositif comporte une cellule de confinement (1) pour contenir ledit échantillon (S) saturé avec un premier fluide, appelé fluide déplacé, et fermée à une extrémité par un preinier embout cylindrique (16) relié à des moyens d'injection sous pression (5a, 5b) d'un second fluide, appelé fluide injecté, et fermée à l'autre extrémité par un second embout cylindrique (15) comportant un système de circulation des fluides (4a, 4b, 4c). Le dispositif selon l'invention est caractérisé en ce qu'il comporte :
- une plaque poreuse semi-perméable (C) placée entre ledit échantillon et ledit premier embout ;
- des moyens (6) de mesures continues de résistivité électrique dudit échantillon ;
- des moyens (3) de mesures du volume de production du fluide déplacé pour chaque palier de stabilisation ;
- des moyens (23, 22) de mesures continues de la pression différentielle d'un des fluides de part et d'autre de l'échantillon, entre l'échantillon et la plaque poreuse semi-perméable.
- on estime des paramètres d'écoulements à partir de ladite mesure continue de pression différentielle ( AP' (t) ) et de ladite courbe continue (V(t)) à
l'aide d'un simulateur numérique d'écoulement.
La relation entre la résistivité électrique et le volume de la production du fluide déplacé peut être estimée en outre à partir d'une courbe de la résistivité en fonction de la saturation de l'échantillon en fin de palier.
Selon l'invention, pour mesurer la résistivité électrique, il est possible de disposer radialement des électrodes de mesure de la résistivité électrique autour de l'échantillon et d'utiliser un échantillon de faible longueur, typiquement de l'ordre de 3 cm.
L'invention concerne également un dispositif pour réaliser sur un échantillon solide poreux (S) des expériences de déplacements multi-débits. Ce dispositif comporte une cellule de confinement (1) pour contenir ledit échantillon (S) saturé avec un premier fluide, appelé fluide déplacé, et fermée à une extrémité par un preinier embout cylindrique (16) relié à des moyens d'injection sous pression (5a, 5b) d'un second fluide, appelé fluide injecté, et fermée à l'autre extrémité par un second embout cylindrique (15) comportant un système de circulation des fluides (4a, 4b, 4c). Le dispositif selon l'invention est caractérisé en ce qu'il comporte :
- une plaque poreuse semi-perméable (C) placée entre ledit échantillon et ledit premier embout ;
- des moyens (6) de mesures continues de résistivité électrique dudit échantillon ;
- des moyens (3) de mesures du volume de production du fluide déplacé pour chaque palier de stabilisation ;
- des moyens (23, 22) de mesures continues de la pression différentielle d'un des fluides de part et d'autre de l'échantillon, entre l'échantillon et la plaque poreuse semi-perméable.
7 Selon l'invention, le premier embout cylindrique peut comporter un premier cylindre placé (18) dans un second cylindre concentrique (17).
Les moyens de mesures continues de la pression différentielle peuvent comporter au moins un conduit de prise de pression (22) relié à un capteur différentiel de pression (23), ce conduit pouvant comporter un ensemble de rainures longitudinales (27) usiné
sur la surface interne dudit second cylindre (17) ou sur la surface externe dudit preznier cylindre (18).
Enfin, selon l'invention, le conduit de prise de pression peut être adapté à
éviter le piégeage d'air lors de l'installation du premier embout.
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 représente schématiquement un dispositif permettant de réaliser des expériences de déplacement multi-débits.
- La figure 2 illustre en détail la cellule de confinement (1) du dispositif selon l'invention.
- La figure 3 illustre en détail l'embout d'injection (16) du dispositif selon l'invention.
- La figure 4 illustre un mode de réalisation du dispositif selon l'invention dans lequel le second cylindre (18) possède des rainures et des canaux (24-27) permettant de réaliser une prise de pression entre la céramique et l'échantillon.
- La figure 5 présente les points de saturation (M) obtenus en fin de période d'injection en fonction des mesures de résistance correspondantes (Re), ainsi que la courbe de résistance (CE) en fonction de la saturation en saumure (Sw).
- La figure 6 illustre, pour chacun des trois débits, les courbes de volume de production (OP) et de pression différentielle (PD) en fonction du temps.
Les moyens de mesures continues de la pression différentielle peuvent comporter au moins un conduit de prise de pression (22) relié à un capteur différentiel de pression (23), ce conduit pouvant comporter un ensemble de rainures longitudinales (27) usiné
sur la surface interne dudit second cylindre (17) ou sur la surface externe dudit preznier cylindre (18).
Enfin, selon l'invention, le conduit de prise de pression peut être adapté à
éviter le piégeage d'air lors de l'installation du premier embout.
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 représente schématiquement un dispositif permettant de réaliser des expériences de déplacement multi-débits.
- La figure 2 illustre en détail la cellule de confinement (1) du dispositif selon l'invention.
- La figure 3 illustre en détail l'embout d'injection (16) du dispositif selon l'invention.
- La figure 4 illustre un mode de réalisation du dispositif selon l'invention dans lequel le second cylindre (18) possède des rainures et des canaux (24-27) permettant de réaliser une prise de pression entre la céramique et l'échantillon.
- La figure 5 présente les points de saturation (M) obtenus en fin de période d'injection en fonction des mesures de résistance correspondantes (Re), ainsi que la courbe de résistance (CE) en fonction de la saturation en saumure (Sw).
- La figure 6 illustre, pour chacun des trois débits, les courbes de volume de production (OP) et de pression différentielle (PD) en fonction du temps.
8 - La figure 7 présente trois points de pression capillaire obtenus de façon analytique (PCA) en fonction de la saturation (SW) ainsi qu'une loi analytique (PCHM) de pression capillaire de type puissance.
- La figure 8A montre trois couples expérimentaux (CE) de valeurs de débit d'injection (q) et de pression différeiitielle (DP) ainsi que le résultat d'une expression analytique de type exponentielle (LA).
- La figure 8B illustre les résultats obtenus en terme de 1w;,, analytiques ( Iwa,nj) en fonction de la saturation (Sti,,.
- La figure 9 illustre le calage des données expérimentales en présentant, pour chacun des trois débits, les courbes de volume de production simulées (OPS) en fonction du temps et les courbes de pression différentielle siinulées (PDS).
- La figure 10 illustre les courbes de perméabilité relatives (bs) obtenues après simulation.
- La figure 11 illustre les profils de saturation obtenus par simulation.
Description détaillée du dispositif et de la méthode 1. Le dispositif Le dispositif selon l'invention permet de réaliser, sur un même échantillon issu d'un milieux poreux, des expériences de déplacements multi-débits au cours desquelles il est possible d'effectuer des mesures pétrophysiques (résistivité, saturation, pression différentielle,...) permettant d'évaluer conjointement les paramètres d'écoulements polyphasiques (perméabilités relatives) et la courbe de résistivité.
Ce dispositif est représenté schématiquement sur la figure 1. Ce dispositif comporte une cellule de confinement allongée 1, qui peut être cylindrique, et qui contient l'échantillon à tester S. Cette cellule est placée à l'intérieur d'une enceinte thermostatée 2 dont la température peut monter généralement jusqu'à 150 C. La face supérieure de la cellule 1 est équipée d'un système permettant de réaliser une circulation de fluide conduisant à un balayage de la face supérieure de l'échantillon. En effet, par soucis de précision, la face de sortie de l'échantillon est balayée à fort débit par de la
- La figure 8A montre trois couples expérimentaux (CE) de valeurs de débit d'injection (q) et de pression différeiitielle (DP) ainsi que le résultat d'une expression analytique de type exponentielle (LA).
- La figure 8B illustre les résultats obtenus en terme de 1w;,, analytiques ( Iwa,nj) en fonction de la saturation (Sti,,.
- La figure 9 illustre le calage des données expérimentales en présentant, pour chacun des trois débits, les courbes de volume de production simulées (OPS) en fonction du temps et les courbes de pression différentielle siinulées (PDS).
- La figure 10 illustre les courbes de perméabilité relatives (bs) obtenues après simulation.
- La figure 11 illustre les profils de saturation obtenus par simulation.
Description détaillée du dispositif et de la méthode 1. Le dispositif Le dispositif selon l'invention permet de réaliser, sur un même échantillon issu d'un milieux poreux, des expériences de déplacements multi-débits au cours desquelles il est possible d'effectuer des mesures pétrophysiques (résistivité, saturation, pression différentielle,...) permettant d'évaluer conjointement les paramètres d'écoulements polyphasiques (perméabilités relatives) et la courbe de résistivité.
Ce dispositif est représenté schématiquement sur la figure 1. Ce dispositif comporte une cellule de confinement allongée 1, qui peut être cylindrique, et qui contient l'échantillon à tester S. Cette cellule est placée à l'intérieur d'une enceinte thermostatée 2 dont la température peut monter généralement jusqu'à 150 C. La face supérieure de la cellule 1 est équipée d'un système permettant de réaliser une circulation de fluide conduisant à un balayage de la face supérieure de l'échantillon. En effet, par soucis de précision, la face de sortie de l'échantillon est balayée à fort débit par de la
9 PCT/FR2006/000911 saumure en circuit fermé de manière à faire circuler rapidement l'intégrité de l'huile produite de l'échantillon vers un système 3 permettant de déterminer le volume de fluide expulsé et comportant un système de séparation et de détection du niveau d'huile tel que des burettes graduées. Ce système de circulation se compose d'une poznpe à
haut débit HFR ("Higl1 Flow Rate") 4a qui injecte de la sauinure vers la surface supérieure de l'échantillon à travers un canal 4b, ainsi qûe d'un canal 4c qui évacue la saumure et l'huile de la face'de l'échantillon vers le système de séparation et de détection du niveau d'huile 3. Au niveau de la face opposée, le dispositif comporte une pompe à
faible débit LFR (Low Flow rate) 5a qui injecte de la saumure vers la surface inférieure de l'échantillon à travers un canal 5b. Le dispositif comporte également un appareil de mesure de la résistivité d'un type connu (impédancemètre "RLC meter") 6, un "isolator"
8 ainsi que des électrodes 9 en contact direct avec l'échantillon S.
La figure 2 illustre plus en détail le dispositif au niveau de la cellule de confinement. L'échantillon à tester S, est placé à l'intérieur d'une gaine déformable 10, et l'enseinble est disposé dans l'enceinte 1. L'espace annulaire 11 autour de la gaine 10 est mis en communication avec une source de fluide sous pression 12 délivrant du fluide par l'intermédiaire d'une ouverture ménagée 13 dans la paroi extérieure de l'enceinte 1, et d'une canalisation 14 raccordée à cette ouverture. Il peut s'agir d'une bouteille de gaz comprimé, d'un vérin ou encore d'une pompe hydraulique. L'application de ce fluide, permet de mettre sous une pression déterminée l'échantillon à tester.
Une plaque poreuse semi perméable C(imperméable au fluide déplacé contenu dans l'échantillon mais perméable au fluide injecté), par exemple une céramique, est placée à la base de l'échantillon S, c'est-à-dire au niveau de la face d'injection, pour obtenir des niveaux de saturation extrêmement bas comme décrit dans la demande de brevet français 03/15.199 du demandeur.
Le dispositif selon l'invention comporte deux embouts cylindriques 15 et 16 permettant de fermer la cellule de confineinent : un embout supérieur 15 (embout de balayage) situé au sommet de l'échantillon, du côté circulation de fluide, ainsi qu'un embout inférieur 16 (embout d'injection) situé à la base de l'échantillon, du côté de l'injection de fluide.
L'einbout de balayage 15 est constitué d'un cylindre en acier très inoxydable tel que de l'hastelloy. Au travers de ce cylindre, les deux canaux 4b et 4c sont orientés parallèlement à l'axe de l'enceinte.
L'embout d'injection comporte deux cylindres concentriques 17, 18. Le premier 5 cylindre 17, creux, est placé de façon permanente en contact avec la gaine.
Il permet de retirer le second cylindre 18, plein, sans modifier les conditions de pression exercées par la gaine sur l'échantillon. Ce premier cylindre a donc un diamètre sensibleinent équivalent à celui de l'échantillon. La taille de sa paroi doit être suffisamment importante pour ne pas se déformer sous l'effet de la pression exercée par la gaine, mais
haut débit HFR ("Higl1 Flow Rate") 4a qui injecte de la sauinure vers la surface supérieure de l'échantillon à travers un canal 4b, ainsi qûe d'un canal 4c qui évacue la saumure et l'huile de la face'de l'échantillon vers le système de séparation et de détection du niveau d'huile 3. Au niveau de la face opposée, le dispositif comporte une pompe à
faible débit LFR (Low Flow rate) 5a qui injecte de la saumure vers la surface inférieure de l'échantillon à travers un canal 5b. Le dispositif comporte également un appareil de mesure de la résistivité d'un type connu (impédancemètre "RLC meter") 6, un "isolator"
8 ainsi que des électrodes 9 en contact direct avec l'échantillon S.
La figure 2 illustre plus en détail le dispositif au niveau de la cellule de confinement. L'échantillon à tester S, est placé à l'intérieur d'une gaine déformable 10, et l'enseinble est disposé dans l'enceinte 1. L'espace annulaire 11 autour de la gaine 10 est mis en communication avec une source de fluide sous pression 12 délivrant du fluide par l'intermédiaire d'une ouverture ménagée 13 dans la paroi extérieure de l'enceinte 1, et d'une canalisation 14 raccordée à cette ouverture. Il peut s'agir d'une bouteille de gaz comprimé, d'un vérin ou encore d'une pompe hydraulique. L'application de ce fluide, permet de mettre sous une pression déterminée l'échantillon à tester.
Une plaque poreuse semi perméable C(imperméable au fluide déplacé contenu dans l'échantillon mais perméable au fluide injecté), par exemple une céramique, est placée à la base de l'échantillon S, c'est-à-dire au niveau de la face d'injection, pour obtenir des niveaux de saturation extrêmement bas comme décrit dans la demande de brevet français 03/15.199 du demandeur.
Le dispositif selon l'invention comporte deux embouts cylindriques 15 et 16 permettant de fermer la cellule de confineinent : un embout supérieur 15 (embout de balayage) situé au sommet de l'échantillon, du côté circulation de fluide, ainsi qu'un embout inférieur 16 (embout d'injection) situé à la base de l'échantillon, du côté de l'injection de fluide.
L'einbout de balayage 15 est constitué d'un cylindre en acier très inoxydable tel que de l'hastelloy. Au travers de ce cylindre, les deux canaux 4b et 4c sont orientés parallèlement à l'axe de l'enceinte.
L'embout d'injection comporte deux cylindres concentriques 17, 18. Le premier 5 cylindre 17, creux, est placé de façon permanente en contact avec la gaine.
Il permet de retirer le second cylindre 18, plein, sans modifier les conditions de pression exercées par la gaine sur l'échantillon. Ce premier cylindre a donc un diamètre sensibleinent équivalent à celui de l'échantillon. La taille de sa paroi doit être suffisamment importante pour ne pas se déformer sous l'effet de la pression exercée par la gaine, mais
10 suffisamment réduite pour que le second cylindre, dont la surface de sa base est supérieure à celle de la céramique, puisse y être encastré. Ce second cylindre 18 est décrit en détail par la figure 3. Il est également en acier très inoxydable tel que de l'hastelloy. Il possède un canal 19, parallèle à son axe, le traversant de part et d'autre, et dans lequel un fluide peut être injecté. Du côté en contact avec l'échantillon, le cylindre possède une cavité 20 permettant de contenir la plaque de céramique C, tout en assurant le contact entre la plaque et l'échantillon.
L'étanchéité entre les deux cylindres de l'embout d'injection 16 est assurée par des joints "O-Ring" 21 disposés entre les deux cylindres.
La présence de la céramique semi perméable permet d'obtenir des niveaux de saturation extrêmement bas. Mais en revanche, elle constitue un inconvénient majeur pour la détermination des perméabilités relatives. En effet, du fait de sa faible perméabilité, la pression différentielle engendrée sur l'ensemble "céramique +
échantillon" est si forte que la mesure de la pression différentielle globale ne permet pas de détecter la réponse de l'échantillon. La pression différentielle produite sur l'échantillon lui-même étant très faible en comparaison. Cela rend la correction d'une mesure globale (céramique et échantillon) de pression différentielle difficile et peu précise en pratique.
C'est pourquoi le dispositif selon l'invention possède, au niveau de l'embout d'injection, une prise de pression 22 (figure 2) entre l'échantillon et la céramique, permettant, à l'aide d'un capteur différentiel de pression 23, de mesurer la pression
L'étanchéité entre les deux cylindres de l'embout d'injection 16 est assurée par des joints "O-Ring" 21 disposés entre les deux cylindres.
La présence de la céramique semi perméable permet d'obtenir des niveaux de saturation extrêmement bas. Mais en revanche, elle constitue un inconvénient majeur pour la détermination des perméabilités relatives. En effet, du fait de sa faible perméabilité, la pression différentielle engendrée sur l'ensemble "céramique +
échantillon" est si forte que la mesure de la pression différentielle globale ne permet pas de détecter la réponse de l'échantillon. La pression différentielle produite sur l'échantillon lui-même étant très faible en comparaison. Cela rend la correction d'une mesure globale (céramique et échantillon) de pression différentielle difficile et peu précise en pratique.
C'est pourquoi le dispositif selon l'invention possède, au niveau de l'embout d'injection, une prise de pression 22 (figure 2) entre l'échantillon et la céramique, permettant, à l'aide d'un capteur différentiel de pression 23, de mesurer la pression
11 différentielle produite sur l'échantillon lui-mêine et non sur l'enseinble "céra.inique +
échantillon".
Selon un premier mode de réalisation, le second cylindre comporte cette prise de pression 22 entre l'échantillon et la céramique pennettant de mesurer la pression différentielle produite sur l'échantillon lui-même.
Selon un second mode de réalisation (figure 4), ce second cylindre possède également une rainure circulaire 24 perpendiculaire à l'axe du cylindre. Cette rainure est placée à une hauteur comprise entre la céramique C et les joints 21. Ce second cylindre possède au moins un second canal 25 parallèle à son axe, le traversant depuis sa base jusqu'à une hauteur coinprise entre la rainure circulaire 24 et les joints 21.
Ce second canal 25 est en communication avec la rainure circulaire 24. Enfin, à
l'intérieur de ce cylindre au niveau où ce deuxième canal s'arrête, au moins un troisième canal 26, perpendiculaire à l'axe du cylindre, permet de relier le second canal 25 à la périphérie du cylindre, entre l'échantillon et le joint. Au moins un conduit comportant des rainures verticales 27 (parallèles à l'axe du cylindre) partant de l'échantillon (en périphérie de la plaque poreuse) et s'arrêtant à hauteur du troisième canal 26, est disposé
entre les deux cylindres. Au moins l'une de ces rainures est en coiumunication avec le troisième canal.
Selon le mode de réalisation de la figure 4, ces rainures sont au nombre de 4 et sont disposées de façon régulière autour de la périphérie des cylindres. Selon les modes de réalisation, ces rainures verticales sont soit sur la surface (nécessairement externe car le cylindre est plein) du cylindre central 18, soit sur la surface interne du cylindre périphérique 17.
Dans ce second mode de réalisation (figure 4), la prise de pression 22 est constituée du conduit comportant des rainures verticales 27. En effet, ce conduit permet de réaliser une prise de pression entre la céramique et l'échantillon, pression qui peut être mesurée à la sortie du canal 25.
Par ailleurs, grâce à cet embout d'injection 16 possédant des rainures 27, il est possible de monter le cylindre mobile 18 sans piéger de l'air.
Ce dispositif permet également de mesurer la pression différentielle sur l'écha.ntillon au cours de l'expérience sans tenir compte de la céramique, en utilisant le
échantillon".
Selon un premier mode de réalisation, le second cylindre comporte cette prise de pression 22 entre l'échantillon et la céramique pennettant de mesurer la pression différentielle produite sur l'échantillon lui-même.
Selon un second mode de réalisation (figure 4), ce second cylindre possède également une rainure circulaire 24 perpendiculaire à l'axe du cylindre. Cette rainure est placée à une hauteur comprise entre la céramique C et les joints 21. Ce second cylindre possède au moins un second canal 25 parallèle à son axe, le traversant depuis sa base jusqu'à une hauteur coinprise entre la rainure circulaire 24 et les joints 21.
Ce second canal 25 est en communication avec la rainure circulaire 24. Enfin, à
l'intérieur de ce cylindre au niveau où ce deuxième canal s'arrête, au moins un troisième canal 26, perpendiculaire à l'axe du cylindre, permet de relier le second canal 25 à la périphérie du cylindre, entre l'échantillon et le joint. Au moins un conduit comportant des rainures verticales 27 (parallèles à l'axe du cylindre) partant de l'échantillon (en périphérie de la plaque poreuse) et s'arrêtant à hauteur du troisième canal 26, est disposé
entre les deux cylindres. Au moins l'une de ces rainures est en coiumunication avec le troisième canal.
Selon le mode de réalisation de la figure 4, ces rainures sont au nombre de 4 et sont disposées de façon régulière autour de la périphérie des cylindres. Selon les modes de réalisation, ces rainures verticales sont soit sur la surface (nécessairement externe car le cylindre est plein) du cylindre central 18, soit sur la surface interne du cylindre périphérique 17.
Dans ce second mode de réalisation (figure 4), la prise de pression 22 est constituée du conduit comportant des rainures verticales 27. En effet, ce conduit permet de réaliser une prise de pression entre la céramique et l'échantillon, pression qui peut être mesurée à la sortie du canal 25.
Par ailleurs, grâce à cet embout d'injection 16 possédant des rainures 27, il est possible de monter le cylindre mobile 18 sans piéger de l'air.
Ce dispositif permet également de mesurer la pression différentielle sur l'écha.ntillon au cours de l'expérience sans tenir compte de la céramique, en utilisant le
12 moyen de mesure de pression.23 disposé à la sortie du second cana125 parallèle à l'axe du cylindre.
Ainsi, le dispositif perinet de réaliser des expériences de déplacement inulti-débits en conditions de réservoir tout en effectuant :
- des mesures continues de résistivité électrique de l'échantillon, noté R(t) à
l'aide de l'appareil 6 ;
- des mesures ponctuelles du volume de la production du fluide déplacé pour cllaque palier d'injection, notée V(tl ), V(t2 ),...), à l'aide du système 3;
- des mesures continues de la pression différentielle, pour un des fluides, de part et d'autre de l'échantillon, notée AP'(t), grâce au capteur de pression différentielle 23 et à la prise de pression 22 placée entre la céramique et l'échantillon.
Ce dispositif fonctionne de la façon suivante :
On considère le cas d'un échantillon initialement saturé en huile par exemple, qui est placé dans l'enceinte 1 et mis sous pression par injection de fluide dans l'espace annulaire 11 autour de la gaine 10. Au moyen de la pompe 5a, par le canal 5b, on injecte de l'eau sous pression avec un débit fixé à une première extréznité de l'échantillon. A
l'extrémité opposée, on établit une circulation d'eau (4a, 4b, 4c) à une pression basse constante (la pression atmosphérique par exemple ou bien une pression fixe imposée par un réservoir de fluide à une pression de référence) et avec un débit faible.
Cette eau est prélevée par le système 3. L'eau est injectée par la pompe 4a dans le canal 4b, et elle ressort par le canal 4c entraînant les fluides (eau et huile) au fur et à
mesure qu'ils sont expulsés hors de l'échantillon S. De cette manière on assure à tout instant la continuité
entre le fluide circulant et l'huile expulsée de l'échantillon.
2. La méthode d'inter~rétation des mesures La méthode selon l'invention permet d'exploiter ce dispositif pour évaluer des courbes de perméabilités relatives malgré la présence de la céramique et la longueur
Ainsi, le dispositif perinet de réaliser des expériences de déplacement inulti-débits en conditions de réservoir tout en effectuant :
- des mesures continues de résistivité électrique de l'échantillon, noté R(t) à
l'aide de l'appareil 6 ;
- des mesures ponctuelles du volume de la production du fluide déplacé pour cllaque palier d'injection, notée V(tl ), V(t2 ),...), à l'aide du système 3;
- des mesures continues de la pression différentielle, pour un des fluides, de part et d'autre de l'échantillon, notée AP'(t), grâce au capteur de pression différentielle 23 et à la prise de pression 22 placée entre la céramique et l'échantillon.
Ce dispositif fonctionne de la façon suivante :
On considère le cas d'un échantillon initialement saturé en huile par exemple, qui est placé dans l'enceinte 1 et mis sous pression par injection de fluide dans l'espace annulaire 11 autour de la gaine 10. Au moyen de la pompe 5a, par le canal 5b, on injecte de l'eau sous pression avec un débit fixé à une première extréznité de l'échantillon. A
l'extrémité opposée, on établit une circulation d'eau (4a, 4b, 4c) à une pression basse constante (la pression atmosphérique par exemple ou bien une pression fixe imposée par un réservoir de fluide à une pression de référence) et avec un débit faible.
Cette eau est prélevée par le système 3. L'eau est injectée par la pompe 4a dans le canal 4b, et elle ressort par le canal 4c entraînant les fluides (eau et huile) au fur et à
mesure qu'ils sont expulsés hors de l'échantillon S. De cette manière on assure à tout instant la continuité
entre le fluide circulant et l'huile expulsée de l'échantillon.
2. La méthode d'inter~rétation des mesures La méthode selon l'invention permet d'exploiter ce dispositif pour évaluer des courbes de perméabilités relatives malgré la présence de la céramique et la longueur
13 courte de l'échantillon utilisé (typiquement de l'ordre de 3 cm) et avec une prise en compte complète et précise des effets capillaires. La méthode permet dans le même temps de déterminer la quantité du fluide déplacé en place dans le inilieu.
Dans le cadre de la description, on va se placer dans une application particulière qui consiste à injecter de la saumure dans un échantillon saturé
principalement en huile avec une saturation réduite en saumure. Cette méthode est appliquée dans le cadre de l'exploration pétrolière pour détenniner les courbes de lcrs ou pour faire varier la saturation en eau dans un échantillon. On ne sortirait pas du cadre de l'invention, si l'on inversait les fluides ou si l'on utilisait deux fluides quelconque non iniscibles.
Étape 1= expérience de déplacement de type multi-débits et mesures pétrophysiques On réalise tout d'abord une expérience de déplacement de type multi-débits, bien connue de l'homme du métier et décrite par exemple dans Egermaiul P. and Lenormand R., "A new methodology to evaluate the inipact of the local heterogeneity on petr oplzysical parameters ~(lw, Pc) :
application on carbonate rocks", SCA n 2004-18, Sciety of Core Analysts, Abu Dhabi, 2004.
Le procédé consiste essentiellement à placer un échantillon solide (S) saturé
avec un premier fluide (de l'huile par exemple) dans une cellule de confinement close (1), d'injecter à un débit fixé par des moyens de pression (5a, 5b) un autre fluide sous pression (une saumure par exemple) à une première extrémité de l'enceinte et de balayer l'extrémité opposée par un fluide à basse pression mis en circulation par des moyens de pompage (4a, 4b, 4c), qui entraîne au dehors le fluide drainé.
Ensuite, une nouvelle injection est réalisée avec un débit différent. Le nombre d'injection est fonction de la qualité et du teinps souhaité pour obtenir les mesures.
Au cours de cette expérience, on réalise :
- des mesures continues de résistivité électrique de l'échantillon, noté R(t) ;
- des mesures ponctuelles du volume de la production du fluide déplacé pour chaque palier d'injection, notée V(t, ), V(tz ),...) ;
Dans le cadre de la description, on va se placer dans une application particulière qui consiste à injecter de la saumure dans un échantillon saturé
principalement en huile avec une saturation réduite en saumure. Cette méthode est appliquée dans le cadre de l'exploration pétrolière pour détenniner les courbes de lcrs ou pour faire varier la saturation en eau dans un échantillon. On ne sortirait pas du cadre de l'invention, si l'on inversait les fluides ou si l'on utilisait deux fluides quelconque non iniscibles.
Étape 1= expérience de déplacement de type multi-débits et mesures pétrophysiques On réalise tout d'abord une expérience de déplacement de type multi-débits, bien connue de l'homme du métier et décrite par exemple dans Egermaiul P. and Lenormand R., "A new methodology to evaluate the inipact of the local heterogeneity on petr oplzysical parameters ~(lw, Pc) :
application on carbonate rocks", SCA n 2004-18, Sciety of Core Analysts, Abu Dhabi, 2004.
Le procédé consiste essentiellement à placer un échantillon solide (S) saturé
avec un premier fluide (de l'huile par exemple) dans une cellule de confinement close (1), d'injecter à un débit fixé par des moyens de pression (5a, 5b) un autre fluide sous pression (une saumure par exemple) à une première extrémité de l'enceinte et de balayer l'extrémité opposée par un fluide à basse pression mis en circulation par des moyens de pompage (4a, 4b, 4c), qui entraîne au dehors le fluide drainé.
Ensuite, une nouvelle injection est réalisée avec un débit différent. Le nombre d'injection est fonction de la qualité et du teinps souhaité pour obtenir les mesures.
Au cours de cette expérience, on réalise :
- des mesures continues de résistivité électrique de l'échantillon, noté R(t) ;
- des mesures ponctuelles du volume de la production du fluide déplacé pour chaque palier d'injection, notée V(t, ), V(tz ),...) ;
14 - des mesures continues de la pression différentielle, pour un des fluides, de part et d'autre de l'échantillon, notée AP' (t) .
Étape 2: détennination d'une courbe continue du volume de production du fluide déplacé en fonction du temM
On établit une courbe continue du volume de production du fluide déplacé en fonction du temps V(t). Le dispositif selon l'invention permet d'effectuer des mesures ponctuelles du volume de la production de fluide déplacé V(tl ), V(tz ), ...
Il serait alors possible d'estimer directement à partir de ces mesures ponctuelles une courbe continue V(t). Cependant, on ne dispose d'aucune information a priori sur la forme de la courbe V(t). De plus n'ayant à disposition que peu de points de mesures du volume de la production du fluide déplacé V(tl ), V(t2 ), ... l'estimation d'une courbe complexe n'est pas envisageable. C'est pourquoi la méthode selon l'invention tire profit des mesures continues de la résistivité électrique de l'échantillon qui donne "l'allure"
de la courbe V(t). A partir de couples connus pour ces deux types de valeurs, il est possible de déterminer une loi reliant la résistivité au volume de la production de fluide déplacé.
Cette détermination peut se faire, par exemple, à l'aide d'une régression car cette relation est linéaire ou quasi linéaire, ce qui permet de coininettre moins d'erreur sur cette estimation que sur l'estimation directe de la courbe V(t). Ensuite, on applique cette loi pour transformer les mesures continues de résistivité en courbe continue de volume de production du fluide déplacé en fonction du temps V(t).
Il est classique dans le domaine pétrophysique d'exprimer les différents paramètres en fonction de la saturation d'un des fluides. Ainsi, que ce soit pour représenter les pressions capillaires, les perméabilités relatives, les mesures de résistivité, on représente ces grandeurs en fonction de la saturation car c'est une grandeur qui ne dépend pas du volume poreux. Cela permet ainsi de comparer facilement les résultats obtenus sur différents échantillons dont les propriétés macroscopiques ne sont jamais exactement les mêmes.
C'est pourquoi, selon un mode de réalisation, on établit au préalable de la courbe continue de volume dë production du fluide déplacé en fonction du temps V(t) une courbe de la résistivité en fonction de la saturation.
Pour ce faire, à partir des valeurs de volume de production de fluide déplacé
5 V(tl ), V(t2 ), ... mesurées à la fin de chaque période d'injection, et à
partir du volume de fluide inj ecté ( V,, (tl ), Vw (t2 ), ... ), on déduit des valeurs de saturation moyenne en fluide injecté ( Sw) et en fluide déplacé ( S, ) :
SN,=1-So=Sti,,;+VP
avec : S,,,, la saturation initiale en fluide injecté
10 VP le volume poreux de l'échantillon (mesuré au préalable, par exemple par une technique de RMN) A partir de couples connus de saturation et de résistivité il est possible de déterminer, une loi reliant la résistivité à la saturation. Cette détermination peut se faire, par exemple, à l'aide d'une régression linéaire ou d'une régression non linéaire. Ensuite,
Étape 2: détennination d'une courbe continue du volume de production du fluide déplacé en fonction du temM
On établit une courbe continue du volume de production du fluide déplacé en fonction du temps V(t). Le dispositif selon l'invention permet d'effectuer des mesures ponctuelles du volume de la production de fluide déplacé V(tl ), V(tz ), ...
Il serait alors possible d'estimer directement à partir de ces mesures ponctuelles une courbe continue V(t). Cependant, on ne dispose d'aucune information a priori sur la forme de la courbe V(t). De plus n'ayant à disposition que peu de points de mesures du volume de la production du fluide déplacé V(tl ), V(t2 ), ... l'estimation d'une courbe complexe n'est pas envisageable. C'est pourquoi la méthode selon l'invention tire profit des mesures continues de la résistivité électrique de l'échantillon qui donne "l'allure"
de la courbe V(t). A partir de couples connus pour ces deux types de valeurs, il est possible de déterminer une loi reliant la résistivité au volume de la production de fluide déplacé.
Cette détermination peut se faire, par exemple, à l'aide d'une régression car cette relation est linéaire ou quasi linéaire, ce qui permet de coininettre moins d'erreur sur cette estimation que sur l'estimation directe de la courbe V(t). Ensuite, on applique cette loi pour transformer les mesures continues de résistivité en courbe continue de volume de production du fluide déplacé en fonction du temps V(t).
Il est classique dans le domaine pétrophysique d'exprimer les différents paramètres en fonction de la saturation d'un des fluides. Ainsi, que ce soit pour représenter les pressions capillaires, les perméabilités relatives, les mesures de résistivité, on représente ces grandeurs en fonction de la saturation car c'est une grandeur qui ne dépend pas du volume poreux. Cela permet ainsi de comparer facilement les résultats obtenus sur différents échantillons dont les propriétés macroscopiques ne sont jamais exactement les mêmes.
C'est pourquoi, selon un mode de réalisation, on établit au préalable de la courbe continue de volume dë production du fluide déplacé en fonction du temps V(t) une courbe de la résistivité en fonction de la saturation.
Pour ce faire, à partir des valeurs de volume de production de fluide déplacé
5 V(tl ), V(t2 ), ... mesurées à la fin de chaque période d'injection, et à
partir du volume de fluide inj ecté ( V,, (tl ), Vw (t2 ), ... ), on déduit des valeurs de saturation moyenne en fluide injecté ( Sw) et en fluide déplacé ( S, ) :
SN,=1-So=Sti,,;+VP
avec : S,,,, la saturation initiale en fluide injecté
10 VP le volume poreux de l'échantillon (mesuré au préalable, par exemple par une technique de RMN) A partir de couples connus de saturation et de résistivité il est possible de déterminer, une loi reliant la résistivité à la saturation. Cette détermination peut se faire, par exemple, à l'aide d'une régression linéaire ou d'une régression non linéaire. Ensuite,
15 on applique cette loi pour transformer les mesures continues de résistivité
en courbe continue de saturation en fonction du temps.
Puis à partir de cette courbe, on en déduit logiquement la courbe continue de volume de produçtion du fluide déplacé en fonction du temps.
Il faut noter que pour établir de telles courbes, il est nécessaire que la valeur ponctuelle de résistivité électrique utilisée soit la moins sensible possible aux profils de saturation de l'échantillon. C'est pourquoi ces relations sont obtenues de façon préférentielle avec une géométrie radiale d'électrodes de mesures et un échantillon court.
Étape 3= estimation des paramètres d'écoulements polyphasigues On estilne des paramètres d'écoulements polyphasiques à partir des mesures continues de pression différentielle et de la courbe continue de volume de production du fluide déplacé en fonction du temps. De nombreuses techniques, bien connues de
en courbe continue de saturation en fonction du temps.
Puis à partir de cette courbe, on en déduit logiquement la courbe continue de volume de produçtion du fluide déplacé en fonction du temps.
Il faut noter que pour établir de telles courbes, il est nécessaire que la valeur ponctuelle de résistivité électrique utilisée soit la moins sensible possible aux profils de saturation de l'échantillon. C'est pourquoi ces relations sont obtenues de façon préférentielle avec une géométrie radiale d'électrodes de mesures et un échantillon court.
Étape 3= estimation des paramètres d'écoulements polyphasigues On estilne des paramètres d'écoulements polyphasiques à partir des mesures continues de pression différentielle et de la courbe continue de volume de production du fluide déplacé en fonction du temps. De nombreuses techniques, bien connues de
16 l'homme du métier, permettent de réaliser cette étape. On peut par exemple utiliser la méthode décrite par Egermann et Lenormand :
Egermann P. and Lenormand R., "A new methodology to evaluate the ifnpact of the local heterogeneity on petrophysical parameters (lu-, Pc) : application on carbonate rocks", SCA n 2004-18, Sciety of Core Analysts, Abu Dhabi, 2004.
Cependant, au niveau de l'interprétation des mesures obtenues à partir du dispositif selon l'invention, le point délicat relève de la longueur courte de l'échantillon.
Les échantillons classiquement utilisés pour déterminer les paramètres d'écoulements polyphasiques sont généraleinent plus longs (6 - 10 cin). Cette courte longueur requiert une prise en compte complète des effets capillaires pour obtenir des résultats cohérents et représentatifs en terme de perméabilités relatives et permettant de prendre en coinpte la non uniformité du profil de saturation. La procédure d'injection la plus adaptée repose sur une injection multi-débit avec attente des stabilisations des volumes de production de fluide déplacé, qui perniet de collecter des informations sur la pression capillaire et les perméabilités relatives de la phase injectée de manière analytique avant le passage sur simulateur numérique.
Interprétation analytique des résultats :
L'étape d'interprétation analytique des données expérimentales est décrite dans la demande de brevet 04/09.273 du demandeur. L'état de stabilisation (du volume de production du fluide déplacé) correspondant à un équilibre entre les forces visqueuses et les forces capillaires, il devient possible d'évaluer, pour une valeur de saturation donnée, à la fois la valeur moyenne de la pression capillaire à l'échelle de l'échantillon analysé, notée Pc, la valeur moyenne de la perméabilité relative au fluide injecté, notée kr;n,, ainsi que la valeur moyenne de la perméabilité relative au, fluide déplacé, notée kr .eP .
Pour ce faire, on peut par exemple utiliser les expressions développées dans le cadre de la méthode semi-dynamique décrite dans le document FR 2 708 742 (US 5 679 885).
Dans cette méthode, on détermine la'saturation en entrée de l'échantillon (S) à l'aide de la formule suivante :
S'(t~P')=S+q ~S (1) q
Egermann P. and Lenormand R., "A new methodology to evaluate the ifnpact of the local heterogeneity on petrophysical parameters (lu-, Pc) : application on carbonate rocks", SCA n 2004-18, Sciety of Core Analysts, Abu Dhabi, 2004.
Cependant, au niveau de l'interprétation des mesures obtenues à partir du dispositif selon l'invention, le point délicat relève de la longueur courte de l'échantillon.
Les échantillons classiquement utilisés pour déterminer les paramètres d'écoulements polyphasiques sont généraleinent plus longs (6 - 10 cin). Cette courte longueur requiert une prise en compte complète des effets capillaires pour obtenir des résultats cohérents et représentatifs en terme de perméabilités relatives et permettant de prendre en coinpte la non uniformité du profil de saturation. La procédure d'injection la plus adaptée repose sur une injection multi-débit avec attente des stabilisations des volumes de production de fluide déplacé, qui perniet de collecter des informations sur la pression capillaire et les perméabilités relatives de la phase injectée de manière analytique avant le passage sur simulateur numérique.
Interprétation analytique des résultats :
L'étape d'interprétation analytique des données expérimentales est décrite dans la demande de brevet 04/09.273 du demandeur. L'état de stabilisation (du volume de production du fluide déplacé) correspondant à un équilibre entre les forces visqueuses et les forces capillaires, il devient possible d'évaluer, pour une valeur de saturation donnée, à la fois la valeur moyenne de la pression capillaire à l'échelle de l'échantillon analysé, notée Pc, la valeur moyenne de la perméabilité relative au fluide injecté, notée kr;n,, ainsi que la valeur moyenne de la perméabilité relative au, fluide déplacé, notée kr .eP .
Pour ce faire, on peut par exemple utiliser les expressions développées dans le cadre de la méthode semi-dynamique décrite dans le document FR 2 708 742 (US 5 679 885).
Dans cette méthode, on détermine la'saturation en entrée de l'échantillon (S) à l'aide de la formule suivante :
S'(t~P')=S+q ~S (1) q
17 où :- S' est la saturation à l'entrée - q est le débit du fluide injecté dans l'échantillon - AP' est la pression différentielle mesurée à l'entrée de l'échantillon - S est la saturation moyenne d'un des deux fluides. S peut être calculée à partir du volume de la production du fluide déplacé mesurée à la sortie de l' échantillon (V(0, V (tZ ), ... ).
A la fin de chacune des périodes d'injection, lorsque la stabilisationest atteinte, il est important de noter que la pression au niveau du fluide déplacé est uniforme et égale à la pression de sortie (P,,*,P). Ainsi, la pression locale du fluide injecté
correspond directement à la pression capillaire plus la pression de sortie. Dans le cas où la stabilisation du volume de production du fluide déplacé n'est pas tout à fait atteinte (pression dans le fluide déplacé non uniforme), il convient de noter que les profils de pression dans les deux fluides sont lisses et qu'il est aussi possible d'utiliser directement les profils de pression capillaire simulés pour déduire les propriétés locales.
Ainsi, à la sortie de l'échantillon, on considère que la pression du fluide injecté est égale à la pression du fluide déplacé. Dans la mesure où la pression capillaire se définit comme la différence entre ces deux pressions, la pression capillaire en sortie (P.') est nulle : P,' = 0.
A l'entrée de l'échantillon, si la pression en fluide déplacée est nulle (P,,',p = 0), la pression différentielle à l'entrée, qui correspond à la pression en entrée du fluide injecté, est égale à la pression capillaire en entrée ( P,): P,' = AP' (P,' = AP' +
Pdep si Pep ~0 ).
Ainsi, à chaque palier, on calcule une saturation en entrée S', et on mesure la pression différentielle en entrée qui correspond à la pression capillaire en entrée. On obtient de cette façon un ensemble de couple, un pour chaque palier, (S' ;
P~i).
A la fin de chacune des périodes d'injection, lorsque la stabilisationest atteinte, il est important de noter que la pression au niveau du fluide déplacé est uniforme et égale à la pression de sortie (P,,*,P). Ainsi, la pression locale du fluide injecté
correspond directement à la pression capillaire plus la pression de sortie. Dans le cas où la stabilisation du volume de production du fluide déplacé n'est pas tout à fait atteinte (pression dans le fluide déplacé non uniforme), il convient de noter que les profils de pression dans les deux fluides sont lisses et qu'il est aussi possible d'utiliser directement les profils de pression capillaire simulés pour déduire les propriétés locales.
Ainsi, à la sortie de l'échantillon, on considère que la pression du fluide injecté est égale à la pression du fluide déplacé. Dans la mesure où la pression capillaire se définit comme la différence entre ces deux pressions, la pression capillaire en sortie (P.') est nulle : P,' = 0.
A l'entrée de l'échantillon, si la pression en fluide déplacée est nulle (P,,',p = 0), la pression différentielle à l'entrée, qui correspond à la pression en entrée du fluide injecté, est égale à la pression capillaire en entrée ( P,): P,' = AP' (P,' = AP' +
Pdep si Pep ~0 ).
Ainsi, à chaque palier, on calcule une saturation en entrée S', et on mesure la pression différentielle en entrée qui correspond à la pression capillaire en entrée. On obtient de cette façon un ensemble de couple, un pour chaque palier, (S' ;
P~i).
18 On détermine kr;,,j à l'aide de la formule suivante, pour une saturation d'entrée donnée et donc un AP' donné (pour chaque palier) :
7crõ, = (2) KA dOP' où :,u est la viscosité du fluide injecté
L est la longueur de l'échantillon K est la perméabilité
A est la surface d'injection Pour déterminer krdep , on peut utiliser une interprétation de type JBN bien connue de l'homine du métier, et décrite par exemple dans le document suivant :
Johnson E.F., Bossler D.P. et Nauinann V.O. :" Calculation of relative perineability fi=om displacement experiments ", Trans. AIME, 216, 370-372, 1959..
Cette métliode nécessite de mesurer la saturation en sortie S .
On obtient ainsi pour chaque palier d'injection, une saturation S' , Pci(S' ), Icr;,, (S' ) et Iwdep (S ). L'ensemble des perméabilités relatives au fluide déplacé et au fluide injecté sera noté krs.
Interprétation numérique des résultats :
Après avoir évalué à chaque palier, à la fois une valeur de saturation ( S' ), une valeur moyenne de pression capillaire (P,'(S' )), une valeur de perméabilité
relative au fluide injecté (kr;,,(S' )) et une valeur de penuéabilité relative au fluide déplacé
lkrdep (S )), on réalise une étape d'interprétation numérique des mesures expérimentales.
On procède alors de la façon suivante :
7crõ, = (2) KA dOP' où :,u est la viscosité du fluide injecté
L est la longueur de l'échantillon K est la perméabilité
A est la surface d'injection Pour déterminer krdep , on peut utiliser une interprétation de type JBN bien connue de l'homine du métier, et décrite par exemple dans le document suivant :
Johnson E.F., Bossler D.P. et Nauinann V.O. :" Calculation of relative perineability fi=om displacement experiments ", Trans. AIME, 216, 370-372, 1959..
Cette métliode nécessite de mesurer la saturation en sortie S .
On obtient ainsi pour chaque palier d'injection, une saturation S' , Pci(S' ), Icr;,, (S' ) et Iwdep (S ). L'ensemble des perméabilités relatives au fluide déplacé et au fluide injecté sera noté krs.
Interprétation numérique des résultats :
Après avoir évalué à chaque palier, à la fois une valeur de saturation ( S' ), une valeur moyenne de pression capillaire (P,'(S' )), une valeur de perméabilité
relative au fluide injecté (kr;,,(S' )) et une valeur de penuéabilité relative au fluide déplacé
lkrdep (S )), on réalise une étape d'interprétation numérique des mesures expérimentales.
On procède alors de la façon suivante :
19 - on détennine, par exeznple à l'aide d'expériences de centrifugation bien connues de l'hoinme du métier, une courbe de valeurs de pression capillaire moyenne en fonction de la saturation, sur un échantillon voisin de l'échantillon étudié
et de propriétés siinilaires. On note Pc"' (S) cette courbe ;
- on affine la courbe P." (S) . Une première méthode consiste à prendre en compte les variations de perméabilité, de porosité et de la tension interfaciale par le biais de la fonction de Leverett. Cette fonction est décrite dans les documents suivants :
- Leverett M.C. :"Capillary belaavior in porous solids' ; Trans. AIME, 142, 152-168, 1941.
- Leverett M. C., Lewis W.B. et True M.E. :"Dimensional Model Studies of Oil-field behavior ", Trans. AIME, 146, 175-193, 1942.
Une seconde méthode consiste à utiliser les résultats de l'étape d'interprétation analytique de l'expérience de déplacement multi-débits. En effet, cette étape permettant de calculer des valeurs de pression capillaire ( Pc(S' )) pour un certain nombre de saturation (pression capillaire d'entrée et saturation d'entrée pour chaque palier d'injection), il est possible de normaliser la courbe de pression capillaire P,'~'(S) à partir de ces points. C'est-à-dire que l'on déforme la courbe P.7 (S) pour qu'elle passe par ces points précisément calculés sur l'échantillon étudié. Cette courbe normalisée est notée PC (S) - Ensuite, à l'aide de P. (S) et du jeu initial de krs issu de l'interprétation analytique des données expérimentales, un simulateur d'écoulement est utilisé
pour caler les données expérimentales (AP'(t) et V(t)) en déformant les courbes de krs jusqu'à
obtenir le minimum d'une fonction objectif.
A l'issue de cette phase d'interprétation numérique, on dispose donc d'un jeu de lrs correspondant au calage de la pression différentielle ( OP' (t) ), des courbes de volume de production (V(t)) et qui prend en compte les effets capillaires d'extrémité.
Exemple d'application à un cas réel Données générales L'exemple traité a été conduit sur une roche carbonatée de faible perméabilité
en conditions de réservoir (101 C). Les caractéristiques de l'échantillon et des fluides 5 (huile et saumure) sont rassemblées dans le Tableau 1.
Longueur (cm) 3.12 Diamètre (cm) 3.98 Porosité (%) 20.6 VP (cm3) 8 Saturation initiale en saumure S; (%) 46 Perméabilité absolue (md) 0.13 Viscosité de l'huile à 101 C (cP) 1.8 Viscosité de la saumure à 101 C (cP) 0.45 Tableau 1 : caractéristiques de l'échantillon et des fluides Étape 1= expérience de déplacement de type multi-débits et mesures électrigues Trois périodes avec des débits d'injection de saumure distincts ont été
envisagées 10 au cours de l'expérience : 0.1, 0.2 et 1 cm3/h. Ces débits assez faibles ont été choisis de manière à monter progressivement en pression de manière à explorer de manière graduelle les propriétés capillaires de l'échantillon.
Étape 2= détermination d'une courbe continue de volume de production du fluide 15 déplacé en fonction du temps Dans cette exemple, la courbe continue de volume de production du fluide déplacé en fonction du temps a été obtenue à partir d'une courbe de la résistance Re (issue des mesures de résistivité) en fonction de la saturation.
La figure 5 présente une comparaison entre les points de saturation (déduits du voluine de la production du fluide déplacé) obtenus en fin de période d'injection pour un débit donné ainsi que les mesures de résistance correspondantes (1lI). On observe une très bonne corrélation entre ces deux mesures qui permet de proposer une loi de type puissance pour la courbe de la résistance Re en fonction de la saturation :
Re = 33.741U
La figure 5 représente la résistance (Re) en fonction de la saturation en saumure (Sw). Les cercles représentent les points de mesure (M) et la ligne continue représente la loi de type puissance choisie pour la courbe de la résistance Re en fonction de la saturation (CE).
La mesure continue de la résistance permet alors de calculer l'évolution de la saturation moyenne au cours du temps qui est ensuite convertie en terme de volume de production pour obtenir une représentation plus standard en terme de perméabilité
relative. La figure 6 illustre, pour chacun des trois débits, les courbes de volume de production (OP) en fonction du temps (T) ainsi obtenues et les courbes de pression différentielle (PD) obtenues par mesure continues. Les mesures ponctuelles du volume de la production de fluide déplacé V(tl ), V(t2 ) sont également indiquées au niveau des deux premiers paliers de stabilisation.
Les courbes obtenues sont tout à fait classiques d'une expérience d'injection multi-débit avec une production d'huile additionnelle à chaque incrément de débit ("bumps").
Après une augmentation au moment du changement de débit, la pression différentielle baisse ensuite de manière progressive à mesure que l'huile est produite. Les légères fluctuations de la pression différentielle sont dues à des effets de variation de température dans le laboratoire qui ont été identifiés et corrigés.
Étape 3= estimation des paramètres d'écouleinents bolyphasigues 1- Interprétation analytique A chaque palier, on calcule une saturation moyenne en saumure S, et on mesure la pression différentielle en entrée qui correspond à la pression capillaire en entrée (P.').
On obtient de cette façon un ensemble de couple, un pour chaque palier, (S.,;
P,'). Dans le cadre de cette expérience, on obtient donc trois points analytiques correspondant aux trois débits de la partie forcée de l'imbibition. La figure 7 représente ces trois points de pression capillaire obtenus de façon analytique (PCA) en fonction de la saturation (S,õ).
L'arrondi suggéré par les points analytiques montre que l'échantillon possède au moins une mouillabilité partielle à l'huile. Les points ont servi de base au calage d'une loi analytique (PCHM) de pression capillaire de type puissance de manière à
introduire une courbe continue dans un simulateur d'écoulement pour conduire la simulation.
Cette courbe est également illustrée sur la figure 7.
Le calcul de kr;,, de façon analytique ( Iwa,õ,) requiert les valeurs de la dérivée du débit d'injection q par rapport à la pression différentielle. La figure 8A
montre les trois couples expérimentaux (CE) de valeurs de débit d'injection (q) et de pression différentielle (DP) ainsi . que le résultat d'une expression analytique de type exponentielle (LA) pour éviter de procéder à des calculs manuels de la dérivée. Les résultats obtenus en terme de krõ, analytiques ( kra;,,) montrent sur la figure 8B que le point limite de kra;,, à Sorw (saturation en huile résiduelle) reste faible, inférieur à 0.1.
2- Interprétation numérique L'ensemble des parties précédentes a servi de base pour compléter un jeu de données pour le simulateur d'écoulement en vue d'une inversion des perméabilités relatives (la~s). L'inversion permet d'obtenir un jeu de pennéabilités relative simulées (krs) en bon accord avec les données expérimentales. La figure 9 illustre ce calage des données expérimentales. Cette figure reprend la figure 6 sur laquelle ont été
représentées, pour chacun des trois débits, les courbes de volume de production simulées (OPS) en fonction du temps et les courbes de pression différentielle simulées (PDS).
Les courbes de perméabilité relatives (k7 s) obtenues après simulation ont une allure tout à fait classique comme on peut l'observer sur la figure 10. Le point initiale de la courbe de Irsdep est assez bas (0.14) du fait de la présence d'une saturation importante en saumure en début d'expérience (0.46). On observe par ailleurs une bonne cohérence avec les points de ki s;,, calculés de manière analytique (1ra;,,,).
Validation On représente les profils de saturation simulés (P01, 02 et P1) pour évaluer l'importance de l'erreur de calcul sur la saturation moyenne dans l'échantillon à partir de la mesure de résistance. En effet, tout est basé sur une calibration RI(Sw) ;
or la courbe RI(Sw) est établie à partir de la mesure de la résistance moyenne et de la saturation moyenne en présence d'un profil de saturation. On souhaite évaluer si la présence de ce profil de saturation, qui change de forme suivant le débit d'injection, joue beaucoup sur la calibration. Ces profils obtenus par simulation sont représentés sur la figure 11, où L
représente la distance par rapport à la face d'injection de l'échantillon, et P01 représente le profil simulé pour un débit de 0,1 cm3/h, P02 représente le profil simulé
pour un débit de 0,2 cm3/h etP 1 représente le profil siinulé pour un débit de 1 cm3/h.
Le tableau 2 ci-après montre la comparaison entre les saturations moyennes mesurées expérimentalement à la fin de chaque période de débit et les valeurs de saturation obtenues en moyennant le profil siinulé avec un systèlne en parallèle, c'est-à-dire tel que le voient les électrodes placées en radial. On compare donc ce que l'on mesure dans la burette en sortie (variation totale de la saturation) avec ce que mesure réellement les électrodes électriques au milieu de l'échantillon (entre le premier et le second centimètre). On constate que l'écart entre les deux mesures n'excède pas 1% ce qui reste tout à fait acceptable par rapport à la précision du calage et ce qui valide l'approche retenue pour évaluer la courbe de volume de production à partir de la mesure de résistance.
Sw Sw simulée sur Erreur expérimentale (1-2 cm) Fin d'injection à 0.1 cm3/h 0.66 0.658 -0.2 %
Fin d'injection à 0.2 cm3/h 0.695 0.69 -0.5 %
Fin d'injection à 1 cm3/h 0.744 0.748 0.4 %
Tableau 2 : évaluation de la pertinence de la courbe de la résistance en fonction de la saturation utilisée
et de propriétés siinilaires. On note Pc"' (S) cette courbe ;
- on affine la courbe P." (S) . Une première méthode consiste à prendre en compte les variations de perméabilité, de porosité et de la tension interfaciale par le biais de la fonction de Leverett. Cette fonction est décrite dans les documents suivants :
- Leverett M.C. :"Capillary belaavior in porous solids' ; Trans. AIME, 142, 152-168, 1941.
- Leverett M. C., Lewis W.B. et True M.E. :"Dimensional Model Studies of Oil-field behavior ", Trans. AIME, 146, 175-193, 1942.
Une seconde méthode consiste à utiliser les résultats de l'étape d'interprétation analytique de l'expérience de déplacement multi-débits. En effet, cette étape permettant de calculer des valeurs de pression capillaire ( Pc(S' )) pour un certain nombre de saturation (pression capillaire d'entrée et saturation d'entrée pour chaque palier d'injection), il est possible de normaliser la courbe de pression capillaire P,'~'(S) à partir de ces points. C'est-à-dire que l'on déforme la courbe P.7 (S) pour qu'elle passe par ces points précisément calculés sur l'échantillon étudié. Cette courbe normalisée est notée PC (S) - Ensuite, à l'aide de P. (S) et du jeu initial de krs issu de l'interprétation analytique des données expérimentales, un simulateur d'écoulement est utilisé
pour caler les données expérimentales (AP'(t) et V(t)) en déformant les courbes de krs jusqu'à
obtenir le minimum d'une fonction objectif.
A l'issue de cette phase d'interprétation numérique, on dispose donc d'un jeu de lrs correspondant au calage de la pression différentielle ( OP' (t) ), des courbes de volume de production (V(t)) et qui prend en compte les effets capillaires d'extrémité.
Exemple d'application à un cas réel Données générales L'exemple traité a été conduit sur une roche carbonatée de faible perméabilité
en conditions de réservoir (101 C). Les caractéristiques de l'échantillon et des fluides 5 (huile et saumure) sont rassemblées dans le Tableau 1.
Longueur (cm) 3.12 Diamètre (cm) 3.98 Porosité (%) 20.6 VP (cm3) 8 Saturation initiale en saumure S; (%) 46 Perméabilité absolue (md) 0.13 Viscosité de l'huile à 101 C (cP) 1.8 Viscosité de la saumure à 101 C (cP) 0.45 Tableau 1 : caractéristiques de l'échantillon et des fluides Étape 1= expérience de déplacement de type multi-débits et mesures électrigues Trois périodes avec des débits d'injection de saumure distincts ont été
envisagées 10 au cours de l'expérience : 0.1, 0.2 et 1 cm3/h. Ces débits assez faibles ont été choisis de manière à monter progressivement en pression de manière à explorer de manière graduelle les propriétés capillaires de l'échantillon.
Étape 2= détermination d'une courbe continue de volume de production du fluide 15 déplacé en fonction du temps Dans cette exemple, la courbe continue de volume de production du fluide déplacé en fonction du temps a été obtenue à partir d'une courbe de la résistance Re (issue des mesures de résistivité) en fonction de la saturation.
La figure 5 présente une comparaison entre les points de saturation (déduits du voluine de la production du fluide déplacé) obtenus en fin de période d'injection pour un débit donné ainsi que les mesures de résistance correspondantes (1lI). On observe une très bonne corrélation entre ces deux mesures qui permet de proposer une loi de type puissance pour la courbe de la résistance Re en fonction de la saturation :
Re = 33.741U
La figure 5 représente la résistance (Re) en fonction de la saturation en saumure (Sw). Les cercles représentent les points de mesure (M) et la ligne continue représente la loi de type puissance choisie pour la courbe de la résistance Re en fonction de la saturation (CE).
La mesure continue de la résistance permet alors de calculer l'évolution de la saturation moyenne au cours du temps qui est ensuite convertie en terme de volume de production pour obtenir une représentation plus standard en terme de perméabilité
relative. La figure 6 illustre, pour chacun des trois débits, les courbes de volume de production (OP) en fonction du temps (T) ainsi obtenues et les courbes de pression différentielle (PD) obtenues par mesure continues. Les mesures ponctuelles du volume de la production de fluide déplacé V(tl ), V(t2 ) sont également indiquées au niveau des deux premiers paliers de stabilisation.
Les courbes obtenues sont tout à fait classiques d'une expérience d'injection multi-débit avec une production d'huile additionnelle à chaque incrément de débit ("bumps").
Après une augmentation au moment du changement de débit, la pression différentielle baisse ensuite de manière progressive à mesure que l'huile est produite. Les légères fluctuations de la pression différentielle sont dues à des effets de variation de température dans le laboratoire qui ont été identifiés et corrigés.
Étape 3= estimation des paramètres d'écouleinents bolyphasigues 1- Interprétation analytique A chaque palier, on calcule une saturation moyenne en saumure S, et on mesure la pression différentielle en entrée qui correspond à la pression capillaire en entrée (P.').
On obtient de cette façon un ensemble de couple, un pour chaque palier, (S.,;
P,'). Dans le cadre de cette expérience, on obtient donc trois points analytiques correspondant aux trois débits de la partie forcée de l'imbibition. La figure 7 représente ces trois points de pression capillaire obtenus de façon analytique (PCA) en fonction de la saturation (S,õ).
L'arrondi suggéré par les points analytiques montre que l'échantillon possède au moins une mouillabilité partielle à l'huile. Les points ont servi de base au calage d'une loi analytique (PCHM) de pression capillaire de type puissance de manière à
introduire une courbe continue dans un simulateur d'écoulement pour conduire la simulation.
Cette courbe est également illustrée sur la figure 7.
Le calcul de kr;,, de façon analytique ( Iwa,õ,) requiert les valeurs de la dérivée du débit d'injection q par rapport à la pression différentielle. La figure 8A
montre les trois couples expérimentaux (CE) de valeurs de débit d'injection (q) et de pression différentielle (DP) ainsi . que le résultat d'une expression analytique de type exponentielle (LA) pour éviter de procéder à des calculs manuels de la dérivée. Les résultats obtenus en terme de krõ, analytiques ( kra;,,) montrent sur la figure 8B que le point limite de kra;,, à Sorw (saturation en huile résiduelle) reste faible, inférieur à 0.1.
2- Interprétation numérique L'ensemble des parties précédentes a servi de base pour compléter un jeu de données pour le simulateur d'écoulement en vue d'une inversion des perméabilités relatives (la~s). L'inversion permet d'obtenir un jeu de pennéabilités relative simulées (krs) en bon accord avec les données expérimentales. La figure 9 illustre ce calage des données expérimentales. Cette figure reprend la figure 6 sur laquelle ont été
représentées, pour chacun des trois débits, les courbes de volume de production simulées (OPS) en fonction du temps et les courbes de pression différentielle simulées (PDS).
Les courbes de perméabilité relatives (k7 s) obtenues après simulation ont une allure tout à fait classique comme on peut l'observer sur la figure 10. Le point initiale de la courbe de Irsdep est assez bas (0.14) du fait de la présence d'une saturation importante en saumure en début d'expérience (0.46). On observe par ailleurs une bonne cohérence avec les points de ki s;,, calculés de manière analytique (1ra;,,,).
Validation On représente les profils de saturation simulés (P01, 02 et P1) pour évaluer l'importance de l'erreur de calcul sur la saturation moyenne dans l'échantillon à partir de la mesure de résistance. En effet, tout est basé sur une calibration RI(Sw) ;
or la courbe RI(Sw) est établie à partir de la mesure de la résistance moyenne et de la saturation moyenne en présence d'un profil de saturation. On souhaite évaluer si la présence de ce profil de saturation, qui change de forme suivant le débit d'injection, joue beaucoup sur la calibration. Ces profils obtenus par simulation sont représentés sur la figure 11, où L
représente la distance par rapport à la face d'injection de l'échantillon, et P01 représente le profil simulé pour un débit de 0,1 cm3/h, P02 représente le profil simulé
pour un débit de 0,2 cm3/h etP 1 représente le profil siinulé pour un débit de 1 cm3/h.
Le tableau 2 ci-après montre la comparaison entre les saturations moyennes mesurées expérimentalement à la fin de chaque période de débit et les valeurs de saturation obtenues en moyennant le profil siinulé avec un systèlne en parallèle, c'est-à-dire tel que le voient les électrodes placées en radial. On compare donc ce que l'on mesure dans la burette en sortie (variation totale de la saturation) avec ce que mesure réellement les électrodes électriques au milieu de l'échantillon (entre le premier et le second centimètre). On constate que l'écart entre les deux mesures n'excède pas 1% ce qui reste tout à fait acceptable par rapport à la précision du calage et ce qui valide l'approche retenue pour évaluer la courbe de volume de production à partir de la mesure de résistance.
Sw Sw simulée sur Erreur expérimentale (1-2 cm) Fin d'injection à 0.1 cm3/h 0.66 0.658 -0.2 %
Fin d'injection à 0.2 cm3/h 0.695 0.69 -0.5 %
Fin d'injection à 1 cm3/h 0.744 0.748 0.4 %
Tableau 2 : évaluation de la pertinence de la courbe de la résistance en fonction de la saturation utilisée
Claims (10)
1) Méthode pour évaluer conjointement la résistivité électrique et des paramètres d'écoulement d'un milieu poreux, à partir d'au moins une expérience de déplacement de type multi-débits consistant à injecter un premier fluide, appelé fluide injecté, dans un échantillon dudit milieu initialement saturé avec un second fluide, appelé fluide déplacé, dans laquelle le débit d'injection du fluide injecté
est augmenté
au moins une fois après un palier de stabilisation du volume de la production du fluide déplacé, ladite méthode comporte les étapes suivantes :
- on mesure de façon continue en fonction du temps la pression différentielle entre deux extrémités dudit échantillon (.DELTA.P i(t)) et la résistivité
électrique dudit échantillon (R(t)) ;
- on mesure le volume de production du fluide déplacé à chaque palier de stabilisation (V(t1), V(t2));
- on estime une relation entre la résistivité électrique de l'échantillon et le volume de la production du fluide déplacé à partir desdites mesures de résistivités (R(t)) et desdites valeurs de volume de production du fluide déplacé (V(t1), V(t2));
- on établit une courbe continue de volume de production du fluide déplacé en fonction du temps (V(t)) à partir de ladite relation ;
- on estime des paramètres d'écoulements à partir de ladite mesure continue de pression différentielle (.DELTA.P i(t)) et de ladite courbe continue (V(t)) à l'aide d'un simulateur numérique d'écoulement.
est augmenté
au moins une fois après un palier de stabilisation du volume de la production du fluide déplacé, ladite méthode comporte les étapes suivantes :
- on mesure de façon continue en fonction du temps la pression différentielle entre deux extrémités dudit échantillon (.DELTA.P i(t)) et la résistivité
électrique dudit échantillon (R(t)) ;
- on mesure le volume de production du fluide déplacé à chaque palier de stabilisation (V(t1), V(t2));
- on estime une relation entre la résistivité électrique de l'échantillon et le volume de la production du fluide déplacé à partir desdites mesures de résistivités (R(t)) et desdites valeurs de volume de production du fluide déplacé (V(t1), V(t2));
- on établit une courbe continue de volume de production du fluide déplacé en fonction du temps (V(t)) à partir de ladite relation ;
- on estime des paramètres d'écoulements à partir de ladite mesure continue de pression différentielle (.DELTA.P i(t)) et de ladite courbe continue (V(t)) à l'aide d'un simulateur numérique d'écoulement.
2) Méthode selon la revendication 1, dans laquelle ladite relation est estimée en outre à partir d'une courbe de la résistivité en fonction de la saturation de l'échantillon en fin de palier.
3) Méthode selon l'une des revendications précédentes, dans laquelle on dispose radialement des électrodes de mesure de la résistivité électrique autour dudit échantillon et l'on utilise un échantillon de faible longueur, typiquement de l'ordre de 3 cm, pour mesurer ladite résistivité électrique.
4) Dispositif pour réaliser sur un échantillon solide poreux (S) des expériences de déplacements multi-débits, comportant une cellule de confinement (1) pour contenir ledit échantillon (S) saturé avec un premier fluide, appelé fluide déplacé, et fermée à
une extrémité par un premier embout cylindrique (16) relié à des moyens d'injection sous pression (5a, 5b) d'un second fluide, appelé fluide injecté, et fermée à
l'autre extrémité par un second embout cylindrique (15) comportant un système de circulation des fluides (4a, 4b, 4c), ledit dispositif comportant également une plaque poreuse semi-perméable (C) placée entre ledit échantillon et ledit premier embout, des moyens (6) de mesures continues de résistivité électrique dudit échantillon, des moyens (3) de mesures du volume de production du fluide déplacé pour chaque palier de stabilisation, ledit dispositif étant caractérisé en ce qu'il comporte des moyens (23, 22) de mesures continues de la pression différentielle d'un des fluides de part et d'autre de l'échantillon, entre l'échantillon et la plaque poreuse semi-perméable.
une extrémité par un premier embout cylindrique (16) relié à des moyens d'injection sous pression (5a, 5b) d'un second fluide, appelé fluide injecté, et fermée à
l'autre extrémité par un second embout cylindrique (15) comportant un système de circulation des fluides (4a, 4b, 4c), ledit dispositif comportant également une plaque poreuse semi-perméable (C) placée entre ledit échantillon et ledit premier embout, des moyens (6) de mesures continues de résistivité électrique dudit échantillon, des moyens (3) de mesures du volume de production du fluide déplacé pour chaque palier de stabilisation, ledit dispositif étant caractérisé en ce qu'il comporte des moyens (23, 22) de mesures continues de la pression différentielle d'un des fluides de part et d'autre de l'échantillon, entre l'échantillon et la plaque poreuse semi-perméable.
5) Dispositif selon la revendication 4, dans lequel ledit premier embout cylindrique comporte un premier cylindre placé (18) dans un second cylindre concentrique (17).
6) Dispositif selon l'une des revendications 4 et 5, dans lequel lesdits moyens de mesures continues de la pression différentielle comporte au moins un conduit de prise de pression (22) relié à un capteur différentiel de pression (23).
7) Dispositif selon la revendication 6, dans lequel ledit conduit comporte un ensemble de rainures longitudinales (27).
8) Dispositif selon la revendication 7, dans lequel ledit ensemble de rainures longitudinales (27) est usiné sur la surface externe dudit premier cylindre (18).
9) Dispositif selon la revendication 7, dans lequel ledit ensemble de rainures longitudinales (27) est usiné sur la surface interne dudit second cylindre (17).
10) Dispositif selon l'une des revendications 6 et 7, dans lequel ledit conduit de prise de pression est adapté à éviter le piégeage d'air lors de l'installation dudit premier embout.
Applications Claiming Priority (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| FR0504272A FR2884923B1 (fr) | 2005-04-26 | 2005-04-26 | Methode et dispositif pour evaluer des parametres d'ecoulement et des parametres electriques d'un milieu poreux. |
| FR0504272 | 2005-04-26 | ||
| PCT/FR2006/000911 WO2006114509A1 (fr) | 2005-04-26 | 2006-04-24 | Methode et dispositif pour evaluer des parametres d'ecoulement et des parametres electriques d'un milieu poreux |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CA2605741A1 true CA2605741A1 (fr) | 2006-11-02 |
Family
ID=35500978
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CA002605741A Abandoned CA2605741A1 (fr) | 2005-04-26 | 2006-04-24 | Methode et dispositif pour evaluer des parametres d'ecoulement et des parametres electriques d'un milieu poreux |
Country Status (5)
| Country | Link |
|---|---|
| US (1) | US8024960B2 (fr) |
| EP (1) | EP1877755A1 (fr) |
| CA (1) | CA2605741A1 (fr) |
| FR (1) | FR2884923B1 (fr) |
| WO (1) | WO2006114509A1 (fr) |
Families Citing this family (27)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| GB0716120D0 (en) * | 2007-08-17 | 2007-09-26 | Imp Innovations Ltd | Apparatus and method to measure properties of porous media |
| CA2821297C (fr) * | 2010-12-13 | 2019-11-12 | University Of Wyoming | Appareil et procede d'inondation de carottes a contre-pression de recirculation constante |
| US20120156787A1 (en) * | 2010-12-15 | 2012-06-21 | Saudi Arabian Oil Company | Laboratory Testing Procedure to Select Acid or Proppant Fracturing Stimulation Treatment for a Given Carbonate Formation |
| US20130002258A1 (en) * | 2011-06-30 | 2013-01-03 | Schlumberger Technology Corporation | Device for dielectric permittivity and resistivity high temperature measurement of rock samples |
| GB2492795B (en) * | 2011-07-12 | 2013-11-20 | Schlumberger Holdings | Testing imbibition of liquid |
| KR101214376B1 (ko) * | 2011-10-31 | 2012-12-24 | 한국지질자원연구원 | 미고결 시료 투과율 산정용 홀더 |
| CN104122181B (zh) * | 2013-04-26 | 2016-09-07 | 中国石油天然气集团公司 | 入井液对储层渗透率损害评价装置 |
| CN103267836B (zh) * | 2013-05-16 | 2014-12-10 | 西南石油大学 | 一种岩心窜流系数测试实验装置 |
| US10884084B2 (en) | 2013-09-05 | 2021-01-05 | Saudi Arabian Oil Company | Systems and methods for tri-axial NMR testing |
| US11241701B2 (en) | 2013-10-21 | 2022-02-08 | Saudi Arabian Oil Company | Tri-axial centrifuge apparatus with electrical sensor, acoustic sensor, and x-ray instrument |
| US10900945B2 (en) | 2013-10-21 | 2021-01-26 | Saudi Arabian Oil Company | Tri-axial centrifuge apparatus with electrical sensor, acoustic sensor, and X-ray instrument |
| CN103558109B (zh) * | 2013-10-28 | 2015-08-26 | 东北石油大学 | 一种测试电阻-含油饱和度标准关系曲线的方法 |
| RU2580177C1 (ru) * | 2014-12-30 | 2016-04-10 | Шлюмберже Текнолоджи Б.В. | Способ определения изменений параметров пористой среды под действием загрязнителя |
| CN105225583A (zh) * | 2015-11-17 | 2016-01-06 | 中国地质大学（武汉） | 一种孔隙介质中的渗流实验方法及模具 |
| CN106761511B (zh) * | 2016-12-02 | 2019-03-22 | 中国石油大学(华东) | 双重介质渗流窜流量表征方法及实验装置 |
| CN107023282A (zh) * | 2017-06-16 | 2017-08-08 | 北京永瑞达科贸有限公司 | 油气水三相驱替可视化实验装置 |
| CN109254134B (zh) * | 2017-07-12 | 2021-12-31 | 北京康普瑞基石油工程技术有限公司 | 一种测量岩石电阻指示油气突破压力的新方法和设备 |
| EP3743718A1 (fr) * | 2018-01-23 | 2020-12-02 | Total Se | Procédé et appareil d'analyse d'un échantillon de roche |
| FR3088118B1 (fr) * | 2018-11-05 | 2020-10-23 | Ifp Energies Now | Systeme et procede de tracage d’especes ioniques par mesures electriques multiples |
| US11353385B2 (en) | 2019-04-24 | 2022-06-07 | Saudi Arabian Oil Company | Testing petro-physical properties using a tri-axial pressure centrifuge apparatus |
| US11249065B2 (en) | 2019-04-24 | 2022-02-15 | Saudi Arabian Oil Company | Testing petro-physical properties using a tri-axial pressure centrifuge apparatus |
| CN112504936B (zh) * | 2020-11-30 | 2021-12-03 | 中国地质大学（北京） | 一种模拟研究深部煤层气渗透率试验装置及其试验方法 |
| DE112021000127T5 (de) * | 2021-03-30 | 2022-12-29 | Northwest Engineering Corporation Limited | Intelligentes mehrzweck-sickerverformungstestsystem und sein betriebsverfahren |
| CN113866053A (zh) * | 2021-09-15 | 2021-12-31 | 中国石油大学(北京) | 一种岩样自发渗吸驱油率的测定方法和装置 |
| WO2023214193A1 (fr) * | 2022-05-06 | 2023-11-09 | Totalenergies Onetech | Cellule de mesure de propriétés de roche et système de mesure de propriétés de roche associé |
| CN116539300B (zh) * | 2023-05-16 | 2024-03-19 | 中国石油大学(北京) | 控水阀性能评价方法 |
| CN116858749B (zh) * | 2023-07-26 | 2024-04-23 | 北京建工环境修复股份有限公司 | 一种阻隔帷幕复电阻率信号与渗透性关系的确定方法 |
Family Cites Families (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4679422A (en) * | 1986-08-28 | 1987-07-14 | The United States Of America As Represented By The Secretary Of The Interior | Method and apparatus for steady-state measurement of liquid conductivity in porous media |
| FR2652162B1 (fr) * | 1989-09-15 | 1991-10-31 | Inst Francais Du Petrole | Dispositif et methode d'evaluation de l'aptitude qu'a un corps a expulser un produit et leur application a un echantillon de roche. |
| US5105154A (en) * | 1991-03-19 | 1992-04-14 | Mobil Oil Corporation | Apparatus for measuring radial resistivities in cylindrical core samples of porous rock |
| US5209104A (en) * | 1992-01-23 | 1993-05-11 | Mobil Oil Corporation | Method for desaturating a porous rock for electrical resistivity measurements |
| FR2699282B1 (fr) * | 1992-12-15 | 1995-01-13 | Inst Francais Du Petrole | Procédé et dispositif perfectionnés pour l'étude de propriétés d'un matériau perméable. |
| US5297420A (en) * | 1993-05-19 | 1994-03-29 | Mobil Oil Corporation | Apparatus and method for measuring relative permeability and capillary pressure of porous rock |
| FR2708742B1 (fr) * | 1993-07-29 | 1995-09-01 | Inst Francais Du Petrole | Procédé et dispositiphi pour mesurer des paramètres physiques d'échantillons poreux mouillables par des fluides. |
| US5493226A (en) * | 1994-04-08 | 1996-02-20 | Mobile Oil Corporation | Method and apparatus for measuring properties of core samples including heating and pressurizing the core sample and measuring the dynamic and static capillary pressure of water in the core sample |
| FR2724460B1 (fr) * | 1994-09-09 | 1997-01-17 | Inst Francais Du Petrole | Dispositif de mesure petrophysique et mehode de mise en oeuvre |
| FR2731073B1 (fr) * | 1995-02-27 | 1997-04-04 | Inst Francais Du Petrole | Methode et dispositif pour determiner differents parametres physiques d'echantillons de materiaux poreux, en presence de fluides di- ou triphasiques |
| FR2758881B1 (fr) * | 1997-01-30 | 1999-02-26 | Inst Francais Du Petrole | Dispositif pour faire des mesures de la resistivite electrique d'un echantillon solide |
-
2005
- 2005-04-26 FR FR0504272A patent/FR2884923B1/fr not_active Expired - Fee Related
-
2006
- 2006-04-24 EP EP06743730A patent/EP1877755A1/fr not_active Withdrawn
- 2006-04-24 WO PCT/FR2006/000911 patent/WO2006114509A1/fr active Application Filing
- 2006-04-24 CA CA002605741A patent/CA2605741A1/fr not_active Abandoned
- 2006-04-24 US US11/912,512 patent/US8024960B2/en not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| US8024960B2 (en) | 2011-09-27 |
| FR2884923B1 (fr) | 2007-10-19 |
| US20090126462A1 (en) | 2009-05-21 |
| FR2884923A1 (fr) | 2006-10-27 |
| EP1877755A1 (fr) | 2008-01-16 |
| WO2006114509A1 (fr) | 2006-11-02 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CA2605741A1 (fr) | Methode et dispositif pour evaluer des parametres d'ecoulement et des parametres electriques d'un milieu poreux | |
| Gupta et al. | Intercept method—A novel technique to correct steady-state relative permeability data for capillary end effects | |
| CA2383764C (fr) | Methode pour modeliser des deplacements de fluides dans un milieu poreux tenant compte d'effets d'hysteresis | |
| FR2909448A1 (fr) | Methode pour caracteriser la distribution de la permeabilite absolue d'un echantillon heterogene | |
| CN107991335B (zh) | 致密砂岩水锁损害评价试验方法 | |
| CA2517244A1 (fr) | Methode de determination de parametres d'ecoulement multiphasique d'un milieu poreux avec prise en compte de l'heterogeneite locale | |
| CA2594086A1 (fr) | Methodes et appareillage permettant de caracteriser le fond de trou des fluides de formations | |
| FR2853071A1 (fr) | Methode et dispositif pour evaluer des parametres physiques d'un gisement souterrain a partir de debris de roche qui y sont preleves | |
| WO2017167567A1 (fr) | Dispositif pour la determination de parametres petrophysiques d'une formation souterraine | |
| FR2767919A1 (fr) | Procede et dispositif de debitmetrie pour effluents petroliers | |
| CN110296931B (zh) | 一种致密砂岩油水相对渗透率信息的表征方法及系统 | |
| CA2278449C (fr) | Methode de mesure rapide de l'indice de resistivite d'echantillons solides tels que des roches | |
| FR2890445A1 (fr) | Methode pour determiner la pression capillaire d'entree d'un milieu poreux | |
| FR2931189A1 (fr) | Procede d'estimation de parametres physiques d'une formation geologique | |
| FR2907555A1 (fr) | Methode et dispositif pour mesurer la pression minimale de miscibilite de deux phases | |
| EP1522850A2 (fr) | Méthode et dispositif pour mesurer l'anisotropie de résistivité d'échantillons de roche présentant des litages | |
| FR2864238A1 (fr) | Methode pour determiner la permeabilite d'un milieu souterrain a partir de mesures par rmn de la permeabilite de fragments de roche issus du milieu | |
| FR2855267A1 (fr) | Affinement de signaux pour l'analyse optique de fluides | |
| EP2034308A1 (fr) | Méthode de mesure rapide de la saturation et de la résistivité d'un milieu poreux | |
| FR2569762A1 (fr) | Procede d'essai de puits d'hydrocarbures | |
| EP2264355A1 (fr) | Méthode pour réduire la perte de charge d'un liquide en écoulement dans une conduite en tenant compte de la dégradation d'agents réducteurs de traînée | |
| Dick et al. | Measuring relative permeability with NMR | |
| WO2020094396A1 (fr) | Systeme et procede de traçage d'especes ioniques par mesures electriques multiples | |
| CN114577677B (zh) | 基于储层非均质性的致密油藏非线性渗流表征方法 | |
| FR3039275B1 (fr) | Dispositif de surveillance d'un produit liquide generant un gaz |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| EEER | Examination request | ||
| FZDE | Discontinued |
Effective date: 20140424 |