Spatially resolved dynamic longitudinal relaxometry in single-sided NMR

doi:10.21203/rs.3.rs-3104684/v1

Download PDF

Research Article

Spatially resolved dynamic longitudinal relaxometry in single-sided NMR

https://doi.org/10.21203/rs.3.rs-3104684/v1

This work is licensed under a CC BY 4.0 License

Journal Publication

published 01 Aug, 2023

Read the published version in Applied Magnetic Resonance →

You are reading this latest preprint version

Relaxation in nuclear magnetic resonance (NMR), both transverse and longitudinal, provide information on microscopic features of a wide variety of systems and may be used to monitor dynamic processes such as cementation, chemical reactions, gelatinization, and evaporation. Dynamic relaxometry, in combination with spatial resolution, is a useful technique that provides deep insight into complex systems evolution. In this work, we explore the range of applicability of single-sided NMR to determine the evaporation kinetics of fluid from porous media. We show that, due to technical experimental restrictions, the determination of the time-dependent amount of fluid in different voids as a function of the position is in general not feasible with transverse relaxation experiments. However, as opposed to common intuition, longitudinal relaxation experiments provide reliable and fast acquisition, compatible with the requirements needed to monitor a water evaporation process from a model oil-reservoir rock sample.

Profiles

Single-Sided NMR

Relaxation times

Porous Media

Evaporation Kinetics

Porous materials with voids that may be filled with fluids, liquid, or gas, are present in a great variety of systems, such as soils, concrete, sponges, foodstuff, or oil-bearing rocks among others. Some of the most relevant parameters that define the exploitability of a porous media are the pore sizes, the number of pores, their connectivity, or the ability of a fluid to flow through the medium. The geometrical structure can be precisely determined by a variety of techniques such as gas sorption-desorption, mercury porosimetry, X-rays, neutron dispersion, and electronic microscopy. However, to fully understand the performance of a porous medium, the confined fluid dynamics must be determined. In this regard, nuclear magnetic resonance (NMR) stands out as it directly detects the spin-bearing particles from the gas or liquid, without the need for tracer particles that may alter the flow or diffusion. It can detect fluids non-invasively in opaque media and may provide spatial and microscopic information. Spatial information is often obtained by magnetic resonance imaging (MRI) whilst microscopic information is obtained by the determination either of restricted diffusion coefficients or relaxation times. Transverse and longitudinal relaxation times will probe the system geometry due to magnetization loss by liquid/surface interactions [1, 2]. For systems changing with time, for instance, during an evaporation process, the term dynamic relaxometry is often used [3]. A great variety of systems have been studied by probing changes in relaxation times, for instance, the cementation process [4–7], porosity, spore size distributions, and permeability in rock from oil wells among others [8–11]. This approach can also be combined with MRI, often in one dimension. In this way not only the total amount of fluid in a slice can be determined, but also its distribution between different environments [12–15].

Single-sided NMR englobes a family of devices where the magnetic field used for NMR purposes ex-situ, with a strong magnetic field gradient for the main magnetic field as for the radiofrequency field. Examples are well-logging tools [16], the stray field of superconducting magnets [17, 18] or tailored designed benchtop equipment [19, 20] as the MOLE (Mobile Lateral Explorer) [21, 22] and the most developed by the industry, the NMR-MOUSE (Mobile Universal Surface Explorer) [23–25]. The NMR-MOUSE presents the advantage of mobility and non-stringent sample sizes, being the penetration depth as the main geometrical limitation. The technique has found a great variety of applications in the study of solid materials, polymers, cultural heritage [26–28], and liquids in oil rocks [29, 30]. Signals are acquired as a function of the penetration depth, namely a 1D image or more commonly referred to as a profile. As the signals are detected with a train of echoes, different relaxation times will reflect changes in polymer cross-linking, fluid-restricted mobility, or changes in material composition. By acquiring multiple echoes and selecting different time intervals for signal addition, the contrast between different components is readily achieved [31]. The main drawback with this strategy arises when the detected signals have relaxation times on the order of hundreds of milliseconds. As the influence of diffusion is in general to be avoided, short echo times are used leading to the application of thousands of radiofrequency pulses. This induces heating of the excitation/detection coil, which will eventually detune the resonant circuit, decreasing the performance of the system.

