Two-dimensional T2 distribution mapping in porous solids with phase encode MRI

Controlling susceptibility mismatch effects, signal lifetimes, and SNR through variation of B0 in MRI of rock core plugs

¹H relaxometry measurements of petroleum core plugs are commonly performed on low field magnets (<0.5 Tesla) to reduce the influence of magnetic susceptibility mismatch on measurements of the spin-spin relaxation time, T₂. The Signal to Noise Ratio (SNR) of the MR signal, however, generally decreases with lower magnetic fields. Higher magnetic fields (>3 Tesla) are typically employed in small animal MRI studies to improve SNR and image resolution. For many rock core plug samples, susceptibility mismatch effects can be severe at these higher fields leading to decreased T₂ and T₂*. In this work we seek an answer to the general question of what is the best field for MRI of rock core plugs, anticipating that it will be both sample and measurement method dependent. Free Induction Decay (FID) relaxation time measurements were undertaken to investigate the conditions under which the SNR in Centric Scan SPRITE (Single Point Ramped Imaging with T₁ Enhancement) MRI measurements is maximized. The image SNR benefits from greater signal at higher fields, but is negatively impacted by the correspondingly shorter signal lifetimes. Depending on the noise regime of the sample, the maximum SNR may be predicted for Centric Scan SPRITE MRI with T₂* being B₀ field dependent. In this work we describe a series of simple experimental considerations to determine the optimal B₀ field for SPRITE MRI. Selection of the best field is aided by a new generation of superconducting magnets which allows the experimentalist to readily vary the field strength. Such magnets allow one to experimentally control sample magnetization for high sensitivity MRI measurements of core plug samples, while controlling the effect of susceptibility mismatch on the signal lifetimes.

Spatially-resolved 1H NMR relaxation-exchange measurements in heterogeneous media

In the last decades, the ¹H NMR T₂-T₂ relaxation-exchange (REXSY) technique has become an essential tool for the molecular investigation of simple and complex fluids in heterogeneous porous solids and soft matter, where the mixing-time-evolution of cross-correlated T₂-T₂ peaks enables a quantitative study of diffusive exchange kinetics in multi-component systems. Here, we present a spatially-resolved implementation of the T₂-T₂ correlation technique, named z-T₂-T₂, based on one-dimensional spatial mapping along z using a rapid frequency-encode imaging scheme. Compared to other phase-encoding methods, the adopted MRI technique has two distinct advantages: (i) is has the same experimental duration of a standard (bulk) T₂-T₂ measurement, and (ii) it provides a high spatial resolution. The proposed z-T₂-T₂ method is first validated against bulk T₂-T₂ measurements on homogeneous phantom consisting of cyclohexane uniformly imbibed in finely-sized α-Al₂O₃ particles at a spatial resolution of 0.47 mm; thereafter, its performance is demonstrated, on a layered bed of multi-sized α-Al₂O₃ particles, for revealing spatially-dependent molecular exchange kinetics properties of intra- and inter-particle cyclohexane as a function of particle size. It is found that localised z-T₂-T₂ spectra provide well resolved cross peaks whilst such resolution is lost in standard bulk T₂-T₂ data. Future prospective applications of the method lie, in particular, in the local characterisation of mass transport phenomena in multi-component porous media, such as rock cores and heterogeneous catalysts.

Accelerating the estimation of 3D spatially resolved T2 distributions

Obtaining quantitative, 3D spatially-resolved T₂ distributions (T₂ maps) from magnetic resonance data is of importance in both medical and porous media applications. Due to the long acquisition time, there is considerable interest in accelerating the experiments by applying undersampling schemes during the acquisition and developing reconstruction techniques for obtaining the 3D T₂ maps from the undersampled data. A multi-echo spin echo pulse sequence is used in this work to acquire the undersampled data according to two different sampling patterns: a conventional coherent sampling pattern where the same set of lines in k-space is sampled for all equally-spaced echoes in the echo train, and a proposed incoherent sampling pattern where an independent set of k-space lines is sampled for each echo. The conventional reconstruction technique of total variation regularization is compared to the more recent techniques of nuclear norm regularization and Nuclear Total Generalized Variation (NTGV) regularization. It is shown that best reconstructions are obtained when the data acquired using an incoherent sampling scheme are processed using NTGV regularization. Using an incoherent sampling pattern and NTGV regularization as the reconstruction technique, quantitative results are obtained at sampling percentages as low as 3.1% of k-space, corresponding to a 32-fold decrease in the acquisition time, compared to a fully sampled dataset.

In Situ Chemically-Selective Monitoring of Multiphase Displacement Processes in a Carbonate Rock Using 3D Magnetic Resonance Imaging

