Fast spin echo MRI of reservoir core plugs with a variable field magnet.

Fast Spin Echo MRI is now widely employed in biomedicine for proton density and T2 contrast imaging. Fast Spin Echo methods provide rapid data acquisition by employing multiple echoes to determine multiple k-space lines with single excitations. Due to the multi-exponential behavior of T2 in typical porous media, and the strong dependence of T2 on the details of the experiment, acquiring a proton density image with Fast Spin Echo methods requires favorable sample and acquisition parameters. In recent years, we have shown the value of pure phase encode Free Induction Decay based methods such as SPRITE. However, in a reservoir rock, a typical T2* is hundreds of µs, whereas a typical T2 is hundreds of ms. Hence, there is merit in considering spin echo-based MRI measurements such as the Fast Spin Echo for rock core plug studies. A variable field superconducting magnet was employed in this study. This is a new class of magnet for MR/MRI. These magnets have the flexibility of operation in the field range of 0.01 Tesla to 3 Tesla. This is advantageous when working with rock core plugs, as it allows one to maximize sample magnetization, by increasing the static field while controlling magnetic susceptibility mismatch effects, and thereby T2 and T2*, through reducing the static field. The magnetic fields employed in the study were 0.79, 1.5, and 3 Tesla. Measurements were undertaken on five brine-saturated reservoir rock core plugs (Bentheimer, Berea, Buff Berea, Nugget, and Wallace). The results show that Fast Spin Echo measurements are more sensitive than SPRITE methods in amenable samples and usually feature higher resolution. Quantification of saturation with Fast Spin Echo methods requires correction for T2 attenuation. The results also show that 3 Tesla is too high a static field in general for rock core MRI studies with either method. While the current study is focused on five representative reservoir rock cores, the conclusions which result are general for MRI of fluids in porous media.

T1-T2* relaxation correlation measurements.

The majority of low field Magnetic Resonance (MR) analyses rely on T2 lifetime measurements. Modification of the T2 measurement to include a T1 dimension has made the T1-T2 measurement a very powerful analytical technique. The T1-T2 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 T1-T2 measurement is challenging or impossible if the T2 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 T2 will not however, in general, preclude observation of a free induction decay with signal decay governed by T2*. As outlined in this paper a T1-T2* measurement is a useful analog to the T1-T2 experiment. T1-T2* measurement enables one to differentiate species as a function of T2* in one dimension and T1 in the other dimension. Monitoring changes of the T1-T2* 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 T1-T2* measurement is very general in application, but in this paper we focus on cement-based mortars to develop and illustrate the essential ideas. T1-T2* results show a multi-modal behaviour of the MR signal lifetimes, T1 and T2*, in mortar samples under study, indicating at least two different water populations. The short T2* lifetime was assigned to interlayer water (water between C-S-H layers) where the associated T1 is also short lived. The longer T2* lifetime was assigned to water in the pore space, where T1 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.

Mapping three-dimensional oil distribution with π-EPI MRI measurements at low magnetic field.

Magnetic resonance imaging (MRI) is a robust tool to image oil saturation distribution in rock cores during oil displacement processes. However, a lengthy measurement time for 3D measurements at low magnetic field can hinder monitoring the displacement. 1D and 2D MRI measurements are instead often undertaken to monitor the oil displacement since they are faster. However, 1D and 2D images may not completely reflect the oil distribution in heterogeneous rock cores. In this work, a high-speed 3D MRI technique, π Echo Planar Imaging (π-EPI), was employed at 0.2T to monitor oil displacement. Centric scan interleaved sampling with view sharing in k-t space was employed to improve the temporal resolution of the π-EPI measurements. A D2O brine was employed to distinguish the hydrocarbon and water phases. A relatively homogenous glass bead pack and a heterogeneous Spynie core plug were employed to show different oil displacement behaviors. High quality 3D images were acquired with π-EPI MRI measurements. Fluid quantification with π-EPI compared favorably with FID, CPMG, 1D-DHK-SPRITE, 3D Fast Spin Echo (FSE) and 3D Conical SPRITE measurements. π-EPI greatly reduced the gradient duty cycle and improved sensitivity, compared to FSE and Conical SPRITE measurements, enabling dynamic monitoring of oil displacement processes. For core plug samples with sufficiently long lived T2, T2(∗), π-EPI is an ideal method for rapid 3D saturation imaging.

FID-SPI pulse sequence for quantitative MRI of fluids in porous media.

