基于COMSO仿真软件的阵列感应测井偏心响应计算方法

引言

阵列感应测井是重要的电阻率测井方法之一，其优点是分辨率高、对侵入反映明显、探测深度深、测量信息丰富。 由于阵列感应测井仪有多个提供高分辨率的短子阵列，仪器偏心严重影响测量结果，尤其是大井眼和低电阻率泥浆情形，导致后续信号合成处理出错，产生错误的解释结论，因此必须进行偏心影响校正。

1、井眼模型建立、网格剖分和求解器选取

1.1建立偏心模型阵列 感应测井仪MIT的偏心模型见图1。

地层由井眼和地层组成，带有偏心块的仪器轴心偏离井眼中心。建模时井眼和地层是等高的圆柱，根据趋肤效应计算地层的大小，由电磁场的传播特性考虑地层为球形，仪器是处在井眼偏心状态的有限高度圆柱。图1给出了偏心模型剖面图和横截面图。

整体剖分效果图
1.3求解COMSOL仿真软件在求解时提供了多种迭代方式，
如广义最小余量法（GMRES）、灵活的广义最小余量法（FGMRES）、稳定的双共轭梯度法（BiCG－Stab）和共轭梯度法（CG）。
求解中发现默认迭代方式GMRES在泥浆和地层电导率对比度大时迭代收敛速度慢或者不收敛，会出现停滞现象。这里选用其他3种迭代方式。
实验发现当仪器内电导率取０时，保证一定精度情况下，CG优于其他２种方法，比BiCG－Stab所用内存空间小，比FGMRES迭代次数少。当不考虑仪器电导率，即仪器电导率和泥浆电导率相同时，FGMRES优于其他2种方法，比CG稳定，比BiCG－Stab所用内存空间小。
因此，考虑仪器响应时用CG计算，不考虑仪器响应时用FGMRES计算。
2、计算数据的比较及分析为准确验证模型建立的正确性以及网格剖分合理性，首先计算均匀地层模型，将结果与解析解比较。取均匀地层电导率为 0.1S/m，地层半径40m，计算结果见表1。

（1）子阵列6的偏心视电导率整体小于居中视电导率，地层电导率不变，随着泥浆电导率增大，偏心与居中的视电导率差异增大。
（2）子阵列7当地层电导率大于泥浆电导率时，偏心视电导率大于居中视电导率，反之偏心视电导率小于居中视电导率，随着泥浆电导率增大差异变大。
（3）比较2图，子阵列6和7受偏心影响完全相反，前者是负差值，后者是正差值。2个子阵列受偏心影响的共同点是随泥浆电导率与地层电导率的比值增大，差值增大，即偏心影响严重。这些现象表明，偏心时仪器测量到的信号与居中时有较大差别。

结 论
（1）通过建立合理模型，针对具体问题实现网格剖分，选择有效的迭代方式使得解稳定有效，将数值解的相对误差控制在一定范围内，达到所需的精度要求，用COMSOL软件可以准确计算阵列感应测井仪MIT的偏心响应。
（2）对阵列感应测井仪MIT建立偏心响应数据库，对偏心严重的响应绘制图版进行分析，得出不同子阵列偏心下的影响，可以有效分析偏心时对视电导率的影响。
（3）利用本文研究的方法可以将COMSOL仿真软件推广到感应类仪器其他问题的分析与解决，为解决此类数值计算问题提供有效途径。

来源：YHJX57