方位电阻率成像仪的三维有限元模拟.pdfVIP

方位电阻率成像仪的三维有限元模拟’ 杨铧 文亚凤 (石油大学，北京102249) (华北电力大学北京10呻B5) 擅蔓：为了研究方位电阻率戚像仪的定量解释方法，我们开发了一套麓横担方位电阻率成 像仅在三维非均匀地层摸型中响应的三维有限元程序。在简化的地层模型中，将谈程序与鲁折 方法和二维有限差分方法进行了对比，结果表明此程序的误差不超过5％。在此程序得到验证 后，我们用该程序考查了方位电阻率成悼仪在一些典型三维章均匀地层模型中的响应，结皋表 明此程序能为方位电阻率成像仅的定量解释提供支持。 美t胃：有限元方位电阻事 成像仅 非均匀地层 方位电阻率成像仪(ARI)是美国sclll一∞b瞎er公司90年代初推出的新一代的侧向测井 仪。它围绕井眼进行定向测量，利用12个方位电极可测得12条方位电阻事曲线，这些曲线 被显示成一个方位电阻率成像图。尽管这个全覆盖成像的分辨率比声波成像和徽电阻率扫描 成像低，但是对它们是很好的补充。利用这些图像可以形象地发现裂缝、溶洞等，但是，要 进行定量评价，就必须知道方位电阻率成像仪(ARI)在这些地层模型下的精确响应。因 此，有必要利用三维的数值模拟技术分析方位电阻率成像仪(ARI)的精确理论响应问题， 而这方面的研究至今设有见公开报道。目前可以见剜的文献有三种：第一种是叙述如何用三 维的有限元或有限差分方法模拟别的仪器如双翻向、阵列感应等在三维非均匀地层中的响 应，但相当简略u“1；第二种是在某些特殊的情况下将ARI的三维计算简化为二维计算[73； 第三种是提到用三维有限元计算ARI的响应．但只字不提是如何计算的，所给出的只是 ARI的井眼校正图版等I”…。在本文中，我们较详尽的描述了用三维有限元方法模拟方位 电阻率成像仪(ARI)响应的三维有限元方法。 1三维有■元方法 在三维的直角坐标系中，测井仪所激发的电场可以用下面的偏微分方程表示： c·， 未(吉髦)+毒(吉考)十麦(吉髦)=，*acz—z。，y—y。，z一砧●|}f．；11●I■， 式中z，v，z——直角坐标系中的三个坐标系变量； zo，，D，20——供电点电极的坐标； 卜函数； p——地层中的电阻率。 其中电阻率p为分块函数： ·国家自然科学基金资助项目(49帅4003)，中国石油天魅气集团公司剖新基盒赍助项目及中国石油天然气集团公司 。九五’嘎日。 f．；}：。；■{，} ·74· pm (z，，，z∈nl Px (z，y，g∈n2) Pt (z，，，z∈n3) D= (2) p叫 (z，，，z∈n4) p耐 (z，y，。∈n5) 其中nl，n2，n3，n4，如代表井筒、侵入带、原状地层、上围岩、下围岩所占的区域， P。，Pl，Pt，P辄，～，分别代表区域0l，n2，n3，m，口5内的电阻率。 利用变分原理可将上述的徽分方程边值问题转化为电位”的变分： cs， 圣c“，={』，[(髦)2+(考)2+(髦)2]如。如一莓咖。 其中口为电导率。IE，“E分别代表电极E的电流及电位。n为体积分区域。实际计算时， 由于此问题计算区域的特殊性，宜采用曲边六面体，因此需要采用等参元，于是上述变分问 题可转化为： ～虬，浏叫学学+学学+学警] xl I