Tortuosity measurement and the effects of finite pulse widths on xenon gas diffusion NMR studies of porous media.

We have extended the utility of NMR as a technique to probe porous media structure over length scales of approximately 100-2000 microm by using the spin 1/2 noble gas 129Xe imbibed into the system's pore space. Such length scales are much greater than can be probed with NMR diffusion studies of water-saturated porous media. We utilized Pulsed Gradient Spin Echo NMR measurements of the time-dependent diffusion coefficient, D(t), of the xenon gas filling the pore space to study further the measurements of both the pore surface-area-to-volume ratio, S/V(p), and the tortuosity (pore connectivity) of the medium. In uniform-size glass bead packs, we observed D(t) decreasing with increasing t, reaching an observed asymptote of approximately 0.62-0.65D(0), that could be measured over diffusion distances extending over multiple bead diameters. Measurements of D(t)/D(0) at differing gas pressures showed this tortuosity limit was not affected by changing the characteristic diffusion length of the spins during the diffusion encoding gradient pulse. This was not the case at the short time limit, where D(t)/D(0) was noticeably affected by the gas pressure in the sample. Increasing the gas pressure, and hence reducing D(0) and the diffusion during the gradient pulse served to reduce the previously observed deviation of D(t)/D(0) from the S/V(p) relation. The Pade approximation is used to interpolate between the long and short time limits in D(t). While the short time D(t) points lay above the interpolation line in the case of small beads, due to diffusion during the gradient pulse on the order of the pore size, it was also noted that the experimental D(t) data fell below the Pade line in the case of large beads, most likely due to finite size effects.     

The narrow pulse approximation and long length scale determination in xenon gas diffusion NMR studies of model porous media

We report a systematic study of xenon gas diffusion NMR in simple model porous media, random packs of mono-sized glass beads, and focus on three specific areas peculiar to gas-phase diffusion. These topics are: (i) diffusion of spins on the order of the pore dimensions during the application of the diffusion encoding gradient pulses in a PGSE experiment (breakdown of the narrow pulse approximation and imperfect background gradient cancellation), (ii) the ability to derive long length scale structural information, and (iii) effects of finite sample size. We find that the time-dependent diffusion coefficient, D(t), of the imbibed xenon gas at short diffusion times in small beads is significantly affected by the gas pressure. In particular, as expected, we find smaller deviations between measured D(t) and theoretical predictions as the gas pressure is increased, resulting from reduced diffusion during the application of the gradient pulse. The deviations are then completely removed when water D(t) is observed in the same samples. The use of gas also allows us to probe D(t) over a wide range of length scales and observe the long time asymptotic limit which is proportional to the inverse tortuosity of the sample, as well as the diffusion distance where this limit takes effect (approximately 1-1.5 bead diameters). The Padé approximation can be used as a reference for expected xenon D(t) data between the short and the long time limits, allowing us to explore deviations from the expected behavior at intermediate times as a result of finite sample size effects. Finally, the application of the Padé interpolation between the long and the short time asymptotic limits yields a fitted length scale (the Padé length), which is found to be approximately 0.13b for all bead packs, where b is the bead diameter.     

Probing four orders of magnitude of the diffusion time in porous silica glass with unconventional NMR techniques

The combined use of two unconventional NMR diffusometry techniques permits measurements of the self-diffusion coefficient of fluids confined in porous media in the time range from 100 microseconds to seconds. The fringe field stimulated echo technique (FFStE) exploits the strong steady gradient in the fringe field of a superconducting magnet. Using a standard 9.4 T (400 MHz) wide-bore magnet, for example, the gradient is 22 T/m at 375 MHz proton resonance and reaches 60 T/m at 200 MHz. Extremely short diffusion times can be probed on this basis. The magnetization grid rotating frame imaging technique (MAGROFI) is based on gradients of the radio frequency (RF) field. The RF gradients not necessarily need be constant since the data are acquired with spatial resolution along the RF gradient direction. MAGROFI is also well suited for unilateral NMR applications where all fields are intrinsically inhomogeneous. The RF gradients reached depend largely on the RF coil diameter and geometry. Using a conic shape, a value of at least 0.3 T/m can be reached which is suitable for long-time diffusion measurements. Both techniques do not require any special hardware and can be implemented on standard high RF power NMR spectrometers. As an application, the influence of the tortuosity increasing with the diffusion time is examined in a saturated porous silica glass.     

