A modeling and inversion method of spin echoes to measure magnetic resonance sounding transverse relaxation time in surface applications
Authors: Yang Yu-Jing , Ye Rui , Zhao Han-Qing , Wan Ling , Lin Ting-Ting ,   College of Instrumentation and Electrical Engineering, Jilin University, Changchun 130061, China,   Key Laboratory of Geophysical Exploration Equipment, Ministry of Education (Jilin University), Changchun 130061, China
Published: March 2021

Surface magnetic resonance sounding (MRS) has generally been considered to be an efficient tool for hydrological investigations. As is well known, the effective relaxation time <inline-formula><tex-math id="M1">\begin{document}$ T_2^*$\end{document}</tex-math><alternatives><graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="6-20201427_M1.jpg"/><graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="6-20201427_M1.png"/></alternatives></inline-formula> which characterizes the decay rate of MRS free-decay-induction (FID) signal and is used to measure pore-scale properties, is particularly limited for several special cases (e.g. areas with magnetic rock subsurfaces). Recent years, the transverse relaxation time <inline-formula><tex-math id="M2">\begin{document}$ T_2$\end{document}</tex-math><alternatives><graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="6-20201427_M2.jpg"/><graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="6-20201427_M2.png"/></alternatives></inline-formula> obtained from spin-echo signal was adopted to implement the surface MRS, and showed great potentials for estimating the porosity and permeability. However, owning to the short period of development, the related modeling and inversion strategies are rarely introduced and summarized. Actually, the general practice for surface MRS <inline-formula><tex-math id="M3">\begin{document}$ T_2$\end{document}</tex-math><alternatives><graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="6-20201427_M3.jpg"/><graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="6-20201427_M3.png"/></alternatives></inline-formula> measurement fits the spin-echo by the exponential function and the fitting line was directly used as the FID signal for inversion. This scheme not only limits the precision of interpretation, but also loses part of valid information about original field data. Aiming at these problems, in this paper, we introduce the calculation of forward model and thus a two-stage framework with singular value decomposition (SVD) linear inversion involved is derived to quantify the <inline-formula><tex-math id="M4">\begin{document}$ T_2$\end{document}</tex-math><alternatives><graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="6-20201427_M4.jpg"/><graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="6-20201427_M4.png"/></alternatives></inline-formula> distributed with depth. Considering the fact that the inversion result of SVD is always strongly affected by the noise level, an improved method which combines the simultaneous iterative reconstruction technology (SIRT) with SVD is proposed. To be specific, we compare the measurement schemes with kernel functions between <inline-formula><tex-math id="M5">\begin{document}$ T_2$\end{document}</tex-math><alternatives><graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="6-20201427_M5.jpg"/><graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="6-20201427_M5.png"/></alternatives></inline-formula> and the original theory in MRS, and then provide the forward and inversion formulations. In order to substantiate the effectiveness of this method, we conduct the synthetic experiments for Carr-Purcell-Meiboom-Gill sequence and explain the dataset with the mentioned strategies. As expected, the combined approach possesses a better performance in shallow layer with an error of 1.5% and 0.02 s for water content and <inline-formula><tex-math id="M6">\begin{document}$ T_2$\end{document}</tex-math><alternatives><graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="6-20201427_M6.jpg"/><graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="6-20201427_M6.png"/></alternatives></inline-formula> for the contaminated data, respectively. With these advantages, it is expected to realize the adoption of the SVD with SIRT in field applications and further investigate the aquifer characterizations in the future.