地面核磁共振理论模拟与反演方法初步研究.pdfVIP

F，=l～“K(Q，z皿． (6) 3地面核磁共振反演 31反演公式 对激发电流脉冲强度Q和深度z离散化后， (6)式有 Eo，=KⅣ胛，， (7) 其中nj为第．，含水层的含水量，并且0≤门，≤1。进一步可将上式写成如下的矩阵形式 Kn=E． (8) 上式等价于下面的优化问题 Wx 2， rain：0y-Ax旷+口2 S．，． (9) XI≥0， 式中a2称为正则化因子，W是体积为K×N的模型约束条件矩阵，取决于模型限制方式。 32提高积分核计算精度 、 在相对深部，地面核磁共振数据的反演精度与分辨力下降。究其原因，部分来自目前计算方法中，反演 含水层厚度随深度增加逐步增大，导致积分核函数计算精度在深部下降所致。为此，我们利用数值积分计算核 函数积分，其中积分节点处的核函数值由尽可能少的深度剖分节点上的核函数值通过插值插值方法求取。由于 插值方法的精度对积分有重要影响，我们选用精度较高的连分式插值。然后，基于插值结果，利用高斯求积公 式计算核函数的积分。 33利用投影迭代实现反演 优化问题(9)的数值实现有多种方法，最常用的是正则化反演技术。但由于正则化反演有时会得到没有意 义的负含水率，因此，这里采用负解迭代消去法，具体过程如下： (1)．利用正则化方法求解(9)式，获得谱分布的一个可能解； (2)．对每一个解分量，判断是否非负。若否，则删除对应于该分量的矩阵与右手边向量的行与元素。转(1)； (3)．输出反演结果。 34反演中变正则化参数技术 在常规的正则化反演中，正则化因子口‘是固定的，这显然没有考虑不同数据对模型参数的分辨能力。变 正则化地面核磁共振数据反演方法的基本原理是，假设正则化因子是各个含水层位置的空间交量，并且由模型 参数分辨率矩阵及其展布函数决定。对于高分辨率的参数，要求正则化因子值小，而分辨率较低的参数，其相 应的正则化因子被赋以较大的值。这样，在反演过程中，将根据参数灵敏度对优化问题中的含水量模型与数据 误差进行平衡，达到同时提高了反演的分辨率和反演过程的稳定性。具体的实现过程大致如下：首先由(9) 式用小的正则化参数定义模型参数分辨率矩阵足，并求出每个模型参数的空间展布函数SPi；接着由展布函数 计算变正则化参数，其公式为 五=砧in+瓦月m。a。x一--／％]mini。(SP／一SPrain)· (10) 式中九m、九。。、IS％in、．SRax分别为正则化参数和展布函数的上下限。最后，利用获得的变正则化参数 重新求解(9)式，获得最终的含水层模型。 本文的研究由国家自然科学基金和吉林省科技发展计划联合资助。 参考文献 l翁爱华等．层状导电介质中地面核磁共振响应特征理论研究，地球物理学报，44，1，2004 2 AVetai．InversionofsurfaceNMR Legchenk data，Geophysics，63，l，1998 · ·33 地面核磁共振理论模拟与反演方法初步研究 作者： 翁爱华 作者单位： 吉林大学,应用地球物理系,长春,130026;浙江大学,理学院,杭州,310027 引用本文格式：翁爱华 地面核磁共振理论模拟与反演方法初步研究[会议论文] 2005