Studies of porous media by thermally polarized gas NMR: current status

Three examples of thermally polarized gas NMR performed at New Mexico Resonance are presented to demonstrate its unique advantages in porous media studies. 1) In-vivo animal lung imaging by Kuethe et al., in which useful quality 3D images of rat lungs were obtained in 30 min. It is conjectured that comparable human lung images would take much less time to make, possibly by the ratio of body weights, a factor of several hundred. 2) The success of the lung imaging suggested other porous media as candidates for thermally polarized gas NMR. Caprihan and coworkers obtained excellent images from partially sintered ceramics and Vycor glass. Since then, Beyea has developed the technique of spatially resolved BET curves for ceramics and other nanoporous solids. In this way, surface area, pore size, and porosity, averaged over an image voxel, can be spatially resolved. This greatly aids in the characterization of such materials, especially with regards to spatial heterogeneities. 3) Finally, we describe Codd's propagator experiments on propane gas flowing through a packed bed of 300 microm beads. In order to increase signal-to-noise ratio, the flowing gas was pressurized to 170 kPa. Excellent quality propagators showing the discrete nature of the bead pack were obtained. This type of information is not available in comparable liquid studies because most spins will not diffuse far enough to sample the walls in the time available.     

Ion transport in porous media studied by NMR.

Moisture and salt transport in masonry can give rise to damages. Therefore a detailed knowledge of the moisture and salt transport is essential for understanding the durability of masonry. A special NMR apparatus has been made allowing quasi-simultaneous measurements of both moisture and Na profiles in porous building materials. Using this apparatus both the absorption of a 4 M NaCl solution in a calcium silicate brick and the drying of a 3 M NaCl capillary saturated fired-clay brick have been studied. It was found that during the absorption process the Na ions clearly stay behind, which this is caused by adsorption of these ions to the pore surface. For the drying it was found that at the beginning of the drying process the ions accumulate near the surface. As the drying rate decreases, diffusion becomes dominant and the ion profile levels off again.     

Retardation of diffusion in porous media with stagnation zones

The problem of diffusion in a porous medium with stagnation zones, which reduce the rate of transport of particles, is considered. A method based on the concept of entropy barriers associated with the small size of inlets into the stagnation zone is used to solve the diffusion problem. This method made it possible to reduce the diffusion problem in 3D structures to a set of one-dimensional equations that can be solved analytically. The calculations yielded an exact value of the effective diffusion coefficient, which is reduced due to random delays in the stagnation zones, and the transition time required to attain this reduction.     

Research on the effective gas diffusion coefficient in dry porous media embedded with a fractal-like tree network

A theoretical model for the relative gas diffusion coefficient in dry porous media embedded with a Y-shaped fractal-like tree network is presented under the combination of bulk diffusion and Knudsen diffusion. The proposed model is expressed as a function of the length ratio, the diameter ratio, the branching level, the branching angle and the relative areal porosity. The effect of the structural parameters of the medium and the tortuosity on gas diffusion is analyzed in detail. Model predictions are compared with available experimental data, and a fair agreement between them is found.     

Novel NMR techniques for porous media research

Diffusion in porous media has been used as a probe of pore geometry in various NMR techniques. We will examine the effect of time-dependent diffusion in CPMG by showing that the diffusion time in CPMG is approximately the echo time, even in grossly inhomogeneous magnetic fields. Extension of the diffusion time in modified CPMG sequences is discussed. Diffusion in the susceptibility-contrast induced internal field is discussed as a means to probe pore size and pore shape. Finally, we present the general concept of two-dimensional relaxation-type experiments for study of molecules, fluids, materials and their dynamics that are characterized by spin relaxation and diffusion.     

Selective imaging of biofilms in porous media by NMR relaxation.

Nuclear magnetic resonance imaging (NMRI) techniques were employed to identify and selectively image biological films (biofilm) growing in aqueous systems. Biofilms are shown to affect both the longitudinal (T1) and transverse (T2) NMR relaxation time values of proximal water hydrogens. Results are shown for biofilm growth experiments performed in a transparent parallel-plate reactor. A comparison of biofilm distributions by both NMR and optical imaging yielded general agreement for both an open-flow system and an idealized porous system (the reactor without and with packed glass beads, respectively). The selective imaging of biofilm by relaxation NMRI is dependent upon the resolution of relaxation times for the fluid phases, dynamic range, and signal-to-noise ratio. For open-flow systems, the use of a rapid and quantitative T2-sorted NMRI technique was preferred. For porous systems where T2 values are generally more similar, a T1-weighted technique was preferred.     

