过套管电阻率测井正演数值模拟方法研讨.pdf

过套管电阻率测井正演数值模拟方法研究 高杰‘、刘福平2、包德洲’，程希3，彭斐1 (1．中国石油大擘责舜与信息学院，北京102249；2．北京印刷学院，北京102600； l 3。中固石油集匝测井有限鸯司技术中心。陕西西安710021 摘要：本文对用积分方程法、有限差分法及传输线方法计算过套管电阻率问题进行了研究和考查，用三 种方法编程进行了计算．并对三种方法计算结果进行了比较．为利用传输线方法解决复杂模型问题．对经典 的传输线方程进行修正．得到所谓的修正的传输线方程法，并计算了复杂模型的测井响应问题。结果表明修 正的传输线方程法计算结果良好，最为适台计算过套管电阻率响应问题． 关键词：过套管电阻率测井；积分方程法；有限差分法；修正传输线方程法 一、引言 进展“…。1988年，美国的顺磁测井公司(PML)试验成功第一套过套管电阻率测井样机 用效果，得到油公司的认可。但是过套管电阻率测井仍有许多理论问题需要进一步研究和解 决，如通过漏失电流大小的数值模拟计算解决过套管电阻率实现的可能性和测量的可靠性、 正确性问题：通过讨论接触电阻、水泥环和套管等环境影响因素，解决测量数据的影响和校 正问题；建立针对过套管电阻率测井资料的解释模型．解决实际应用等问题。积分方程法己 被用于过套管电阻率数值模拟计算，并给出了三层层状地层的测井响应f8,14]，有限差分和有 限元法还未见到相关报道，目前所见到的文献多为基于传输线方程的相关讨论，并来涉及更 复杂的地层模型。我们对用积分方程法、有限差分法及传输线方法计算过套管电阻率问题进 行了研究和考查，用三种方法编程进行了计算，著对计算结果进行了比较．在此基础上，对 经典的传输线方程进行修正，得到了可信的结果。 二、有限长套管分层地层模型条件下的积分方程法(IEM) 1．含任意曲面的电势积分方程 使用曲面积分方程来描述半空间井眼流体、套管和纵向分层的非均质地层的套管井电阻 339 方程及其边界条件上，得到关于电位函数的电势积分方程18’121，即 u每，=鲁西，+莓警P，警，纠∞A㈤ 式中 G(只扔—Gmen函数；厶一电源电流；【，(局一空间尸点电势；吼—编号为f的地层 电导率；k一积分边界面。 2．电势积分方程的数值模拟计算 为使问题简化只考虑轴对称条件下的积分方程解，此条件下Green函数为∞，161 毋力=丢p∞_冰棚一”卅盯啦卸坳c：) 入式(1)可将二维积分转变为一维积分，并离散化方程，求解方程组可得电势空间分布。 三、复杂地层模型条件下的电势计算的有限差分法(FDM) 积分方程法计算问题有其独特的优势，即所涉及的未知量只分布在不同介质的边界上， 其未知量较其他方法大为减少。因此积分方程法一直受到大家的关注和应用。但从我们近年 来的计算发现，积分方程在解决某些问题时，特别是当矩阵元变化比较大时．边界元积分方 程组是病态的，这对问题的求解精度会有一定的限制，我们试图利用有限差分法(FDM)计 U—c■=0，所以可认为o_．v=0；在截面上，电流连续，在轴对称坐标系中， ’、 ’ 百02U+土罢+尝：三砸√)砸√)(3) 毋2 2勿 如2 盯、 不等间距网格差分方程为 ‰ [去+等卜南2 +i翻u，，一+i：五i2：r丽：+hru+·，，+j褊配．t．，+孑I 。 (4) ，“ 式中 虹wr方向的步长：尼厂-z方向的步长；％一点G力的径向半径· 四，修正的传输线方程法 若沿z轴方向套管是均匀或分段均匀的，地层