一种核磁共振二维谱反演的方法.pdf

(19)中华人民共和国国家知识产权局*CN103116148A*(12)发明专利申请(10)申请公布号(10)申请公布号CNCN103116148103116148A(43)申请公布日2013.05.22(21)申请号201310035140.8(22)申请日2013.01.30(71)申请人上海理工大学地址200093上海市杨浦区军工路516号(72)发明人周小龙聂生东王远军张英力杨培强(74)专利代理机构上海申汇专利代理有限公司31001代理人吴宝根(51)Int.Cl.G01R33/56(2006.01)权权利要求书1页利要求书1页说明书7页说明书7页附图3页附图3页(54)发明名称一种核磁共振二维谱反演的方法(57)摘要本发明涉及一种核磁共振二维谱反演的方法，包括第一步，噪声提取与估计：使用小波变换对采集数据CPMG回波串中的噪声进行提取并估计其标准差；第二步，数据压缩：生成反演核，并利用核矩阵的秩进行截断奇异值分解与重构来完成数据压缩；第三步，数据拟合：对压缩后数据的拟合问题进行正则化处理，并用结合了非精确一维搜索的牛顿法进行正则化因子和反演谱的迭代求解，得到反演谱。大幅提高了二维反演算法的执行效率与二维谱的分辨率，同时还具有很好的鲁棒性。CN103116148ACN103648ACN103116148A权利要求书1/1页1.一种核磁共振二维谱反演的方法，其特征在于，包括如下步骤：1）噪声提取与估计：使用小波变换对采集数据CPMG回波串中的噪声进行提取并估计其标准差；2）数据压缩：生成反演核，并利用核矩阵的秩进行截断奇异值分解与重构来完成数据压缩；3）数据拟合：对压缩后数据的拟合问题进行正则化处理，并用结合了非精确一维搜索的牛顿法进行正则化因子和反演谱的迭代求解，得到反演谱。2.根据权利要求1所述核磁共振二维谱反演的方法，其特征在于，所述步骤1）中使用小波去噪的方法对CPMG数据串进行去噪，通过比较去噪后数据与滤波前数据的幅度差异，进行“噪声”提取，首先，使用sym8小波基对原始数据进行4级分解，然后使用全局阈值对小波系数进行软阈值化，最后进行信号重构，得到滤波后相对光滑的CPMG数据串；然后计算原始数据与滤波后数据的差值，使用高斯白噪声的幅度模型对差值数据进行拟合，得到对噪声标准差的估计。3.根据权利要求1所述核磁共振二维谱反演的方法，其特征在于，所述步骤3）数据拟合首先对核心矩阵各元素平方后得到一个新矩阵，将新矩阵的迹(trace)除以对角元素个数的值作为初始的正则化因子；其次对于给定的正则因子，使用牛顿法（NewtonMethod）对极小值问题进行求解，在牛顿方向上使用的步长由Wolfe-Powell非精确一维搜索确定。4.根据权利要求1或3所述核磁共振二维谱反演的方法，其特征在于，所述步骤3）正则化因子更新方式为：根据离差原理计算正则化因子的理想值，如果理想值大于原来的值，那么将原值的1/2作为新的正则化因子的值；否则，直接使用理想值作为正则化因子的新值。5.根据权利要求1所述核磁共振二维谱反演的方法，其特征在于，所述步骤3）反演谱的迭代求解的终值条件为：①正则化因子下一轮循环的新值已经不大于一个预设的小值；正则化因子在连续两次迭代的过程中变化幅度小于预设阈值；拟合残差在连续两次迭代的过程中变化幅度小于预设阈值；拟合残差不大于第一步估计的噪声水平。2CN103116148A说明书1/7页一种核磁共振二维谱反演的方法技术领域[0001]本发明涉及一种核磁共振领域中信号处理，特别涉及一种核磁共振二维谱反演的方法。背景技术[0002]核磁共振技术在能源勘探领域已经得到了广泛应用。NMR测井能够快速、无损地为储层流体识别、岩石物性评价和产能评估等常规测井工作流中的各个步骤提供精准的信息。相比电法测井、声波测井等传统方法，NMR测井法是惟一一种不使流体流动就能对流体类型进行识别的方法。基于一维测井实验的解决方案在效率和精度上都存在明显的不足，二维NMR测井技术应运而生。使用D-T2（扩散-横向弛豫）二维谱可以快速、直观地对水、油等不同成分进行区分，通过计算T1/T2（纵向弛豫/横向弛豫）这一比值还能准确地对烃的类型进行判断，二维NMR测井方法在定性和定量分析中都有着得天独厚的优势。二维NMR检测方法的实现还拓宽了常规NMR检测实验的应用范围，为食品、农业、生物材料等领域提供了更精确、可靠的解决方案。[0003]NMR检测方法采集的数据并不能直接使用，真正需要的“谱”信息需要通过这些实验数据进行反演后才能得到。目前世界上只有极少数NMR测井设备制造商有能力生产二维测井设备，而且这些国际巨头只向各国能源公司提供仪器租赁和有偿的解释服务，并不出售设备和反演软件，二维反演算法几近处于被垄断状态。另一方面，二维实验采集的数据量远