Fractal geometry of surface areas of sand grains probed by pulsed field gradient NMR

Pore Structure Characterization of Eocene Low-Permeability Sandstones via Fractal Analysis and Machine Learning: An Example from the Dongying Depression, Bohai Bay Basin, China

Poroperm analysis, mercury injection capillary pressure (MICP), and nuclear magnetic resonance (NMR) measurements were performed to delineate the pore structures and fractal behaviors of the Eocene low-permeability sandstones in the Dongying Depression, Bohai Bay Basin, China. Three types of pore structures (I, II, and III) have been classified by applying the self-organizing map (SOM) clustering model. Comparative analysis of three different fractal models indicates that the MICP tubular model and NMR model are quite effective for pore structure characterization. The results show that the reservoirs generally exhibit high fractal dimensions, indicative of complex pore structures. The presence of small pore throats is primarily responsible for the heterogeneities and complexities in the Eocene low-permeability sandstones. A modified Winland model was established for the permeability estimation using MICP data. Different from high-permeability reservoirs or unconventional (e.g., shale and tight formation) reservoirs, r ₁₀ is the best parameter for permeability estimation, indicating that the permeability of the Eocene low-permeability sandstones is largely controlled by the large pore systems. Additionally, a porosity model derived from movable fluids using NMR data has been established and provided better prediction effect compared with the classic Coates and Schlumberger Doll Research (SDR) models. Fractal analysis and permeability estimation are shown to be quite effective for investigating microscopic behaviors and in predicting the reservoir quality of low-permeability sandstone reservoirs.

On modeling

Mapping tissue microstructure with MRI holds great promise as a noninvasive window into tissue organization at the cellular level. Having originated within the realm of diffusion NMR in the late 1970s, this field is experiencing an exponential growth in the number of publications. At the same time, model-based approaches are also increasingly incorporated into advanced MRI acquisition and reconstruction techniques. However, after about two decades of intellectual and financial investment, microstructural mapping has yet to find a single commonly accepted clinical application. Here, we suggest that slow progress in clinical translation may signify unresolved fundamental problems. We outline such problems and related practical pitfalls, as well as review strategies for developing and validating tissue microstructure models, to provoke a discussion on how to bridge the gap between our scientific aspirations and the clinical reality. We argue for recalibrating the efforts of our community toward a more systematic focus on fundamental research aimed at identifying relevant degrees of freedom affecting the measured MR signal. Such a focus is essential for realizing the truly revolutionary potential of noninvasive three-dimensional in vivo microstructural mapping.

Measurement of interfacial area from NMR time dependent diffusion and relaxation measurements

The interfacial area between two immiscible phases in porous media is an important parameter for describing and predicting 2 phase flow. Although present in several models, experimental investigations are sparse due to the lack of appropriate measurement techniques. We propose two NMR techniques for the measurement of oil-water interfacial area: (i) a time dependent NMR diffusion technique applicable in static conditions, similar to those used for the measurement of the solid specific surface of a porous media, and (ii) a fast relaxation technique applicable in dynamic conditions while flowing, based on an interfacial relaxation mechanism induced by the inclusion of paramagnetic salts in the water phase. For dodecane relaxing on doped water, we found an oil interfacial relaxivity of 1.8μm/s, large enough to permit the measurement of specific interfacial surface as small as 1000cm²/cm³. We demonstrate both NMR techniques in drainage followed by imbibition, in a model porous media with a narrow pore size distribution. While flowing, we observe that the interfacial area is larger in imbibition than in drainage, implying a different organization of the oil phase. In a carbonate sample with a wide pore size distribution, we evidence the gradual invasion of the smallest pores as the oil-water pressure difference is increased.

Grain sizing in porous media using Bayesian magnetic resonance

