FR2844355A1 - Methode de mesure de la mouillabilite de roches par resonnance magnetique nucleaire - Google Patents
Methode de mesure de la mouillabilite de roches par resonnance magnetique nucleaire Download PDFInfo
- Publication number
- FR2844355A1 FR2844355A1 FR0211283A FR0211283A FR2844355A1 FR 2844355 A1 FR2844355 A1 FR 2844355A1 FR 0211283 A FR0211283 A FR 0211283A FR 0211283 A FR0211283 A FR 0211283A FR 2844355 A1 FR2844355 A1 FR 2844355A1
- Authority
- FR
- France
- Prior art keywords
- water
- oil
- sample
- wettability
- relaxation
- 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.)
- Granted
Links
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 68
- 238000005259 measurement Methods 0.000 title claims abstract description 20
- 239000011435 rock Substances 0.000 title claims abstract description 18
- 229920006395 saturated elastomer Polymers 0.000 title claims abstract description 18
- 238000005481 NMR spectroscopy Methods 0.000 title claims abstract description 13
- 239000012530 fluid Substances 0.000 claims abstract description 41
- 238000000034 method Methods 0.000 claims abstract description 31
- 238000006073 displacement reaction Methods 0.000 claims abstract description 7
- 238000005119 centrifugation Methods 0.000 claims description 6
- 239000003921 oil Substances 0.000 description 52
- 239000000523 sample Substances 0.000 description 34
- 238000009826 distribution Methods 0.000 description 20
- 239000011148 porous material Substances 0.000 description 11
- SNRUBQQJIBEYMU-UHFFFAOYSA-N dodecane Chemical compound CCCCCCCCCCCC SNRUBQQJIBEYMU-UHFFFAOYSA-N 0.000 description 10
- 238000005213 imbibition Methods 0.000 description 5
- 230000003993 interaction Effects 0.000 description 5
- 238000004458 analytical method Methods 0.000 description 4
- 239000012267 brine Substances 0.000 description 4
- 238000002347 injection Methods 0.000 description 4
- 239000007924 injection Substances 0.000 description 4
- 239000007788 liquid Substances 0.000 description 4
- 239000012528 membrane Substances 0.000 description 4
- HPALAKNZSZLMCH-UHFFFAOYSA-M sodium;chloride;hydrate Chemical compound O.[Na+].[Cl-] HPALAKNZSZLMCH-UHFFFAOYSA-M 0.000 description 4
- 238000002474 experimental method Methods 0.000 description 3
- 239000003208 petroleum Substances 0.000 description 3
- 238000012360 testing method Methods 0.000 description 3
- 230000015572 biosynthetic process Effects 0.000 description 2
- 238000009792 diffusion process Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 239000011159 matrix material Substances 0.000 description 2
- 238000013421 nuclear magnetic resonance imaging Methods 0.000 description 2
- LFQSCWFLJHTTHZ-UHFFFAOYSA-N Ethanol Chemical compound CCO LFQSCWFLJHTTHZ-UHFFFAOYSA-N 0.000 description 1
- 230000032683 aging Effects 0.000 description 1
- 230000004888 barrier function Effects 0.000 description 1
- 239000004566 building material Substances 0.000 description 1
- 239000000919 ceramic Substances 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000004140 cleaning Methods 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 238000000151 deposition Methods 0.000 description 1
- 238000002592 echocardiography Methods 0.000 description 1
- 230000005684 electric field Effects 0.000 description 1
- 230000002964 excitative effect Effects 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 239000007789 gas Substances 0.000 description 1
- 238000010438 heat treatment Methods 0.000 description 1
- 229930195733 hydrocarbon Natural products 0.000 description 1
- 150000002430 hydrocarbons Chemical class 0.000 description 1
- 125000004435 hydrogen atom Chemical group [H]* 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 239000002184 metal Substances 0.000 description 1
- 230000010287 polarization Effects 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 238000005086 pumping Methods 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 230000006833 reintegration Effects 0.000 description 1
- 238000000679 relaxometry Methods 0.000 description 1
- 239000005871 repellent Substances 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 239000002904 solvent Substances 0.000 description 1
- 238000010561 standard procedure Methods 0.000 description 1
- 238000011282 treatment Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
-
- 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
- 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
- G01N24/00—Investigating or analyzing materials by the use of nuclear magnetic resonance, electron paramagnetic resonance or other spin effects
- G01N24/08—Investigating or analyzing materials by the use of nuclear magnetic resonance, electron paramagnetic resonance or other spin effects by using nuclear magnetic resonance
- G01N24/082—Measurement of solid, liquid or gas content
Landscapes
- Physics & Mathematics (AREA)
- Chemical & Material Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Immunology (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Pathology (AREA)
- Analytical Chemistry (AREA)
- High Energy & Nuclear Physics (AREA)
- Dispersion Chemistry (AREA)
- Fluid Mechanics (AREA)
- Engineering & Computer Science (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Remote Sensing (AREA)
- Food Science & Technology (AREA)
- Medicinal Chemistry (AREA)
- Geology (AREA)
- Environmental & Geological Engineering (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
- Production Of Liquid Hydrocarbon Mixture For Refining Petroleum (AREA)
Abstract
- Méthode pour la mesure de la mouillabilité de roches par résonance magnétique nucléaire à faible champ.- La méthode comporte essentiellement la détermination de la surface poreuse mouillée par l'eau et de la surface poreuse mouillée par l'huile quand l'échantillon est saturé en eau et en huile, par le biais de mesures de temps de relaxation (T1, T2) de l'échantillon placé dans un appareil à résonance magnétique nucléaire, amené préalablement à différents états de saturation en eau ou en huile, et le calcul de l'indice de mouillabilité par combinaison des valeurs obtenues des dites surfaces.- Applications à l'étude ou l'exploitation de gisements d'hydrocarbures ou en génie civil, etc.
Description
Désignation du domaine technique La présente invention concerne une
méthode de mesure de la mouillabilité d'échantillons de roches par la mesure du temps de relaxation par résonance magnétique
nucléaire à faible champ.
La méthode trouve des applications notamment pour l'analyse de roches prélevées dans une formation souterraine recelant ou susceptible de receler des hydrocarbures La connaissance de différents paramètres et notamment de la mouillabilité des
roches, est aussi utile notamment quand on doit procéder à une récupération assistée d'une formation, par une injection d'un fluide sous pression, et que l'on doit déterminer par des 10 tests préalables le fluide (liquide ou gaz) qui convient le mieux pour déplacer les effluents.
L'invention trouve aussi des applications en génie civil pour faire de l'hydrologie de terrains pour évaluer leur degré de pollution par exemple, ou encore dans le bâtiment pour tester des matériaux de construction afin notamment de décider de traitements hydrofuges
par exemple.
Présentation de l'art antérieur On connaît des méthodes pour déterminer la mouillabilité des roches vis à vis de l'eau et de l'huile qui peuvent y être contenues, comportant la réalisation de cycles de drainage de la roche c'est-à-dire de déplacement des fluides visant à diminuer la saturation en eau, suivi d'une imbibition, en désignant par ce terme un déplacement des fluides 20 permettant d'augmenter la saturation en eau (Sw) de la roche. La pression capillaire Pc en un point se définit comme la différence à l'équilibre entre la pression P(huile) de l'huile et celle P(eau) de l'eau. Ce paramètre n'a de sens que si les deux fluides sont en phase continue dans le milieu poreux. Pour un milieu mouillable à l'eau, seules les valeurs positives ont un sens. Par contre, lorsque le milieu a une mouillabilité mixte, les fluides 25 peuvent rester en phase continue aussi bien pour les pressions capillaires (Pc) positives que négatives. Pour une application de ce type, un cycle complet de mesure de la pression capillaire comporte généralement (Fig. 1):
2 2844355
a) un drainage primaire positif d'un échantillon saturé initialement en eau à 100% (courbe 1); b) une imbibition positive (courbe 2); c) une imbibition négative (courbe 3); d) un drainage négatif (courbe 4); et
e) un drainage secondaire positif (courbe 5).
l existe différents types de dispositif capables de dresser les courbes de la Fig. 1.
Dans un dispositif dit à "plaque poreuse", décrit notamment par le brevet US 4 506 542, l'échantillon de roche poreuse contenant deux fluides en phases continues, est 10 placé dans une cellule allongée terminée à ses deux extrémités opposées par des barrières
capillaires perméables à un premier fluide. Ce premier fluide est injecté sous pression au travers de la première membrane et on mesure la différence de pression entre la pression d'injection et la pression du fluide évacué à l'autre extrémité. Les pressions des deux fluides ainsi que la pression capillaire Pc sont constantes le long de l'échantillon, et la 15 saturation est supposée uniforme.
Il est également connu de procéder par centrifugation à vitesse progressive au moyen de dispositifs de centrifugation tels que ceux décrits par exemple par les brevets FR 2 772 477 (US 6 185 985) ou FR 2 763 690 ou encore dans les brevets EP 603 040
(US 5 463 894) ou FR 2 798 734) du demandeur.
L'échantillon initialement saturé avec un premier liquide (de la saumure par exemple) est placé (Fig.10) dans un récipient ou godet allongé A contenant un deuxième fluide de densité différente de l'huile par exemple) . Le godet A est fixé au bout d'un bras tournant B, et on lui applique une force centrifuge de façon à étudier les déplacements des fluides dans l'échantillon au cours d'au moins deux phases distinctes. Durant une première 25 phase de drainage, on soumet alors l'ensemble à une force centrifuge dirigée suivant la longueur du récipient de manière à exercer sur lui une force d'expulsion qui tend à faire sortir une partie du premier fluide. Le champ de pression créé par la centrifugation s'exprime en fonction de la densité r, du rayon R et de la vitesse angulaire w, par la relation: 1/2 w2**-. p (Rmax2- R2), pour chaque fluide. On impose que la pression des 30 deux fluides à la sortie de l'échantillon soit la même et qu'elle s'annule à la sortie. Dans le
3 2844355
même temps, le deuxième fluide pénètre à l'intérieur de l'échantillon. Les deux fluides se déplacent à l'intérieur de l'échantillon jusqu'à une position d'équilibre o la force due à la pression capillaire dans les pores, compense la force centrifuge exercée. Une sonde de mesure, est placée dans le godet, sur le côté de l'échantillon. La sonde peut être de type 5 capacitif par exemple et comporter une tige métallique, isolée par une gaine en céramique.
On mesure la capacité entre la tige et le fluide conducteur (la saumure), qui est proportionnelle à la hauteur immergée. Avec ce moyen de mesure, la précision de mesure est de 1,5% du volume poreux. La sonde détecte la position de l'interface entre les deux liquides dans le godet et transmet les signaux de mesure à un automate E de commande et 10 d'acquisition de signaux de mesure regroupant des moyens hydrauliques de circulation des liquides et un appareil d'acquisition Dans la phase de ré-imbibition, on diminue la vitesse de façon à étudier la
réintégration du fluide initial dans l'échantillon. Les saturations locales mesurées avec ce type de dispositif, sont calculées par un programme d'inversion à partir de la quantité totale 15 d'eau expulsée hors de l'échantillon.
Suivant une autre méthode dite "dynamique", on place un échantillon dans une cellule allongée terminée à ses deux extrémités par des membranes perméables à l'eau. A une première extrémité, on injecte directement de l'huile sous pression dans l'enceinte. De l'eau y est également injectée mais cette injection est effectuée au travers de la membrane 20 et à une pression inférieure. A l'extrémité opposée, l'huile est évacuée directement alors que l'eau sort au travers de la membrane terminale. En réglant les débits d'injection d'huile et d'eau, on s'arrange pour que la pression capillaire soit la même à l'entrée de l'enceinte qu'à sa sortie, ce qui entraîne une saturation uniforme pouvant être déduite du bilan des fluides. La pression capillaire est obtenue par exemple en mesurant la différence entre les 25 pressions de l'huile et de l'eau à la sortie de l'enceinte. Une telle méthode est décrite notamment par Brown H.W. in "Capillary Pressure Investigations" Petroleum Transaction AIME, vol. 192, 1951. Des exemples de mise en oeuvre sont décrits par exemple dans les
brevets EP 729 022 (US 5 698 772) ou EP-A- 974 839 (US 6 229 312) du demandeur.
On connaît également une méthode dite semi dynamique dans laquelle on confine 30 un échantillon de roche imbibé avec un premier fluide dans une cellule close, on injecte un autre fluide sous pression à une première extrémité de l'enceinte et on balaie l'extrémité opposée par un fluide à basse pression mis en circulation par des moyens de pompage qui
4 2844355
entraîne au dehors le fluide drainé. Le dispositif comporte des moyens de mesure de la pression et la saturation de l'échantillon différents de sa longueur, la quantité de fluide évacué et la résistivité électrique de l'échantillon. Cette méthode est mise en oeuvre par
exemple dans le brevet FR 2 708 742 (US 5 679 885) du demandeur.
Ces courbes de drainage et d'imbibition ayant été établies, il est connu de calculer les indices de mouillabilité WI à partir des aires Ad et Ai délimitées par les courbes de pression capillaire positives et négatives, comme le montre la Fig.1, par la relation WI(USBM) = Lg Ad Ai La méthode selon l'invention La méthode de mesure de la mouillabilité d'un échantillon de roche poreuse en présence d'eau et d'huile, selon l'invention comporte la détermination de la surface poreuse mouillée *par l'en eau et de la surface poreuse mouillée par l'huile quand l'échantillon est saturé en eau et en huile, et le calcul de l'indice de mouillabilité par combinaison des
valeurs obtenues des dites surfaces.
La détermination de la surface poreuse mouillée en eau et de la surface poreuse mouillée par l'huile quand l'échantillon est saturé en eau et en huile, est obtenue par exemple par le biais de mesures de temps de relaxation de l'échantillon placé dans un
appareil à résonance magnétique nucléaire.
On détermine par exemple l'indice de mouillabilité par la relation: IMN = SM -SM ou la relation IN = îocr SMw SMNS + SM, Mv C10 SMo dans lesquelles SM. est la surface poreuse mouillée par l'eau et SM, est la surface poreuse
mouillée par l'huile lorsque le milieu poreux est saturé en eau et huile.
Suivant un mode de mise en oeuvre, on détermine l'indice de mouillabilité par les opérations suivantes: a) on mesure les temps caractéristiques de relaxation de l'échantillon saturé en eau;
2844355
b) on mesure les temps caractéristiques de relaxation de l'huile dans l'échantillon en présence d'eau, dans une zone proche de la saturation (Swir) de l'échantillon; c) on mesure les temps caractéristiques de relaxation de de l'eau dans l'échantillon en présence d'huile, dans une zone proche de la saturation résiduelle (Sor); d) on mesure les temps de relaxation de l'échantillon dans un état o son point de saturation à 100% en huile est atteint; et
e) on combine les mesures des différentes temps de relaxation pour obtenir le dit indice.
Suivant un mode préféré de mise en oeuvre, on détermine les temps de relaxation
des étapes a) à c) après application à l'échantillon d'une centrifugation.
Suivant un mode préféré de mise en oeuvre, on détermine les temps de relaxation de l'étape d) après déplacement forcé des fluides dans l'échantillon placé dans une cellule de confinement. De préférence, on choisit une huile dont le temps de relaxation intrinsèque (TB) est
le plus grand possible et le plus proche de l'eau, du dodécane par exemple.
Les temps caractéristiques de relaxation sont par exemple ceux qui correspondent
soit au maxima des courbes de saturation soit à des valeurs moyennes des dites courbes.
On vérifie à l'expérience que la mesure de IRMN obtenue par la méthode est aussi sensible mais nécessite bien moins de temps à établir et qu'elle est applicable à un grand
nombre d'échantillons.
Présentation succincte des figures Les caractéristiques et avantages de la méthode et du dispositif selon l'invention,
apparaîtront plus clairement à la lecture de la description ci-après d'un exemple non
limitatif de réalisation, en se référant aux dessins annexés o: - la figure 1 montre des exemples de courbes de pression capillaire nécessaires à la 25 détermination de l'indice de mouillabilité classique USBM (United States Bureau of Mines);
6 2844355
- les figure 2a et 2b montrent respectivement la distribution des temps de relaxation T2 pour une roche mouillable à l'eau saturée en eau et en huile, et une représentation de la distribution de phase (matrice en hachuré, eau en gris clair et huile en gris plus foncé); - les figure 3a et 3b montrent respectivement la distribution des temps de relaxation T2 5 pour une roche de mouillabilité intermédiaire saturée en eau et en huile (même partie centrale que pour la figure 2), et une représentation de la distribution de phase (matrice en hachuré, eau en gris clair et huile en gris plus foncé); - la figure 4 montre sur des courbes de pression capillaire, les valeurs de saturation o les mesures de relaxation RMN sont faites - la figure 5 montre un exemple de distribution des temps de relaxation T2 aux valeurs de saturation indiquées à la figure 4 (de bas en haut: 100% eau, Swi, Sor et 100% dodecane); et - la figure 6 montre la relation entre l'indice IRMN et l'indice IUSBM standard; - la figure 7 montre schématiquement une cellule de confinement avec un circuit 15 permanent de chauffage par circulation de fluide; et - la figure 8 montre schématiquement l'agencement d'un dispositif de centrifugation
utilisé pour la mise en oeuvre de la méthode.
Description détaillée de la méthode
On rappelle que la technique d'analyse RNIN consiste essentiellement on le 20 rappelle, à appliquer à un objet à tester un premier champ magnétique statique Bo de polarisation destiné à aligner les noyaux des protons d'hydrogène initialement orientés au hasard, suivant la direction du champ puis à un deuxième champ magnétique impulsionnel oscillant à la fréquence de Larmor perpendiculaire au premier créé par des bobinages excités par un signal de commande pour réaliser une expérience de résonance magnétique 25 nucléaire. Quand ce champ impulsionnel cesse, le retour des noyaux vers leur état initial ou
relaxation, engendre des signaux électromagnétiques (échos) qui sont détectés et analysés.
A partir des caractéristiques d'amplitude de ces signaux, on détermine la présence de
paramètres physiques de l'objet.
7 2844355
On propose de définir un nouvel indice de mouillabilité construit en combinant des valeurs de la surface poreuse SM, mouillée par l'eau et de la surface poreuse SM, mouillée par l'huile lorsque le milieu poreux est saturé en eau et huile. L'indice peut être calculé par exemple par la relation suivante: RMN SMW + SMo (la) ou encore par la relation: IRPMN =logo SM (lb) Les quantitées SM, et SMo sont obtenues par exemple par la mesure du temps de
relaxation dominant dans une expérience de résonance magnétique nucléaire à faible 10 champ.
En effet, dans une telle expérience, la relation fondamentale liant le temps de relaxation longitudinal T1 ou transversal T2 à la surface S et au volume de ce pore est la suivante:
1 1 S (
T,2 TIB,2B V
ou T1B,2B est le temps de relaxation du fluide hors du milieu poreux. Fondamentalement, cette relation provient du fait que les molécules dans le pore subissent des mouvements effets de diffusion et interagissent avec la paroi poreuse pendant le temps de mesure (le temps de relaxation typique est de l'ordre de 100 ms). Les interactions de surface sont rassemblées dans le coefficient Pî,2 appelé relaxivité de surface. Nous avons négligé par 20 ailleurs un terme de gradient qui est important lorsque la mesure est faite à haut champ magnétique (>0. lT). La relation 2 est strictement valable pour un pore saturé d'un seul fluide. En général, les milieux poreux naturels présentent une distribution de tailles de pore et donc du rapport S/V. On observe ainsi en général une distribution des temps de
relaxation mas ceci ne modifie pas la méthode décrite ici.
Lorsque deux fluides sont présents dans un pore au sein du milieu poreux, le même mécanisme physique de diffusion est valable, mais la distribution des deux fluides vis à vis
8 2844355
de la surface poreuse a une importance cruciale. Par exemple, lorsque le milieu est mouillable à l'eau, l'eau se trouve en surface et interagit avec elle, tandis que l'huile se trouve à l'intérieur du pore et n'interagit pas avec la surface. Lorsque l'on mesure la distribution des temps de relaxation pour un tel système, on observe le résultat de la figure 5 2a. L'eau est marquée par des temps de relaxation plus courts qu'à la saturation 100% en eau, car le volume d'eau V a diminué (voir équation 1). L'huile est caractérisée par un temps de relaxation qui est le même que si la mesure était faite hors milieu poreux (TIB,2B) car il n'y a pas d'interaction avec la surface poreuse. fl est par ailleurs connu que les interactions à l'interface eau/huile ne donnent pas lieu à une relaxation différente de celle 10 de l'huile hors milieu poreux, à condition que l'huile considérée soit légère (par une huile raffinée tel le dodecane). Schématiquement, la distribution de phase est représentée à la figure 2b. Lorsque le même milieu poreux a une mouillabilité différente, non franche à l'eau, la distribution des temps de relaxation est modifiée (figure 3a). Essentiellement, le temps caractéristique de l'huile sera plus faible du fait d'une interaction de l'huile à la 15 surface poreuse car l'eau n'est plus en premier voisin de la surface. On représente la distribution possible des deux fluides à la figure 3b. Un tel effet a été obtenu en utilisant les procédures standards permettant de reproduire les propriétés de surface de roches provenant d'un réservoir pétrolier: nettoyage puis vieillissement dans l'huile réservoir à la
température du réservoir à la saturation irréductible en eau.
Ainsi, on généralise la relation (2) à un système diphasique eau-huile. Pour l'eau, nous avons: i i Sm T = p1 V (3) Tl,2w TBw,2Bw W et pour l'huile:
1 1 + SM, (4)
T1T,20 lo,2Bo 10,20 VO On peut donc déterminer les quantités SM, et SM, par la mesure des temps de relaxation Tl ou T2 en milieux poreux, des temps de relaxation TIB ou T2B des fluides hors milieux poreux, des volumes liquides VO ou V, présents dans le milieux poreux. Les constantes d'interaction p peuvent être déterminées séparement, mais nous verrons que
seul le rapport de ces constantes est nécessaire.
9 2844355
En introduisant les relations 3 et 4 dans la relation de base 1, nous obtenons: Sw( - T) - CpSo( - T) r _ w bw o b IRMN T Tb (5) sw( 1- -) + CpSo(- - -)
T T T T
I w lbw so lbo o Sw=Vw/Vp est la saturation en eau, So=Vo/Vp est la saturation en huile Vp le volume poreux Cp= P- est le rapport des relaxivités pour les surface mouillées par l'eau (w) et
par l'huile (o). Pour alléger la notation, les indices 1 et 2 ont été supprimés.
Il reste donc à déterminer les saturations les plus appropriées pour mesurer les quantités SM, et SM0. Pour des raisons liées au calcul des distributions de temps de 10 relaxation, on choisira la saturation irréductible en eau (Swi, figure 4) pour déterminer la surface mouillée par l'huile en présence d'eau (SMw), et la saturation résiduelle en huile (Sor, figure 4) pour déterminer la surface mouillée par l'eau en présence d'huile. On choisit également une huile dont le temps de relaxation intrinsèque (TB) est le plus grand possible et le plus proche de l'eau. En effet, si TB est trop faible, la méthode sera limitée à 15 des milieux poreux dont le rapport S/V est grand (des tailles de pore petites), ou à des milieux dont la relaxivité de surface est grande. Le dodecane par exemple est une huile raffinée qui est adaptée à la mesure proposée car son temps de relaxation intrinsèque (TBO) est de 1 s, proche du temps de relaxation de l'eau (TBW, environ 2.7 s). Il faut en général exclure les huiles pétrolières brutes car elles ont des temps de relaxation trop faibles et 20 possèdent également une distribution de temps de relaxation qui gêne considérablement
l'analyse. Les détails pratiques du calcul sont détaillés ci-dessous.
A partir de la distribution des temps de relaxation aux différentes saturations (voir l'exemple de la figure 5), on peut facilement déterminer le temps de relaxation correspondant au pic dominant (Tw à Sor, To à Swi), et ce sont ces valeurs que l'on utilise 25 dans la relation 5. Il est connu que la détermination du temps de relaxation dominant est robuste et dépend peu du procédé de calcul des distributions. On note que le pic dominant correspond à celui de l'huile à Swi, et à l'eau à Sor. La présence du deuxième fluide à ces
2844355
différentes saturations est importante du point de vue physique mais perturbe peu le calcul.
On peut montrer que quelque soit le temps de relaxation du deuxième fluide présent en faible quantité (l'eau à Swi, l'huile à Sor), le temps de relaxation du fluide dominant fluctue peu. Cet aspect est important pour la robustesse du calcul. La figure 5 présente la 5 distribution des temps de relaxation transversal T2. On peut également utiliser le temps de relaxation longitudinal Ti mais cette mesure est moins favorable en général car T1>T2 mais
T1B=T2B pour un fluide donné.
Pour déterminer le rapport de relaxivité de surface Cp, on utilise les temps de relaxation dominants Twloo et T.Ioo respectivement aux deux saturations Sw=100% et 10 So=100%. La formule utilisée est la suivante: C p= 1/TW100 -1TBW pO 1/TT100 -1/TBO (6) La figure 5 illustre les deux mesures effectuées. Ainsi, au total, on doit donc mesurer les distributions de temps de relaxation à quatre saturations: 100% saumure, 100% huile et les saturations extrêmes Swi et Sor. Pour cela on pourra utiliser la technique 15 de la centrifugation pour atteindre, en partant de Sw=100%, Swi puis Sor. Un palier de rotation unique à vitesse maximale permet de les atteindre rapidement. Pour la mesure à So=100%, on utilise une succession de déplacements miscibles en injectant un alcool puis
l'huile raffinée.
Mise en oeuvre Pour la mise en oeuvre de la méthode, on utilise par exemple un appareil de mesure RMN qui comporte classiquement (Fig.7) des aimants 1 avec un entrefer dans lequel sont placés des bobinages 2. Les bobinages sont connectés à un boîtier électronique pour engendrer les signaux excitateurs (signaux dans la gamme des radiofréquences par exemple) créant un champ électrique oscillant et acquérir la réponse de l'échantillon S aux 25 signaux excitateurs. Un dispositif de ce type est mis en oeuvre par exemple dans la
demande de brevet FR-A-2 823 308.
il 2844355 Exemple pratique de mise en oeuvre L'indice de mouillabilité est obtenu en réalisant par exemple la succession d'étapes suivantes avec un échantillon de réservoir dont la mouillabilité est intermédiaire: 1 On sature l'échantillon à 100% avec de la saumure de gisement et on le place dans 5 l'appareil de mesure RMN (celui de la Fig.7-8 par exemple) pour mesurer les distributions des temps de relaxation TI, T2 et en déduire le paramètre TwlOO de l'équation 6; 2 On le place ensuite dans le dispositif de centrifugation (tel que celui schématisé à la Fig.10) et on le centrifuge à la vitesse maximale dans un godet rempli d'huile jusqu'à atteindre la saturation irréductible en eau Swir; après quoi on mesure les distributions des 10 temps de relaxation T1, T2 dans l'appareil de RMN, de façon à en déduire To de l'équation ; 3 Placé de nouveau dans le dispositif de centrifugation, on centrifuge l'échantillon à vitesse maximale jusqu'à atteindre le point de saturation résiduelle en huile Sor; puis on le transfère dans l'appareil RMN de nouveau pour mesurer les temps de relaxation TI, T2 et en 15 déduire T. 4 Ceci fait, on place l'échantillon dans une cellule de confinement telle que celles
décrites dans les brevets précités et on injecte un solvant miscible à l'eau et à l'huile puis de l'huile (par exemple du dodécane) jusqu'à atteindre la saturation 100% en huile. Ce point étant atteint, on mesure de nouveau les temps de relaxation TI, T2 dans l'appareil RMN 20 d'o l'on déduit la valeur de Toloo.nécessaire à la détermination de Cp.
On connaît ainsi tous les paramètres pour calculer le coefficient de mouillabilité IRMN. Comparaison avec l 'indice IuSBM
La comparaison entre l'indice de mouillabilité lUSBM et le nouvel indice IRMN 25 montre une bonne corrélation entre ces deux quantités (figure 6) pour des roches réservoir.
On peut donc avantageusement remplacer la mesure de IUSBM par une mesure de IRMN aussi sensible mais nécessitant bien moins de temps à établir et applicable à un grand nombre d'échantillons.
12 2844355
On a -décrit un exemple de mise en oeuvre o l'on utilise une méthode de relaxométrie de type RMIN pour mesurer la surface et le volume des pores. Cette méthode n'a cependant rien de limitatif. On peut utiliser toute autre méthode d'analyse et notamment la cryomicroscopie.
13 2844355
Claims (5)
1) Méthode de mesure de la mouillabilité d'un échantillon de roche poreuse en
présence d'eau et d'huile, caractérisée en ce qu'elle comporte la détermination de la surface poreuse mouillée en eau et de la surface poreuse mouillée par l'huile quand l'échantillon est 5 saturé en eau et en huile, et le calcul de l'indice de mouillabilité par combinaison des valeurs obtenues des dites surfaces.
2 Méthode selon la revendication 1, caractérisée en ce que la détermination de la surface poreuse mouillée en eau et de la surface poreuse mouillée par l'huile quand l'échantillon est saturé en eau et en huile, est obtenue par le biais de mesures de temps de 10 relaxation (TI, T2) de l'échantillon placé dans un appareil à résonance magnétique nucléaire 3) Méthode selon la revendication 1 ou 2, caractérisée en ce que l'indice de mouillabilité est obtenue par la relation:
= SM" -SM"
RMVSMW + SM0
o SM, est la surface poreuse mouillée par l'eau et SM, est la surface poreuse mouillée 15 par l'huile lorsque le milieu poreux est saturé en eau et huile.
4) Méthode selon la revendication 1 ou 2, caractérisée en ce que l'indice de mouillabilité est obtenu par la relation: RMN = 1og1O SM0
o SM, est la surface poreuse mouillée par l'eau et SM, est la surface poreuse mouillée 20 par l'huile lorsque le milieu poreux est saturé en eau et huile.
) Méthode selon l'une des revendications précédentes, caractérisée en ce que l'on
détermine l'indice de mouillabilité par les opérations suivantes: a) on mesure les temps caractéristiques de relaxation de l'échantillon saturé en eau; b) on mesure les temps caractéristiques de relaxation de l'huile dans l'échantillon en 25 présence d'eau, dans une zone proche de la saturation (Swir) de l'échantillon;
- 14 2844355
c) on mesure les temps caractéristiques de relaxation de de l'eau dans l'échantillon en présence d'huile, dans une zone proche de la saturation résiduelle (Sor) d) on mesure les temps de relaxation de l'échantillon dans un état o son point de saturation à 100% en huile est atteint; et e) on combine les mesures des différentes temps de relaxation pour obtenir le dit indice. 6) Méthode selon la revendication 5, caractérisée en ce que l'on mesure les temps
de relaxation des étapes a) à c) après application à l'échantillon d'une centrifugation.
7) Méthode selon la revendication 5, caractérisée en ce que l'on mesure les temps
de relaxation de l'étapes d) après déplacement forcé des fluides dans l'échantillon placé 10 dans une cellule de confinement.
8) Méthode selon l'une des revendications précédentes, caractérisée en ce que l'on
choisit une huile dont le temps de relaxation intrinsèque (TB) est le plus grand possible et le
plus proche de l'eau.
9) Méthode selon la revendication 5, caractérisée en ce que les temps 15 caractéristiques de relaxation sont ceux qui correspondent soit au maxima des courbes de
saturation soit à des valeurs moyennes des dites courbes.
Priority Applications (5)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
FR0211283A FR2844355B1 (fr) | 2002-09-11 | 2002-09-11 | Methode de mesure de la mouillabilite de roches par resonnance magnetique nucleaire |
PCT/FR2003/002544 WO2004025317A1 (fr) | 2002-09-11 | 2003-08-18 | Methode de mesure de la mouillabilite de roches par resonance magnetique nucleaire |
AU2003274261A AU2003274261A1 (en) | 2002-09-11 | 2003-08-18 | Method of measuring rock wettability by means of nuclear magnetic resonance |
US10/527,088 US7397240B2 (en) | 2002-09-11 | 2003-08-18 | Method of measuring rock wettability by means of nuclear magnetic resonance |
EP03758245A EP1540363A1 (fr) | 2002-09-11 | 2003-08-18 | Methode de mesure de la mouillabilite de roches par resonance magnetique nucleaire |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
FR0211283A FR2844355B1 (fr) | 2002-09-11 | 2002-09-11 | Methode de mesure de la mouillabilite de roches par resonnance magnetique nucleaire |
Publications (2)
Publication Number | Publication Date |
---|---|
FR2844355A1 true FR2844355A1 (fr) | 2004-03-12 |
FR2844355B1 FR2844355B1 (fr) | 2005-03-25 |
Family
ID=31726034
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
FR0211283A Expired - Fee Related FR2844355B1 (fr) | 2002-09-11 | 2002-09-11 | Methode de mesure de la mouillabilite de roches par resonnance magnetique nucleaire |
Country Status (5)
Country | Link |
---|---|
US (1) | US7397240B2 (fr) |
EP (1) | EP1540363A1 (fr) |
AU (1) | AU2003274261A1 (fr) |
FR (1) | FR2844355B1 (fr) |
WO (1) | WO2004025317A1 (fr) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2034308A1 (fr) * | 2007-09-07 | 2009-03-11 | Ifp | Méthode de mesure rapide de la saturation et de la résistivité d'un milieu poreux |
CN112858364A (zh) * | 2020-07-27 | 2021-05-28 | 苏州泰纽测试服务有限公司 | 一种利用核磁共振测量岩心物性的方法 |
CN113933333A (zh) * | 2021-09-09 | 2022-01-14 | 西安石油大学 | 致密储层co2驱沥青质沉积对岩石润湿性变化特征评价方法 |
Families Citing this family (48)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2008088335A1 (fr) * | 2007-01-18 | 2008-07-24 | Halliburton Energy Services, Inc. | Inversion simultanee de temps de relaxation |
US7567079B2 (en) | 2007-06-08 | 2009-07-28 | University Of New Brunswick | Methods suitable for measuring capillary pressure and relative permeability curves of porous rocks |
EP2144053A1 (fr) * | 2008-07-08 | 2010-01-13 | Services Pétroliers Schlumberger | Détermination de paramètres de formations géologiques par mesures RMN de dispersion du temps de relaxation longitudinal T1 |
US8278922B2 (en) * | 2009-03-23 | 2012-10-02 | Schlumberger Technology Corporation | Continuous wettability logging based on NMR measurements |
US8076933B2 (en) | 2009-04-29 | 2011-12-13 | Schlumberger Technology Corporation | Method for determining wettability of an oil reservoir |
US8201439B2 (en) * | 2009-11-02 | 2012-06-19 | Schlumberger Technology Corporation | Material wettability characterization and chemical additive evaluation |
WO2011073608A1 (fr) * | 2009-12-16 | 2011-06-23 | Bp Exploraton Operating Company Limited | Procédé permettant de mesurer la mouillabilité d'une roche |
EP2341372A1 (fr) * | 2009-12-16 | 2011-07-06 | BP Exploration Operating Company Limited | Procédé pour mésurer la mouillabilité de roche |
US20130057277A1 (en) * | 2010-01-22 | 2013-03-07 | Lukasz Zielinski | Method for Determining Rock Formation Fluid Interaction Properties Using Nuclear Magnetic Resonance Well Logging Measurements |
US8768628B2 (en) * | 2010-10-20 | 2014-07-01 | Shawket Ghedan | Rise in core wettability characterization method |
US9678185B2 (en) | 2013-03-15 | 2017-06-13 | Pepsico, Inc. | Method and apparatus for measuring physico-chemical properties using a nuclear magnetic resonance spectrometer |
CN104316554B (zh) * | 2014-10-27 | 2016-07-06 | 中国石油天然气股份有限公司 | 储层孔隙水可动性测试方法 |
CN105891248A (zh) * | 2015-04-17 | 2016-08-24 | 北京大学 | 一种高温高压岩石物性及渗流机理核磁共振在线测试装置 |
US9599581B2 (en) | 2015-04-22 | 2017-03-21 | Saudi Arabian Oil Company | Nuclear magnetic resonance gas isotherm technique to evaluate reservoir rock wettability |
US10533933B2 (en) | 2015-04-22 | 2020-01-14 | Saudi Arabian Oil Company | Nuclear magnetic resonance gas isotherm technique to evaluate reservoir rock wettability |
US10718701B2 (en) * | 2015-05-12 | 2020-07-21 | Schlumberger Technology Corporation | NMR based reservoir wettability measurements |
WO2016210151A1 (fr) * | 2015-06-24 | 2016-12-29 | Conocophillips Company | Déterminations de la mouillabilité d'une roche |
CN107315024A (zh) * | 2016-04-26 | 2017-11-03 | 中国石油化工股份有限公司 | 一种识别致密砂岩储层油水层的方法 |
CN106525888B (zh) * | 2016-09-26 | 2018-10-16 | 中国石油天然气股份有限公司 | 一种测试致密油藏润湿性的方法及装置 |
CN106814018B (zh) * | 2017-01-18 | 2023-04-28 | 西南石油大学 | 一种致密岩石气相相对渗透率测量装置及方法 |
AU2018212812A1 (en) | 2017-01-26 | 2019-08-15 | Dassault Systemes Simulia Corp. | Multi-phase flow visualizations based on fluid occupation time |
US10488352B2 (en) * | 2017-01-27 | 2019-11-26 | Saudi Arabian Oil Company | High spatial resolution nuclear magnetic resonance logging |
FR3067811B1 (fr) * | 2017-06-19 | 2019-06-21 | IFP Energies Nouvelles | Procede de mesure du coefficient de diffusion de l'eau au sein d'un milieu poreux par une methode de resonance magnetique nucleaire |
FR3069643B1 (fr) * | 2017-07-27 | 2021-10-01 | Ifp Energies Now | Procede de determination par centrifugation de la saturation residuelle d'un premier fluide dans un milieu poreux suite a l'injection d'un deuxieme fluide |
CN107561109A (zh) * | 2017-08-03 | 2018-01-09 | 中国石油大学(华东) | 基于核磁共振横向弛豫谱技术的鲜花生含水率检测方法 |
CN107831186B (zh) * | 2017-09-28 | 2019-09-20 | 同济大学 | 岩芯毛细压力与饱和度曲线、油气藏参数的一维1h核磁共振成像测量表征方法 |
CN107505350A (zh) * | 2017-10-20 | 2017-12-22 | 河南工业大学 | 一种基于低场核磁的葡萄籽油掺伪快速鉴别方法 |
CN108020488B (zh) * | 2017-11-21 | 2019-11-19 | 中国石油大学(北京) | 核磁共振评价致密砂岩油气储层润湿性的方法与装置 |
US10451571B2 (en) | 2017-12-11 | 2019-10-22 | Saudi Arabian Oil Company | Measuring the wettability of porous media based on the temperature sensitivity of nuclear magnetic resonance relaxation time |
US11714040B2 (en) * | 2018-01-10 | 2023-08-01 | Dassault Systemes Simulia Corp. | Determining fluid flow characteristics of porous mediums |
CN108414560B (zh) * | 2018-03-06 | 2020-07-07 | 中国石油大学(华东) | 一种核磁-驱替联用装置评价致密油充注过程的方法 |
US10495589B2 (en) | 2018-04-17 | 2019-12-03 | Saudi Arabian Oil Company | Determining permeability of porous media based on nuclear magnetic resonance measurement |
CN108827853B (zh) * | 2018-05-22 | 2024-02-06 | 西南石油大学 | 基于核磁共振的致密储层岩电测量装置及测量方法 |
CN109100384B (zh) * | 2018-08-21 | 2021-07-30 | 中国石油天然气股份有限公司 | 致密油藏水驱油微观机理信息确定方法、装置及系统 |
CN109030292B (zh) * | 2018-09-26 | 2019-12-17 | 西南石油大学 | 一种致密岩石润湿性确定的新方法 |
CN110646331B (zh) | 2019-09-10 | 2020-08-11 | 中国石油天然气股份有限公司 | 高含粘土岩心的有效孔隙度确定方法及装置 |
CN110687612B (zh) * | 2019-09-17 | 2020-09-08 | 中国石油天然气股份有限公司 | 吸附油和游离油含量连续表征的页岩油分析方法及装置 |
CN110672495A (zh) * | 2019-10-30 | 2020-01-10 | 哈尔滨工业大学 | 一种基于低场磁共振技术的水泥基材料水分渗透率的预测方法 |
CN111022009B (zh) * | 2019-12-27 | 2024-07-12 | 延长油田股份有限公司志丹采油厂 | 一种脉冲作用下渗吸实验装置及实验方法 |
CN111537543B (zh) * | 2020-06-03 | 2021-04-16 | 中国矿业大学 | 低场核磁共振测定页岩黏土与脆性矿物相对含量的方法 |
US11847391B2 (en) | 2020-06-29 | 2023-12-19 | Dassault Systemes Simulia Corp. | Computer system for simulating physical processes using surface algorithm |
US11187766B1 (en) | 2020-07-09 | 2021-11-30 | Saudi Arabian Oil Company | Methods and systems for determining fluid content in formation samples using low field nuclear magnetic resonance |
US11907625B2 (en) | 2020-12-29 | 2024-02-20 | Dassault Systemes Americas Corp. | Computer simulation of multi-phase and multi-component fluid flows including physics of under-resolved porous structures |
CN112924356A (zh) * | 2021-01-28 | 2021-06-08 | 中国石油大学(北京) | 一种储层动态渗流特征获取方法及装置 |
US11821861B2 (en) * | 2021-04-22 | 2023-11-21 | Baker Hughes Oilfield Operations Llc | Wettability estimation using magnetic resonance |
US11493461B1 (en) * | 2021-06-28 | 2022-11-08 | Baker Hughes Oilfield Operations Llc | Wettability estimation using T2 distributions of water in wetting and non-wetting phases |
US11614417B2 (en) | 2021-07-06 | 2023-03-28 | Saudi Arabian Oil Company | Determining saturation in low resistivity pay zones |
CN113834840B (zh) * | 2021-09-24 | 2024-05-14 | 西安工程大学 | 一种测试岩心渗吸效率的方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4543821A (en) * | 1983-12-14 | 1985-10-01 | Texaco Inc. | Method and apparatus for measuring relative permeability and water saturation of a core |
EP0974839A1 (fr) * | 1998-07-24 | 2000-01-26 | Institut Francais Du Petrole | Méthode de mesure rapide de l'indice de résistivité d'échantillons solides tels que des roches |
Family Cites Families (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4506542A (en) | 1983-04-22 | 1985-03-26 | Chandler Engineering Company | Apparatus and procedure for relative permeability measurements |
US5162733A (en) * | 1991-02-26 | 1992-11-10 | Phillips Petroleum Company | Method for determining relative wettability |
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. |
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. |
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 |
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 |
FR2798734B1 (fr) | 1999-09-21 | 2001-10-12 | Inst Francais Du Petrole | Methode optimisee pour determiner des parametres physiques d'un echantillon soumis a centrifugation |
US6522136B1 (en) * | 1999-12-10 | 2003-02-18 | Schlumberger Technology Corporation | Well logging technique and apparatus for determining pore characteristics of earth formations using magnetic resonance |
AU2001280678A1 (en) * | 2000-07-21 | 2002-02-05 | Schlumberger Holdings Limited | Nuclear magnetic resonance measurements and methods of analyzing nuclear magnetic resonance data |
FR2823308B1 (fr) | 2001-04-05 | 2003-12-19 | Inst Francais Du Petrole | Cellule de confinement thermostatee pour echantillon destine a des mesures par rmn et une methode pour sa mise en oeuvre |
US6765380B2 (en) * | 2002-05-23 | 2004-07-20 | Schlumberger Technology Corporation | Determining wettability of an oil reservoir using borehole NMR measurements |
-
2002
- 2002-09-11 FR FR0211283A patent/FR2844355B1/fr not_active Expired - Fee Related
-
2003
- 2003-08-18 AU AU2003274261A patent/AU2003274261A1/en not_active Abandoned
- 2003-08-18 US US10/527,088 patent/US7397240B2/en not_active Expired - Fee Related
- 2003-08-18 WO PCT/FR2003/002544 patent/WO2004025317A1/fr not_active Application Discontinuation
- 2003-08-18 EP EP03758245A patent/EP1540363A1/fr not_active Withdrawn
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4543821A (en) * | 1983-12-14 | 1985-10-01 | Texaco Inc. | Method and apparatus for measuring relative permeability and water saturation of a core |
EP0974839A1 (fr) * | 1998-07-24 | 2000-01-26 | Institut Francais Du Petrole | Méthode de mesure rapide de l'indice de résistivité d'échantillons solides tels que des roches |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2034308A1 (fr) * | 2007-09-07 | 2009-03-11 | Ifp | Méthode de mesure rapide de la saturation et de la résistivité d'un milieu poreux |
FR2920876A1 (fr) * | 2007-09-07 | 2009-03-13 | Inst Francais Du Petrole | Methode de mesure rapide de la saturation et de la resistivite d'un milieu poreux. |
US7642774B2 (en) | 2007-09-07 | 2010-01-05 | Institut Francais Du Petrole | Method for fast measurement of the saturation and the resistivity of a porous medium |
CN112858364A (zh) * | 2020-07-27 | 2021-05-28 | 苏州泰纽测试服务有限公司 | 一种利用核磁共振测量岩心物性的方法 |
CN113933333A (zh) * | 2021-09-09 | 2022-01-14 | 西安石油大学 | 致密储层co2驱沥青质沉积对岩石润湿性变化特征评价方法 |
CN113933333B (zh) * | 2021-09-09 | 2024-06-07 | 西安石油大学 | 致密储层co2驱沥青质沉积对岩石润湿性变化特征评价方法 |
Also Published As
Publication number | Publication date |
---|---|
WO2004025317A1 (fr) | 2004-03-25 |
US20060132131A1 (en) | 2006-06-22 |
EP1540363A1 (fr) | 2005-06-15 |
US7397240B2 (en) | 2008-07-08 |
AU2003274261A1 (en) | 2004-04-30 |
FR2844355B1 (fr) | 2005-03-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
FR2844355A1 (fr) | Methode de mesure de la mouillabilite de roches par resonnance magnetique nucleaire | |
EP1167948B1 (fr) | Méthode pour évaluer des paramètres physiques d'un gisement souterrain à partir de débris de roche qui y sont prélevés | |
CA2474712C (fr) | Methode et dispositif pour evaluer des parametres physiques d'un gisement souterrain a partir de debris de roche qui y sont preleves | |
FR2853071A1 (fr) | Methode et dispositif pour evaluer des parametres physiques d'un gisement souterrain a partir de debris de roche qui y sont preleves | |
CA2482751C (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 | |
FR2869694A1 (fr) | Methode de determination des proprietes des fluides de formation | |
NO335578B1 (no) | Bestemmelse av fuktningsevne i et oljereservoar ved bruk av NMR-borehullsmålinger | |
FR2871235A1 (fr) | Methode et appareil pour utilisation de mesures rmn de gradient de champ pulse pour determiner les proprietes des fluides dans un outil de diagraphie de puits par echantillonnage des fluides | |
NO312266B1 (no) | Fremgangsmåte for bestemmelse av permeabilitet i poröst medium, spesielt et fluidummettet medium | |
EP2034308B1 (fr) | Méthode de mesure rapide de la saturation et de la résistivité d'un milieu poreux | |
FR2724460A1 (fr) | Dispositif de mesure petrophysique et mehode de mise en oeuvre | |
EP1548455A1 (fr) | Détermination de la perméabilité d'un milieu souterrain à partir de mesures par RMN de la perméabilité de fragments de roche issu du milieu | |
EP1398630B1 (fr) | Méthode pour déterminer l'indice de résistivité en fonction de la saturation en eau, de certaines roches de porosité complexe | |
FR3079931A1 (fr) | Procede de mesure de l'imbibition spontanee et/ou de la diffusion d'une phase aqueuse dans un milieu poreux par une methode de resonance magnetique nucleaire. | |
FR2772477A1 (fr) | Procede et dispositif pour mesurer des caracteristiques physiques d'un echantillon poreux en y deplacant des fluides par centrifugation | |
Dick et al. | Spatially resolved wettability measurements using NMR wettability index | |
EP0603040A1 (fr) | Procédé et dispositif perfectionnés pour l'étude des propriétés d'un matériau perméable | |
FR2864244A1 (fr) | Methode et dispositif pour mesurer des caracteristiques physiques d'un echantillon solide poreux | |
FR2846365A1 (fr) | Procede et dispositif de localisation d'une interface par rapport a un trou fore | |
FR3069643A1 (fr) | Procede de determination par centrifugation de la saturation residuelle d'un premier fluide dans un milieu poreux suite a l'injection d'un deuxieme fluide | |
CA2472087A1 (fr) | Methode de detection et de suivi par resonance magnetique nucleaire de la cinetique de floculation des fractions lourdes d'un fluide complexe | |
FR2836227A1 (fr) | Methode pour evaluer des parametres physiques d'un gisement souterrain a partir de debris de roche qui y sont preleves | |
Johannesen et al. | Evaluation of wettability distributions in experimentally aged core | |
Siqveland et al. | Aging time control by NMR relaxation | |
FR3106207A1 (fr) | procédé de détermination de la tension interfaciale entre deux fluides au sein d’un milieu poreux |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
CD | Change of name or company name | ||
ST | Notification of lapse |
Effective date: 20150529 |