基于三次样条插值的地面核磁共振反演方法研究

作     者：冷舒喆 

作者单位：吉林大学 

学位级别：硕士

导师姓名：林婷婷

授予年度：2017年

学科分类：0809[工学-电子科学与技术（可授工学、理学学位）] 07[理学] 070205[理学-凝聚态物理] 08[工学] 0805[工学-材料科学与工程（可授工学、理学学位）] 080502[工学-材料学] 0702[理学-物理学] 

主      题：地面核磁共振 核函数 三次样条插值 含水量反演 二分法 

摘      要：地面核磁共振技术利用地下水中氢质子的自旋性质,通过在地面铺设探测线圈,对地下水进行直接的无损的探测,凭借其方法的高效性和对结果评估的准确性,该技术成为主要的现代地球物理找水方法。近年来,随着磁共振探测方法研究的不断深入,已经由一维探测向二维、三维探测发展。目前常用的二维磁共振探测通常采用阵列式线圈,在地面同一条测线上铺设一个发射线圈和多个接收线圈,接收线圈采用边对边或者半覆盖方式内嵌在发射线圈内部,然后利用多通道磁共振仪器同时进行探测。该方法可以对地下水进行二维探测成像,但是针对地下复杂含水构造仍然不能进行准确的探测,由于地面测点有限,其反演结果精度和成像平滑度都有待提高。若铺设多个探测线圈,可以提高结果精度,但是会增加探测工作量,降低工作效率。针对上述缺点,本文提出基于三次样条插值的地面核磁共振反演方法,该方法可以在提高反演精度和成像平滑度的同时,减少野外探测铺设线圈的工作量。本文首先对三维空间进行网格剖分,计算接收和发射线圈周围空间分布的磁场强度,计算正反演所需的核函数。建立地下含水体模型,正演计算得到E0-q曲线。然后对初始振幅E0进行三次样条插值计算,得到新的数据矩阵,再进行二维含水量反演成像,利用二分法搜索最优正则化因子,提高反演算法效率。最后,针对不同的地下含水体模型进行了正演仿真和反演计算,验证该方法的适用性。本文在进行正演时,提出利用三次样条插值方法对初始振幅E0进行插值计算,扩大数据矩阵,然后重新计算核函数使之与新的数据矩阵对应,通过对核函数进行SVD分解得到分辨率矩阵,对比插值前后的模型分辨率信息,验证基于三次样条插值的正演方法对浅层探测区域的高分辨能力。针对经过三次样条插值方法对各探测线圈的初始振幅0E进行插值后,数据矩阵和核函数矩阵扩大,导致反演计算量增大的问题,本文提出利用二分法搜索最优正则化因子?,不用计算每一个?对应的目标函数的值,可以大幅减少计算量,提高反演算法的效率。为了检验该方法得到的反演结果的精确度,本文进行了初始振幅插值结果和含水量反演结果的误差分析。在理想情况下,初始振幅插值结果均方根误差均小于7 n V,在加入20 n V噪声后,插值结果最大均方根误差为13.475 n V。在反演结果方面,基于三次样条插值的反演方法可以降低含水量反演结果整体均方根误差4个百分点。最后,使用该方法对长春市烧锅镇和安徽省黄山市阳台村的野外实测数据进行了反演解释,验证该方法的有效性。通过上述研究,本文得出基于三次样条插值的地面核磁共振反演方法能够在减少探测线圈数量的同时,提高反演精度和成像平滑度,实际运用中可以有效降低工作量,提高工作效率。此外,本文提出的方法运用在三维探测和数据解释方面更具实际意义,三维探测常用的测量模式是利用5个线圈在四个水平方向进行半覆盖形式铺设,在测量区域的对角线上移动测量,其探测线圈的增加数量是按平方形式增长,因此基于三次样条插值的反演方法可以大幅减少三维探测的线圈铺设工作量,并且得到平滑精确的地下水成像结果。