Probing the Pore Size of Porous Ceramics with Controlled Amount of Magnetic Impurities via Diffusion Effects on the CPMG Technique

In our work, we will explore the possibility of implementing the well-known Carr–Purcell–Meiboom–Gill pulse sequence to determine the pore size of porous ceramics with magnetic impurities. The proposed approach exploits the diffusion dependence of the spin-echo signal in the presence of internal gradients occurring as a result of susceptibility contrast between the porous matrix and the confined liquid. For calibrating the technique, a comparison of the pore size data with those extracted from the so-called DDIF technique (DDIF, decay due to diffusion in the internal fields) is performed. This approach can be applied for nondestructive in situ characterization of soils, concrete, biological tissues or other structures with micrometer pore size.     

Pore sizes and pore connectivity in rocks using the effect of internal field.

A novel nuclear magnetic resonance method has been applied to several sandstone rocks to measure the pore size distribution using the magnetization decay due to diffusion in the internal magnetic field (DDIF). By comparing the results of the DDIF and Hg porosimetry experiments, a clear picture of pore connectivity emerges. The pore body diameter can be defined using the DDIF data and is found to have a clear trend as a function of porosity.     

The internal magnetic field distribution, and single exponential magnetic resonance free induction decay, in rocks

When fluid saturated porous media are subjected to an applied uniform magnetic field, an internal magnetic field, inside the pore space, is induced due to magnetic susceptibility differences between the pore-filling fluid and the solid matrix. The microscopic distribution of the internal magnetic field, and its gradients, was simulated based on the thin-section pore structure of a sedimentary rock. The simulation results were verified experimentally. We show that the 'decay due to diffusion in internal field' magnetic resonance technique may be applied to measure the pore size distribution in partially saturated porous media. For the first time, we have observed that the internal magnetic field and its gradients in porous rocks have a Lorentzian distribution, with an average gradient value of zero. The Lorentzian distribution of internal magnetic field arises from the large susceptibility contrast and an intrinsic disordered pore structure in these porous media. We confirm that the single exponential magnetic resonance free induction decay commonly observed in fluid saturated porous media arises from a Lorentzian internal field distribution. A linear relationship between the magnetic resonance linewidth, and the product of the susceptibility difference in the porous media and the applied magnetic field, is observed through simulation and experiment.     

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.     

Quantitative measurement and imaging of transport processes in plants and porous media by 1H NMR.

NMR and MRI have been applied to transport processes, that is, net flow and diffusion/perfusion, of water in whole plants, cells, and porous materials. By choosing proper time windows and pulse sequences, magnetic resonance imaging can be made selective for each of the two transport processes. For porous media and plant cells the evolution of the spatial distribution of excited spins has been determined by q-space imaging, using a 20 MHz pulsed 1H NMR imager. The results of these experiments are explained by including spin-relaxation and exchange at boundaries. A 10 MHz portable 1H NMR spectrometer is described, particularly suitable to study the response of net flow in plants and canopies to changing external conditions.     

NMR Relaxation and Diffusion Measurements on Iron(III)-Doped Kaolin Clay

Kaolin clay samples were mixed with various amounts of Fe₂O₃ powder. The influence of this magnetic impurity on NMR relaxation and diffusion measurements on the water in this porous material was investigated. The NMR relaxation measurements showed a nearly mono-exponential decay, leading to the conclusion that the pore size distribution of the clay samples is either narrow and/or that the pores are interconnected very well. Both the longitudinal and the transverse relaxation rate depend linearly on the concentration of the Fe₂O₃ impurity. The NMR diffusion measurements revealed that the Fe₂O₃ causes internal magnetic field gradients that largely exceed the maximum external gradient that could be applied by our NMR apparatus (0.3 T/m). Additional SQUID measurements yielded the magnetization and magnetic susceptibility of the samples at the magnetic field strength used in the NMR measurements (0.8 T). A theoretical estimate of the internal magnetic field gradients leads to the conclusion that the water in the porous clay samples cannot be described by the commonly observed motional averaging regime. Probably an intermediate or a localization regime is induced by the large internal gradients, which are estimated to be on the order of 1 to 10 T/m in the pore volume and may exceed 1000 T/m at the pore surface.     

