A Robust Algorithm for 2-D NMR Diffusion–Relaxation Spectra Inversion
- 31 July 2018
- journal article
- research article
- Published by Institute of Electrical and Electronics Engineers (IEEE) in IEEE Geoscience and Remote Sensing Letters
- Vol. 15 (10) , 1545-1549
Abstract
A two-step iterative shrinkage/thresholding (TIST) algorithm was proposed to solve the objective function of the L1 regularization for the 2-D nuclear magnetic resonance diffusion-relaxation (D-T 2 ) spectra inversion. Moreover, a simple and effective method combining the S-curve method and the Chi-square method was proposed to select the optimal L1 regularization parameter. We compared the inverted D-T 2 spectra derived from the TIST algorithm with those from the truncated singular value decomposition (TSVD), Butler-Reeds- Dawson (BRD), original IST (OIST), and IST algorithms using numerical experiments. The experimental results indicated that the inverted D-T 2 spectra derived from the OIST, IST, and TIST algorithms were affected by noise to a smaller degree, and the robustness was higher for the OIST, IST, and TIST algorithms than for the TSVD and BRD algorithms. In addition, we compared the accuracy and operating time of the TIST, OIST, and IST algorithms for D-T 2 spectra inversion and concluded that the three algorithms had an equivalent accuracy but the TIST algorithm had a faster convergence rate.Keywords
Funding Information
- National Natural Science Foundation of China (41674126)
This publication has 23 references indexed in Scilit:
- Joint inversion ofspectrum combining the iterative truncated singular value decomposition and the parallel particle swarm optimization algorithmsComputer Physics Communications, 2016
- NMR petrophysical interpretation method of gas shale based on core NMR experimentJournal of Petroleum Science and Engineering, 2015
- A novel method for NMR data compressionComputational Geosciences, 2015
- A new inversion method for (T2, D) 2D NMR logging and fluid typingComputers & Geosciences, 2013
- A bound optimization approach to wavelet-based image deconvolution2005
- A global inversion method for multi-dimensional NMR loggingJournal of Magnetic Resonance, 2005
- Signal Recovery by Proximal Forward-Backward SplittingMultiscale Modeling & Simulation, 2005
- An iterative thresholding algorithm for linear inverse problems with a sparsity constraintCommunications on Pure and Applied Mathematics, 2004
- An EM algorithm for wavelet-based image restorationIEEE Transactions on Image Processing, 2003
- Estimating Solutions of First Kind Integral Equations with Nonnegative Constraints and Optimal SmoothingSIAM Journal on Numerical Analysis, 1981