Accurate monitoring of multiphase displacement processes is essential for the development, validation and benchmarking of numerical models used for reservoir simulation and for asset characterization. Here we demonstrate the first application of a chemically-selective 3D magnetic resonance imaging (MRI) technique which provides high-temporal resolution, quantitative, spatially resolved information of oil and water saturations during a dynamic imbibition core flood experiment in an Estaillades carbonate rock. Firstly, the relative saturations of dodecane (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$S_{\mathrm{o}})$$\end{document}So) and water (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$S_{\mathrm{w}})$$\end{document}Sw), as determined from the MRI measurements, have been benchmarked against those obtained from nuclear magnetic resonance (NMR) spectroscopy and volumetric analysis of the core flood effluent. Excellent agreement between both the NMR and MRI determinations of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$S_{\mathrm{o}}$$\end{document}So and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$S_{\mathrm{w}}$$\end{document}Sw was obtained. These values were in agreement to 4 and 9% of the values determined by volumetric analysis, with absolute errors in the measurement of saturation determined by NMR and MRI being 0.04 or less over the range of relative saturations investigated. The chemically-selective 3D MRI method was subsequently applied to monitor the displacement of dodecane in the core plug sample by water under continuous flow conditions at an interstitial velocity of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$1.27\times 10^{-6}\,\hbox {m}\,\hbox {s}^{-1}$$\end{document}1.27×10-6ms-1 (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$0.4\,\hbox {ft}\,\hbox {day}^{-1})$$\end{document}0.4ftday-1). During the core flood, independent images of water and oil distributions within the rock core plug at a spatial resolution of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$0.31\,\hbox {mm}\times 0.39\,\hbox {mm} \times 0.39\,\hbox {mm}$$\end{document}0.31mm×0.39mm×0.39mm were acquired on a timescale of 16 min per image. Using this technique the spatial and temporal dynamics of the displacement process have been monitored. This MRI technique will provide insights to structure–transport relationships associated with multiphase displacement processes in complex porous materials, such as those encountered in petrophysics research.

Fast imaging of laboratory core floods using 3D compressed sensing RARE MRI

Three-dimensional (3D) imaging of the fluid distributions within the rock is essential to enable the unambiguous interpretation of core flooding data. Magnetic resonance imaging (MRI) has been widely used to image fluid saturation in rock cores; however, conventional acquisition strategies are typically too slow to capture the dynamic nature of the displacement processes that are of interest. Using Compressed Sensing (CS), it is possible to reconstruct a near-perfect image from significantly fewer measurements than was previously thought necessary, and this can result in a significant reduction in the image acquisition times. In the present study, a method using the Rapid Acquisition with Relaxation Enhancement (RARE) pulse sequence with CS to provide 3D images of the fluid saturation in rock core samples during laboratory core floods is demonstrated. An objective method using image quality metrics for the determination of the most suitable regularisation functional to be used in the CS reconstructions is reported. It is shown that for the present application, Total Variation outperforms the Haar and Daubechies3 wavelet families in terms of the agreement of their respective CS reconstructions with a fully-sampled reference image. Using the CS-RARE approach, 3D images of the fluid saturation in the rock core have been acquired in 16min. The CS-RARE technique has been applied to image the residual water saturation in the rock during a water-water displacement core flood. With a flow rate corresponding to an interstitial velocity of vi=1.89±0.03ftday(-1), 0.1 pore volumes were injected over the course of each image acquisition, a four-fold reduction when compared to a fully-sampled RARE acquisition. Finally, the 3D CS-RARE technique has been used to image the drainage of dodecane into the water-saturated rock in which the dynamics of the coalescence of discrete clusters of the non-wetting phase are clearly observed. The enhancement in the temporal resolution that has been achieved using the CS-RARE approach enables dynamic transport processes pertinent to laboratory core floods to be investigated in 3D on a time-scale and with a spatial resolution that, until now, has not been possible.

Contributed review: nuclear magnetic resonance core analysis at 0.3 T

Nuclear magnetic resonance (NMR) provides a powerful toolbox for petrophysical characterization of reservoir core plugs and fluids in the laboratory. Previously, there has been considerable focus on low field magnet technology for well log calibration. Now there is renewed interest in the study of reservoir samples using stronger magnets to complement these standard NMR measurements. Here, the capabilities of an imaging magnet with a field strength of 0.3 T (corresponding to 12.9 MHz for proton) are reviewed in the context of reservoir core analysis. Quantitative estimates of porosity (saturation) and pore size distributions are obtained under favorable conditions (e.g., in carbonates), with the added advantage of multidimensional imaging, detection of lower gyromagnetic ratio nuclei, and short probe recovery times that make the system suitable for shale studies. Intermediate field instruments provide quantitative porosity maps of rock plugs that cannot be obtained using high field medical scanners due to the field-dependent susceptibility contrast in the porous medium. Example data are presented that highlight the potential applications of an intermediate field imaging instrument as a complement to low field instruments in core analysis and for materials science studies in general.

Spatially resolved measurements of mean spin-spin relaxation time constants

