Retaining both discrete and smooth features in 1D and 2D NMR relaxation and diffusion experiments

Simple MATLAB and Python scripts for multi-exponential analysis

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.

Automated component analysis in DOSY NMR using information criteria

This study introduces a model selection technique based on Bayesian information criteria for estimating the number of components in a mixture during Diffusion-Ordered Spectroscopy (DOSY) Nuclear Magnetic Resonance (NMR) data analysis. As the accuracy of this technique is dependent on the efficiency of parameter estimators, we further investigate the performance of the Weighted Least Squares (WLS) and Maximum a Posteriori (MAP) estimators. The WLS method, enhanced with meticulously tuned L2-regularization, effectively detects components when the difference in self-diffusion coefficients is more than two-fold, especially when the component with the smaller coefficient has a larger weight ratio. The MAP method, strengthened by a substantial database of prior information, exhibits outstanding precision, decreasing this threshold to 1.5 times. Both estimators provide weight ratio estimates with standard deviations of approximately around 1 percentage point, although the MAP method tends to overestimate the component with a larger self-diffusion coefficient. Deviations from the expected values can exceed 10 percentage points, often due to inaccuracies in component detection. The error estimates are determined using data resampling techniques derived from a large-scale 1000-point experiment and an additional five measurements from a single-component mixture. This approach allowed us to thoroughly examine data distribution characteristics, thereby laying a robust groundwork for future refinement efforts.

In Situ Determination of Carbon Number Distributions of Mixtures of Linear Hydrocarbons Confined within Porous Media Using Pulsed Field Gradient NMR

Pulsed field gradient (PFG) NMR measurements, combined with a novel optimization method, are used to determine the composition of hydrocarbon mixtures of linear alkanes (C7-C16) in both the bulk liquid state and when imbibed within a porous medium of mean pore diameter 28.6 nm. The method predicts the average carbon number of a given mixture to an accuracy of ±1 carbon number and the mole fraction of a mixture component to within an average root-mean-square error of ±0.036 with just three calibration mixtures. Given that the method can be applied at any conditions of temperature and pressure at which the PFG NMR measurements are made, the method has the potential for application in characterizing hydrocarbon liquid mixtures inside porous media and at the operating conditions relevant to, for example, hydrocarbon recovery and heterogeneous catalysis.

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.