In this work, we analyze the limitation imposed by coil heating during the acquisition of transverse relaxation for liquids confined in porous media. We find that correctly determining the distribution of fluids in different void sizes requires an experimental time that may very well exceed the system's characteristic time in an evaporation process. We then analyze the approach of using longitudinal relaxation times to probe the fluid distribution. This counterintuitive approach is based on the boost in signal-to-noise ratio that may be achieved by adding several echoes during a multi-pulse sequence used as a detection block. We have recently applied this approach to show that diffusion-relaxation maps may be acquired in a shorter time using D-T₁ rather than D-T₂ [32]. We determine a set of parameters to be applied to a particular configuration of a PM25 NMR-MOUSE developed by Magritek GmbH and apply these conditions to the determination of spatial and microscopic information during the evaporation process of water from a model oil rock.

2.1 Profiles with transverse relaxation times

Experiments were carried out in a single-sided NMR device PM25, from Magritek GmbH, which has a static gradient of 7 T/m and operates at a frequency of 12.99 MHz for ¹H, with detection volume 26.7 mm away from the magnets, and a sensitive area of approximately (4 × 4) cm². The sample is positioned over a flat holder and the magnet position is moved with a precision lift, in this way, control of the position of the selective slice is achieved. The system may be set to five different configurations setting the excitation/detection coil further from the magnet and closer to the sample, defining different penetration depths ranging from 5.4 to 25 mm. In this work, without loss of generality, we restrict our discussion to a penetration depth of 10.6 mm. The length of the radiofrequency (RF) pulses was of ${t}_{RF}=$9.5 µs. The excited bandwidth for this configuration defines a maximum slice thickness of 350 µm [25, 31]. Due to the large static gradient, the magnetization is rapidly dephased after a 90° excitation pulse, in a time shorter than the typical receiver dead times, therefore the free induction decay (FID) cannot be acquired. This leads to the use of multipulse refocusing sequences that capture the evolution of multiple echoes. For liquids, a Carr-Purcell-Meiboom-Gill (CPMG)[33, 34] pulse sequence is generally used to detect the magnetization evolution.

For the acquisition of proton densities as a function of the sample height, namely a profile, the addition of several echoes is used to enhance the signal-to-noise ratio (SNR). However, if the sample consists of different proton pools with varying transverse relaxation times, the acquisition of the whole echo train will provide information on the different proton pools as a function of the sensor´s position. Relaxation rates may be determined by fitting with multiple exponential decay functions or, more commonly, via a numerical data inversion, the so-called Inverse Laplace Transform (ILT) where the most favored on-dimensional algorithm is CONTIN [35]. The data acquired with a CPMG pulse sequence depend on molecular diffusion which, in the presence of the strong magnetic field gradient, will render an apparent relaxation time, T_2D. The magnetization of a diffusing spin-bearing molecule in the presence of a constant magnetic field gradient during a CPMG sequence with echo time ${t}_{E}$ at the top of the $m-th$ echo, decays as [36–38]:

$$\text{ln}\left(\frac{S\left(m{t}_{E}\right)}{{S}_{0}}\right)=-\left(\frac{1}{{T}_{2S}}+\frac{1}{12}{\left(\gamma G{t}_{E}\right)}^{2}D\right)m{t}_{E}=-\left(\frac{1}{{T}_{2D}}\right)m{t}_{E},$$

where $\gamma$is the gyromagnetic ratio of ¹H and $D$ the diffusion coefficient. ${T}_{2D}$ is the apparent relaxation time driven both by diffusion and surface relaxation, ${T}_{2S},$ which, in the surface-limited relaxation regime is expressed as [1, 5, 8]:

$$\frac{1}{{T}_{2S}}= \frac{1}{{T}_{2B}}+\frac{6{\rho }_{2}}{d}$$

where the bulk relaxation time of the fluid is ${T}_{2B}$, ${\rho }_{2}$ is the surface relaxivity and $d$ is the pore diameter. To reduce as much as possible the influence of diffusion in Eq. (1), echo times must be set as short as possible. This implies that thousands of radiofrequency pulses must be applied, leading to the heating of the surface coil. The change in the coil’s resistance produces a shift of the signal’s phase and eventually a pulse drop. A straightforward solution to avoid coil heating beyond a critical level is to set the recycle delay between experiments to longer times than those usually applied in NMR, namely 3–5 T₁. For the PM25 NMR MOUSE, considering the pulses of the CPMG as an on-period, and the recycle delay ($rd$) as an off-period, a duty cycle of 1.6 provides is good enough performance, with less than 3° dephasing. This sets the waiting time between experiments as $rd={t}_{RF}*NE/1.6$, where $NE$ is the total number of echoes to be acquired. It is then clear the total experimental time set by the relaxation time ${T}_{2D}$, is defined by the diffusion coefficient and the echo time.

