Random-walk technique for simulating NMR measurements and 2D NMR maps of porous media with relaxing and permeable boundaries

We revisit random-walk methods to simulate the NMR response of fluids in porous media. Simulations reproduce the effects of diffusion within external inhomogeneous background magnetic fields, imperfect and finite-duration B1 pulses, T1/T2 contrasts, and relaxing or permeable boundaries. The simulation approach consolidates existing NMR numerical methods used in biology and engineering into a single formulation that expands on the magnetic-dipole equivalent of spin packets. When fluids exhibit low T1/T2 contrasts and when CPMG pulse sequences are used to acquire NMR measurements, we verify that classical NMR numerical models that neglect T1 effects accurately reproduce surface magnetization decays of saturated granular porous media regardless of the diffusion/relaxation regime. Currently, analytical expressions exist only for the case of arbitrary pore shapes within the fast-diffusion limit. However, when fluids include several components or when magnetic fields are strongly inhomogeneous, we show that simulations results obtained using the complete set of Bloch's equations differ substantially from those of classical NMR models. In addition, our random-walk formulation accurately reproduces magnetization echoes stemming from coherent-pathway calculations. We show that the random-walk approach is especially suited to generate parametric multi-dimensional T1/T2/D NMR maps to improve the characterization of pore structures and saturating fluids.