Method for determining the properties of porous materials

Abstract

A sample of porous material is placed in a calorimeter cell and a pressure in the cell is increased starting from a pressure value of a first step by filling the cell with a wetting fluid. Measurements are taken of a heat flow to the cell and a fluid volume at each step. Then, the pressure in the cell is decreased to the pressure value of a first step with continued measurements of the heat flow to the cell. Increase and following decrease of the fluid pressure in the cell are repeated at least once. Then a temperature in the cell is decreased below a wetting fluid crystallization point. Once the fluid has been fully crystallized in sample pores, the temperature in the cell is increased above a wetting fluid melting point. Wetting limiting angle of the pores filled with fluid, and pore sizes are determined based on the results of heat flow measurements with due consideration of heat effect of fluid compression.

Claims

1. A method for determining properties of porous materials comprising: placing a sample of a porous material into a calorimetric cell, increasing step by step a hydrostatic pressure in the calorimetric cell with the sample starting from a pressure value of a first step by filling the cell with a wetting fluid and keeping until stabilization of a heat flow at each step, measuring the heat flow to the calorimetric cell and a volume of fluid at each step, decreasing the hydrostatic pressure of the wetting fluid in the calorimetric cell down to the pressure value of the first step with continuous measurement of the heat flow to the cell, at least once repeating step by step increase and following decrease of the fluid pressure in the calorimetric cell down to the pressure value of the first step, decreasing a temperature of the calorimetric cell at constant hydrostatic pressure below a crystallization point of the wetting fluid with continuous measurement of the heat flow and the fluid volume, after completion of fluid crystallization in pores of the sample increasing the temperature of the calorimetric cell above a melting point of the wetting fluid with continuous measurement of the heat flow and the fluid volume, calculating contact angles and pore size distribution of the pores of the sample filled with the wetting fluid based on results of heat flow measurement and taking into account a heat effect from fluid compressibility.

2. The method of claim 1 wherein the wetting fluid is preliminary filled into the calorimetric cell without the sample, the hydrostatic pressure in the calorimetric cell without the sample is increased step by step starting from a pressure value of a first step and then decreased to the value of the first step with measurement of heat flow into the cell.

3. The method of claim 1 wherein the pressure value of the first step is 0.1 MPa (1 bar).

4. The method of claim 1 wherein the wetting fluid is oil.

5. The method of claim 1 wherein the wetting fluid is water.

6. The method of claim 1 wherein the wetting fluid is a water based solution of salt.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

(1) The invention is illustrated by drawings, where

(2) FIG. 1 shows a heat flow as a time function while brine is penetrating into sandstone;

(3) FIG. 2 shows ice/water phase transition in a sample with known pore sizes.

DETAILED DESCRIPTION

(4) The proposed method for determining properties of porous materials is based on heat flow measurements by micro-calorimetric methods in the process of fluid (water, solution or oil) penetration into a porous rock structure (sandstone, limestone, etc.), and subsequent crystallization/melting in a pore volume.

(5) The novelty of the method is based on determining a value of temperature shift in a heat effect associated with phase transition (for instance, ice/water) in sample pores, and a heat effect from wetting associated with fluid (for instance, water) injection.

(6) In case of the simplest systema drop of fluid on a homogeneous, flat and inert surfacea limiting angle is determined from the Young equation:
.sub.sv=.sub.sl+.sub.lv cos ,(1)
where a contact angle, .sub.sva surface energy at the solid/vapor interface, .sub.sla surface energy at the solid/liquid interface, .sub.lva surface energy at the liquid/vapor interface.

(7) It is known that a temperature shift associated with substance melting/solidification phase transition depends on size of particles. If a given phase transition takes place in a porous body where the size of particles is limited by size of pores, the temperature shift will characterize pores of a rock sample. So, fluid melting temperature in bulk T.sub.0 and in a porous structure T.sub.m can be measured by a micro-calorimeter. The temperature shift associated with phase transition and dependent upon the pore size may be calculated by the Gibson-Thompson formula. Notably, a non-freezing liquid layer should be considered in calculations (0.5-2 nm). That correction is especially important in analyzing samples with nanopores:

(8) $\begin{array}{cc}{T}_m=T_0-T_m=\frac{2T_0{}_{\mathrm{sl}}{}_l}{\left(R-t\right).Math.H}=\frac{2T_0{}_{\mathrm{sl}}{}_l}{r_{\mathrm{eff}}.Math.H}& \left(2\right)\end{array}$

(9) .sub.sla surface energy at a solid/liquid (ice/water) interface, .sub.la liquid molar volume, Rpore size, ta non-freezing liquid layer, r.sub.effan effective pore radius, Ha phase transition enthalpy.

(10) With fluid penetrating into a porous medium, the process of surface/liquid interface starts under controlled pressure. Variation in free energy values (F of area unit) can be described by the following equations (3, 4):

(11) $\begin{array}{cc}U=F-T\frac{F}{T}& \left(3\right)\\ F={}_{\mathrm{sl}}-{}_{\mathrm{sv}}& \left(4\right)\end{array}$
where Uchange in the system inner energy in response to fluid penetration.

(12) Using the Young equation (1) and energy variation, the limiting angle can be expressed as:

(13) $\begin{array}{cc}U=-{}_{\mathrm{lv}}\cos +T\frac{{}_{\mathrm{lv}}\cos }{T}=\left(-{}_{\mathrm{lv}}+T\frac{{}_{\mathrm{lv}}}{T}\right)\cos +T\frac{{}_{\mathrm{lv}}\cos }{T}& \left(5\right)\end{array}$

(14) Equation (5) can be approximated to (6) when the limiting angle is not controlled by temperature.

(15) $\begin{array}{cc}\cos =\frac{-U}{\left({}_{\mathrm{lv}}-T\frac{{}_{\mathrm{lv}}}{T}\right)}& \left(6\right)\end{array}$

(16) The advantage of calorimetric methods is an option for conducting experiments in which the system initial and finite states are well-defined, which is not always the case with other methods, like conventional methods for determining contact angle.