A second factor to consider is the spatial resolution of the profile. As all radiofrequency pulses are applied in the presence of the magnetic field gradient, the length of the pulse can be used to select a desired slice thickness. However, multiple points can be acquired at the top of each echo and the acquisition time can be set to acquire the desired field of view of a rectangular slice in the direction of the gradient, or in this case the desired resolution (Res) [25, 39]. Setting the acquisition bandwidth constant, in this case, with a dwell-time of $dw=$ 0.5 µs, the number of points to acquire in each echo are:

$$n= \text{round}\left(\frac{1}{G\bullet dw\bullet \text{Res}}\right)$$

For the used configuration this corresponds to $n=$ 19 for a resolution of 350 µm and defines a minimum echo time of ${t}_{E}=$ 67.5 µs. As the slice thickness decreases, the number of acquired data points will increase. For instance, for a slice of 100 µm, $n=$ 67 and ${t}_{E}=$ 91 µs. Figure 1A shows the transverse relaxation time values calculated with Eq. (1) as a function of the resolution for an arbitrary value of ${T}_{2S}$ = 500 ms and $D=$ 1.10⁻⁹ m²/s. As the slice becomes thicker echo times are shorter, therefore the effective relaxation time increases its value. If the acquisition is to ${t}_{acq}=3{T}_{2D}$ the number of echoes increases from c.a. 4000 to 8500 for resolutions of 100 µm and 350 µm respectively. Therefore, the recycle delay for each of these resolutions changes from 24 s to 50 s respectively. Considering a minimum of 4 phase cycled scans, which are required to remove unwanted coherence pathways [40–43], the total experimental time calculated to acquire a single profile is shown in Fig. 1B. The number of slices required to map the complete proton density of a 1 cm sample decreases with increasing slice thickness. Even though relaxation times for thinner slices are shorter, the experimental time is much longer than that required for thicker slices. However, for the worst resolution achievable in the configuration used with the NMR-MOUSE, more than 90 hs are required for a single profile, a prohibitive time for monitoring dynamic processes.

2.2 Profiles with longitudinal relaxation times

Implementation of a saturation-recovery [44] sequence for T₁ measurements is straightforward in single-side NMR as the presence of the strong magnetic field gradient is used to destroy the magnetization at the beginning of each experiment by applying an odd number of multiple 90° pulses separated by increasing time delays [25]. A given number of echoes, NE, can be accumulated to improve the SNR and are a key feature for obtaining shorter measurement times [32]. For liquids confined in a porous media with a distribution of pore sizes, a distribution of transverse relaxation times will be present, following Eq. (1). Addition of echoes inevitably introduces a contrast between fluid confined in small and large cavities that could hinder a correct quantification of the smaller voids. A threshold value must then be set, recently we proposed a penalty of 7% in the quantification of the populations associated with shorter relaxation times [32]. As echoes are coadded, the sum of NE echoes should satisfy the relation:

$$S\left(NE\right)={S}_{0}\left[\sum _{i=1}^{NE}\text{exp}\left(-\frac{i{t}_{E}}{{T}_{2D,Sh}}\right)\right]/NE={0.93 S}_{0}$$

where ${T}_{2D,Sh}$ is the shortest transverse relaxation time that depends both on the liquid/surface interactions and on diffusion.

Figure 2A ) Number of echoes that can be added in a T ₁ measurement considering that the contrast of the shortest relaxation time, ${T}_{2D,S}$, is 93% from its zero-time value. As low resolution has a shorter echo time, a larger number of echoes can be acquired for a given relaxation time. B) Total experimental time for the acquisition of a saturation-recovery experiment for a single slice as a function of the longest relaxation time in a sample and the number of encoding points.

Figure 2 shows the number of echoes that can be added as a function of the shortest surface-defined transverse relaxation time,${T}_{2S,Sh}$, and of the selected resolution. As shown in Fig. 1A, for each value of ${T}_{2S,Sh}$ the diffusion relaxation time changes as a function of the selected resolution, this is the reason that for a given surface-defined relaxation time, the number of echoes to add decreases with increasing resolution. For instance, for ${T}_{2S,Sh}=$ 60 ms, NE = 110 for a slice thickness of 350 µm and NE = 68 for a slice thickness of 100 µm.

Once a resolution is chosen, and the number of echos is set to define the maximum SNR achievable, the encoding strategy of T₁ must be chosen. It is customary to sample the recovery curve to values of 3–5 T_1S,L, where T_1S,L stands for the longest surface-defined longitudinal relaxation time present in the sample. For the sake of speed, we use 3 T_1S,L in the present work. Several encoding data points, Np_T1 must be chosen with a compromise between numerical inversion accuracy and experimental time. Equally spaced recovery times in a logarithmic scale are considered, and the stability of the inversion procedure is discussed in the next section. Figure 2B shows the total experimental time for a single slice considering a fixed recycle delay of rd = 200 ms and 4 phase cycling steps which are needed to minimize the effect of unwanted coherence pathways. A sampling of 5 to 50 recovery times is considered. As an example, for T_1S,L = 800 ms, 93 s are required for Np_T1 = 50, while 15 s are required for Np_T1 = 5.

