[发明专利]基于动态轮廓模型的地面核磁共振反演方法

说明书

技术领域

本发明涉及地面核磁共振(Surface Nuclear Magnetic Resonance，简称SNMR)领域，具体涉及一种基于动态轮廓模型的地面核磁共振反演方法。

背景技术

地面核磁共振(Surface Nuclear Magnetic Resonance，简称SNMR)技术是目前世界上唯一的一种直接找水的物探方法，该项技术已在探测地下水、考古、地下水污染检测等领域得到了一定的应用。

近年来，随着专家和学者们的逐渐深入研究，SNMR技术得到了进一步的完善。反演计算含水率是该技术研究过程中的关键环节，而反演结果的准确度和反演过程的运算效率是衡量反演算法性能的关键指标。其中，一维正反演理论较为成熟，已经相继刊登出多种有效算法，如：改进的模拟退火算法反演，提高了现有反演算法的稳定度和收敛速度；QT反演算法，利用各个激发脉冲矩对应的全部采样点数据进行反演，充分挖掘了接收信号信息，在一定程度上提高了反演精度，但是，由于接收信号呈现近似指数衰减，晚期信号信噪比很低，该方法只适用于高信噪比环境；后来也有学者采用了积分门技术接收信号，提高了各个采样点数据的精度，并进行全衰减反演，是对QT反演的一种改进。

SNMR技术在二维反演方面，Boucher、Girard和Legchenko等研究了在二维剖面方向上E0-q曲线随地下含水构造的变化趋势，但他们只对二维反演做了定性研究，没有给出具体的二维反演公式。Legchenko等对三维反演做了一定的研究，虽然能在三维空间反演出模型的含水构造，但是由于在三维空间设定的网格尺寸较大，只能粗略的估计出地下含水构造，其反演分辨率有待提高。由于二维、三维反演算法存在运算量大、待求解变量数多、非线性等问题，目前世界上唯一商业版反演软件NUMISPLUS仍采用一维反演。

地面核磁共振反演能抽象为一个求解病态矩阵方程的问题，现有的文献中多采用正则化方法对其进行求解，但是，现有的标准正则化方法以及在其基础上的改进方法存在两方面的不足，分别是：正则化参数选取困难；将反演问题的近似解强制添加光滑性约束，弱化了相邻地质层之间含水量的差异。

发明内容

本发明所要解决的技术问题是提供一种基于动态轮廓模型的地面核磁共振反演方法，其能够实现对地面核磁共振反演问题的高精度求解。

为解决上述问题，本发明是通过以下技术方案实现的：

一种基于动态轮廓模型的地面核磁共振反演方法，包括如下步骤：

步骤1，建立的地面核磁共振反演的动态轮廓模型，即

An＝E    ①

式中，A为核函数矩阵；n为待求解的含水量向量，M_N为模型中含水层的微元数；E为地面核磁共振信号初始振幅值向量，E=(E(q1),E(q2),...,E(qIN))T,]]>I_N为激发脉冲矩个数；

步骤2，初始化探测区域内的含水层数为L，各层的厚度和含水量值分别为H_L＝(h₁,h₂,...,h_L)^T和N_L＝(n₁,n₂,...,n_L)^T；最大迭代次数为N_max；初始迭代次数为N_I＝0；反演精度阈值设定为θ；