(17) The proposed method is used as follows. A cell containing a sample is placed in a calorimeter, for instance BT2.15 (SETARAM, France, http://www.setaram.ru/BT-2.15-ru.htm). Then the cell is filled with a wetting fluid (oil, water, salt solution) at a constant flow rate until a hydrostatic pressure exceeds atmospheric pressure, for instance, 0.1 MPa (1 bar). A time needed for heat flow stabilization is a function of a flow rate and temperature of the wetting fluid filling the calorimeter cell (usually, several hours). Then, the hydrostatic pressure is increased step-by-step. A size of a pressure step depends on specifics of the sample pore volume, an amount of macro- and micro-pores. In our case, the pressure change steps were as follows: 0.2 MPa (2 bar), 0.4 MPa (4 bar), 0.8 MPa (8 bar), 1.0 MPa (10 bar), and 2.0 MPa (20 bar). The sample is held at each pressure step for several hours for stabilization of the heat flow. While increasing the pressure the heat flow to the cell is measured, i.e., the number of mW/sec (FIG. 1), and a volume of the fluid injected into the sample-containing cell.

(18) Then, the fluid pressure is reduced to 0.1 MPa (1 bar), and the system is kept until the heat flow is stabilized. Then, at least once step by step increase and following decrease of the fluid pressure in the calorimetric cell down to value of the first step are repeated.

(19) At the next stage, fluid pressure is stabilized at a certain value, for instance, 0.1 MPa (1 bar). Meanwhile, the sample remains filled with fluid. Then, the temperature of the cell with the sample is decreased at constant pressure. An important parameter for determining pore size distribution in the sample is a phase transition temperature at which fluid would crystallize. For instance, crystallization temperature for water solutions in bulk is somewhere around 0 C., depending on salts dissolved in them. If water solution is injected in the experiment at temperature above crystallization point, than, at that stage, the temperature is decreased below water solution crystallization point. In the solid/liquid phase transition, change in the fluid volume depending on the sample temperature is measured, as well as heat flow.

(20) Then, after full fluid crystallization, the temperature is raised above the melting point of that fluid. While raising the temperature, the heat flow and change in the fluid volume are measured.

(21) At each step of increasing and decreasing the pressure, a heat effect of sample wetting by the fluid is measured, and it is necessary to consider the heat effect of the injected fluid compression. Further, a wetting contact angle is calculated, for instance, by formula (6). Formula (6) may be simplified for the purpose of wetting angle calculation assuming that surface energy (.sub.lv) does not change in the temperature range in question. Then the formula is simplified and the limiting angle can be calculated from the ratio of the system inner energy change (U) to surface energy

(22) $\begin{array}{cc}\cos =\frac{-U}{{}_{\mathrm{lv}}}& \left(7\right)\end{array}$

(23) By integrating the heat flow over time at a given hydrostatic pressure (FIG. 1) and deducting the heat effect of fluid compression (FIG. 1, Peak 1) we obtain the value of the system inner energy change from wetting (U). Surface energy value .sub.lv can be calculated by multiplying tabular reference value of fluid surface tension (A. P. Babichev, N. A. Babushkina, A. M. Bratkovsky et al., Physical Values: Manual, Moscow, Energoatomizdat, 1991, p. 1232) and the sample specific surface. Sample specific surface can be measured, for instance, by gas adsorption method (Stephen Brunauer, P. H. Emmett, Edward Teller, J. Am. Chem. Soc., 1938, 60 (2), pp 309-319).

(24) Measured temperature shifts in sample pores associated with phase transition (liquid/solid) relative to a typical phase transition temperature in a void space (FIG. 2, peaks around 6 C. and 1 C.) may be used to determine pore sizes, for instance by formula (2). Fluid melting temperature in bulk is T.sub.0, and T.sub.m is the difference between fluid melting temperature in pores and in bulk, i.e., between peak values (FIG. 2). Both values are determined from experiment. Surface tension at solid/liquid interface (.sub.sl, ice/water in our case), fluid molar volume (.sub.l), and phase transition enthalpy (H) are tabular reference values (A. P. Babichev, N. A. Babushkina, A. M. Bratkovsky et al., Physical Values: Manual, Moscow, Energoatomizdat, 1991, p. 1232). Using said values, pore sizes can be determined for each temperature shift (difference) Tm. Then, pore volumes and typical pore sizes as measured are used to determine pore size distribution for a given sample.

(25) FIG. 1 relates to brine injection into a sandstone sample under a pressure of 1.0 MPa (10 bar); heat flow functional connection with time is shown under fluid pressure of 10 bars: reference experiment 1, fluid injection2, and repeated fluid injection3.

(26) The major narrow peak represents heat emission from fluid compression, and heat flow value is a function of pressure change value. Additional heat effects observed after pressure stabilization, for instance, the heat flow wide peak (FIG. 1, indicated by arrow) are associated with brine injection into the sample pore volume. Heat effects associated with brine injection are observed under other pressure values as well, for instance, 0.8 MPa (8 bar). Heat effect value, i.e., wetting heat value is used to determine the limiting angle (45) by formula (6). The resulting contact angle is in good agreement with the Ammot-Harvey index measured for that sample independently.

(27) Samples made by Asahi Company and used in earlier experimental calorimetric studies of ice/water phase transition (CPG100A and CPG300A) were taken for reference as samples with known pore sizes. Water melting temperature in pores was measured in compliance with the ISO 11357-1 Standard. Additional heat flow peaks occurring below water freezing point were observed in both CPG samples (FIG. 2). Temperature shift in phase transition depends on the size of pores. Considering that .sub.l, and H values for water and surface tension at the ice/water interface (.sub.sl=60.5 mJ/m2) are provided in reference tables, pore sizes can be calculated for the CPG samples (Formula 2). As follows from the calculation, the effective pore radius for CPG 100A (r.sub.eff=120 ); for CPG300A, (r.sub.eff=380 ). Those data show good correlation with pore sizes indicated by manufacturer.