CA2319962C - Methode optimisee pour determiner des parametres physiques d`un echantillon soumis a centrifugation - Google Patents
Methode optimisee pour determiner des parametres physiques d`un echantillon soumis a centrifugation Download PDFInfo
- Publication number
- CA2319962C CA2319962C CA2319962A CA2319962A CA2319962C CA 2319962 C CA2319962 C CA 2319962C CA 2319962 A CA2319962 A CA 2319962A CA 2319962 A CA2319962 A CA 2319962A CA 2319962 C CA2319962 C CA 2319962C
- Authority
- CA
- Canada
- Prior art keywords
- saturation
- sample
- centrifugation
- curve
- parameters
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 32
- 238000005119 centrifugation Methods 0.000 title claims abstract description 20
- 239000012530 fluid Substances 0.000 claims abstract description 37
- 230000035699 permeability Effects 0.000 claims abstract description 23
- 229920006395 saturated elastomer Polymers 0.000 claims abstract description 7
- 238000004519 manufacturing process Methods 0.000 abstract description 16
- 230000006641 stabilisation Effects 0.000 abstract description 12
- 238000011105 stabilization Methods 0.000 abstract description 12
- 238000005259 measurement Methods 0.000 abstract description 10
- 238000004088 simulation Methods 0.000 abstract description 5
- 239000000523 sample Substances 0.000 description 39
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 6
- 238000012360 testing method Methods 0.000 description 5
- 239000012267 brine Substances 0.000 description 4
- 230000000052 comparative effect Effects 0.000 description 4
- 238000002474 experimental method Methods 0.000 description 4
- 239000011435 rock Substances 0.000 description 4
- HPALAKNZSZLMCH-UHFFFAOYSA-M sodium;chloride;hydrate Chemical compound O.[Na+].[Cl-] HPALAKNZSZLMCH-UHFFFAOYSA-M 0.000 description 4
- 230000001052 transient effect Effects 0.000 description 4
- 230000000694 effects Effects 0.000 description 3
- 239000011148 porous material Substances 0.000 description 3
- 230000001133 acceleration Effects 0.000 description 2
- 239000004020 conductor Substances 0.000 description 2
- 239000007788 liquid Substances 0.000 description 2
- JEYCTXHKTXCGPB-UHFFFAOYSA-N Methaqualone Chemical compound CC1=CC=CC=C1N1C(=O)C2=CC=CC=C2N=C1C JEYCTXHKTXCGPB-UHFFFAOYSA-N 0.000 description 1
- 239000000919 ceramic Substances 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 229930195733 hydrocarbon Natural products 0.000 description 1
- 150000002430 hydrocarbons Chemical class 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 239000002184 metal Substances 0.000 description 1
- 230000008520 organization Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000010200 validation analysis Methods 0.000 description 1
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
-
- 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
-
- 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
- G01N2015/0813—Measuring intrusion, e.g. of mercury
-
- 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
- G01N2015/0853—Investigating permeability, pore-volume, or surface area of porous materials by electrical capacitance measurement
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Health & Medical Sciences (AREA)
- Physics & Mathematics (AREA)
- Immunology (AREA)
- Analytical Chemistry (AREA)
- Pathology (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- Biochemistry (AREA)
- Medicinal Chemistry (AREA)
- Environmental & Geological Engineering (AREA)
- Remote Sensing (AREA)
- Food Science & Technology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geology (AREA)
- Fluid Mechanics (AREA)
- Dispersion Chemistry (AREA)
- Centrifugal Separators (AREA)
- Sampling And Sample Adjustment (AREA)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
- Investigating Or Analysing Biological Materials (AREA)
Abstract
- Méthode pour déterminer plus rapidement certains paramètres physiques (courbes de pression capillaire et perméabilités relatives par exemple) d'un échantillon géologique initialement saturé avec un premier fluide placé dans un récipient et soumis à centrifugation en présence d'un deuxième fluide. - Par le biais d'une modélisation analytique mettant en jeu des fonctions exponentielles, s'ajustant au mieux à une courbe de production effective, pour chaque palier de vitesse, on détermine les saturations d'équilibre bien avant d'atteindre la complète stabilisation. On utilise une base de données obtenues par diverses simulations, pour obtenir les courbes de perméabilités relatives correspondant à l'échantillon centrifugé. - Applications à des mesures de caractéristiques pétrophysiques.
Description
MÉTHODE OPTIMISÉE POUR DÉTERMINER DES PARAMETRES PHYSIQUES
D'UN ÉCHANTILLON SOUMIS A CENTRIFUGATION
La présente invention concerne une méthode de modélisation permettant d'obtenir plus rapidement certains paramètres physiques d'un échantillon géologique à partir de résultats expérimentaux obtenus en le soumettant à
centrifugation.
La méthode trouve des applications importantes dans le domaine de la pétrophysique où la connaissance de paramètres tels que la pression capillaire, les courbes de perméabilité relative des roches poreuses, permet aux ingénieurs de réservoirs, de caractériser les conditions d'écoulement d'hydrocarbures au sein d'un gisement, de faire des prévisions de production et d'élaborer des stratégies d'exploitation.
Une méthode bien connue pour déterminer rapidement à la fois les courbes de pression capillaires et les courbes de perméabilité relative d'un échantillon de roche poreuse initialement saturé par un premier fluide consiste essentiellement à la soumettre à centrifugation pour déplacer le fluide en place par un autre fluide.
On utilise par exemple des dispositifs de centrifugation tels que décrits dans les brevets FR 2 699 282 (US 5,463,894), FR 2 763 690, FR 2 772 477 du demandeur.
L'échantillon initialement saturé avec un liquide A (de la saumure par exemple) est placé dans un récipient ou godet allongé contenant un autre fluide B de densité différente (de l'huile par exemple. Le godet est fixé au bout d'un bras tournant, et on lui applique une force centrifuge de façon à étudier les déplacements des fluides dans l'échantillon au cours d'au moins deux phases distinctes.
Durant une première phase de drainage, on soumet alors l'ensemble à une force centrifuge dirigée suivant la longueur du récipient de manière à exercer sur lui une force la d'expulsion qui tend à faire sortir une partie du premier fluide A. Dans le même temps, du fluide A pénètre à l'intérieur de l'échantillon. Les deux fluides se déplacent à l'intérieur de l'échantillon jusqu'à une position d'équilibre où
la force due à la pression capillaire dans les pores, compense la force centrifuge exercée.
Une sonde de mesure, est placée dans le godet, sur le côté de l'échantillon.
La sonde
D'UN ÉCHANTILLON SOUMIS A CENTRIFUGATION
La présente invention concerne une méthode de modélisation permettant d'obtenir plus rapidement certains paramètres physiques d'un échantillon géologique à partir de résultats expérimentaux obtenus en le soumettant à
centrifugation.
La méthode trouve des applications importantes dans le domaine de la pétrophysique où la connaissance de paramètres tels que la pression capillaire, les courbes de perméabilité relative des roches poreuses, permet aux ingénieurs de réservoirs, de caractériser les conditions d'écoulement d'hydrocarbures au sein d'un gisement, de faire des prévisions de production et d'élaborer des stratégies d'exploitation.
Une méthode bien connue pour déterminer rapidement à la fois les courbes de pression capillaires et les courbes de perméabilité relative d'un échantillon de roche poreuse initialement saturé par un premier fluide consiste essentiellement à la soumettre à centrifugation pour déplacer le fluide en place par un autre fluide.
On utilise par exemple des dispositifs de centrifugation tels que décrits dans les brevets FR 2 699 282 (US 5,463,894), FR 2 763 690, FR 2 772 477 du demandeur.
L'échantillon initialement saturé avec un liquide A (de la saumure par exemple) est placé dans un récipient ou godet allongé contenant un autre fluide B de densité différente (de l'huile par exemple. Le godet est fixé au bout d'un bras tournant, et on lui applique une force centrifuge de façon à étudier les déplacements des fluides dans l'échantillon au cours d'au moins deux phases distinctes.
Durant une première phase de drainage, on soumet alors l'ensemble à une force centrifuge dirigée suivant la longueur du récipient de manière à exercer sur lui une force la d'expulsion qui tend à faire sortir une partie du premier fluide A. Dans le même temps, du fluide A pénètre à l'intérieur de l'échantillon. Les deux fluides se déplacent à l'intérieur de l'échantillon jusqu'à une position d'équilibre où
la force due à la pression capillaire dans les pores, compense la force centrifuge exercée.
Une sonde de mesure, est placée dans le godet, sur le côté de l'échantillon.
La sonde
2 peut être de type capacitif par exemple et comporter une tige métallique, isolée par une gaine en céramique. On mesure la capacité entre la tige et le fluide conducteur (la saumure), qui est proportionnelle à la hauteur immergée. Avec ce système, la précision de mesure est de 1,5% du volume poreux. La sonde détecte la position de l'interface entre les deux liquides dans le godet et transmet les signaux de mesure à
un automate de commande et d'acquisition de signaux de mesure regroupant des moyens hydrauliques de circulation des liquides et un appareil d'acquisition Les différents points de la courbe de pression capillaire (Pc) sont obtenus à
partir de la mesure de la saturation moyenne de l'échantillon à l'équilibre pour des io vitesses de rotation données. Pour la détermination des courbes de perméabilités relatives (Kr), on s'intéresse à la variation de la saturation au cours du temps entre deux paliers de vitesse (transitoire). Aucune de ces courbes n'est obtenue directement ; la première doit être calculée et la deuxième nécessite un simulateur numérique d'écoulement comportant (Fig. 11) un module de simulation M qui reproduit les effets de l'expérience grâce à un modèle numérique, et un module inversion I permettant la détermination des courbes de pression capillaire et perméabilité relatives par ajustement itératif d'un modèle numérique initial aux résultats expérimentaux.
Parmi les simulateurs d'écoulement. connus, on peut citer par exemple le simulateur d'écoulement ACEX basé sur les équations généralisées de Darcy, que l'on trouve décrit dans :
Chardaire-Rivière, C. ; 1992, Improving the Centrifuge Technique by Measuring Local Saturations ; SPE 24882.
Bien que la technique de centrifugation soit largement utilisée dans l'industrie pétrolière, du fait de sa rapidité relative si on la compare aux autres méthodes existantes, certaines difficultés de mise en oeuvre demeurent tant en ce qui concerne la détermination de la courbe de pression capillaire que celle des courbes de perméabilité relative.
un automate de commande et d'acquisition de signaux de mesure regroupant des moyens hydrauliques de circulation des liquides et un appareil d'acquisition Les différents points de la courbe de pression capillaire (Pc) sont obtenus à
partir de la mesure de la saturation moyenne de l'échantillon à l'équilibre pour des io vitesses de rotation données. Pour la détermination des courbes de perméabilités relatives (Kr), on s'intéresse à la variation de la saturation au cours du temps entre deux paliers de vitesse (transitoire). Aucune de ces courbes n'est obtenue directement ; la première doit être calculée et la deuxième nécessite un simulateur numérique d'écoulement comportant (Fig. 11) un module de simulation M qui reproduit les effets de l'expérience grâce à un modèle numérique, et un module inversion I permettant la détermination des courbes de pression capillaire et perméabilité relatives par ajustement itératif d'un modèle numérique initial aux résultats expérimentaux.
Parmi les simulateurs d'écoulement. connus, on peut citer par exemple le simulateur d'écoulement ACEX basé sur les équations généralisées de Darcy, que l'on trouve décrit dans :
Chardaire-Rivière, C. ; 1992, Improving the Centrifuge Technique by Measuring Local Saturations ; SPE 24882.
Bien que la technique de centrifugation soit largement utilisée dans l'industrie pétrolière, du fait de sa rapidité relative si on la compare aux autres méthodes existantes, certaines difficultés de mise en oeuvre demeurent tant en ce qui concerne la détermination de la courbe de pression capillaire que celle des courbes de perméabilité relative.
3 Pour un palier de vitesse donné, l'état stationnaire est difficile à
déterminer puisque l'on observe une stabilisation très lente de la production. La recherche de cet équilibre conduit à réaliser des manipulations de longue durée, qui s'étalent sur plusieurs semaines, dont la mise en oeuvre présente certaines difficultés pratiques surtout si l'on considère la vitesse de rotation souvent élevée atteinte par la centrifugeuse. La durée requise est en outre allongée si l'on cherche à
augmenter le volume poreux de l'échantillon (en augmentant sa longueur) pour améliorer la précision de la mesure de saturation.
L'inversion, c'est à dire l'ajustement d'une modélisation numérique aux résultats expérimentaux, est difficile à réaliser car l'information relative aux perméabilités relatives est concentrée sur la partie transitoire de la courbe de production, c'est à dire dans une gamme étroite de valeurs expérimentales.
L'accélération appliquée à la centrifugeuse exerce une influence prépondérante sur la forme de cette partie transitoire, notamment sur sa durée. La détermination, en parallèle, de la courbe de pression capillaire conduit à appliquer plusieurs paliers de vitesse constante plutôt qu'une longue accélération, ce qui réduit les courbes transitoires.
La présente invention vise une méthode pour déterminer les paramètres physiques relatifs à une circulation de fluides à travers un échantillon géologique poreux, comprenant:
a) placer dans un récipient l'échantillon initialement saturé avec un premier fluide;
b) effectuer une centrifugation de l'échantillon en présence d'un second fluide à
des vitesses de rotation variant progressivement par paliers;
c) effectuer pour chaque vitesse de rotation une mesure d'une quantité des premiers et seconds fluides progressivement déplacés par la centrifugation et en
déterminer puisque l'on observe une stabilisation très lente de la production. La recherche de cet équilibre conduit à réaliser des manipulations de longue durée, qui s'étalent sur plusieurs semaines, dont la mise en oeuvre présente certaines difficultés pratiques surtout si l'on considère la vitesse de rotation souvent élevée atteinte par la centrifugeuse. La durée requise est en outre allongée si l'on cherche à
augmenter le volume poreux de l'échantillon (en augmentant sa longueur) pour améliorer la précision de la mesure de saturation.
L'inversion, c'est à dire l'ajustement d'une modélisation numérique aux résultats expérimentaux, est difficile à réaliser car l'information relative aux perméabilités relatives est concentrée sur la partie transitoire de la courbe de production, c'est à dire dans une gamme étroite de valeurs expérimentales.
L'accélération appliquée à la centrifugeuse exerce une influence prépondérante sur la forme de cette partie transitoire, notamment sur sa durée. La détermination, en parallèle, de la courbe de pression capillaire conduit à appliquer plusieurs paliers de vitesse constante plutôt qu'une longue accélération, ce qui réduit les courbes transitoires.
La présente invention vise une méthode pour déterminer les paramètres physiques relatifs à une circulation de fluides à travers un échantillon géologique poreux, comprenant:
a) placer dans un récipient l'échantillon initialement saturé avec un premier fluide;
b) effectuer une centrifugation de l'échantillon en présence d'un second fluide à
des vitesses de rotation variant progressivement par paliers;
c) effectuer pour chaque vitesse de rotation une mesure d'une quantité des premiers et seconds fluides progressivement déplacés par la centrifugation et en
4 déduire une courbe de saturation par rapport au temps pour chaque vitesse de rotation;
d) établir une courbe de saturation par modélisation en sélectionnant des paramètres d'ajustement pour chaque palier de vitesse de rotation, qui correspond au mieux avec une portion mesurée de la courbe de saturation déjà obtenue, et en déduire une valeur de saturation asymptotique; et e) pour chaque palier de vitesse de rotation, prolonger la portion mesurée de la courbe de saturation déjà obtenue en fonction du temps jusqu'à ce que la saturation asymptotique soit stabilisée de façon à ainsi obtenir des valeurs de saturation asymptotique stabilisées pour tous les paliers de vitesses de rotation et en déterminer les paramètres physiques.
De préférence, la méthode de modélisation selon l'invention permet d'obtenir avec un gain de temps considérable les valeurs de paramètres physiques d'un échantillon géologique initialement saturé avec un premier fluide placé dans un récipient et soumis à centrifugation en présence d'un deuxième fluide. Elle comprend une mesure de la quantité des fluides déplacée par la centrifugation de la quantité des fluides déplacée par la centrifugation et l'établissement de courbes de saturation, à des intervalles de temps fixés et elle est caractérisée en ce qu'elle comporte, pour chaque palier de vitesse de rotation, la modélisation d'une courbe de saturation s'ajustant au mieux à
chaque pas de temps avec la portion de courbe déjà obtenue correspondant aux différentes mesures déjà effectuées, la valeur de la saturation à l'équilibre étant déduite de caractéristiques de la courbe modélisée.
De préférence, pour chaque palier de vitesse, on modélise, une courbe de saturation de type exponentiel par un ajustement de paramètres, l'un de ces paramètres correspondant à la valeur de la saturation à l'équilibre. La courbe de 4a saturation est modélisée par exemple au moyen par exemple d'une somme de fonctions exponentielles dépendant desdits paramètres d'ajustement.
De préférence, on choisit les vitesses de centrifugation permettant d'obtenir une distribution régulière de points sur chaque courbe de saturation.
De préférence, suivant un mode de mise en oeuvre, la méthode comporte en outre - l'utilisation d'un simulateur d'écoulement pour modéliser une relation entre. les valeurs des perméabilités relatives respectivement aux deux fluides d'un échantillon de dimensions et de perméabilité définies et les valeurs de saturation respectives à ces mêmes fluides en fonction de différents coefficients et établir une base de données reliant des paramètres d'ajustement aux dits coefficients, pour plusieurs valeurs de saturation à au moins un des fluides ; et - la détermination des perméabilités relatives de l'échantillon soumis à
centrifugation par référence à des valeurs de correspondance extraites de la base de données.
D'autres caractéristiques et avantages de la méthode selon l'invention, apparaîtront à la lecture de la description ci-après d'un exemple non limitatif de réalisation, en se référant aux dessins annexés où :
- la Fig. 1 montre un dispositif de centrifugation permettant une mise en oeuvre de la méthode ;
- la Fig.2 montre un exemple d'ajustement d'une courbe de modélisation sur une courbe de production expérimentale ;
- la Fig.3 montre à titre de comparaison différentes courbes de modélisation obtenues à des instants différents d'un même palier de vitesse ;
- la Fig.4a, 4b, 4c montrent les évolutions comparées respectivement d'une courbe de production, d'une courbe de variation de la constante de temps T2 intervenant
d) établir une courbe de saturation par modélisation en sélectionnant des paramètres d'ajustement pour chaque palier de vitesse de rotation, qui correspond au mieux avec une portion mesurée de la courbe de saturation déjà obtenue, et en déduire une valeur de saturation asymptotique; et e) pour chaque palier de vitesse de rotation, prolonger la portion mesurée de la courbe de saturation déjà obtenue en fonction du temps jusqu'à ce que la saturation asymptotique soit stabilisée de façon à ainsi obtenir des valeurs de saturation asymptotique stabilisées pour tous les paliers de vitesses de rotation et en déterminer les paramètres physiques.
De préférence, la méthode de modélisation selon l'invention permet d'obtenir avec un gain de temps considérable les valeurs de paramètres physiques d'un échantillon géologique initialement saturé avec un premier fluide placé dans un récipient et soumis à centrifugation en présence d'un deuxième fluide. Elle comprend une mesure de la quantité des fluides déplacée par la centrifugation de la quantité des fluides déplacée par la centrifugation et l'établissement de courbes de saturation, à des intervalles de temps fixés et elle est caractérisée en ce qu'elle comporte, pour chaque palier de vitesse de rotation, la modélisation d'une courbe de saturation s'ajustant au mieux à
chaque pas de temps avec la portion de courbe déjà obtenue correspondant aux différentes mesures déjà effectuées, la valeur de la saturation à l'équilibre étant déduite de caractéristiques de la courbe modélisée.
De préférence, pour chaque palier de vitesse, on modélise, une courbe de saturation de type exponentiel par un ajustement de paramètres, l'un de ces paramètres correspondant à la valeur de la saturation à l'équilibre. La courbe de 4a saturation est modélisée par exemple au moyen par exemple d'une somme de fonctions exponentielles dépendant desdits paramètres d'ajustement.
De préférence, on choisit les vitesses de centrifugation permettant d'obtenir une distribution régulière de points sur chaque courbe de saturation.
De préférence, suivant un mode de mise en oeuvre, la méthode comporte en outre - l'utilisation d'un simulateur d'écoulement pour modéliser une relation entre. les valeurs des perméabilités relatives respectivement aux deux fluides d'un échantillon de dimensions et de perméabilité définies et les valeurs de saturation respectives à ces mêmes fluides en fonction de différents coefficients et établir une base de données reliant des paramètres d'ajustement aux dits coefficients, pour plusieurs valeurs de saturation à au moins un des fluides ; et - la détermination des perméabilités relatives de l'échantillon soumis à
centrifugation par référence à des valeurs de correspondance extraites de la base de données.
D'autres caractéristiques et avantages de la méthode selon l'invention, apparaîtront à la lecture de la description ci-après d'un exemple non limitatif de réalisation, en se référant aux dessins annexés où :
- la Fig. 1 montre un dispositif de centrifugation permettant une mise en oeuvre de la méthode ;
- la Fig.2 montre un exemple d'ajustement d'une courbe de modélisation sur une courbe de production expérimentale ;
- la Fig.3 montre à titre de comparaison différentes courbes de modélisation obtenues à des instants différents d'un même palier de vitesse ;
- la Fig.4a, 4b, 4c montrent les évolutions comparées respectivement d'une courbe de production, d'une courbe de variation de la constante de temps T2 intervenant
5 dans la modélisation de la courbe de simulation, et la courbe de variation de la valeur finale S, ,,.,que l'on peut prédire par la simulation ;
- la Fig.5 est un tableau montrant un exemple d'évolutions comparées de plusieurs paramètres au cours d'un même palier de vitesse ;
- les Fig.6a, 6b, 6c, 6d montrent respectivement les évolutions au cours du temps de la valeur de saturation Sfm,,d'un échantillon, de coefficients Tl,T2 intervenant comme paramètres dans la courbe simulée et la courbe de production correspondante ;
- la Fig.7 est un tableau montrant un deuxième exemple d'évolutions comparées de plusieurs paramètres au cours d'un même palier de vitesse ;
- la Fig.8a montre schématiquement une courbe de variation de saturation Sw et les Fig.8b à 8c montrent respectivement des abaques que l'on a formé, reliant différents coefficients, abaques qui correspondent à deux zones de variation de la saturation ;
- la Fig.9 montre les évolutions comparées de courbes de saturation obtenues suivant la méthode classique (Dl i, D21) et suivant la méthode rapide selon l'invention (D12, D22) ;
- la Fig.10 est un tableau comparatif de valeurs numériques pour trois paliers de vitesse, selon que l'on opère par la méthode classique ou la méthode accélérée selon l'invention ; et - la Fig. 11 montre schématiquement l'organisation d'un simulateur d'écoulement.
- la Fig.5 est un tableau montrant un exemple d'évolutions comparées de plusieurs paramètres au cours d'un même palier de vitesse ;
- les Fig.6a, 6b, 6c, 6d montrent respectivement les évolutions au cours du temps de la valeur de saturation Sfm,,d'un échantillon, de coefficients Tl,T2 intervenant comme paramètres dans la courbe simulée et la courbe de production correspondante ;
- la Fig.7 est un tableau montrant un deuxième exemple d'évolutions comparées de plusieurs paramètres au cours d'un même palier de vitesse ;
- la Fig.8a montre schématiquement une courbe de variation de saturation Sw et les Fig.8b à 8c montrent respectivement des abaques que l'on a formé, reliant différents coefficients, abaques qui correspondent à deux zones de variation de la saturation ;
- la Fig.9 montre les évolutions comparées de courbes de saturation obtenues suivant la méthode classique (Dl i, D21) et suivant la méthode rapide selon l'invention (D12, D22) ;
- la Fig.10 est un tableau comparatif de valeurs numériques pour trois paliers de vitesse, selon que l'on opère par la méthode classique ou la méthode accélérée selon l'invention ; et - la Fig. 11 montre schématiquement l'organisation d'un simulateur d'écoulement.
6 Description détaillée Pour la mise en oeuvre de la méthode, on utilise un dispositif de centrifugation tel que décrit ci-après.
Il comporte (Fig.11) un récipient ou godet allongé 13 pour un échantillon Une sonde détectrice de niveau de type capacitif (non représentée) est disposée dans une chambre latérale tubulaire parallèle à la direction d'allongement de l'échantillon et du godet. Toute variation dans le niveau du fluide conducteur dans le godet, se traduit par une variation de la capacité de la sonde. La centrifugeuse Il comporte une cuve 9, un moteur électrique 10 dont l'axe entraîne en rotation un moyeu 11.
io Deux bras 12 identiques sont montés en opposition sur le moyeu 11. Deux récipients ou godets 13, 14 sont montés pivotants aux extrémités des deux bras de façon à s'aligner spontanément avec la direction de la force centrifuge appliquée, et ils s'équilibrent l'un l'autre en rotation. Dans le godet 13 est placé un échantillon de roche poreuse à mesurer initialement saturé en un fluide électriquement conducteur A, tel que de la saumure par exemple 1 Le godet 13 est rempli initialement d'un autre fluide B électriquement isolant tel que de l'huile.
Deux canalisations 15, 16 fixées le long d'un des bras 9, relient le godet 13 à
un connecteur électro-hydraulique tournant 17 porté par le moyeu 11. La première canalisation 15 sert à injecter de l'huile dans le godet 13. La deuxième, 16, sert à
collecter la saumure qui est drainée hors de l'échantillon par la rotation. Un câble 18 relie en outre la sonde de mesure au connecteur tournant 17. Par l'intermédiaire du connecteur 17, les conducteurs du câble 18 sont connectés électriquement à un autre câble. Ces deux canalisations et cet autre câble forment un élément de liaison entre le godet 13 et un ensemble extérieur E adapté à collecter les signaux émanant de la sonde de mesure et à traiter les données acquises, et aussi à commander les transferts de fluide et à actionner le moteur d'entraînement 10.
La méthode comporte au préalable la sélection d'un certain nombre de paliers de vitesse de rotation choisis pour obtenir des points bien répartis sur la courbe de saturation. L'opérateur peut les choisir en se référant à des valeurs
Il comporte (Fig.11) un récipient ou godet allongé 13 pour un échantillon Une sonde détectrice de niveau de type capacitif (non représentée) est disposée dans une chambre latérale tubulaire parallèle à la direction d'allongement de l'échantillon et du godet. Toute variation dans le niveau du fluide conducteur dans le godet, se traduit par une variation de la capacité de la sonde. La centrifugeuse Il comporte une cuve 9, un moteur électrique 10 dont l'axe entraîne en rotation un moyeu 11.
io Deux bras 12 identiques sont montés en opposition sur le moyeu 11. Deux récipients ou godets 13, 14 sont montés pivotants aux extrémités des deux bras de façon à s'aligner spontanément avec la direction de la force centrifuge appliquée, et ils s'équilibrent l'un l'autre en rotation. Dans le godet 13 est placé un échantillon de roche poreuse à mesurer initialement saturé en un fluide électriquement conducteur A, tel que de la saumure par exemple 1 Le godet 13 est rempli initialement d'un autre fluide B électriquement isolant tel que de l'huile.
Deux canalisations 15, 16 fixées le long d'un des bras 9, relient le godet 13 à
un connecteur électro-hydraulique tournant 17 porté par le moyeu 11. La première canalisation 15 sert à injecter de l'huile dans le godet 13. La deuxième, 16, sert à
collecter la saumure qui est drainée hors de l'échantillon par la rotation. Un câble 18 relie en outre la sonde de mesure au connecteur tournant 17. Par l'intermédiaire du connecteur 17, les conducteurs du câble 18 sont connectés électriquement à un autre câble. Ces deux canalisations et cet autre câble forment un élément de liaison entre le godet 13 et un ensemble extérieur E adapté à collecter les signaux émanant de la sonde de mesure et à traiter les données acquises, et aussi à commander les transferts de fluide et à actionner le moteur d'entraînement 10.
La méthode comporte au préalable la sélection d'un certain nombre de paliers de vitesse de rotation choisis pour obtenir des points bien répartis sur la courbe de saturation. L'opérateur peut les choisir en se référant à des valeurs
7 obtenues dans des essais antérieurs pour des échantillons de même type. Il peut aussi à défaut, réaliser rapidement un essai sur échantillon par une méthode de type Purcell bien connue des spécialistes de façon à établir une courbe de pression capillaire correspondante et selon sa forme, estimer les paliers de vitesse qui conviennent.
L'accès rapide recherché aux valeurs de saturation obtenu par la méthode est basé sur une modélisation analytique des paliers de production de fluide par l'échantillon soumis à centrifugation. On utilise à cet effet une équation de modélisation de type exponentiel.
On a trouvé que la plupart des courbes de production pouvaient être modélisées par l'équation :
1 S(t) = Ség 1- P - e-"' + (1- P) e-`'TZ + Sini OU Séq est la saturation à l'équilibre, S;,,i , la saturation initiale du palier de vitesse précédent, t, le temps, P, un facteur de pondération, et T1, T2 sont des constantes de temps.
La constante de temps la plus faible Ti décrit la partie transitoire C1 de la courbe de saturation (Fig.2), c'est à dire les effets dynamiques. La constante de temps la plus T2 grande rend compte de la zone de stabilisation C2, c'est-à-dire des effets capillaires.
Pour estimer la saturation à l'équilibre, bien avant le stade où ses variations deviennent effectivement inférieures à la précision de la mesure, on recherche les coefficients Tl, T2 qui, appliqués dans l'équation ci-dessus, génère une courbe de saturation s'ajustant à la partie déjà connue de la courbe de production.
A chaque nouveau palier de vitesse de rotation, l'automate procède à une acquisition de mesures et on obtient un nouveau couple de données (temps, production). Le programme ajuste alors la modélisation mathématique bi-exponentielle à la courbe de points dont il dispose déjà (Fig.3). On obtient ainsi
L'accès rapide recherché aux valeurs de saturation obtenu par la méthode est basé sur une modélisation analytique des paliers de production de fluide par l'échantillon soumis à centrifugation. On utilise à cet effet une équation de modélisation de type exponentiel.
On a trouvé que la plupart des courbes de production pouvaient être modélisées par l'équation :
1 S(t) = Ség 1- P - e-"' + (1- P) e-`'TZ + Sini OU Séq est la saturation à l'équilibre, S;,,i , la saturation initiale du palier de vitesse précédent, t, le temps, P, un facteur de pondération, et T1, T2 sont des constantes de temps.
La constante de temps la plus faible Ti décrit la partie transitoire C1 de la courbe de saturation (Fig.2), c'est à dire les effets dynamiques. La constante de temps la plus T2 grande rend compte de la zone de stabilisation C2, c'est-à-dire des effets capillaires.
Pour estimer la saturation à l'équilibre, bien avant le stade où ses variations deviennent effectivement inférieures à la précision de la mesure, on recherche les coefficients Tl, T2 qui, appliqués dans l'équation ci-dessus, génère une courbe de saturation s'ajustant à la partie déjà connue de la courbe de production.
A chaque nouveau palier de vitesse de rotation, l'automate procède à une acquisition de mesures et on obtient un nouveau couple de données (temps, production). Le programme ajuste alors la modélisation mathématique bi-exponentielle à la courbe de points dont il dispose déjà (Fig.3). On obtient ainsi
8 l'évolution de la constante de temps T2 qui détermine la stabilisation et de la saturation finale Sfnai au cours du temps. Dans ce programme, on fixe un écart maximal entre deux valeurs consécutives de T2 et de Sfna1. correspondant par exemple à la précision pratique de l'appareil de mesure Lorsque les écarts calculés sont inférieurs à ces valeurs, on peut considérer que T2 et Sfnal sont suffisamment stabilisés et donc que la courbe de production est entièrement déterminée et on passe alors au palier de vitesse suivant.
Leur évolution est représentée sur les Fig 4a à 4c. On remarque une très nette stabilisation de la constante de temps T2 et de la saturation finale. Cette stabilisation 1o intervient dès que l'écoulement atteint le plateau final, c'est-à-dire dès que les forces capillaires deviennent prépondérantes. La constante T1 se stabilise beaucoup plus vite puisqu'elle ne concerne que la première partie de la courbe.
Le tableau de la Fig.5 montre l'évolution au cours du temps de quatre jeux de paramètres P, T1, T2, Sfinal, calculés à quatre instants différents de l'expérience.
Si la précision de la mesure de la saturation S par la sonde capacitive dans le godet contenant l'échantillon est de l'ordre de 1,5% du volume poreux par exemple, soit 1,501o en saturation, on peut dans un tel cas considérer que l'on connaît la valeur finale de la saturation au bout de 3 heures d'expérience. Cela implique une erreur de 2,6% sur T2 et une erreur de 0,9% sur Sonal, inférieure à la précision de la sonde. Si l'on attendait la stabilisation effective de la saturation où les écarts entre mesures successives deviennent inférieurs à la précision de la sonde, le palier de vitesse devrait durer une dizaine d'heures. On divise donc la durée de manipulation par trois.
A partir d'une courbe expérimentale, on retrouve pour les paramètres T2 et Sfnai, le pic puis la stabilisation observés sur les courbes simulées. On obtient les évolutions montrées aux Fig.6a à 6d.
Le tableau de la Fig.7 montre des résultats analogues pour autre cas expérimental.
Leur évolution est représentée sur les Fig 4a à 4c. On remarque une très nette stabilisation de la constante de temps T2 et de la saturation finale. Cette stabilisation 1o intervient dès que l'écoulement atteint le plateau final, c'est-à-dire dès que les forces capillaires deviennent prépondérantes. La constante T1 se stabilise beaucoup plus vite puisqu'elle ne concerne que la première partie de la courbe.
Le tableau de la Fig.5 montre l'évolution au cours du temps de quatre jeux de paramètres P, T1, T2, Sfinal, calculés à quatre instants différents de l'expérience.
Si la précision de la mesure de la saturation S par la sonde capacitive dans le godet contenant l'échantillon est de l'ordre de 1,5% du volume poreux par exemple, soit 1,501o en saturation, on peut dans un tel cas considérer que l'on connaît la valeur finale de la saturation au bout de 3 heures d'expérience. Cela implique une erreur de 2,6% sur T2 et une erreur de 0,9% sur Sonal, inférieure à la précision de la sonde. Si l'on attendait la stabilisation effective de la saturation où les écarts entre mesures successives deviennent inférieurs à la précision de la sonde, le palier de vitesse devrait durer une dizaine d'heures. On divise donc la durée de manipulation par trois.
A partir d'une courbe expérimentale, on retrouve pour les paramètres T2 et Sfnai, le pic puis la stabilisation observés sur les courbes simulées. On obtient les évolutions montrées aux Fig.6a à 6d.
Le tableau de la Fig.7 montre des résultats analogues pour autre cas expérimental.
9 Au bout de 3,5 heures, on connaît la valeur finale de la saturation à 0,66%
près ce qui est inférieur à la précision de la mesure. Si on avait attendu la stabilisation, l'expérience aurait duré plus de 7h. On peut donc réduire de moitié la durée de l'essai d'un facteur de l'ordre de 2 ou 3. On peut donc estimer, compte-tenu de ces deux exemples, qu'à partir du moment où la valeur de Tl est stable à 5%
près, et la valeur de Sf;n,É à 1% près, les paramètres de la modélisation mathématique sont suffisamment déterminés pour pouvoir stopper l'expérience et passer au palier de vitesse suivant.
Détermination des courbes de perméabilité relative Pour déterminer les courbes de perméabilité relative relatives à l'huile et à
l'eau respectivement Kr,nr et Kr,r à partir des saturations respectives Snr et Sr obtenues par l'expérience, on utilise l'équation de Corey bien connue des spécialistes :
Kr,r = S*ni et Kr,nr = S-n2 (2) Ko r Ko nr r,r r,nr où nl, n2 désignent des coefficients.
A cet effet, on constitue une base de données que l'on peut présenter sous forme d'abaques par exemple, permettant de relier les coefficients la constante T2 avec les exposants nl et n2, pour différentes valeurs de ces exposants. On considère un échantillon de longueur déterminée (6 cm par exemple) avec une perméabilité
donnée (100 mD par exemple) et l'on fait varier les coefficients nl, n2 chacun entre 1 et 5.
Au moyen du simulateur, on détermine les courbes S(t) de production pour différents paliers de vitesse. Comme indiqué précédemment, en appliquant la relation 1, on détermine les valeurs de P, Ti, T2 permettant d'ajuster au mieux la courbe modélisée à la variation de courbe de saturation. Pour tous les triplets de valeurs (T2, nl, n2), on forme des tableaux de valeurs ou de façon plus pratique, des abaques comme celles montrées aux Fig. 8a, 8b, correspondant à des paliers différentes de chaque courbe de saturation (Fig.8c).
Pour un palier de vitesse VI correspondant à une forte saturation Sw à l'eau, on constate effectivement que l'exposant n2 sur l'huile a une influence 5 prédominante. Plus il est grand, c'est-à-dire plus la perméabilité relative à l'huile est incurvée, plus la constante de temps T2 est grande. On remarque toutefois que pour ce palier l'exposant à l'eau joue aussi un rôle, même s'il est moins marqué.
Pour le palier de vitesse V3 correspondant à une saturation à l'eau Sw faible, la prédominance de l'exposant nl sur l'eau est très nette. L'exposant sur l'huile n2 io n'a plus aucune influence.
Pour les paliers intermédiaires, on constate que les deux exposants jouent un rôle non négligeable. A ce stade, l'écoulement est sensible à la faible perméabilité
de chacun d'une perméabilité relative faible à chacun des fluides.
Les tableaux ou abaques formés peuvent être utilisés pour des échantillons de longueur et de perméabilité différentes, moyennant l'application de certains facteurs correctifs.
Utiliser des échantillons courts permet de réduire la durée d'un palier de production. Cependant, l'utilisation d'un échantillon plus long permet d'étaler les parties transitoires sur les courbes et donc d'optimiser la détermination des perméabilités relatives. On constate que pour des échantillons d'une même roche mais de longueurs différentes, on peut dégager une relation entre le temps de stabilisation et la longueur de type polynomial du second degré. Pour passer d'un échantillon de longueur L1 (celui ayant servi à l'établissement des abaques) à
un échantillon de longueur L2 différente, on peut appliquer la relation : T2 (L2) _ (L2) 2 T2 (L1) L1 On constate une évolution quasi-linéaire des constantes de temps T1 et T2 en fonction, de la perméabilité absolue de l'échantillon. On applique donc un facteur correctif : TZ (Kz) = (Ki ) T2 (K1) K1 Validation de la méthode On a centrifugé deux échantillons, le premier mouillable à l'eau, le second plutôt mouillable à l'huile, une première fois suivant la méthode classique, la deuxième suivant la méthode décrite. On a appliqué 3 paliers de vitesse successifs :
470, 600 et 1050 t/min, déterminés à la suite d'essais par simulation et l'on a obtenu les courbes de production montrées aux Fig.9a, 9b.
En ce qui concerne les saturations d'équilibres, on obtient les résultats suivant regroupées dans le tableau de la Fig. 10.
On remarque que la précision de la saturation finale est largement tributaire de la durée du palier.
Ainsi, le drainage rapide pour les paliers 1 et 3 donne une saturation finale légèrement différente de la stabilisation effective. Cependant, la durée de ces deux paliers est fortement réduite. Le palier 3 donne par contre des résultats tout à fait satisfaisants, mais sa durée a été beaucoup moins réduite.
Il existe donc un compromis entre précision de la mesure et durée. Après différents essais, on suggère de fixer les critères de variation des paramètres T2 et Sfinal à, respectivement, 2% et 0.5 %. Cependant ces valeurs sont à estimer au cas par cas. On remarque en effet, qu'il existe des échantillons pour lesquels, la constante de temps T2 se stabilise très difficilement, quand la saturation finale est déjà bien déterminée.
Par ailleurs, on constate que le fait de ne pas attendre la stabilisation des fluides avant de lancer un nouveau palier ne perturbe pas la forme de la courbe de production. On peut, en effet, toujours y associer une modélisation bi-exponentielle.
La précision de l'ajustement de la modélisation aux courbes expérimentales reste très acceptable (de l'ordre de 10-5).
On a décrit un type de modélisation faisant appel à des fonctions exponentielles. On ne sortirait pas du cadre de l'invention toutefois, en adoptant s d'autres types de fonctions (des splines par exemple) que l'on peut relier à
des paramètres physiques significatifs de l'échantillon.
près ce qui est inférieur à la précision de la mesure. Si on avait attendu la stabilisation, l'expérience aurait duré plus de 7h. On peut donc réduire de moitié la durée de l'essai d'un facteur de l'ordre de 2 ou 3. On peut donc estimer, compte-tenu de ces deux exemples, qu'à partir du moment où la valeur de Tl est stable à 5%
près, et la valeur de Sf;n,É à 1% près, les paramètres de la modélisation mathématique sont suffisamment déterminés pour pouvoir stopper l'expérience et passer au palier de vitesse suivant.
Détermination des courbes de perméabilité relative Pour déterminer les courbes de perméabilité relative relatives à l'huile et à
l'eau respectivement Kr,nr et Kr,r à partir des saturations respectives Snr et Sr obtenues par l'expérience, on utilise l'équation de Corey bien connue des spécialistes :
Kr,r = S*ni et Kr,nr = S-n2 (2) Ko r Ko nr r,r r,nr où nl, n2 désignent des coefficients.
A cet effet, on constitue une base de données que l'on peut présenter sous forme d'abaques par exemple, permettant de relier les coefficients la constante T2 avec les exposants nl et n2, pour différentes valeurs de ces exposants. On considère un échantillon de longueur déterminée (6 cm par exemple) avec une perméabilité
donnée (100 mD par exemple) et l'on fait varier les coefficients nl, n2 chacun entre 1 et 5.
Au moyen du simulateur, on détermine les courbes S(t) de production pour différents paliers de vitesse. Comme indiqué précédemment, en appliquant la relation 1, on détermine les valeurs de P, Ti, T2 permettant d'ajuster au mieux la courbe modélisée à la variation de courbe de saturation. Pour tous les triplets de valeurs (T2, nl, n2), on forme des tableaux de valeurs ou de façon plus pratique, des abaques comme celles montrées aux Fig. 8a, 8b, correspondant à des paliers différentes de chaque courbe de saturation (Fig.8c).
Pour un palier de vitesse VI correspondant à une forte saturation Sw à l'eau, on constate effectivement que l'exposant n2 sur l'huile a une influence 5 prédominante. Plus il est grand, c'est-à-dire plus la perméabilité relative à l'huile est incurvée, plus la constante de temps T2 est grande. On remarque toutefois que pour ce palier l'exposant à l'eau joue aussi un rôle, même s'il est moins marqué.
Pour le palier de vitesse V3 correspondant à une saturation à l'eau Sw faible, la prédominance de l'exposant nl sur l'eau est très nette. L'exposant sur l'huile n2 io n'a plus aucune influence.
Pour les paliers intermédiaires, on constate que les deux exposants jouent un rôle non négligeable. A ce stade, l'écoulement est sensible à la faible perméabilité
de chacun d'une perméabilité relative faible à chacun des fluides.
Les tableaux ou abaques formés peuvent être utilisés pour des échantillons de longueur et de perméabilité différentes, moyennant l'application de certains facteurs correctifs.
Utiliser des échantillons courts permet de réduire la durée d'un palier de production. Cependant, l'utilisation d'un échantillon plus long permet d'étaler les parties transitoires sur les courbes et donc d'optimiser la détermination des perméabilités relatives. On constate que pour des échantillons d'une même roche mais de longueurs différentes, on peut dégager une relation entre le temps de stabilisation et la longueur de type polynomial du second degré. Pour passer d'un échantillon de longueur L1 (celui ayant servi à l'établissement des abaques) à
un échantillon de longueur L2 différente, on peut appliquer la relation : T2 (L2) _ (L2) 2 T2 (L1) L1 On constate une évolution quasi-linéaire des constantes de temps T1 et T2 en fonction, de la perméabilité absolue de l'échantillon. On applique donc un facteur correctif : TZ (Kz) = (Ki ) T2 (K1) K1 Validation de la méthode On a centrifugé deux échantillons, le premier mouillable à l'eau, le second plutôt mouillable à l'huile, une première fois suivant la méthode classique, la deuxième suivant la méthode décrite. On a appliqué 3 paliers de vitesse successifs :
470, 600 et 1050 t/min, déterminés à la suite d'essais par simulation et l'on a obtenu les courbes de production montrées aux Fig.9a, 9b.
En ce qui concerne les saturations d'équilibres, on obtient les résultats suivant regroupées dans le tableau de la Fig. 10.
On remarque que la précision de la saturation finale est largement tributaire de la durée du palier.
Ainsi, le drainage rapide pour les paliers 1 et 3 donne une saturation finale légèrement différente de la stabilisation effective. Cependant, la durée de ces deux paliers est fortement réduite. Le palier 3 donne par contre des résultats tout à fait satisfaisants, mais sa durée a été beaucoup moins réduite.
Il existe donc un compromis entre précision de la mesure et durée. Après différents essais, on suggère de fixer les critères de variation des paramètres T2 et Sfinal à, respectivement, 2% et 0.5 %. Cependant ces valeurs sont à estimer au cas par cas. On remarque en effet, qu'il existe des échantillons pour lesquels, la constante de temps T2 se stabilise très difficilement, quand la saturation finale est déjà bien déterminée.
Par ailleurs, on constate que le fait de ne pas attendre la stabilisation des fluides avant de lancer un nouveau palier ne perturbe pas la forme de la courbe de production. On peut, en effet, toujours y associer une modélisation bi-exponentielle.
La précision de l'ajustement de la modélisation aux courbes expérimentales reste très acceptable (de l'ordre de 10-5).
On a décrit un type de modélisation faisant appel à des fonctions exponentielles. On ne sortirait pas du cadre de l'invention toutefois, en adoptant s d'autres types de fonctions (des splines par exemple) que l'on peut relier à
des paramètres physiques significatifs de l'échantillon.
Claims (6)
1. Une méthode pour déterminer les paramètres physiques relatifs à une circulation de fluides à travers un échantillon géologique poreux, comprenant:
a) placer dans un récipient l'échantillon initialement saturé avec un premier fluide;
b) effectuer une centrifugation de l'échantillon en présence d'un second fluide à
des vitesses de rotation variant progressivement par paliers;
c) effectuer pour chaque vitesse de rotation une mesure d'une quantité des premiers et seconds fluides progressivement déplacés par la centrifugation et en déduire une courbe de saturation par rapport au temps pour chaque vitesse de rotation;
d) établir une courbe de saturation par modélisation en sélectionnant des paramètres d'ajustement pour chaque palier de vitesse de rotation, qui correspond au mieux avec une portion mesurée de la courbe de saturation déjà obtenue, et en déduire une valeur de saturation asymptotique; et e) pour chaque palier de vitesse de rotation, prolonger la portion mesurée de la courbe de saturation déjà obtenue en fonction du temps jusqu'à ce que la saturation asymptotique soit stabilisée de façon à ainsi obtenir des valeurs de saturation asymptotique stabilisées pour tous les paliers de vitesses de rotation et en déterminer les paramètres physiques.
a) placer dans un récipient l'échantillon initialement saturé avec un premier fluide;
b) effectuer une centrifugation de l'échantillon en présence d'un second fluide à
des vitesses de rotation variant progressivement par paliers;
c) effectuer pour chaque vitesse de rotation une mesure d'une quantité des premiers et seconds fluides progressivement déplacés par la centrifugation et en déduire une courbe de saturation par rapport au temps pour chaque vitesse de rotation;
d) établir une courbe de saturation par modélisation en sélectionnant des paramètres d'ajustement pour chaque palier de vitesse de rotation, qui correspond au mieux avec une portion mesurée de la courbe de saturation déjà obtenue, et en déduire une valeur de saturation asymptotique; et e) pour chaque palier de vitesse de rotation, prolonger la portion mesurée de la courbe de saturation déjà obtenue en fonction du temps jusqu'à ce que la saturation asymptotique soit stabilisée de façon à ainsi obtenir des valeurs de saturation asymptotique stabilisées pour tous les paliers de vitesses de rotation et en déterminer les paramètres physiques.
2. Une méthode selon la revendication 1, comprenant l'étape d'effectuer une modélisation de la courbe de saturation au moyen d'une somme de fonctions exponentielles dépendant des paramètres d'ajustement.
3. Une méthode selon la revendication 1, comprenant l'étape d'effectuer une sélection de vitesses de centrifugation permettant d'obtenir une distribution régulière de points sur chaque courbe de saturation pour chaque vitesse de rotation.
4. Une méthode selon la revendication 1, dans laquelle les paramètres comprennent au moins un des paramètres suivants :
S éq qui représente une saturation à l'équilibre, P qui représente un facteur de pondération; et Ti et T2 qui représentent des constances de temps.
S éq qui représente une saturation à l'équilibre, P qui représente un facteur de pondération; et Ti et T2 qui représentent des constances de temps.
5. Une méthode selon une quelconque des revendications 1 à 4, dans laquelle une modélisation de chacune des courbes de saturation comprend les étapes suivantes:
une utilisation d'un simulateur d'écoulement pour modéliser une relation entre des valeurs de perméabilité relatives respectivement aux premier et deuxième fluides d'un échantillon de dimensions et de perméabilité définies, et des valeurs de saturation respectives à ces premier et deuxième fluides en fonction de différents coefficients, et établir une base de données reliant des paramètres d'ajustement auxdits coefficients, pour plusieurs valeurs de saturation relatives à au moins un desdits premier et second fluides; et déterminer des perméabilités relatives de l'échantillon soumis à la centrifugation par référence à des valeurs correspondantes extraites de la base de données.
une utilisation d'un simulateur d'écoulement pour modéliser une relation entre des valeurs de perméabilité relatives respectivement aux premier et deuxième fluides d'un échantillon de dimensions et de perméabilité définies, et des valeurs de saturation respectives à ces premier et deuxième fluides en fonction de différents coefficients, et établir une base de données reliant des paramètres d'ajustement auxdits coefficients, pour plusieurs valeurs de saturation relatives à au moins un desdits premier et second fluides; et déterminer des perméabilités relatives de l'échantillon soumis à la centrifugation par référence à des valeurs correspondantes extraites de la base de données.
6. Une méthode selon la revendication 5, dans laquelle les coefficients sont une paire d'exposants.
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
FR9912266A FR2798734B1 (fr) | 1999-09-21 | 1999-09-21 | Methode optimisee pour determiner des parametres physiques d'un echantillon soumis a centrifugation |
FR9912266 | 1999-09-21 |
Publications (2)
Publication Number | Publication Date |
---|---|
CA2319962A1 CA2319962A1 (fr) | 2001-03-21 |
CA2319962C true CA2319962C (fr) | 2010-11-23 |
Family
ID=9550478
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CA2319962A Expired - Fee Related CA2319962C (fr) | 1999-09-21 | 2000-09-19 | Methode optimisee pour determiner des parametres physiques d`un echantillon soumis a centrifugation |
Country Status (7)
Country | Link |
---|---|
US (1) | US6490531B1 (fr) |
BR (1) | BR0004842A (fr) |
CA (1) | CA2319962C (fr) |
FR (1) | FR2798734B1 (fr) |
GB (1) | GB2354590B (fr) |
NL (1) | NL1016223C2 (fr) |
NO (1) | NO318617B1 (fr) |
Families Citing this family (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20020060778A (ko) * | 1999-12-06 | 2002-07-18 | 홀던 데이비드 | 프로세스에서 적어도 하나의 유체 흐름 변수를 예측하는온라인 장치 |
US7126332B2 (en) * | 2001-07-20 | 2006-10-24 | Baker Hughes Incorporated | Downhole high resolution NMR spectroscopy with polarization enhancement |
US7032661B2 (en) * | 2001-07-20 | 2006-04-25 | Baker Hughes Incorporated | Method and apparatus for combined NMR and formation testing for assessing relative permeability with formation testing and nuclear magnetic resonance testing |
FR2837572B1 (fr) * | 2002-03-20 | 2004-05-28 | Inst Francais Du Petrole | Methode pour modeliser la production d'hydrocarbures par un gisement souterrain soumis a une depletion |
FR2839782B1 (fr) | 2002-05-15 | 2004-06-18 | Inst Francais Du Petrole | Methode d'evalution de la courbe de pression capillaire des roches d'un gisement souterrain a partir de mesures sur des debris de roche |
FR2844355B1 (fr) | 2002-09-11 | 2005-03-25 | Inst Francais Du Petrole | Methode de mesure de la mouillabilite de roches par resonnance magnetique nucleaire |
FR2864244B1 (fr) * | 2003-12-19 | 2006-08-04 | Inst Francais Du Petrole | Methode et dispositif pour mesurer des caracteristiques physiques d'un echantillon solide poreux |
US6971260B2 (en) * | 2004-01-13 | 2005-12-06 | Coretest Systems, Inc. | Overburden rock core sample containment system |
US9020793B2 (en) * | 2005-12-22 | 2015-04-28 | Chevron U.S.A. Inc. | Method, system and program storage device for reservoir simulation utilizing heavy oil solution gas drive |
BE1020054A3 (nl) * | 2011-07-06 | 2013-04-02 | Peiffer Herman | Inrichting voor het analyseren van bodemmonsters. |
BE1021504B1 (nl) * | 2012-10-19 | 2015-12-03 | PEIFFER, Herman Hugo Jan | Inrichting voor het analyseren van bodemmonsters |
US10552553B2 (en) | 2015-08-17 | 2020-02-04 | Saudi Arabian Oil Company | Capillary pressure analysis for petrophysical statistical modeling |
CN106493003B (zh) * | 2017-01-09 | 2018-10-02 | 中国工程物理研究院总体工程研究所 | 超大型土工离心机负载可分离结构 |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4740077A (en) * | 1985-06-21 | 1988-04-26 | Atlantic Richfield Company | Centrifugal measurement of core samples |
IT1255750B (it) * | 1992-08-27 | 1995-11-15 | Agip Spa | Procedimento per ricavare sperimentalmente la curva di pressione capillare in un mezzo poroso |
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. |
US5373727A (en) * | 1993-04-16 | 1994-12-20 | New Mexico Tech Research Foundation | Miniporopermeameter |
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. |
FR2763690B1 (fr) * | 1997-05-23 | 1999-07-02 | Inst Francais Du Petrole | Dispositif perfectionne pour faire des mesures de caracteristiques physiques d'un echantillon poreux |
FR2772477B1 (fr) * | 1997-12-11 | 2000-03-24 | Inst Francais Du Petrole | Procede et dispositif pour mesurer des caracteristiques physiques d'un echantillon poreux en y deplacant des fluides par centrifugation |
-
1999
- 1999-09-21 FR FR9912266A patent/FR2798734B1/fr not_active Expired - Fee Related
-
2000
- 2000-09-19 CA CA2319962A patent/CA2319962C/fr not_active Expired - Fee Related
- 2000-09-19 NL NL1016223A patent/NL1016223C2/nl not_active IP Right Cessation
- 2000-09-19 GB GB0022911A patent/GB2354590B/en not_active Expired - Fee Related
- 2000-09-20 BR BR0004842-9A patent/BR0004842A/pt not_active IP Right Cessation
- 2000-09-20 US US09/666,568 patent/US6490531B1/en not_active Expired - Fee Related
- 2000-09-20 NO NO20004699A patent/NO318617B1/no not_active IP Right Cessation
Also Published As
Publication number | Publication date |
---|---|
FR2798734B1 (fr) | 2001-10-12 |
GB2354590B (en) | 2003-03-12 |
FR2798734A1 (fr) | 2001-03-23 |
GB0022911D0 (en) | 2000-11-01 |
NL1016223C2 (nl) | 2001-03-23 |
BR0004842A (pt) | 2001-06-12 |
NO318617B1 (no) | 2005-04-18 |
NO20004699D0 (no) | 2000-09-20 |
NO20004699L (no) | 2001-03-22 |
GB2354590A (en) | 2001-03-28 |
US6490531B1 (en) | 2002-12-03 |
CA2319962A1 (fr) | 2001-03-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CA2319962C (fr) | Methode optimisee pour determiner des parametres physiques d`un echantillon soumis a centrifugation | |
EP0113285B1 (fr) | Procédé et dispositif pour déterminer les caractéristiques d'écoulement d'un fluide dans un puits à partir de mesures de températures | |
WO2003098196A1 (fr) | Methode d'evaluation de la courbe de pression capillaire des roches d'un gisement souterrain a partir de mesures sur des debris de roche | |
WO2004025317A1 (fr) | Methode de mesure de la mouillabilite de roches par resonance magnetique nucleaire | |
FR2461066A1 (fr) | Procede et appareil d'essai dynamique de pieux | |
EP0816629B1 (fr) | Méthode et système d'estimation en temps reel d'au moins un paramètre lié au déplacement d'un outil de forage | |
FR2815124A1 (fr) | Procede de determination de la saturation en hydrocarbure d'une formation | |
FR2953289A1 (fr) | Procede et dispositif de surveillance de vibrations en torsion d'un arbre rotatif d'une turbomachine. | |
EP0026706B1 (fr) | Procédé et appareil pour la détermination de paramètres de direction d'un puits exploré en continu | |
FR2793032A1 (fr) | Procede et appareil pour determiner la resistivite d'une formation traversee par un puits tube | |
EP0031744A2 (fr) | Dispositif et procédé pour la détermination de caractéristiques d'écoulement d'un fluide | |
CA2080483C (fr) | Methode de surveillance automatique de l'etat vibratoire d'une garniture de forage | |
FR2611804A1 (fr) | Procede de controle des operations de forage d'un puits | |
FR2740170A1 (fr) | Procedes et dispositifs de mesure de la resistivite de la boue dans un puits d'hydrocarbure | |
FR2795521A1 (fr) | Procede et dispositif pour determiner la resistivite d'une formation traversee par un puits tube | |
EP0816630A1 (fr) | Méthode et système d'estimation en temps réel d'au moins un paramètre lié au comportement d'un outil de fond de puits | |
FR2636593A1 (fr) | Pale pour voilure tournante d'aeronef et voilure tournante comportant une telle pale | |
EP0603040B1 (fr) | Procédé et dispositif perfectionnés pour l'étude des propriétés d'un matériau perméable | |
FR2763690A1 (fr) | Dispositif perfectionne pour faire des mesures de caracteristiques physiques d'un echantillon poreux | |
CA2254308A1 (fr) | Procede et dispositif pour mesurer des caracteristiques physiques d'un echantillon poreux en y deplacant des fluides par centrifugation | |
FR2864244A1 (fr) | Methode et dispositif pour mesurer des caracteristiques physiques d'un echantillon solide poreux | |
FR2740169A1 (fr) | Procedes et dispositifs de mesure de caracteristiques d'une formation traversee par un trou de forage | |
FR2844360A1 (fr) | Methode pour determiner l'indice de resistivite en fonction de la saturation en eau, de certaines roches de porosite complexe | |
FR2907836A1 (fr) | Procede pour determiner une repartition de phases fluides et quantifier une retention dans un puits de forage | |
EP0394085A1 (fr) | Procédé et dispositif pour la mesure des écoulements multiphasiques, notamment dans les puits d'hydrocarbures |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
EEER | Examination request | ||
MKLA | Lapsed |
Effective date: 20170919 |