Two-dimensional nuclear magnetic resonance petrophysics

Two-dimensional nuclear magnetic resonance (2D NMR) opens a wide area for exploration in petrophysics and has significant impact to petroleum logging technology. When there are multiple fluids with different diffusion coefficients saturated in a porous medium, this information can be extracted and clearly delineated from CPMG measurements of such a system either using regular pulsing sequences or modified two window sequences. The 2D NMR plot with independent variables of T2 relaxation time and diffusion coefficient allows clear separation of oil and water signals in the rocks. This 2D concept can be extended to general studies of fluid-saturated porous media involving other combinations of two or more independent variables, such as chemical shift and T1/T2 relaxation time (reflecting pore size), proton population and diffusion contrast, etc.     

Multidimensionally resolved pore size distributions

A novel method of determining median pore size and pore size distributions as a function of spatial position inside a porous sample is described. Pore sizes have been measured with 1-, 2- and 3-dimensional spatial resolution, using NMR cryoporometry in conjunction with magnetic resonance imaging techniques. The method is suitable for pore diameters in the range of 30 Å to over 2000 Å pore diameter, and is based on the technique of freezing a liquid in the pores and measuring the melting temperature by nuclear magnetic resonance. Since the melting point is depressed for crystals of small size, the melting point depression gives a measurement of pore size.     

Influential Factors of Internal Magnetic Field Gradient in Reservoir Rock and Its Effects on NMR Response

Nuclear magnetic resonance (NMR) plays a significant role in porous media analysis and petroleum exploration, but its response is significantly influenced by the internal magnetic field gradient in fluid saturated porous medium, which obviously limits the accuracy of rock core analysis and logging interpretation. The influential factors of the internal magnetic field gradient in formation and its influences on NMR response are studied in this paper, based on NMR mechanism through one- and two-dimensional core NMR experiments. The results indicate that the internal magnetic field gradient is positively correlated with the static magnetic field strength and the magnetic susceptibility difference between pore fluid and solid grains, while it presents negative correlation with pore radius. The internal magnetic field gradient produces an additional diffusion relaxation in hydrogen relaxation system and accelerates the attenuation of magnetization vector. As a result, T₂ spectrum shifts to the left and NMR porosity and diffusion coefficient of the fluid could be inaccurate. This research sets a foundation for the NMR porosity correction and fluid distribution on T₂-G maps based on the internal magnetic field gradient correction.     

岩心孔隙介质中流体的核磁共振弛豫

弛豫时间是核磁共振研究中的一个重要参数,岩心孔隙介质流体的弛豫过程是自由流体弛豫机制、表面弛豫机制和流体的扩散弛豫机制共同作用的结果,它包含了丰富的孔隙和流体本身的信息. 弛豫时间和自扩散系数的测量及对弛豫时间的分析是核磁共振技术应用于岩心分析和石油勘测的重要内容.     

Internal Magnetic Field Gradients in Heterogeneous Porous Systems: Comparison Between Spin-Echo and Diffusion Decay Internal Field (DDIF) Method

Two techniques used for evaluating internal magnetic field gradient (G _i), spin-echo (SE) and diffusion decay internal field (DDIF), were investigated at 9.4 T and compared in porous systems characterized by different pores size ranging from 4 to 96 μm with magnetic susceptibility difference between solid and liquid phase, \(\Delta \chi\)  ≈ 1.6 ppm. Since diffusion of a fluid in a solid porous matrix plays a role in both SE and DDIF methods, we investigated these two different methods by highlighting their dependence on characteristic parameters and length scales used to describe diffusion behavior of fluids in porous systems. Therefore, G _i behavior as a function of the dephasing length (l _g), diffusion length (l _d) and pores size (l _s) was obtained. Moreover G _i was evaluated by using both free diffusion and measured apparent diffusion coefficient of water, to quantify diffusion effect in different porous samples. This study gives more insight into the physical dynamics process to explain contrast mechanisms recently exploited by DDIF and SE applications for cancellous bone quality measurements.     

