计算阵列感应测井响应的积分方程法.pdfVIP

46 内蒙古石油化工 2005年第4期 计算阵列感应测井响应的积分方程法 孙晓霞，申辉林 (石油大学(华东)资源与信息学院) 摘要 开发了积分方程法计算阵列感应测井仪器各子阵列在轴对称二维地层模型中的响 应。积分方程法的计算区域只是二维电导率异常剖面而不是整个剖面，可以大大提高计算精度。以 纵向成层地层作为背景地层(即原状地层)，以所推导的纵向成层背景地层的Green函数为基础，考 虑到模型中存在井眼和侵入带的情况，采用积分方程计算出接收线圈在这种地层模型中的感应电 动势和实、虚部视电导率。 关键词 阵列感应；格林函数；积分方程 引言 2，以以，N。每个小单元内地层电导率恒定，设为盯，，并 利用积分方程计算轴对称二维电磁场时，由于 假设在每个小单元内磁矢势恒定，由其中心点的值 数据点只分布在电导率异常区域而不是整个子午 代替。由此得到。内(1)式的离散形式： 露。’珂丽与有限元素法和有限差分法以及混合法相 ；艺数据量相对较小，适合于进行反演计算。 i计算二维积分方程的改进型逐次逼近解法 dr7]A(rj，仰)，(i一1，2，K，Ⅳ) (3) 对于纵向成层原状地层中存在井眼和侵入带的 下面讨论采用改进型逐次逼近解法(MSAM) 轴对称二维模型，可将原状地层作为背景地层，采用 计算积分方程(3)。在异常区域内的散射电流密度定 积分方程计算电导率异常区域内(井眼和侵入带)各 义为： 点的磁矢势以及接收线圈处的磁矢势，从而得到接 (4) 收线圈在这种地层模型中的感应电动势和实、虚部 所以(1)式可写为： 视电导率。电导率异常区域内各处磁矢势(只有分 f A。(r，卯)一胁J 量)用积分方程可表示为[1]： (5) f A(，．，卯)一A6(r，rT)一ioL’／10Jnr(r，r『)Aa(一)A (r’，FT)dr7 (1) 三维电导率异常体内，下列关系式成立：上如 式中：A6为磁矢势的背景值，n为电导率异常 区域即井眼和侵入带，r(，．，r!)为纵向成层背景地层 E4+去忙』饼咖 ㈤ 的Green函数，可由(1)式得到，知(r7)一a(r’)一巩 (，．’)，a(r’)和ab(r’)分别为电导率异常区域和背景地 dz，对角度的积分可以先积出来。(7)式变为： 层的电导率。 由(1)式可形成如下迭代序列： 2《r“l--iooA*q-2f瓦”b2d以z≤2《，．≮≠d以2 f A¨’(r，rT)一A6(r，rT)一itopoJnr(，．，r’)△盯(r7) (8) A‘一1(r7，rT)dr7，(志一1，2，以) (2) (8)式对任意大小的异常区域n均成立，故有： 这就是文献[1]的逐次逼近解法，积分过程可通 过将电导率异常区域力离散实现，参见(3)式。采用 锄卜印+2铲ab2--脞4ab (9)