MRI has great potential for providing quantitative, spatially resolved information about fluids imbibed in porous media. The pure phase encode SPRITE technique has proven to be a very general method for the generation of density images in porous media; however, low flip-angle RF pulses and broad filter widths, required by short encoding times, yield sub-optimal S/N images. A 1-D phase-encoding sequence for T2(∗) mapping, named FID-SPI, is presented and analyzed in terms of image quality and accuracy of fluid content distribution in porous media. Extension to 2-D and 3-D imaging was straightforward and images of heterogeneous samples are presented. The FID-SPI measurement results in a series of individual T2(∗) weighted images acquired following RF excitation and pulsed phase-encoding gradients. Key to the performance of the FID-SPI method is high quality control of the magnetic field gradient pulse to ensure each FID point has identical spatial encoding. FID-SPI is intended for a quantitative determination of the spatially resolved fluid content in heterogeneous porous media, having the ability to determine the T2(∗) decay for each image pixel. T2(∗) mapping aids in estimation of the local fluid content.

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.

Determination of spatially-resolved porosity, tracer distributions and diffusion coefficients in porous media using MRI measurements and numerical simulations.

This work is focused on measuring the concentration distribution of a conservative tracer in a homogeneous synthetic porous material and in heterogeneous natural sandstone using MRI techniques, and on the use of spatially resolved porosity data to define spatially variable diffusion coefficients in heterogeneous media. The measurements are made by employing SPRITE, a fast MRI method that yields quantitative, spatially-resolved tracer concentrations in porous media. Diffusion experiments involving the migration of H(2)O into D(2)O-saturated porous media are conducted. One-dimensional spatial distributions of H(2)O-tracer concentrations acquired from experiments with the homogeneous synthetic calcium silicate are fitted with the one-dimensional analytical solution of Fick's second law to confirm that the experimental method provides results that are consistent with expectations for Fickian diffusion in porous media. The MRI-measured concentration profiles match well with the solution for Fick's second law and provide a pore-water diffusion coefficient of 1.75×10(-9)m(2)s(-1). The experimental approach was then extended to evaluate diffusion in a heterogeneous natural sandstone in three dimensions. The relatively high hydraulic conductivity of the sandstone, and the contrast in fluid density between the H(2)O tracer and the D(2)O pore fluid, lead to solute transport by a combination of diffusion and density-driven advection. The MRI measurements of spatially distributed tracer concentration, combined with numerical simulations allow for the identification of the respective influences of advection and diffusion. The experimental data are interpreted with the aid of MIN3P-D - a multicomponent reactive transport code that includes the coupled processes of diffusion and density-driven advection. The model defines local diffusion coefficients as a function of spatially resolved porosity measurements. The D(e) values calculated for the heterogeneous sandstone and used to simulate diffusive and advective transport range from 5.4×10(-12) to 1.0×10(-10)m(2)s(-1). These methods have broad applicability to studies of contaminant migration in geological materials.

Determining the geometry of strongly hydrogen-bonded silanols in a layered hydrous silicate by solid-state nuclear magnetic resonance.

High-resolution solid-state NMR spectroscopy is exploited to obtain structural constraints involving strongly hydrogen-bonded silanols in octosilicate, a prominent member of the layered hydrous sodium silicates. Proton-silicon cross-polarization dynamics reveals that octosilicate contains two types of Q(3) silicons present in hydrogen-bonded -Si-O-Hcdots, three dots, centeredO-Si- and -Si-O(-)-type sites which can only be distinguished by their different abilities to cross polarize and the different mobilities of neighboring hydrous species. The theoretical analysis of the oscillating components of the polarization transfer buildup curves suggests that the model of heteronuclear pairs is an adequate description of the quantum spin system within hydrogen-bonded -Si-O-Hcdots, three dots, centeredO-Si- fragments. We also show that dipolar modulated, slow speed magic-angle (29)Si NMR spectrum provides unique geometric information on strongly hydrogen-bonded silanols. The dipolar modulated spinning sidebands contain all the information necessary to determine the internuclear Sicdots, three dots, centeredH distances as well as the magnitude and orientation of the principal elements of the (29)Si chemical shielding tensor in the molecular frame. The data provide definite proof of the intralayer character of strongly hydrogen-bonded silanol groups in a bridging, albeit not symmetric, position between neighboring tetrahedra. The approach developed in this work may be useful to obtain structural information on related layered alkali metal silicates, silica gels as well as on other classes of microporous materials.