Search

Kerogen-rich shale reservoirs will play a key role during the energy transition, yet the effects of nano-confinement on the NMR relaxation of hydrocarbons in kerogen are poorly understood. We use atomistic MD simulations to investigate the effects of nano-confinement on the $^1$H NMR relaxation times $T_1$ and $T_2$ of heptane in kerogen. In the case of $T_1$, we discover the important role of confinement in reducing $T_1$ by $\sim$3 orders of magnitude from bulk heptane, in agreement with measurements of heptane dissolved in kerogen from the Kimmeridge Shale, without any models or free parameters. In the case of $T_2$, we discover that confinement breaks spatial isotropy and gives rise to residual dipolar coupling which reduces $T_2$ by $\sim$5 orders of magnitude from bulk heptane. We use the simulated $T_2$ to calibrate the surface relaxivity and thence predict the pore-size distribution of the organic nano-pores in kerogen, without additional experimental data.

Atomistic molecular dynamics simulations are used to investigate $^1$H NMR $T_1$ relaxation of water from paramagnetic Gd$^{3+}$ ions in solution at 25$^{\circ}$C. Simulations of the $T_1$ relaxivity dispersion function $r_1$ computed from the Gd$^{3+}$--$^1$H dipole--dipole autocorrelation function agree within $\simeq 8$\% of measurements in the range $f_0 \simeq $ 5 $\leftrightarrow$ 500 MHz, without any adjustable parameters in the interpretation of the simulations, and without any relaxation models. The simulation results are discussed in the context of the Solomon-Bloembergen-Morgan inner-sphere relaxation model, and the Hwang-Freed outer-sphere relaxation model. Below $f_0 \lesssim $ 5 MHz, the simulation overestimates $r_1$ compared to measurements, which is used to estimate the zero-field electron-spin relaxation time. The simulations show potential for predicting $r_1$ at high frequencies in chelated Gd$^{3+}$ contrast-agents used for clinical MRI.

The mechanism behind the $^1$H NMR frequency dependence of $T_1$ and the viscosity dependence of $T_2$ for polydisperse polymers and bitumen remains elusive. We elucidate the matter through NMR relaxation measurements of polydisperse polymers over an extended range of frequencies ($f_0 = 0.01 \leftrightarrow$ 400 MHz) and viscosities ($\eta = 385 \leftrightarrow 102,000$ cP) using $T_{1}$ and $T_2$ in static fields, $T_{1}$ field-cycling relaxometry, and $T_{1\rho}$ in the rotating frame. We account for the anomalous behavior of the log-mean relaxation times $T_{1LM} \propto f_0$ and $T_{2LM} \propto (\eta/T)^{-1/2}$ with a phenomenological model of $^1$H-$^1$H dipole-dipole relaxation which includes a distribution in molecular correlation times and internal motions of the non-rigid polymer branches. We show that the model also accounts for the anomalous $T_{1LM}$ and $T_{2LM}$ in previously reported bitumen measurements. We find that molecular dynamics (MD) simulations of the $T_{1} \propto f_0$ dispersion and $T_2$ of similar polymers simulated over a range of viscosities ($\eta = 1 \leftrightarrow 1,000$ cP) are in good agreement with measurements and the model. The $T_{1} \propto f_0$ dispersion at high viscosities agrees with previously reported MD simulations of heptane confined in a polymer matrix, which suggests a common NMR relaxation mechanism between viscous polydisperse fluids and fluids under confinement, without the need to invoke paramagnetism.

The role of internal motions and molecular geometry on $^1$H NMR relaxation times $T_{1,2}$ in hydrocarbons is investigated using MD (molecular dynamics) simulations of the autocorrelation functions for in{\it tra}molecular $G_R(t)$ and in{\it ter}molecular $G_T(t)$ $^1$H-$^1$H dipole-dipole interactions arising from rotational ($R$) and translational ($T$) diffusion, respectively. We show that molecules with increased molecular symmetry such as neopentane, benzene, and isooctane show better agreement with traditional hard-sphere models than their corresponding straight-chain $n$-alkane, and furthermore that spherically-symmetric neopentane agrees well with the Stokes-Einstein theory. The influence of internal motions on the dynamics and $T_{1,2}$ relaxation of $n$-alkanes are investigated by simulating rigid $n$-alkanes and comparing with flexible (i.e. non-rigid) $n$-alkanes. Internal motions cause the rotational and translational correlation-times $\tau_{R,T}$ to get significantly shorter and the relaxation times $T_{1,2}$ to get significantly longer, especially for longer-chain $n$-alkanes. Site-by-site simulations of $^1$H's along the chains indicate significant variations in $\tau_{R,T}$ and $T_{1,2}$ across the chain, especially for longer-chain $n$-alkanes. The extent of the stretched (i.e. multi-exponential) decay in the autocorrelation functions $G_{R,T}(t)$ are quantified using inverse Laplace transforms, for both rigid and flexible molecules, and on a site-by-site bases. Comparison of $T_{1,2}$ measurements with the site-by-site simulations indicate that cross-relaxation (partially) averages-out the variations in $\tau_{R,T}$ and $T_{1,2}$ across the chain of long-chain $n$-alkanes. This work also has implications on the role of nano-pore confinement on the NMR relaxation of fluids in the organic-matter pores of kerogen and bitumen.

Quantifying the statistics of occupancy of solvent molecules in the vicinity of solutes is central to our understanding of solvation phenomena. Number fluctuations in small `solvation shells' around solutes cannot be described within the macroscopic grand canonical framework using a single chemical potential that represents the solvent `bath'. In this communication, we hypothesize that molecular-sized observation volumes such as solvation shells are best described by coupling the solvation shell with a mixture of particle baths each with its own chemical potential. We confirm our hypotheses by studying the enhanced fluctuations in the occupancy statistics of hard sphere solvent particles around a distinguished hard sphere solute particle. Connections with established theories of solvation are also discussed.