Magnetic resonance is an important noninvasive technology across life sciences and industry. Free induction decay is the simplest
1
H magnetic resonance measurement method and an important means of ...probing fast-decaying signals in porous materials such as rocks, lung, and bone. It is commonly assumed that the free induction decay in geological porous materials is single-exponential. We experimentally observed two regimes of free induction decay behavior in geological porous materials: single-exponential and non-exponential decay. Numerical simulations that match experimental data highlight the effect of mass diffusion, especially in the single-exponential behavior. These two regimes of free induction decay in porous materials are associated with a bifurcation point in the solutions of the Bloch–Torrey equation for diffusion of fluids in confined domains in the presence of internal magnetic field gradients. This finding facilitates the extraction of absolute internal magnetic field gradient intensities from simple free induction decay measurements in the laboratory and field. This work also warns against common single-exponential assumptions in surface magnetic resonance methods employed in surveying underground water aquifers.
Graphic Abstract
Solid-state magnetic resonance is a unique technique that can reveal the dynamics of complex biological systems with atomic resolution. Longitudinal relaxation is a mechanism that returns ...longitudinal nuclear magnetization to its thermal equilibrium by incoherent processes. The measured longitudinal relaxation rate constant however represents the combination of both incoherent and coherent contributions to the change of nuclear magnetization. This work demonstrates the effect of magic angle spinning rate on the longitudinal relaxation rate constant in two model compounds: L-histidine hydrochloride monohydrate and glycine serving as proxies for isotopically-enriched biological materials. Most notably, it is demonstrated that the longitudinal N15 relaxation of the two nitrogen nuclei in the imidazole ring in histidine is reduced by almost three orders of magnitude at the condition of rotational resonance with the amine, while the amine relaxation rate constant is increased at these conditions. The observed phenomenon may have radical implications for the solid-state magnetic resonance in biophysics and materials, especially in the proper measurement of dynamics and as a selective serial transfer step in dynamic nuclear polarization.
Biological tissue, pharmaceutical tablets, wood, porous rocks, catalytic reactors, concrete, and foams are examples of heterogeneous systems that may contain one or several fluid phases. Fluids in ...such systems carry chemical species that may participate in chemical reactions in the bulk of a fluid, as homogeneous reactions, or at the fluid/fluid or fluid/solid interfaces, as heterogeneous reactions. Magnetic resonance relaxation measures the return of
1
H nuclear magnetization in chemical species of these fluids to an equilibrium state in a static magnetic field. Despite the perceived difference between reaction–diffusion and relaxation–diffusion in heterogeneous systems, similarities between the two are remarkable. This work draws a close parallel between magnetic resonance relaxation–diffusion and chemical reaction–diffusion for elementary unitary reaction
A
→
B
in a dilute solution—both in heterogeneous systems. A striking similarity between the dimensionless numbers that characterize their relevant behavior is observed: the Damköhler number of the second kind
Da
II
for reaction and the Brownstein–Tarr number
BT
i
for relaxation. The new vision of analogy between reaction- and magnetic resonance relaxation–diffusion in heterogeneous systems encourages the exploitation of similarities between reaction and relaxation processes to noninvasively investigate the dynamics of chemical species and reactions. One such example of importance in chemical engineering is provided for solid–fluid reaction in packed beds.
We present a systematic study of crude oil–brine–rock interactions in tight chalk cores at reservoir conditions. Flooding experiments are performed on outcrops (Stevns Klint) as well as on reservoir ...core plugs from Dan field, the Ekofisk and Tor formations. These studies are carried out in core plugs with reduced pore volumes, i.e., short core samples and aged with a dynamic ageing method. The method was evaluated by three different oil compositions. A series of synthetic multicomponent brines and designed fluid injection scenarios are investigated; injection flow rates are optimized to ensure that a capillary-dominant regime is maintained. Changes in brine compositions and fluid distribution in the core plugs are characterized using ion chromatography and X-ray computed tomography, respectively. First, we show that polar components in the oil phase play a major role in wettability alteration during ageing; this controls the oil production behavior. We also show that, compared to seawater, both formation water and ten-times-diluted seawater are better candidates for enhanced oil recovery in the Dan field. Finally, we show that the modified flow zone indicator, a measure of rock quality, is likely the main variable responsible for the higher oil recoveries observed in Tor core samples.
An image processing workflow is presented for the characterization of pore and grain size distributions in porous geological samples from X-ray microcomputed tomography (μCT) and scanning electron ...microscopy (SEM) images. The pore and grain size distributions of five sandstone samples including Berea, Buff Berea, Nugget, Castlegate, and Bentheimer, and one carbonate sample, Indiana limestone, are extracted using the proposed workflow. Two-dimensional size distributions acquired from SEM images were found to be biased toward smaller sizes misrepresenting the actual 3D distributions. Stereological techniques unfolded the measured 2D size distributions from SEM images to 3D distributions comparable with μCT results. While larger pores and grains can easily be detected from μCT and SEM images, the quantification of small-scale heterogeneities is severely influenced by their limits of resolution. We show that microstructural details resolved by SEM can significantly impact the pore and grain size distributions in sandstone and carbonate rock samples. For example, SEM-resolved microporosities in Indiana limestone result in bimodal distributions of pore and grain sizes, whereas μCT observations exhibit unimodal distributions. The acquired images and processed results are openly available and may be used by researchers investigating image processing, magnetic resonance relaxation or fluid flow simulations in natural rocks. The proposed methodology can be implemented to process μCT and SEM images of natural rocks as well as other types of porous materials.
•Image processing workflow computes pore/grain size distributions from images.•Stereological unfolding of pore/grain size from SEM agrees with that of μCT images.•SEM resolves microstructural details that affect pore and grain size distributions.•The agreement between μCT and SEM results was consistent in six sedimentary rocks.•A dataset of μCT and SEM images of six sedimentary rock samples made public.
The Brooks-Corey power-law capillary pressure model is commonly imposed on core analysis data without verifying the validity of its underlying assumptions. The Brooks-Corey model, originally ...developed to model the pressure head during the drainage of soil, is only valid at low wetting phase saturations. However, such models are often applied in petroleum production simulations and may lead to erroneous recovery factors when the saturation range of interest is far from the end points. We demonstrate that exponential models work much better for capillary pressure compared to the Brooks-Corey model over a wide saturation range. Mercury injection porosimetry, petrographic image analysis, and magnetic resonance studies suggest that the pore and throat size distribution in many rocks are log-normally distributed. This fact was previously employed to calculate the capillary pressure function as a function of saturation for pore size distributions described by atruncated log-normal distribution. Employing a Taylor series expansion, we simplify the random fractal capillary pressure model of Hunt to
P
c
= exp(
a
−
bS
), where
S
is the wetting phase saturation, and
a
and
b
characteristic of the porous medium. An extensive dataset of seventeen centrifuge capillary pressure measurements were used in this research to demonstrate the merit of the new method. For both sandstones and carbonates, the logarithm of capillary pressure showed a linear relationship with saturation as observed by magnetic resonance imaging centrifuge capillary pressure measurements over a wide saturation range. This work demonstrates that: (a) in semi-log plots of capillary pressure as a function of saturation, capillary pressurewill vary linearly over a wide saturation range, (b) such a plot as described in (a) will show the uni-or bimodal pore size distribution of the rock, (c) the exponential capillary pressure function simplifies analytical modelsthat use the capillary pressure function, for example oil recovery models for fractured reservoirs.
Methane hydrate deposits world-wide are promising sources of natural gas. Magnetic Resonance Imaging (MRI) has proven useful in previous studies of hydrate formation. In the present work, methane ...hydrate formation in a water saturated sand pack was investigated employing an MRI-compatible metallic core holder at low magnetic field with a suite of advanced MRI methods developed at the UNB MRI Centre. The new MRI methods are intended to permit observation and quantification of residual fluids in the pore space as hydrate forms. Hydrate formation occurred in the water-saturated sand at 1500 psi and 4 °C. The core holder has a maximum working pressure of 4000 psi between -28 and 80 °C. The heat-exchange jacket enclosing the core holder enabled very precise control of the sample temperature. A pure phase encode MRI technique, SPRITE, and a bulk T1-T2 MR method provided high quality measurements of pore fluid saturation. Rapid 1D SPRITE MRI measurements time resolved the disappearance of pore water and hence the growth of hydrate in the sand pack. 3D π-EPI images confirmed that the residual water was inhomogeneously distributed along the sand pack. Bulk T1-T2 measurements discriminated residual water from the pore gas during the hydrate formation. A recently published local T1-T2 method helped discriminate bulk gas from the residual fluids in the sample. Hydrate formation commenced within two hours of gas supply. Hydrate formed throughout the sand pack, but maximum hydrate was observed at the interface between the gas pressure head and the sand pack. This irregular pattern of hydrate formation became more uniform over 24 hours. The rate of hydrate formation was greatest in the first two hours of reaction. An SE-SPI T2 map showed the T2 distribution changed considerably in space and time as hydrate formation continued. Changes in the T2 distribution are interpreted as pore level changes in residual water content and environment.
Multi-exponential decay is prevalent in magnetic resonance spectroscopy, relaxation, and imaging. This paper describes simple MATLAB and Python functions and scripts for regularized multi-exponential ...analysis methods for 1D and 2D data and example test problems and experiments. Regularized least-squares solutions provide production-quality outputs with robust stopping rules in ~5 and ~20 lines of code for 1D and 2D inversions, respectively. The software provides an open-architecture simple solution for transforming exponential decay data to the distribution of their decay lifetimes. Examples from magnetic resonance relaxation of a complex fluid, a Danish North Sea crude oil, and fluid mixtures in porous materials-brine/crude oil mixture in North Sea reservoir chalk-are presented. Developed codes may be incorporated in other software or directly used by other researchers, in magnetic resonance relaxation, diffusion, and imaging or other physical phenomena that require multi-exponential analysis.
The majority of low field Magnetic Resonance (MR) analyses rely on T
lifetime measurements. Modification of the T
measurement to include a T
dimension has made the T
-T
measurement a very powerful ...analytical technique. The T
-T
measurement is uniquely well suited to characterization of different spin populations in porous materials, such as fluid bearing reservoir rocks, and in soft biopolymer materials, for example foods. However, the T
-T
measurement is challenging or impossible if the T
relaxation lifetime, or a component lifetime, is short-lived and the signal unobservable in an echo measurement. This occurs in many important classes of materials. A short lifetime T
will not however, in general, preclude observation of a free induction decay with signal decay governed by T
*. As outlined in this paper a T
-T
* measurement is a useful analog to the T
-T
experiment. T
-T
* measurement enables one to differentiate species as a function of T
* in one dimension and T
in the other dimension. Monitoring changes of the T
-T
* coordinate, and associated signal intensity changes, has the potential to reveal structural changes in materials evolving in time. These methods may also be employed to discriminate and identify solid-like species present in static samples. The T
-T
* measurement is very general in application, but in this paper we focus on cement-based mortars to develop and illustrate the essential ideas. T
-T
* results show a multi-modal behaviour of the MR signal lifetimes, T
and T
*, in mortar samples under study, indicating at least two different water populations. The short T
* lifetime was assigned to interlayer water (water between C-S-H layers) where the associated T
is also short lived. The longer T
* lifetime was assigned to water in the pore space, where T
is also longer lived. In addition to mortar samples we also show application of the method to a crystalline organic species, o-phenylenediamine, which features Sinc Gaussian and exponential decays of transverse magnetization.