物理学报2016，Vol.65Issue(10)：107601-1-107601-12,12.DOI:10.7498/aps.65.107601
基于改进非线性拟合的核磁共振T2谱多指数反演∗
Multi-exp onential inversion of T2 sp ectrum in NMR based on improved nonlinear fitting
摘要
Abstract
Multi-exponential inversion algorithm of nuclear magnetic resonance (NMR) T2 spectrum is an important math-ematical tool for the NMR relaxation study of complicated samples. The popular algorithm usually obtains the T2 spectrum by linear fitting under the prescribed distribution of T2. When the T2 spectrum is dispersed, such a procedure is inaccurate because of the lack of adaptive prescription and the limit of linear method. Nonlinear fitting method does not fix the T2 distribution, and it provides the positions and the weights of T2 simultaneously via the nonlinear fitting of multi-exponential function. In this case, the problem of multi-exponential inversion is transformed into a nonlinear optimization problem with non-negative constraints. The optimization objective function is the residual sum of squares (or residual sum of squares with regularization). The nonlinear optimization problem can usually be solved by Levenberg-Marquardt algorithm and evolutionary algorithm. But the results of Levenberg-Marquardt algorithm are dependent on initial values, and the calculation of evolutionary algorithm is complicated. We provide an optimal model for the nonlinear fitting in the inversion of dispersed T2 spectrum based on the linear regression and least-squares. The key idea is that the optimal weights of T2 can be calculated by least square when the positions of T2 are fixed, although the positions of T2 are adjusted adaptively. So we can relate the positions to weights appropriately to improve the popular nonlinear fitting algorithms. Such an improvement can reduce the searching inversion parameters, speed up its convergence and reduce the dependence on initial value. Incorporating it into the Levenberg-Marquardt algorithm or evolutionary algorithm can improve the inversion accuracy and make the algorithm more robust. The validity of our improvement is demonstrated by the inversions of simulation data and practical NMR data by combining Levenberg-Marquardt algorithm and differential evolution algorithm with our improvement. The inversion results of simulation data show that for dispersed T2 spectrum, the algorithm using this improvement can obtain more accurate T2 spectrum than previous ones, especially in the case of low signal-to-noise ratio (SNR) cases. The inversion results also indicate that the improvement can reduce the dependence on initial value of Levenberg-Marquardt algorithm, and can accelerate the convergence of differential evolution algorithm. The inversion results of practical NMR data show that the algorithm using the improvement can obtain more accurate T2 spectrum than the widely used CONTIN program in the case of low signal-to-noise ratio (SNR). The inversion results of oil-water mixture sample NMR data also demonstrate that the relaxation time T2 is independent of dispersion degree of immiscible system components.关键词
核磁共振/多指数反演/非线性拟合/差分进化Key words
nuclear magnetic resonance/multi-exponential inversion/nonlinear fitting/differential evolution引用本文复制引用
吴量,陈方,黄重阳,丁国辉,丁义明..基于改进非线性拟合的核磁共振T2谱多指数反演∗[J].物理学报,2016,65(10):107601-1-107601-12,12.基金项目
国家重点基础研究发展计划(批准号:2013CB910200)和国家自然科学基金(批准号:11405264)资助的课题 （批准号:2013CB910200）