We introduce a Bayesian inference approach to analyze magnetic resonance data of granular solids. To characterize structure using magnetic resonance, it is usual to acquire data in k space which are then Fourier transformed to obtain an image. An alternative approach, adopted here, is to utilize the Rayleigh distribution observed in the signal intensity for a given k when a random selection of grains is measured in k space, to define a likelihood function for Bayesian analysis. This Bayesian likelihood function is used to noninvasively characterize grains within a porous medium on length scales below the practical resolution of magnetic resonance imaging. A pore size distribution is then calculated from the measured grain size distribution using a Monte Carlo approach. We demonstrate this general technique with specific examples of water-saturated rock cores.

Temporal scaling characteristics of diffusion as a new MRI contrast: findings in rat hippocampus

Features of the diffusion-time dependence of the diffusion-weighted magnetic resonance imaging (MRI) signal provide a new contrast that could be altered by numerous biological processes and pathologies in tissue at microscopic length scales. An anomalous diffusion model, based on the theory of Brownian motion in fractal and disordered media, is used to characterize the temporal scaling (TS) characteristics of diffusion-related quantities, such as moments of the displacement and zero-displacement probabilities, in excised rat hippocampus specimens. To reduce the effect of noise in magnitude-valued MRI data, a novel numerical procedure was employed to yield accurate estimation of these quantities even when the signal falls below the noise floor. The power-law dependencies characterize the TS behavior in all regions of the rat hippocampus, providing unique information about its microscopic architecture. The relationship between the TS characteristics and diffusion anisotropy is investigated by examining the anisotropy of TS, and conversely, the TS of anisotropy. The findings suggest the robustness of the technique as well as the reproducibility of estimates. TS characteristics of the diffusion-weighted signals could be used as a new and useful marker of tissue microstructure.

Effect of pore structure on capillary condensation in a porous medium

The Kelvin equation relates the equilibrium vapor pressure of a fluid to the curvature of the fluid-vapor interface and predicts that vapor condensation will occur in pores or irregularities that are sufficiently small. Past analyses of capillary condensation in porous systems with fractal structure have related the phenomenon to the fractal dimension of the pore volume distribution. Recent work, however, suggests that porous systems can exhibit distinct fractal dimensions that are characteristic of both their pore volume and the surfaces of the pores themselves. We show that both fractal dimensions have an effect on the thermodynamics that governs capillary condensation and that previous analyses can be obtained as limiting cases of a more general formulation.

Capillary pressure in a porous medium with distinct pore surface and pore volume fractal dimensions

The relationship between capillary pressure and saturation in a porous medium often exhibits a power-law dependence. The physical basis for this relation has been substantiated by assuming that capillary pressure is directly related to the pore radius. When the pore space of a medium exhibits fractal structure this approach results in a power-law relation with an exponent of 3-D(v), where D(v) is the pore volume fractal dimension. However, larger values of the exponent than are realistically allowed by this result have long been known to occur. Using a thermodynamic formulation for equilibrium capillary pressure we show that the standard result is a special case of the more general exponent (3-D(v))(3-D(s)) where D(s) is the surface fractal dimension of the pores. The analysis reduces to the standard result when D(s)=2, indicating a Euclidean relationship between a pore's surface area and the volume it encloses, and allows for a larger value for the exponent than the standard result when D(s)>2 .

Proton dynamics in the perchloric acid clathrate hydrate HClO4⋅5.5H2O

In the perchloric acid clathrate hydrate HClO4.5.5H2O, the perchlorate anions are contained inside an aqueous host crystalline matrix, positively charged because of the presence of delocalized acidic protons. Our experimental results demonstrate that the microscopic mechanisms of proton conductivity in this system are effective on a time scale ranging from nanosecond to picosecond. In the present paper, we discuss more specifically on the relaxation processes occurring on a nanosecond time scale by combining high-resolution quasielastic neutron scattering and 1H pulse-field-gradient nuclear magnetic resonance experiments. The combination of these two techniques allows us to probe proton dynamics in both space and time domains. The existence of two types of proton dynamical processes has been identified. The slowest one is associated to long-range translational diffusion of protons between crystallographic oxygen sites and has been precisely characterized with a self-diffusion coefficient of 3.5 x 10(-8) cm2/s at 220 K and an activation energy of 29.2+/-1.4 kJ/mol. The fastest dynamical process is due to water molecules' reorientations occurring every 0.7 ns at 220 K with an activation energy of 17.4+/-1.5 kJ/mol. This powerful multitechnique approach provides important information required to understand the microscopic origin of proton transport in an ionic clathrate hydrate.

