二维井间电磁场的正反演计算

提出了一种改进的逐次逼近解法 ,并用该方法对轴对称二维井间电磁场进行了正反演计算。与逐次逼近解法相比 ,该方法收敛性强 ,应用范围广 ,可适用于高电导率对比地层。与直接求解积分方程相比 ,由于不必进行直接的大型复矩阵求逆运算 ,该方法计算速度快 ,所需内存少。反演中采用了Born迭代算法 ,将成像区域集中于一定范围内。考虑到信息量和计算机内存的限制 ,采用双重面元分割法 ,将反演区域分割为一定数目的大面元 ,每个大面元被分割为几个小面元。属于同一大面元的不同小面元具有不同的磁矢势 ,但具有相同的电导率 ,从而减少了反演过程中未知量的数目。根据第一次成像结果确定出更准确的成像范围 ,并进行第二次成像。数值计算结果表明 ,改进的逐次逼近解法是计算二维井间电磁场的一种有效方法 ,将该方法用于反演过程能够得到较高分辨率的二维井间电导率图像。

用改进的逐次逼近解法计算二维井间电磁场

提出了一种计算积分方程的新方法——改进的逐次逼近解法 (MSAM) ,并用该方法计算了轴对称二维井间电磁场。与传统的逐次逼近解法 (SAM)相比 ,该方法收敛性强 ,应用范围广 ,可适用于高电导率对比地层。由于不必进行直接的大型矩阵求逆运算 ,因此与积分方程的直接解法 (IE)相比 ,该方法计算速度更快 ,所需内存更少。采用矩阵求逆方法计算了成层介质中的二维 Green函数 ,并对含 Green函数的积分进行了简化 ,从而加快了计算速度。数值计算结果显示 ,在地层电导率对比度达到 2个数量级时 MSAM仍收敛 ,且计算结果与直接求解积分方程的结果一致 ,因此