核磁共振T2谱M-P广义逆反演研究

基于广义逆和奇异值分解理论,研究核磁共振T2谱反演的截断法和阻尼法。首先给出反演问题的M-P广义逆解;然后对系数矩阵进行奇异值分解求其M-P广义逆,对奇异值进行截断或者加阻尼,保证解的稳定性兼顾其分辨率;最后利用迭代技术实现解的非负约束。数值仿真实验和岩心T2谱反演表明:对于双峰谱,在信噪比SNR≥10时,截断法和阻尼法反演T2谱的分辨率和稳定性都很高;截断法截断因子选择范围小,在SNR<10时解的分辨率和稳定性差;阻尼法阻尼因子选择范围大,在SNR<10时解的稳定性好,但短组分分辨率低,长组分分辨率高。对于三峰谱,当SNR≥20时,截断法和阻尼法反演的T2谱光滑连续,具有三峰结构;当SNR<20时两者均不能反演出与构造谱接近的三峰谱。因此,2种方法可应用于SNR≥10的双峰谱和SNR≥20的三峰谱反演;从分辨率和稳定性、平滑噪声及截断和阻尼因子的选取上,阻尼法略优于截断法。