2.2 Sample description

Bentheimer samples purchased from Kocurek Industries Inc. were used. Disks of 3.8 cm diameter and 8.7 mm height were cut and imbibed with water by placing the submerged disks in a distiller under vacuum conditions for 48 hs. Samples were tightly fit in a sample holder that enables evaporation through the top of the sample. Measurements were carried out at ambient pressure and temperature. The porosity of the samples is 23%, with a permeability of ∼ 400 mD [11, 45, 46].

3.1 Determination of the minimal number of points for T₁ encoding

A CPMG pulse sequence with t_E = 67.5 µs was used to determine the apparent transverse relaxation times. A fully saturated sample sealed with parafilm to avoid evaporation was used. The relaxation distribution is shown in Fig. 3A, the longest T_2D appears shifted to lower values as shown in Fig. 1A in agreement with recently reported data with this same sensor’s configuration [30], where the shorter relaxation time does not change when measured with a PM25 sensor of in a homogeneous field at 3.9 MHz. A single penetration depth was used and a T₁ measurement was carried out with 50 recovery times and 4 steps of phase cycling of the radiofrequency pulses [25]. The detection of the signal was carried out by adding 50 echoes in a CPMG pulse sequence, as discussed in Fig. 2A, the maximum recovery time was set to 2500 ms. Data were inverted using a one-dimensional version of the FISTA algorithm [47]. The mean value (see Fig. 3B) of the longest T₁ present is T_1S,L = 800 ms while the mean value of the shortest one is 50 ms with a contribution of 90 % and 10 % o the toal signal respectively. By setting T_2S = 28 ms in Eq. (1) a value of T_2D = 26 ms is obtained. From the T_1S distribution, a ratio of T₁/T₂ ∼ 2 is obtained for the smaller pore. If a homogeneous relaxivity is considered throughout the sample, a value of T_2S ∼ 400 ms is expected for the longer relaxation time, giving rise to T_2D = 172 ms calculated from Eq. 1, which is in very good agreement with the data shown in Fig. 3A.

A series of experiments changing the number of recovery times was carried out. Figure 3C shows the percentual contribution of each population, obtained by taking the cumulative of the relaxation time distribution spectra, as a function of the number of recovery steps, which were equally spaced in a logarithmic scale in each case. The regularization parameter and inversion time domain were fixed for all data sets. The population associated with the shortest relaxation time is well described with more than ten recovery times. A profile experiment was then carried on to determine empirically the acquisition time between slices. A value of 30 seconds was found, which is in good agreement with the 24 s obtained from Fig. 2B, considering the additional time consumed by data transferring from the Kea2 console to the PC, and the movement of the lift to the position of the next slice.

3.2 Evaporation monitoring

Evaporation of water from sandstones has been previously reported, under both diffusion and convection-driven evaporation conditions[12]. The underlying physical processes involved and the information on the systems that may be obtained are out of the scope of this work, we will center on a description of the performance of the introduced methodology. A second fully saturated sample was placed tightly in a Teflon holder (no background signal was detected from the sample holder) and left to evaporate through the top face of the disk. Data acquisition was started immediately acquiring a total of 31 slices of 350 µm width. A profile without contrast can be obtained by plotting the total cumulative as a function of the position (black trace in Fig. 4A) or two profiles which describe the amount of water in each mean pore size (blue and red traces in Fig. 4A). Adding data points from 1.8 mm to 9.5 mm is equivalent to acquiring the signal from the whole sample on homogeneous field equipment, where the sample is fully contained within the coil. In this way the overall signal evaporation can be obtained, as shown in black dots in Fig. 4B. The initial period, c.a. one hour, corresponds to a stabilization in the temperature of the system with a low drying rate[48, 49] which is followed by a constant drying rate, depicted by the linear decay in the water content with time. This is a characteristic process of a funicular regime [50], where a hydraulic connection throughout the systems facilitates capillary flow, which moves water over to the evaporation front, which is mainly constant in position over the whole evaporation. The individual contribution from large and small pores can be extracted from the longitudinal relaxation distribution of each slice, added up to obtain the blue and red curves of Fig. 4B. Evaporation is produced only from the large pores during almost all the drying process, whilst water contained in the smaller pores does not evaporate until the latest stages, where the funicular regime gives place to a pendular regime [50]where the hydraulic connection is lost.