Short-time restricted diffusion in a static gradient and the attenuation of individual coherence pathways

We experimentally explore some of the implications of a recent theoretical study [J. Magn. Reson. 64 (2003) 145] for the measurement of restricted diffusion in connected porous media in a static gradient. In particular, we examine how restriction affects the short-time attenuation of different coherence pathways, all excited with the same sequence of slice-selective radiofrequency (RF) pulses, and how the various pathways make the transition to the long-time or tortuosity regime. We confirm that every pathway contains equivalent diffusional information and, for short times, yields the surface-to-volume ratio (S/V) of the confining space. We find also, in agreement with the theoretical predictions, that different pathways are controlled by different time scales and, thus, exhibit different sensitivity to restriction. This property might be exploited when designing optimal sequences to study restricted motion.

Boundary effects of molecular diffusion in nanoporous materials: A pulsed field gradient nuclear magnetic resonance study

The boundary conditions of intraparticle diffusion in nanoporous materials may be chosen to approach the limiting cases of either absorbing or reflecting boundaries, depending on the host-guest system under study and the temperature of measurement. Pulsed field gradient nuclear magnetic resonance is applied to monitor molecular diffusion of n-hexane and of an n-hexane-tetrafluoromethane mixture adsorbed in zeolite crystallites of type NaX under either of these limiting conditions. Taking advantage of the thus-established peculiarities of mass transfer at the interface between the zeolite bulk phase and the surrounding atmosphere, three independent routes for probing the crystal size are compared. These techniques are based on (i) the measurement of the effective diffusivity under complete confinement, (ii) the application of the so-called NMR tracer desorption technique, and (iii) an analysis of the time dependence of the effective diffusivity in the short-time limit where, by an appropriate variation of the adsorbate and the measuring conditions, the limiting cases of reflecting and adsorbing boundaries could be considered. All these techniques are found to yield coinciding results, which are in excellent agreement with the crystal sizes determined by microscopy.

Effect of internal gradients in the nuclear magnetic resonance measurement of the surface-to-volume ratio

We consider a system of spins diffusing in a static inhomogeneous (nonuniform-gradient) magnetic field B in a restricted geometry and in the presence of surface relaxation. We show that the short-time diffusional decay of nuclear magnetization is controlled by the field scattering kernel F(t) identical with [B(t)-B(0)](2), which is a measure of the average field inhomogeneity sampled by the spins in time t and does not depend on the particular sequence of radio-frequency pulses used. Magnetization in arbitrary sequences can be straightforwardly computed by evaluating elementary integrals of F(t). Diffusion takes place while the field is on, so that the spins precess as they diffuse, in contrast to the simpler problem of purely classical diffusion considered in [P. P. Mitra, P. N. Sen, and L. M. Schwartz, Phys. Rev. B 47, 8565 (1993)] which is applicable only to the ideal pulsed-field gradient experiment. We compute the short-time asymptotic form of F(t) and find that it depends on the surface-to-volume ratio (S/V) of the pore space as well as on the average of the gradients over the bounding surface. In a system with nonuniform gradients that vary faster near the surface than in the bulk, as for internal susceptibility fields, this gradient surface average may be much larger than the gradients in the bulk, significantly enhancing the apparent S/V. We discuss the application of our results to the widely used Carr-Purcell-Meiboom-Gill (CPMG) pulse sequence as well as proposing a modification of it, which we term "padded" CPMG, that may be preferable in systems with significant surface relaxation. We indicate how each sequence can be used to probe the internal fields.