Characterization of porous media structure by non linear NMR methods.

In this paper we discuss the possibility of modifying the multiple spin echoes existing theory, developed for a homogeneous system, to describe also an inhomogeneous system such as a porous medium. We report here the first experimental application of MSE methods to materials like travertine. The ratio A(2)/A(1) from water in travertine presents minima for characteristic values of the delay time tau, like what was previously observed in the trabecular bone. By a judicious choice of the delay time tau and of the G gradient strength, the MSE sequence can be made sensitive to a specific length-scale of the sample heterogeneity. Furthermore the MSE image shows a particular new contrast that makes the non linear NMR method very attractive for the assessment of variations of the porous structure in porous systems.     

Using NMR displacement imaging to characterize electroosmotic flow in porous media.

Pulsed field gradient nuclear magnetic resonance (PFG-NMR) and NMR imaging were used to study temporal and spatial domains of an electrokinetically-driven mobile phase through open and packed segments of capillaries. Characteristics like velocity distribution and an asymptotic dispersion are contrasted to viscous flow behavior. We show that electroosmotic flow in microchannel geometries can offer a significant performance advantage over the pressure-driven flows at comparable Peclét numbers, indicating that velocity extremes in the pore space of open tubes and packed beds are drastically reduced. An inherent problem of capillary electrochromatography that we finally address is the existence of wall effects when in the general case the surface zeta-potentials of the capillary inner wall and the adsorbent particles are different. Using dynamic NMR microscopy we were able resolve this systematic velocity inequality of the flow pattern which strongly influences axial dispersion and may be responsible for long time-tails of velocity distribution in the mobile phase.     

2D relaxation/diffusion correlations in porous media

2D correlations between NMR relaxation and/or diffusion have been used to investigate water and oil dynamics in food and micro-emulsion systems. In the case of Mozzarella and Gouda cheese samples, a significant change in D/T2 correlation is appearing with cheese aging. In the case of a water/toluene micro-emulsion, some evidence for coalescence effects is suggested by D/D exchange spectra.     

Quantitative evaluation of porous media wettability using NMR relaxometry

We propose a new method to determine wettability indices from NMR relaxometry. The new method uses the sensitivity of low field NMR relaxometry to the fluid distribution in oil-water saturated porous media. The model is based on the existence of a surface relaxivity for both oil and water, allowing the determination of the amount of surface wetted either by oil or by water. The proposed NMR wettability index requires the measurement of relaxation time distribution at four different saturation states. At the irreducible water saturation, we determine the dominant relaxation time of oil in the presence of a small amount of water, and at the oil residual saturation, we determine the dominant relaxation time of water in the presence of a small amount of oil. At 100% water and 100% oil saturation, we determine the surface relaxivity ratio. The interaction of oil with the surface is also evidenced by the comparison of the spin-lattice (T1) and spin-locking (T1rho) relaxation times. The new NMR index agrees with standard wettability measurements based on drainage-imbibition capillary pressure curves (USBM test) in the range [-0.3-1].     

A fractional diffusion equation for a marker in porous media

To study the properties of the flow of a marker through an irregular packed bed, porosity is described by stochastic processes with exponentially decaying correlation functions. We show that the ensemble average concentration of the marker satisfies a fractional diffusion equation. (c) 2001 American Institute of Physics.     

Direct measurement of porous media local hydrodynamical permeability using gas MRI.

The concept of hydraulic permeability is at the core of modeling single phase or multi-phase flow in heterogeneous porous media, as it is the spatial distribution of the permeability that primarily governs the behavior of fluid flow in the medium. To date, the modeling of fluid flow in porous media has been hampered by poor estimates of local permeability. Magnetic Resonance Imaging is well known for its ability to measure non-invasively the local density and flow rate of different fluids saturating porous media [1,2]. In this paper we demonstrate the first non-invasive method for the direct measurement of a single projection of the local permeability tensor of a porous medium using gas-phase MRI. The potential for three-dimensional imaging of the medium permeability is also discussed. The limitations of the method are listed and results are presented in a model porous medium as well as in a real oil reservoir rock.