Figure 5A shows the total profiles obtained for each evaporation time. There is no visible receding of the evaporation front, that is, the whole sample loses water homogeneously until the last stages of evaporation, between 13–15 h, where a top-to-bottom evaporation is observed, that is, a receding evaporation front. Figure 5B shows that water contained in the larger pores evaporates homogeneously and that when the smaller pores start to become depleted (Fig. 5C) the evaporation starts from the top of the sample to the lower part, this period corresponds to a pendular regime. It is interesting to note that, within the high dispersion of the lower relaxation times (see Fig. 6B), these remain mainly constant, with a medium value of 50 ms, until the last stages of evaporation, where T_1S,Sh is reduced. This is more noticeable in the longer relaxation times shown in Fig. 6A, where at the initial stages T_1S,L ∼ 800 ms while at a time of 8 h is reduced to ∼500 ms. This is an expected behavior that corresponds to the presence of a vapor phase in the central part of the pores, therefore the fluid becomes a film of water covering the pore walls, and the available dimension for diffusion is smaller, thus water molecules sense a smaller void. It is worth noting that similar observation was carried out in porous polymer systems where spatial evaporation was determined with single-sided NMR, while the microscopic behavior was determined with a smaller sample in the homogeneous field of a 60 MHz magnet [51]. The simultaneous determination presented in this work is comparable to those obtained by combining relaxation and MRI in homogeneous magnets, with the advantage of less stringent requirements on sample geometry provided by the NMR-MOUSE. Finally, we note that the methodology is not restricted to monitoring evaporation processes but to any time-dependent process that can be discriminated by longitudinal relaxation, thus circumventing the inherent coil heating in this type of commercial devices by avoiding the acquisition of a full transverse relaxation decay [32].

We showed that the acquisition of profiles in a single-sided NMR sensor, combined with the determination of longitudinal relaxation time distributions can be performed in an experimental time short enough to monitor time-dependent phenomena, in this case, the evaporation of water from a model oil-bearing rock. The acquisition and addition of multiple echoes during a CPMG pulse sequence provides a suitable signal-to-noise ratio, without the imposition of undesirable contrast between the different proton pools. As saturation-recovery sequences are implemented straightforwardly in the presence of large magnetic field gradients, which are used as magnetization crushers, the recycle delay between acquisition depends only on the time necessary to dissipate heat from the superficial radiofrequency coil. Individual slices of 350 mm were acquired in 30 seconds, while the whole profile took 15 minutes plus the time required to reposition the sensor to the first slice, where the complete evaporation of water from the rock took around fifteen hours. Therefore, this approach is suitable for the application of dynamic relaxometry.

Ethical approval

Not applicable

Competing interests

The authors have no competing interests as defined by Springer, or other interests that might be perceived to influence the results and/or discussion reported in this paper.

Authors' contributions

Franco A. Milana: Investigation, Methodology, Software

Manuel I. Velasco: Investigation, Methodology, Software, Writing – review & editing

Gustavo A. Monti: Writing – review & editing, Funding Acquisition, Supervision

Rodolfo H. Acosta: Methodology, Funding Acquisition, Project Management, Writing – original draft, Supervision

Funding

We would like to acknowledge the financial support from CONICET (PIP-1111220130100746CO), SeCyT-UNC (33620180100154CB), and ANPCYT (PICT 2017-0957 and PICT-2019-2802).

Availability of data and materials

The datasets generated during and/or analyzed during the current study are available from the corresponding author on reasonable request.

