Efficient numerical Bloch solutions for multipulse surface NMR
- 12 August 2021
- journal article
- research article
- Published by Oxford University Press (OUP) in Geophysical Journal International
- Vol. 227 (3), 1905-1916
- https://doi.org/10.1093/gji/ggab321
Abstract
Simplified solutions of the Bloch equation can lead to inaccurate estimates of hydrogeological parameters from surface nuclear magnetic resonance measurements. Even for single pulse measurements, using simplified forward models is common practice because of the computational intensity of obtaining the full-Bloch solution. These challenges are exacerbated for multipulse sequences. We show parallelizing the full-Bloch solver on a Graphics Processing Unit reduces the solve time by three orders of magnitude. Further optimizations by numerical, analytical and hybrid solutions yield an additional 3× speed up. We simulate the full-Bloch physics for free-induction decay, spin-echo and pseudo-saturation recovery excitation schemes for an unprecedented range of physical scenarios. We explore the time-dependence and relaxation time sensitivity in these solution spaces. Characterizing the solution spaces with polynomials of the relaxation times, the solutions can be rapidly reproduced; a technique known as fast-mapping. By fitting these higher dimensional solution ensembles with polynomials, the original fast-mapping technique is extended to include T1 at arbitrary times. Accuracy of the 7th order polynomial is such that a minimum 96 per cent of the models are within a ±3 per cent relative error. This permits the rapid reproduction of full-Bloch solutions with a matrix multiplication and opens up surface NMR to time-series based inversion of single and multipulse data.This publication has 14 references indexed in Scilit:
- The inversion of surface-NMR T1 data for improved aquifer characterizationGeophysics, 2013
- Accounting for relaxation processes during the pulse in surface NMR dataGeophysics, 2009
- Pore size distributions and hydraulic conductivities of rocks derived from Magnetic Resonance Sounding relaxation data using multi-exponential decay time inversionJournal of Applied Geophysics, 2008
- Scalable Parallel Programming with CUDAQueue, 2008
- The internal magnetic field distribution, and single exponential magnetic resonance free induction decay, in rocksJournal of Magnetic Resonance, 2005
- Magnetic Resonance Sounding Applied to Aquifer CharacterizationGroundwater, 2004
- Pore-Size Distributions and Tortuosity in Heterogeneous Porous MediaJournal of Magnetic Resonance, Series A, 1995
- Selective pulses in NMR imaging: A reply to criticismJournal of Magnetic Resonance (1969), 1979
- Spin EchoesPhysical Review B, 1950
- Nuclear InductionPhysical Review B, 1946