Random-walk technique for simulating NMR measurements and 2D NMR maps of porous media with relaxing and permeable boundaries

We revisit random-walk methods to simulate the NMR response of fluids in porous media. Simulations reproduce the effects of diffusion within external inhomogeneous background magnetic fields, imperfect and finite-duration B(1) pulses, T(1)/T(2) contrasts, and relaxing or permeable boundaries. The simulation approach consolidates existing NMR numerical methods used in biology and engineering into a single formulation that expands on the magnetic-dipole equivalent of spin packets. When fluids exhibit low T(1)/T(2) contrasts and when CPMG pulse sequences are used to acquire NMR measurements, we verify that classical NMR numerical models that neglect T(1) effects accurately reproduce surface magnetization decays of saturated granular porous media regardless of the diffusion/relaxation regime. Currently, analytical expressions exist only for the case of arbitrary pore shapes within the fast-diffusion limit. However, when fluids include several components or when magnetic fields are strongly inhomogeneous, we show that simulations results obtained using the complete set of Bloch's equations differ substantially from those of classical NMR models. In addition, our random-walk formulation accurately reproduces magnetization echoes stemming from coherent-pathway calculations. We show that the random-walk approach is especially suited to generate parametric multi-dimensional T(1)/T(2)/D NMR maps to improve the characterization of pore structures and saturating fluids.     

Restricted diffusion and exchange of water in porous media: average structure determination and size distribution resolved from the effect of local field gradients on the proton NMR spectrum

NMR Pulsed field gradient measurements of the restrained diffusion of confined fluids constitute an efficient method to probe the local geometry in porous media. In most practical cases, the diffusion decay, when limited to its principal part, can be considered as Gaussian leading to an apparent diffusion coefficient. The evolution of the latter as a function of the diffusion interval yields average information on the surface/volume ratio of porosities and on the tortuosity of the network. In this paper, we investigate porous model systems of packed spheres (polystyrene and glass) with known mean diameter and polydispersity, and, in addition, a real porous polystyrene material. Applying an Inverse Laplace Transformation in the second dimension reveals an evolution of the apparent diffusion coefficient as a function of the resonance frequency. This evolution is related to a similar evolution of the transverse relaxation time T2. These results clearly show that each resonance frequency in the water proton spectrum corresponds to a particular magnetic environment produced by a given pore geometry in the porous media. This is due to the presence of local field gradients induced by magnetic susceptibility differences at the liquid/solid interface and to slow exchange rates between different pores as compared to the frequency differences in the spectrum. This interpretation is nicely confirmed by a series of two-dimensional exchange experiments.     

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.     

Simulation and experimental verification of the diffusion in an anisotropic fiber phantom

Diffusion weighted magnetic resonance imaging enables the visualization of fibrous tissues such as brain white matter. The validation of this non-invasive technique requires phantoms with a well-known structure and diffusion behavior. This paper presents anisotropic diffusion phantoms consisting of parallel fibers. The diffusion properties of the fiber phantoms are measured using diffusion weighted magnetic resonance imaging and bulk NMR measurements. To enable quantitative evaluation of the measurements, the diffusion in the interstitial space between fibers is modeled using Monte Carlo simulations of random walkers. The time-dependent apparent diffusion coefficient and kurtosis, quantifying the deviation from a Gaussian diffusion profile, are simulated in 3D geometries of parallel fibers with varying packing geometries and packing densities. The simulated diffusion coefficients are compared to the theory of diffusion in porous media, showing a good agreement. Based on the correspondence between simulations and experimental measurements, the fiber phantoms are shown to be useful for the quantitative validation of diffusion imaging on clinical MRI-scanners.     

An NMR-Based Analysis of Soil–Water Characteristics