J.P. Korb, Nuclear magnetic relaxation of liquids in porous media, New J Phys. 13 (2011). https://doi.org/10.1088/1367-2630/13/3/035016.
P.T. Callaghan, Translational Dynamics and Magnetic Resonance: Principles of Pulsed Gradient Spin Echo NMR, Oxford University Press, Oxford, 2011.
B. Maillet, R. Sidi-Boulenouar, P. Coussot, Dynamic NMR Relaxometry as a Simple Tool for Measuring Liquid Transfers and Characterizing Surface and Structure Evolution in Porous Media, Langmuir. 38 (2022) 15009–15025. https://doi.org/10.1021/acs.langmuir.2c01918.
M. Van Landeghem, J.B. D’Espinose De Lacaillerie, B. Blümich, J.P. Korb, B. Bresson, The roles of hydration and evaporation during the drying of a cement paste by localized NMR, Cem Concr Res. 48 (2013) 86–96. https://doi.org/10.1016/j.cemconres.2013.01.012.
P. McDonald, J.-P. Korb, J. Mitchell, L. Monteilhet, Surface relaxation and chemical exchange in hydrating cement pastes: A two-dimensional NMR relaxation study, Phys Rev E. 72 (2005). https://doi.org/10.1103/PhysRevE.72.011409.
P.F. Faure, S. Car, J. Magat, T. Chaussadent, Drying effect on cement paste porosity at early age observed by NMR methods, Constr Build Mater. 29 (2012) 496–503. https://doi.org/10.1016/j.conbuildmat.2011.07.012.
L. Monteilhet, J.P. Korb, J. Mitchell, P.J. McDonald, Observation of exchange of micropore water in cement pastes by two-dimensional T2 - T2 nuclear magnetic resonance relaxometry, Phys Rev E Stat Nonlin Soft Matter Phys. 74 (2006) 1–9. https://doi.org/10.1103/PhysRevE.74.061404.
Y.Q. Song, Magnetic resonance of porous media (MRPM): A perspective, Journal of Magnetic Resonance. 229 (2013) 12–24. https://doi.org/10.1016/j.jmr.2012.11.010.
Y.-Q. Song, Recent Progress of Nuclear Magnetic Resonance Applications in Sandstones and Carbonate Rocks, Vadose Zone Journal. 9 (2010) 828. https://doi.org/10.2136/vzj2009.0171.
Y. Zhang, L. Xiao, G. Liao, Y.Q. Song, Direct correlation of diffusion and pore size distributions with low field NMR, Journal of Magnetic Resonance. 269 (2016) 196–202. https://doi.org/10.1016/j.jmr.2016.06.013.
E. Lucas-Oliveira, A.G. Araujo-Ferreira, W.A. Trevizan, B. Coutinho, C. Santos, T.J. Bonagamba, Sandstone surface relaxivity determined by NMR T 2 distribution and digital rock simulation for permeability evaluation, J Pet Sci Eng. 193 (2020) 107400. https://doi.org/10.1016/j.petrol.2020.107400.
Y. Peysson, M. Fleury, V. Blázquez-Pascual, Drying Rate Measurements in Convection- and Diffusion-Driven Conditions on a Shaly Sandstone Using Nuclear Magnetic Resonance, Transp Porous Media. 90 (2011) 1001–1016. https://doi.org/10.1007/s11242-011-9829-3.
S. Ghoshal, C. Mattea, S. Stapf, Inhomogeneity in the drying process of gelatin film formation: NMR microscopy and relaxation study, Chem Phys Lett. 485 (2010) 343–347. https://doi.org/10.1016/j.cplett.2009.12.064.
M. Cocusse, M. Rosales, B. Maillet, R. Sidi-Boulenouar, E. Julien, S. Caré, P. Coussot, Two-step diffusion in cellular hygroscopic (vascular plant-like) materials, Sci Adv. 8 (2023) eabm7830. https://doi.org/10.1126/sciadv.abm7830.
B. Maillet, G. Dittrich, P. Huber, P. Coussot, Diffusionlike Drying of a Nanoporous Solid as Revealed by Magnetic Resonance Imaging, Phys Rev Appl. 18 (2022). https://doi.org/10.1103/PhysRevApplied.18.054027.
G.R. Coates, L. Xiao, M.G. Prammer, NMR Logging Principles and Applications, Halliburton Energy Services, Houston, 2001.
P.J. McDonald, Stray field magnetic resonance imaging, Prog Nucl Magn Reson Spectrosc. 30 (1997) 69–99. https://doi.org/https://doi.org/10.1016/S0079-6565(96)01035-7.
K. Smith, A. Burbidge, D. Apperley, P. Hodgkinson, F.A. Markwell, D. Topgaard, E. Hughes, Stray-field NMR diffusion q-space diffraction imaging of monodisperse coarsening foams, J Colloid Interface Sci. 476 (2016) 20–28. https://doi.org/10.1016/j.jcis.2016.04.053.
R. de Oliveira-Silva, É. Lucas-Oliveira, A.G. de Araújo-Ferreira, W.A. Trevizan, E.L.G. Vidoto, D. Sakellariou, T.J. Bonagamba, A benchtop single-sided magnet with NMR well-logging tool specifications – Examples of application, Journal of Magnetic Resonance. 322 (2021) 106871. https://doi.org/https://doi.org/10.1016/j.jmr.2020.106871.
P. Blümler, F. Casanova, CHAPTER 4: Hardware developments: Single-sided magnets, in: New Developments in NMR, 2016: pp. 110–132. https://doi.org/10.1039/9781782628095-00110.
B. Manz, A. Coy, R. Dykstra, C.D. Eccles, M.W. Hunter, B.J. Parkinson, P.T. Callaghan, A mobile one-sided NMR sensor with a homogeneous magnetic field: The NMR-MOLE, Journal of Magnetic Resonance. 183 (2006) 25–31. https://doi.org/https://doi.org/10.1016/j.jmr.2006.07.017.
M.N. D’Eurydice, P. Galvosas, D-T2 correlation using the inhomogeneity of single sided NMR devices, Microporous and Mesoporous Materials. 205 (2015) 40–43. https://doi.org/10.1016/j.micromeso.2014.08.026.
G. Eidmann, R. Savelsberg, P. Blümler, B. Blümich, The NMR MOUSE, a Mobile Universal Surface Explorer, J Magn Reson A. 122 (1996) 104–109. https://doi.org/https://doi.org/10.1006/jmra.1996.0185.
B. Blümich, P. Blümler, G. Eidmann, A. Guthausen, R. Haken, U. Schmitz, K. Saito, G. Zimmer, The NMR-mouse: construction, excitation, and applications, Magn Reson Imaging. 16 (1998) 479–484. https://doi.org/https://doi.org/10.1016/S0730-725X(98)00069-1.
F. Casanova, J. Perlo, B. Blümich, Single-sided NMR, 2011. https://doi.org/10.1007/978-3-642-16307-4.
B. Blümich, J. Anders, When the MOUSE leaves the house, Magn. Reson. 2 (2021) 149–160. https://doi.org/10.5194/mr-2-149-2021.
B. Blümich, D. Jaschtschuk, C. Rehorn, Advances and Adventures with Mobile NMR, in: Magnetic Resonance Microscopy, 2022: pp. 155–172. https://doi.org/https://doi.org/10.1002/9783527827244.ch7.
B. Blümich, J. Perlo, F. Casanova, Mobile single-sided NMR, Prog Nucl Magn Reson Spectrosc. 52 (2008) 197–269. https://doi.org/10.1016/j.pnmrs.2007.10.002.
B. Blümich, S. Anferova, R. Pechnig, H. Pape, J. Arnold, C. Clauser, Mobile NMR for porosity analysis of drill core sections, 1 (2004) 177–180. https://doi.org/10.1088/1742-2132/1/3/001.
S. Costabel, T. Hiller, R. Dlugosch, S. Kruschwitz, M. Müller-Petke, Evaluation of single-sided nuclear magnetic resonance technology for usage in geosciences, Meas Sci Technol. 34 (2023). https://doi.org/10.1088/1361-6501/ac9800.
J. Perlo, F. Casanova, B. Blümich, Profiles with microscopic resolution by single-sided NMR, Journal of Magnetic Resonance. 176 (2005) 64–70. https://doi.org/10.1016/j.jmr.2005.05.017.
E. V. Silletta, M.I. Velasco, G.A. Monti, R.H. Acosta, Comparison of experimental times in T1-D and D-T2 correlation experiments in single-sided NMR, Journal of Magnetic Resonance. 334 (2022) 107112. https://doi.org/10.1016/j.jmr.2021.107112.
H.Y. Carr, E.M. Purcell, Effects of diffusion on free precession in nuclear magnetic resonance experiments, Physical Review. 94 (1954) 630–638. https://doi.org/10.1103/PhysRev.94.630.
S. Meiboom, D. Gill, Modified Spin-Echo Method for Measuring Nuclear Relaxation Times, Review of Scientific Instruments. 29 (1958) 688. https://doi.org/10.1063/1.1716296.
S.W. Provencher, A constrained regularization method for inverting data represented by linear algebraic or integral equations, Comput Phys Commun. 27 (1982) 213–227. https://doi.org/10.1016/0010-4655(82)90173-4.
M.D. Hürlimann, Diffusion and relaxation effects in general stray field NMR experiments, Journal of Magnetic Resonance. 148 (2001) 367–378. https://doi.org/10.1006/jmre.2000.2263.
R. Kimmich, Principles of Soft-Matter Dynamics, Springer Netherlands, Dordrecht, 2012. https://doi.org/10.1007/978-94-007-5536-9.
P.T. Callaghan, Translational Dynamics and Magnetic Resonance: Principles of Pulsed Gradient Spin Echo NMR, Oxford University Press, Oxford, 2011.
P.J. Mcdonald, B. Newling, Stray field magnetic resonance imaging, 1998.
M.D. Hürlimann, D.D. Griffin, Spin dynamics of Carr-Purcell-Meiboom-Gill-like sequences in grossly inhomogeneous B(0) and B(1) fields and application to NMR well logging., J Magn Reson. 143 (2000) 120–35. https://doi.org/10.1006/jmre.1999.1967.
M.D. Hürlimann, L. Venkataramanan, Quantitative Measurement of Two-Dimensional Distribution Functions of Diffusion and Relaxation in Grossly Inhomogeneous Fields, Journal of Magnetic Resonance. 157 (2002) 31–42. https://doi.org/10.1006/jmre.2002.2567.
G. Goelman, M.G. Prammer, The CPMG Pulse Sequence in Strong Magnetic Field Gradients with Applications to Oil-Well Logging, J Magn Reson A. 113 (1995) 11–18. https://doi.org/10.1006/jmra.1995.1050.
Y.-Q. Song, Categories of Coherence Pathways for the CPMG Sequence, Journal of Magnetic Resonance. 157 (2002) 82–91. https://doi.org/10.1006/jmre.2002.2577.
J.L. Markley, W.J. Horsley, M.P. Klein, Spin-Lattice Relaxation Measurements in Slowly Relaxing Complex Spectra, J Chem Phys. 55 (1971) 3604–3605. https://doi.org/10.1063/1.1676626.
W.E. Kenyon, P.I. Day, C. Straley, J.F. Willemsen, Three-part study of NMR longitudinal relaxation properties of water-saturated sandstones, SPE Formation Evaluation. 3 (1988) 622–636. https://doi.org/10.2118/15643-pa.
M. Esteves Ferreira, M. Rodrigues Del Grande, R. Neumann Barros Ferreira, A. Ferreira da Silva, M. Nogueira Pereira da Silva, J. Tirapu-Azpiroz, E. Lucas-Oliveira, A.G. de Araújo Ferreira, R. Soares, C. B Eckardt, T. J Bonagamba, M. Steiner, Full scale, microscopically resolved tomographies of sandstone and carbonate rocks augmented by experimental porosity and permeability values, Sci Data. 10 (2023) 368. https://doi.org/10.1038/s41597-023-02259-z.
P.D. Teal, C. Eccles, Adaptive truncation of matrix decompositions and efficient estimation of NMR relaxation distributions, Inverse Probl. 31 (2015) 045010. https://doi.org/10.1088/0266-5611/31/4/045010.
A.G. Yiotis, I.N. Tsimpanogiannis, A.K. Stubos, Y.C. Yortsos, Pore-network study of the characteristic periods in the drying of porous materials., J Colloid Interface Sci. 297 (2006) 738–48. https://doi.org/10.1016/j.jcis.2005.11.043.
P. Coussot, Scaling approach of the convective drying of a porous medium, The European Physical Journal B-Condensed Matter and Complex Systems. 15 (2000) 557–566. https://doi.org/10.1007/s100510051160.
M. Kaviany, M. Mittal, Funicular state in drying of a porous slab, Int J Heat Mass Transf. 30 (1987) 1407–1418. https://doi.org/10.1016/0017-9310(87)90172-4.
M.I. Velasco, E.V. Silletta, C.G. Gomez, M.C. Strumia, S. Stapf, G.A. Monti, C. Mattea, R.H. Acosta, Spatially Resolved Monitoring of Drying of Hierarchical Porous Organic Networks, Langmuir. 32 (2016) 2067–2074. https://doi.org/10.1021/acs.langmuir.5b04230.