Magnetic Resonance measurements of the T2 distribution have become very common and they are a powerful way to probe microporous fluid bearing solids. While the structure of the T2 distribution, and changes in the structure, are often very informative, it is common to reduce the T2 distribution to a mean numeric quantity in order to provide a quantitative interpretation of the distribution. Magnetic Resonance Imaging measurements of the T2 distribution have recently been introduced, but they are time consuming, especially for 2 and 3 spatial dimensions. In this paper we explore a direct MRI measurement of the arithmetic mean of 1/T2, characterizing the distribution by using the initial slope of the spatially resolved T2 decay in a CPMG prepared Centric Scan SPRITE experiment. The methodology is explored with a test phantom sample and realistic petroleum reservoir core plug samples. The arithmetic mean of 1/T2 is related to the harmonic mean of T2. The mean obtained from the early decay is explored through measurements of uniform saturated core plug samples and by comparison to other means determined from the complete T2 distribution. Complementary data were obtained using SE-SPI T2 distribution MRI measurements. The utility of the arithmetic mean 1/T2 is explored through measurements of centrifuged core plug samples where the T2 distribution varies spatially. The harmonic mean T2 obtained from the early decay was employed to estimate the irreducible water saturation for core plug samples.

Low-field permanent magnets for industrial process and quality control

In this review we focus on the technology associated with low-field NMR. We present the current state-of-the-art in low-field NMR hardware and experiments, considering general magnet designs, rf performance, data processing and interpretation. We provide guidance on obtaining the optimum results from these instruments, along with an introduction for those new to low-field NMR. The applications of lowfield NMR are now many and diverse. Furthermore, niche applications have spawned unique magnet designs to accommodate the extremes of operating environment or sample geometry. Trying to capture all the applications, methods, and hardware encompassed by low-field NMR would be a daunting task and likely of little interest to researchers or industrialists working in specific subject areas. Instead we discuss only a few applications to highlight uses of the hardware and experiments in an industrial environment. For details on more particular methods and applications, we provide citations to specialized review articles.

B(1) mapping with a pure phase encode approach: quantitative density profiling

In MRI, it is frequently observed that naturally uniform samples do not have uniform image intensities. In many cases this non-uniform image intensity is due to an inhomogeneous B1 field. The 'principle of reciprocity' states that the received signal is proportional to the local magnitude of the applied B1 field per unit current. Inhomogeneity in the B1 field results in signal intensity variations that limit the ability of MRI to yield quantitative information. In this paper a novel method is described for mapping B1 inhomogeneities based on measurement of the B1 field employing centric-scan pure phase encode MRI measurements. The resultant B1 map may be employed to correct related non-uniformities in MR images. The new method is based on acquiring successive images with systematically incremented low flip angle excitation pulses. The local image intensity variation is proportional to B1(2), which ensures high sensitivity to B1 field variations. Pure phase encoding ensures the resultant B1 field maps are free from distortions caused by susceptibility variation, chemical shift and paramagnetic impurities. Hence, the method works well in regions of space that are not accessible to other methods such as in the vicinity of conductive metallic structures, such as the RF probe itself. Quantitative density images result when the centric scan pure phase encode measurement is corrected with a relative or absolute B1 field map. The new technique is simple, reliable and robust.

Restricted k-space sampling in pure phase encode MRI of rock core plugs

In the study of rock core plugs with multidimensional MRI, the samples are of a regular cylindrical shape that yields well defined intensity distributions in reciprocal space. The high intensity k-space points are concentrated in the central region and in specific peripheral regions. A large proportion of the k-space points have signal intensities that are below the noise level. These points can be zero-filled instead of being collected experimentally. k-space sampling patterns that collect regions of high intensity signal while neglecting low intensity regions can be naturally applied to a wide variety of pure phase encoding measurements, such as T2 mapping SESPI, hybrid-SESPI and SPRITE, since all imaging dimensions can be under-sampled. With a shorter acquisition time, as fewer experimental data points are required, the RF and gradient duty cycles are reduced, while the image SNR is improved.

Two-dimensional T2 distribution mapping in rock core plugs with optimal k-space sampling

Spin-echo single point imaging has been employed for 1D T(2) distribution mapping, but a simple extension to 2D is challenging since the time increase is n fold, where n is the number of pixels in the second dimension. Nevertheless 2D T(2) mapping in fluid saturated rock core plugs is highly desirable because the bedding plane structure in rocks often results in different pore properties within the sample. The acquisition time can be improved by undersampling k-space. The cylindrical shape of rock core plugs yields well defined intensity distributions in k-space that may be efficiently determined by new k-space sampling patterns that are developed in this work. These patterns acquire 22.2% and 11.7% of the k-space data points. Companion density images may be employed, in a keyhole imaging sense, to improve image quality. T(2) weighted images are fit to extract T(2) distributions, pixel by pixel, employing an inverse Laplace transform. Images reconstructed with compressed sensing, with similar acceleration factors, are also presented. The results show that restricted k-space sampling, in this application, provides high quality results.