Both direct and indirect methods for determining soil–water characteristic curves rely on determination of some empirical coefficients, which may not necessarily represent real microscopic mechanisms. Proton nuclear magnetic resonance (NMR) is a powerful tool for investigating water content and their interaction with solid particles in porous media. The NMR technique is widely used in food science and petroleum. In the present study, proton NMR spin–spin relaxation time (T ₂) distribution measurement is integrated with a Tempe apparatus to characterize the hydraulic processes of unsaturated soils, shedding insights into the microscopic mechanisms of pore water distribution and migration in the soil during hydraulic cycles. It is revealed that during a drying process the drainage of pore water occurs sequentially from larger pores to smaller pores, whereas in a wetting process the water invades into the soil sequentially from smaller pores to larger pores. A new procedure is developed which can be used to determine the pore size distribution of the soil based on the NMR T ₂ distribution measurements; compared to the traditional methods, the new method is rapid and non-destructive. The new procedure is validated by comparing the new result with the measurement of the mercury intrusion porosimetry.     

Magnetic resonance for fluids in porous media at the University of Bologna

The magnetic resonance in porous media (MRPM) community is now a vast community of scientists from all over the world who recognize magnetic resonance as an instrument of choice for the characterization of pore space and of the distribution, diffusion and flow of fluids inside a vast range of different materials. The MRPM conferences are the occasions in which this community gets together, compares notes and grows. The scene was different in 1990, when this series of conferences was promoted at Bologna. I will go briefly over the history of these events, showing the role played by the University of Bologna and in particular by the intuition, ingenuity and passion of Giulio Cesare Borgia. The MRPM work at Bologna began in the mid-1980s. New correlations were found among parameters from NMR relaxation measurements and oil field parameters such as porosity, permeability to fluid flow, irreducible water saturation, residual oil saturation and pore-system surface-to-volume ratio, and fast algorithms were developed to give the different NMR parameters. Interest in valid interpretation of data led to extensive work also on the inversion of multiexponential relaxation data and the effects of inhomogeneous fields from susceptibility differences on distributions of relaxation times. In the last few years, extensive developments were made of combined magnetic resonance imaging and relaxation measurements in different fields.     

Determining pore sizes using an internal magnetic field

A concept is proposed to measure the pore size length scale by the internal magnetic field (B_i) in porous materials. The spatial distribution of the magnetic field inhomogeneity, a result of the magnetic susceptibility contrast between the porous material and the fluid, reflects the underlying pore geometry. Diffusion in B_i causes the initial decay of magnetization. At long times, the effect of B_i saturates when the diffusion length reaches a characteristic pore size. This method is independent of surface spin relaxation in determining pore sizes. Nuclear magnetic resonance experiments on packed glass beads and sedimentary rock samples will be presented.     

Xenon porometry: a novel method for the derivation of pore size distributions

Xenon porometry is a novel method used for characterizing porous materials by the (129)Xe nuclear magnetic resonance of xenon gas. With the method, the diffusion of gas is slowed down by immersing the material in a medium, which can be in liquid or solid state during measurements. Because of slow diffusion, the signal of a xenon atom is characteristic of the properties of only one pore, and the composite signal of all atoms represents the distribution of properties. The method is especially applicable for determining pore size distribution because the chemical shifts of two different xenon signals (one from liquid and the other from gas pockets in solid) are dependent on pore size. Therefore, the shapes of these signals represent pore size distribution function. In addition, the porosity of the material can be determined by comparing the intensities of two signals. This article focuses on describing xenon signals observed from gas pockets in a solid medium, which has turned out to be most convenient for pore size determination.     

Probing porous media with gas diffusion NMR.

We show that gas diffusion nuclear magnetic resonance (GD-NMR) provides a powerful technique for probing the structure of porous media. In random packs of glass beads, using both laser-polarized and thermally polarized xenon gas, we find that GD-NMR can accurately measure the pore space surface-area-to-volume ratio, S/V rho, and the tortuosity, alpha (the latter quantity being directly related to the system's transport properties). We also show that GD-NMR provides a good measure of the tortuosity of sandstone and complex carbonate rocks.