核磁共振地下水探测理论及反演研究.pdfVIP

摘要     摘 要 核磁共振(Magnetic Resonance Sounding, MRS)地下水探测技术是目前唯一一种直 接探测地下水的方法，以其高分辨力、高效性和无损性等优点被应用于地下淡水资源 的寻找、地下水污染监测及滑坡等地质灾害监测中。本文主要围绕MRS 信号正反演 技术展开研究。近年来虽然 MRS 正反演理论得到快速发展，但是由于 MRS 信号本 身比较微弱(nV 级) ，加上求解含水量的系数矩阵是病态的，造成求解的含水量不精 确。为提高 MRS 地下水探测技术的准确度，本文对 MRS 正反演理论进行了深入研 究，主要的工作为： (1) 为高效的提取特征参数，本文在研究 MRS 勘探地下水原理的基础上，对 MRS 地下水探测的特征参数进行了分析，并采用基于最小二乘法的非线性拟合法提取 MRS 地下水探测数据的特征参数，结果与法国 NUMIS 系统提取的参数一致。 (2) 正演是反演的基础。本文从正演的数学模型入手，推导简化模型下 MRS 正 演计算公式，计算均匀半空间内单个和多个含水层模型的 MRS 响应，并依次从导电 性、含水层参数(含水量，含水层深度，含水层厚度)和平均衰减时间三个方面对影响 MRS 响应的因素进行详细分析，得到以下结论：含水层上方覆有小于100 m 良导 层时，大地电阻率值对初始振幅E (q ) 影响较大；含水层的含水量、深度、厚度的增 0 大，都会引起初始振幅E (q ) 的幅度明显变大，且深度越大需要激发脉冲矩也越大； 0 * 平均衰减时间T 越大，越有利于探测 MRS 信号。 2 (3) 本文中 MRS 反演主要求解含水量，由方程A n = E0 可求得。但是由于A 是 病态的，通常采用奇异值分解方法(SVD)、截断奇异值分解方法(TSVD) 以及Tikhonov 正则化等方法求解。由于采用奇异值分解(SVD)求解的反演结果不稳定且出现负数 解，而采用基于共轭梯度的Tikhonov 正则化反演结果虽然稳定，但是过于依赖初始 模型，容易陷入局部最优，难以获得全局最优。针对这些问题，本文提出了改进的模 拟原子跃迁算法(RSATA)，该算法是一种结合正则化的随机搜索算法，不依赖初始模 型。它将原问题转化为一个带非负约束的最小二乘问题，先用 SVD 法获得局部最优 解，而后通过模拟原子跃迁法获得全局最优。相比改进的模拟退火算法，该方法收敛 速度平均提高 57.1%，准确度提高23.2% ，信噪比的要求降低 50% 。 (4) 根据研究结果，完成样机版 MRS 工程应用软件的设计，经过野外实测数据 验证，该软件满足实际探测需求。 关键词：核磁共振反演；非线性拟合；Tikhonov 正则化；模拟原子跃迁法 -I- 万方数据 Abstract Abstract Magnetic Resonance Sounding (MRS) for underground water detection is the only method for detecting groundwater directly. This method is applied to f