No competing interests reported.

Download PDF

Journal Publication

published 01 Aug, 2023

Read the published version in Applied Magnetic Resonance →

Editorial decision: Major revision
12 Jul, 2023
Reviews received at journal
10 Jul, 2023
Reviewers agreed at journal
02 Jul, 2023
Reviewers invited by journal
28 Jun, 2023
Editor assigned by journal
28 Jun, 2023
Submission checks completed at journal
28 Jun, 2023
First submitted to journal
24 Jun, 2023

You are reading this latest preprint version

Spatially resolved dynamic longitudinal relaxometry in single-sided NMR

Status:

Journal Publication

Version 1

Abstract

Figures

1. Introduction

2. Materials and methods

2.1 Profiles with transverse relaxation times

2.2 Profiles with longitudinal relaxation times

2.2 Sample description

3. Results and discussion

3.1 Determination of the minimal number of points for T₁ encoding

3.2 Evaporation monitoring

4. Conclusions

Declarations

References

Additional Declarations

Status:

Journal Publication

Version 1

Spatially resolved dynamic longitudinal relaxometry in single-sided NMR

Status:

Journal Publication

Version 1

Abstract

Figures

1. Introduction

2. Materials and methods

2.1 Profiles with transverse relaxation times

2.2 Profiles with longitudinal relaxation times

2.2 Sample description

3. Results and discussion

3.1 Determination of the minimal number of points for T1 encoding

3.2 Evaporation monitoring

4. Conclusions

Declarations

References

Additional Declarations

Status:

Journal Publication

Version 1

3.1 Determination of the minimal number of points for T₁ encoding