各向异性介质多分量感应测井三维Born几何因子理论
各向异性地层性质及测井响应
各向异性地层性质及测井响应1. 均匀各向异性介质中的电阻率测井响应均匀各向异性介质中电偶极⼦和磁偶极⼦视电阻率表达式,即普通电阻率测井和感应测井的测量结果为:a R R =gm R =2V H H VR R σλσ==为各向异性系数;θ为相对地层倾⾓;对于层状地层,垂向电阻率总是⼤于⽔平电阻率:V H R R ≥,因此,各向异性系数α通常总是⼤于1的。
下⾯给出两种特殊情况的结果:1)对于θ=0的特殊情况,即直井情况,H a R R =;2)对于θπ=/2(90度)的特殊情况,即⽔平井情况,a R =。
因此,在有倾⾓的各向异性地层中,普通电阻率测井或感应测井仪器反映的是地层垂向电阻率和⽔平电阻率的加权平均:由0度时的H R 变化到90度时的从物理机制看,感应测井在直井中的涡流是⽔平⽅向的,因此仅得到⽔平电阻率,⽽在斜井中由于涡流存在于两个⽅向,所以其读数为垂向和⽔平向电阻率的平均值。
在直井中,感应测井只是反映地层⽔平电阻率,普通电阻率测井或侧向测井也主要反映地层⽔平电阻率,这是常规电测井在反映各向异性⽅⾯的局限性。
2. 各向异性指数系数[2,4,12]在各向异性地层中,电阻率测井的响应还与井下仪器的结构有关,不同测井仪器测量出的视电阻率之间往往存在明显的差异。
如在泥岩层中,感应测井仪测量的视电阻率明显低于0.4m 电位测井仪;在砂泥岩互层,梯度测井仪测量值异常地低。
因⽽利⽤各种电测井仪的响应差异可识别电阻率各向异性地层。
在垂直井眼中,假定地层是⽔平的,砂泥岩薄互层、不同粒度⼤⼩的砂岩层、岩层中薄层的电阻性或电导性条带等都使地层表现为各向异性。
各向异性指数主要与砂泥岩电阻率反差程度和砂泥岩相对厚度有关。
图1是砂泥岩互层⽔平电阻率、垂直电阻率与各向异性指数关系图。
图1中R sh =1.0Ω·m ,Rsd=10.0Ω·m 。
h sh =h sd ⽬处是各向异性指数2λ最⼤的地⽅,约为3,此时R h 约为1.8Ω·m ,⽽R v 约为5.5Ω·m 。
泥浆侵入电阻率各向异性对多分量感应测井影响分析
J u n l fOi a d Ga e h o o y J n 2 1 Vo. 4 No 6 o r a l n s T c n lg u . 0 2 o 13 .
பைடு நூலகம்
21年 6 02 月 第3卷 第6 4 期
・ 5 8 ・
泥 浆 侵 入 电 阻 率 各 向 异 性 对 多 分 量 感 应 测 井
一
二 T ( 一 Z ) 】 x ( i s e p - ih 2 . )2 T r ep - k ) E x ( kR)
( 4)
引 入 Het 势 , rz 与电场 的关 系为 : rz Het 势
( 5)
X、 Z方 向上 的单位 磁偶极 子在 均匀 各 向异 性介 质 中产生 的 Het 势可 分别 表示 为嘲 : Y、 rz
(, r = = = A
——■ 一
( r ) 一 — p( v ) r, e i S x - k
的磁 偶极 子 。 收 线圈 中测得 的磁 场张 量 再 为 : 接
F  ̄ H x H X ] H x y Z
H = x Hy H I = Hy :l y
L z H H x
[ 收稿 日期 ] 2 1 0 2—0 —2 3 o
() 3
H z Z
[ 基金项 目] 国家 自然科学基金项目 ( 0 70 3 。 4 7 4 7 ) [ 作者简介]刘迪仁 ( 9 5一 ,男 ,18 16 ) 96年江汉石油学院毕业 ,博士 ,副教授 ,现主要从 事电法测井正反 演、复杂储层 测井评价 及 光 纤 传 感 技 术 等 方 面 的 理论 和应 用研 究 工 作 。
・
8 6・
石 油 天 然 气 学 报
最新钻井地球物理勘探教案——第四章 感应测井
第四章感应测井感应测井可在井眼不导电的情况下(如油基泥浆井,空气钻井等)测量地层的电导率。
这种方法对低阻层反应灵敏,因此更适合区分低阻油、水层和油水过渡带。
第一节感应测井的基本理论一、基本原理感应测井是利用交变电磁场研究岩石导电性的一种方法。
发射线圈 T ,通以 20kHz 交变电流,该电流在周围介质中产生交变电磁场中。
φ 1 在介质中适应出环形电流 i 1 ,同时在接收线圈 R 中,产生感应电动势 E 1 。
环形电流 i 1 ,在介质中亦将引起二次磁场φ 2 ,φ 1 在 R 中引起适应由动势 E 2 。
φ 1 在 R 中引起的电动势为无用信号,而φ 2 在 R 中引起的感应电动势 E 2 与 i 1 有关, i 1 又与地层导电能力有关,因而,通过测量 E 2 的大小,便可测量介质的导电能力。
在均匀无限介质的条件下,通过求解电磁场的基本方程可得出,接收线圈中,总适应电动势的表达式:该式展开后,可简化为:上式中,虚部是无用信号,实部与σ成正比,是有用信号,二者相位上差 90 °。
这就是感应测井的基础。
上式的得出是忽略了三次方以上的高次项的结果,是忽略了趋肤效应影响的一种近似方法。
这样就可把有用信号看作是介质各部分所引起的感应由动势线性相加的结果,这种方法就是“几何因子”理论。
几何因子理论要点:①认为发射电磁场与每个单元环电磁场之间互不发生作用(即幅度衰减和相位移动场可忽略)②认为电磁波瞬间便可通过地层,(而实际地层中电磁波传播速度仅为自由空间的 0.15% )。
根据几何因子理论,得到的接收线圈中的有用信号为:dE 2 = kg σ· ds几何因子 g 的物理意义:在均匀无限介质中,任意一点上截面积为一个单位的单元介质环对总信号的相对供献。
二、均匀介质双线圈系感应测井的电磁理论1 .传播效应(趋肤效应)2 .麦克斯韦方程组及其解3 .感应测井信号的虚、实分量第二节感应测井线圈系特性空间各部分介质对总的感应电动势贡献大小是由每部分介质的电导率与它的几何因子两部分因素决定的,因此,必须研究几何因子的空间分布,才能研究各部分介质对感应电动势的贡献,而几何因子的空间分布与线圈系结构有关,因而必须研究线圈系的特性。
感应测井几何因子理论计算及其应用
感应测井几何因子理论计算及其应用彭玲;潘和平【摘要】几何因子对于计算地层真实电阻率、了解地层特性、划分储层等起着相当重要的作用.这里建立双线圈系感应测井基本模型,分别得出纵、径向简单多层地层模型的几何因子理论曲线,从中可得知感应测井响应值受到目的层厚度、围岩电阻率、泥浆电阻率、侵入带等多种情况的因素影响,①目的层越厚,感应测井值越接近原状地层电阻率值;②在过渡带厚度不变的情况下,冲洗带深度越大,测量的电阻率越小,它的径向几何因子有下降的趋势,但是整体来说,随着冲洗带的半径增加,视电导率增加幅度变缓,且越来越接近过渡带地层的电导率;③在冲洗带厚度不变的情况下,过渡带的深度越小,越接近原始地层的电导率;④泥浆侵入会导致得到的地层电阻率与实际地层电阻率存在较大偏差.将这些结论应用于计算塔河油田石炭系区域井的真实电阻率中,得到塔河油田石炭系区域地层电阻率真值,最后解释结果与试油结果对比证实,这一过程确实提高了地层电阻率和冲洗带电阻率的精度.【期刊名称】《物探化探计算技术》【年(卷),期】2016(038)001【总页数】11页(P8-18)【关键词】双感应测井;几何因子;塔河油田石炭系区域【作者】彭玲;潘和平【作者单位】中国地质大学地球物理与空间信息学院,武汉 430074;中国地质大学地球物理与空间信息学院,武汉 430074【正文语种】中文【中图分类】P631.8几何因子可以很好地解释周围地层对测量信号的贡献,从而可以计算得到地层的真实电阻率,所以研究几何因子对于了解地层特性、划分储层等起着相当重要的作用。
Doll[1]几何因子是Doll在1949年提出感应测井方法的同时就给出了的,它是从电磁场方程在忽略趋肤效应的影响的条件下导出的,因存在一定不足,此后有学者相继提出新的几何因子, Gianzero 等人[2]提出了考虑趋肤效应的几何因子-Gianzero几何因子; Moran[3]推导出的Moran几何因子,指出其具有对称性和实际物理意义;Howard 等人[4]将量子散射问题中的 Born 近似方法应用于非均匀地层中的感应测井问题,导出感应测井中的Born 近似公式。
岩石物理、井眼环境和地层几何参数对大斜度井多分量感应测井测量的影响
-. 001
0
00 .1
一
Hz 正 交 ( m) z A/ 各向同性 一一各向异性
_ 6
一
_ ‘ 、
' 4 - 2 0 2
H x 正交 x
一
( m) A/
各向同性
一 一各向异性
和侵入对不同测量分量的相对影响。接着完成了围 岩各向异性 、 各种倾角情况下多分量感应测量对侵
商用化 , 可用于提供更多的测量 , 这些测量对于精确
描述 近井眼地层各 向异性性 质是必不 可少 的 。
已有大量文献报告了复杂井眼环境同轴阵列感 应测量的应用 , 但很少讨论环境 和岩石物理性质对 多分量感应井眼测量影响。Tm k s o pi 等人分析了井 n 眼和侵入对大斜度井多分量感应测量 的影响, ag W n 和 Fn 等人研究了大斜度井多分 量感应仪器的响 ag 应 , ag W n 等人还研究 了无限厚地层井眼性质( 泥浆 电导率 和仪器偏 心 ) 的影 响。
各 向异性存 在 的情 况下 多分 量 感应 测量对侵 入 带范 围的灵敏 度 , 们作 了进 一步 的分析 。最后 , 我 我 们研 究 了水平 井 中泥浆滤 液产 生的 非均 匀侵 入情 况下仪 器的响 应。 数值模 拟表 明, 层 两边 围岩 为各 向异性 时 , 影响很 大 , 砂 其 即使 井眼倾 角较 小 ; 砂 岩层 中点进 在
型完成 了一 系列模拟 , 以检 查井 眼 、 围岩层各 向异性
宏观 电各 向异性 的存 在 , 如薄层状砂 一泥 岩和砂
一
砂 岩层序 产生的各 向异性 , 能极 大 地影 响井眼常 可
规感应测量, 使油气饱和度 的评价不正确 。此外 , 大 斜度井和水平井泥浆滤液的非对称和非均匀侵入也
感应测井
二 、感应测井线圈系的选择
双线圈系(一个发射线圈,一个接收线圈)的探测特性. 复合线圈系( 多个发射线圈,多个接收线圈)的探测特性. 感应测井仪器(线圈系)的探测特性(探测深度和纵 向分辨率)是评估感应测井仪器的重要指标,它与地层
参数(如冲洗带电导率、原状地层电导率、围岩电导率)
分布有关。现在流行的横向探测深度和纵向分辨率的定 义是基于道尔(Doll)的几何因子理论发展起来的。
s
围岩
t
s
i
m
侵入带 围岩
原状地层
分区均匀介质地层模型
无用信号Vx(发射线圈和接收线圈间的直接耦合电动势):
VX
2 i nT nR S0 I 3 2 L
将VR与Vx相比,并带入常数得
| VR | 8% | VX |
在测井中,应把VX 消除掉。为消除无用信号,一般采用 下面两个措施: 1) 采用复合线圈系,有意识地压制无用信号。 2) 利用VX与VR间90°的相位差,可在线路中加入相敏检 波器来进一步压制。
2. 复合线圈系
复合线圈系是由串联在一起的多个发射线圈和串联在
一起的多个接收线圈所组成。它们分别用 T0 , T1 , …… , TL 和 R0 , R1 , …… , Rm 代表 , 其匝数分别是 nT0 , ……, nTL和nR0,……,nRm,其中T0和R0称为主发射和主接收线 圈,它们的匝数 nT0和nR0)一定是最大的
围岩影响较大 薄层探测不准
定义纵向积分几何因子等于0.5时对应的层厚为纵向分辨率。
Cs
Ct
Cs
一个h=1m,Ct=200ms/m的水平地层,在Cs=400ms/m的围岩 中,求线圈距为1m的双线圈系在地层中点时的视电导率值和 围岩贡献的大小? Ca=Ct*Gt+Cs*(1-Gt) 通过查图,可知Gt=0.5 Ca=200*0.5+400*0.5=300ms/m 围岩的贡献率: Cs*(1-Gt)/Ca=400*0.5/300=66.7%
地层方位各向异性
地层方位各向异性第l1卷第6期圉虾曲乞勘摞747④谣方位各向异性Bertchenin等着Ch,enlnJB主题词方位各向异性交叉偶极横波分裂裂缝你对钻井过程中的井眼稳定性存在疑问吗?产层具有裂缝渗透率吗?裂缝方向是怎样的?局部应力状态是怎样影响增产计划的?厚殖,衽叶,l寸,所有这些问题都有一个共同点.要回答这些问题,就需要测量地层的方位各向异性自1927年在法国首次采用电缆测井以来,井眼中记录的大多数信息都表明地层特性是随井中的深度或井旁的深度而变化的.这个数据是一维的,因为只有深度是可变的.当今复杂的油藏评估需要更广泛地反映地下情况.不仅要弄清物理性质是怎样随深度变化的,而且还要掌握井眼附近这些特性的方位分布本文介绍了测量地层方位各向异性的创新技术,这就是交叉偶极声波测量法.探测仪探测仪除了具有一个单极发射源和一组单极接收阵列以外,还包括能排成一线或排成互相垂直的两个偶极声渡发射器和一组8个偶极的接收器.偶极发射源及接收器发送或接收与其表面垂直的质点运动.每个偶极接收器记录着称作轴和Y轴的两个水平,正交轴上的数据.实际上,每个X,Y轴方向都有一个发射器和8个接收器.仪器简图如图1所示对于8个当中的每一个接收器,测井仪记录4个分量,这4个分量为:删——方向作为发射源,方向接收器记录的波形皿L_一方向作为发射源,Y方向接收器记录的波形.n._一Y方向作为发射源,Y方向接收器记录的渡形.——Y方向作为发射源,方向接收器记录的波形.图2示出了由8个接收器当中的一个得到的4个单一分量的数据集.取决于岩石的压实作用,探测深度是1R一4n,这个穿透深度是很重要的特征.与之相反,以往研究的大多数技术都是依据探测深度很浅的电阻率数据或井下声波成像仪获得的数据来估算各向异性的.横波和各向异性在各向同性介质中,偶极源产生的挠曲渡(横渡)沿井眼传播,且质点运动与传播方向垂直.让我们假设一个各向异性模型是由南北向的一叠垂向平板建成的(图3).这个简单模型代表裂缝性地层.平行于平板方向的横渡速度()较快,而垂直于平板方向的速度()较慢.当偶极源和接收器顺着平板的方向(南北向)排列时,质点运动方向与平板是平行的,并且*收稿日期:1999年6月1日.原文有图8幅.编时删去其中1幅,图号重新排序.ZL,CPT泽荣宝牛748国外油气勘探10f3图1撂测杖.图中显示了发射器,隔离器厦接收器部分,该仪器直径3.88i~.,可在直径45in一21in的井中使用.谊杖器适台于在400.F的温度,:~0000psi的压力下剐量999缸图28个接收器之一记录的波形.该图显示了4个分量:,,W厦YX图3描述各向异性介质中横渡性能的简单模型图4措波的分裂.产生的横波分裂成两个子波,较?的子波与裂缝平行,较?l蔓的子波与裂缝垂直第11卷第6期地层方位各向异性"裂缝"是看不见的.而测井攸仅测出.然而,当偶极源和接收器与平板方向垂直时,质点运动就会受到每个板之间界面的影响,这时仪器测出慢横波速度(K).当偶极源既不平行也不垂直于"裂缝"方向时,原先的挠曲横波分裂成两个子波.一个子渡以快横波速度()传播,并且质点运动与裂缝方向平行;而另一个子波则以慢横波速度()传播.这种现象称作横渡分裂或横渡双折射(图4).交叉偶极数据处理处理交叉偶极数据的常规方法是利用Afford(1986)描述的分析法将4分量偶极阵列(,皿,,YY,YX)旋转成两个主波阵列,然后单独处理每一个单一的波列,从而获得快,慢横渡速度称作Afford旋转法的该项技术还应用于垂直地震剖面(VSP)这项技术的主要缺点是,当速度分析误差值达到快,慢横波速度的差值时,用:Mford旋转法确定方位有90~的双值性这要求我们开发出新的处理方法Tang和Chunduru(1997)研究了一种能降低方位双值性,又能准确测量各向异性参数的新的反演技术.新方法是根据各接收器之间的波形拟合进行的.根据这种仪器的构型,能在每个深度点对32对接收器进行估值.在对数据集中所有接收器对接收的时间和振幅延迟进行多次对比之后,再应用迭代全局最小值的处理方法,Tang-Chundum方法同时得到快速方向和,值因此,该项新技术不受与Afford旋转法有关的方位双值性的影响.美国中大陆区应用实例我们阐述该探1测器在美国中大陆的应用实例.常规测井分析表明该区有70ft厚的含少许白云岩的钙质砂岩.孔隙度范围4%~8%,层段上部孔隙度较高该层顶部一个小泥饼的形成表明层段上部的渗透率有限.由于该砂层孔隙度和渗透率较低,因此不像是一个好的远景区为了证明该测量结果具有重复性,用探测器在相同层段记录了3次.图5表明了两个轴向分量XX和的基本处理过程.左边记录道的曲线是伽马射线测井曲线2号道(左边第二个道)示出了分量方位该实例中,探测仪在砂层段旋转360~.3号道是XX和IT分量的视横波速度曲线.视速度随探测仪在井下的旋转而变化,这清楚地表明了地层的方位各向异性.然而,在该处理阶段,由于测量值是有限的,故没有指出各向异性值和方位.图6显示了利用Tang和Chundum(1997)的反演技术进行全交叉偶极处理的结果. 所有3次都经过处理,表明了处理的稳定性.第1道显示了伽马射线测井曲线和仪器方位.在所有这3次中,仪器每50fi旋转一圈第2道的曲线表明了地层的慢与快横波速度以及井孔XY方向的直径.第3道的曲线,称作ISO.s1,是一个反演过程的有效性标志.反演过程中,同时考虑了各向同性和各向异性模型.当快横渡分量的模拟和测量的速度拟合较好时,ISO—s1的值为正,否则为负.当IS0一S1为负值时,推测地层是各向同性的.第4道显示了各向异性程度和由多次计算得出的标准偏差值右边的最后一道是各向异性方向(或裂缝走向)以及根据3次测井计算出的标准偏差.该特殊实例中,各向异性位于砂岩段,并且总是东西向.特别强调的是,各向异性方向与仪器的实际位置无关.尽管测井当中探测仪总是快速旋转,但各向异性的方向仍保持不变.这种恒定性证实了该过程的有效性.以更多的"用户友好的"方式显示信息是可能的图7的显示中,图反映的是井眼附近横5'50国外油气勘探999矩■畦船譬科艇薯三uq霹辞匠霉譬田惟智田惦辉譬赳.崃定雌崃譬掣霹蝌鲁丰面譬掣纂寒蓉k哪奄嫜÷窄恬田星星广,<)'∽P-,,一,宕—量一置霉毒昌惜0:}昌鲁:厘趔噬:柽蘑霉厘牝暑虽莩C§ll\,/.,—一窖^-一,r,斟亲羞.●第11卷第6期地层方位各向异性75l图7旋转360~所显示的信息图一图中色彩的变化表示井眼周围横波速度的变化波速度的变化情况.这种情况下,快横渡速度指向东西向各向异性与平行于快横波速度的的裂缝系统有关.裂缝的显示表明该砂层是一很有前景的层段,即使岩石渗透率较低,从该地层借助裂缝也可产出有经济价值的流体.译自TheLeod/ng&,,T(.1.18,No.2,1999译者简介牛宝荣女.助理工程师.1962年生,自学英语大专毕业,一直在吐哈研究院从事情报工作曾发表译文数篇.(本文编辑宋焙)。
《地球物理测井》-课后思考题
思考题第一课自然电位测井SP?*1.分析自然电位的成因,写出扩散电动势、扩散吸附电动势、总电动势表达式。
答:自然电场的产生(原理)扩散电动势、扩散吸附电动势、过滤电动势1.扩散电动势产生原因:泥浆和地层水矿化度不同——电化学过程——电动势——自然电场产生过程:溶液浓度不同——离子扩散——离子迁移率不同——两边分别富集正、负离子 (延缓离子迁移速度)——产生电动势(直到正负离子达到动态平衡为止 ) 公式:2.扩散吸附电动势产生原因:泥浆和地层水矿化度不同——产生阳离子交换——产生电动势——自然电场产生过程:溶液浓度不同——带电离子扩散——阳离子交换——孔隙内溶液阳离子增多——浓度小的一方富集正电荷,浓度大的一方富集负电荷产生电动势(扩散吸附)公式:3.过滤电动势产生原因:泥浆柱与地层之间的压差造成离子的扩散。
一般在近平衡钻井情况下不考虑。
总电动势公式:*2.不同Cw、Cmf情况下自然电位测井曲线有哪些特征?1.当Cw>Cmf:(Rmf>Rw,E<0)负异常(淡水泥浆)2.当Cw<Cmf:(Rmf<Rw,E>0)正异常(咸水泥浆)3.当Cw=Cmf:(Rmf=Rw, E=0)无异常,自然电位测井失效*4.自然电位测井曲线在油田勘探开发中应用于哪些方面?1.划分渗透层(半幅点法,砂泥岩剖面较常用)2.估算泥质含量3.地层对比依据: 1)相同沉积环境下沉积的地层岩性特征相似; 2)同一段地层有相同或相似的沉积韵律组合; 3)由1)和2)决定同层、同沉积(相)的SP曲线特征一致。
4.确定、划分沉积相5.确定油水层及油水界面(△USP油小于△USP水)6.识别水淹层(依据 Cw <或> Cwz) 渗透层水淹后SP基线偏移,偏移量与Cw/Cwz(注入)有关7.确定地层水电阻率Rw3.影响自然电位测井的因素有哪些?1.Cw/Cmf影响(地层水矿化度/泥浆滤液矿化度)当Cw>Cmf:(Rmf>Rw,E<0)负异常(淡水泥浆).当Cw<Cmf:(Rmf<Rw,E>0)正异常(咸水泥浆)当Cw=Cmf:(Rmf=Rw, E=0)无异常,自然电位测井失效2 .岩性影响砂泥岩剖面泥岩(纯泥岩)——基线纯砂岩——SSP(h>4d)当储层Vsh 增大,自然电位幅度△USP(变小)<SSP 靠近泥岩基线3..温度影响温度对离子运动,离子扩散速率有影响不同深度地层温度不同4.地层水、泥浆滤液中含盐性质影响(溶液中离子类型不同,迁移速率不同,直接影响Kd、Kda)5.地层电阻率影响(当地层电阻率较大时,其影响不容忽视。
三维各向异性油藏水平井产能新公式
三维各向异性油藏水平井产能新公式贾晓飞; 雷光伦; 孙召勃; 姚传进【期刊名称】《《油气地质与采收率》》【年(卷),期】2019(026)002【总页数】7页(P113-119)【关键词】水平井; 产能预测; 三维各向异性; 线性流; 新公式【作者】贾晓飞; 雷光伦; 孙召勃; 姚传进【作者单位】中国石油大学(华东)石油工程学院山东青岛 266580; 中海石油(中国)有限公司天津分公司渤海石油研究院天津 300459【正文语种】中文【中图分类】TE32+8水平井能够有效提高生产井段与储层的接触面积,增大泄油半径以及渗流面积,大幅度提高单井产能和采收率,广泛应用于实际油田开发中[1-5]。
水平井产能是水平井部署决策的重要技术指标,而目前常规水平井产能公式只考虑了外部平面径向流和内部垂向径向流,预测结果往往偏高[6-13]。
事实上,水平井的渗流机理较复杂,不仅包括外部平面径向流和内部垂向径向流,还包括中间平面线性流[14-15],传统方法预测结果偏大的重要原因就是对中间平面线性流的忽视。
近年来,诸多学者充分考虑这3种渗流形态,推导建立了均质油藏的综合形式的水平井产能预测公式[16-19],也有学者通过平面和纵向的二维异性对公式进行了校正[20-21],并没有考虑平面的各向异性。
大量研究和开发实践表明,不同沉积环境下,由于物源方向、水动力条件不同,沉积过程中岩石颗粒大小不同且其空间排列方式各异,从而导致油藏不同方向的渗透率不同[22-28]。
在对水平井产能进行计算评价时必须考虑各向异性的影响,尤其是平面各向异性与水平井部署方位的配置关系对水平井的产能影响较为显著。
为此,笔者通过坐标变换将三维各向异性油藏转换为等效的各向同性油藏,在此基础上,利用水电相似原理和等值渗流阻力方法,推导出了考虑油藏三维各向异性和平面线性流的水平井产能计算新公式,可以对三维各向异性油藏不同方位部署的水平井的产能进行预测。
1 三维各向异性等效转换假设一个三维各向异性油藏上、下边界封闭,主渗透率方向与坐标轴方向平行(图1),油藏厚度为h,泄油半径为r e,x,y,z方向的渗透率分别为Kx,Ky,Kz,水平井长度为L,井筒半径为r w,与坐标轴x夹角为α。
三维非对称介质电阻率各向异性反演及应用
三维非对称介质电阻率各向异性反演及应用吴意明;张伟;刘保银;张中庆【期刊名称】《测井技术》【年(卷),期】2016(040)005【摘要】针对海上水平井和大斜度井电阻率测井因受泥浆侵入、地层各向异性等因素影响,视电阻率会偏离真实地层电阻率的现象,提出了三维电阻率反演方法,旨在获取准确的地层电阻率剖面和井眼地层位置关系。
利用哈里伯顿公司ADR方位电阻率测井数据,在正演仿真基础上采用基于M arquardt方法的迭代反演方法计算测井时刻的泥浆侵入半径和侵入带电阻率(有侵入地层)以及原状地层水平电阻率、垂直电阻率、层边界距离(水平井)等参数。
数值模拟结果显示,反演方法精度较高,反演结果准确可靠。
【总页数】8页(P541-548)【作者】吴意明;张伟;刘保银;张中庆【作者单位】中海石油深圳分公司开发部,广东深圳518067;中海油服油田技术事业部,河北燕郊065201;杭州迅美科技有限公司,浙江杭州 310012;杭州迅美科技有限公司,浙江杭州 310012; 浙江大学,浙江杭州310012【正文语种】中文【中图分类】P631.84【相关文献】1.各向异性介质三维纵横波联合叠前反演方法及应用 [J], 李录明;罗省贤;王明春;孔选林;周怀来2.基于VTI各向异性介质的频率域海洋可控源电磁三维约束反演 [J], 赵宁;王绪本;秦策;余年;周军;徐玉聪3.基于各向异性介质模型讨论地电阻率各向异性变化 [J], 马占虎;杜学彬;陈军营;谭大诚4.任意各向异性介质三维非结构谱元法直流电阻率正演模拟研究 [J], 朱姣;殷长春;任秀艳;刘云鹤;惠哲剑;谷宇5.三维主轴各向异性介质中张量CSAMT正反演研究(英文) [J], 王涛;王堃鹏;谭捍东因版权原因,仅展示原文概要,查看原文内容请购买。
层状各向异性地层中三维感应测井响应快速计算及资料处理
层状各向异性地层中三维感应测井响应快速计算及资料处理肖加奇;张国艳;洪德成;杨善德【摘要】Using generalized reflection coefficients, we deduced a complete form of solution for the electromagnetic field in layered anisotropic formations. Those generalized reflection coefficients can be evaluated with recursion formula for each electromagnetic wave mode. A fast forward modeling is achieved through the Fast Hankel Transformation (FHT) method. The responses of a 3D induction logging tool are determined by the horizontal conductivity, the vertical conductivity, the relative dip of the formation and the azimuth angle of the tool, and a single component of the logging curve can't meet the need of data interpretation. Gaining insight into the response characteristics of the measurement components of a 3D induction tool, we developed a quick-look interpretation method by combining the measurement components in different ways. The processing results of numerically simulated data have shown the fact that the combination of the cross-components accurately outlines the boundary location and qualitatively recognize anisotropic layers, while the combination of the primary components enhances resolution and reduces nonlinearity relating to electrical parameters.%本文采用广义反射系数法推导了水平层状各向异性地层中电磁场的积分解析解,并利用快速汉克尔变换技术实现了三维感应仪器测井响应的快速计算.三维感应测井响应与地层水平电导率、垂直电导率和井斜角及仪器方位角同时有关,单一分量的测井曲线不能满足资料解释的需要.通过对仪器测量分量响应特征的考察,本文提出了一种基于组合量测井曲线的资料直观解释方法.数值模拟显示,交叉分量相关组合量可准确划分地层纵向边界,并可直观识别各向异性层;与单独分量相比,主分量相关组合量提高了纵向分辨率、减弱了与地层电导率参数的非线性关系.【期刊名称】《地球物理学报》【年(卷),期】2013(056)002【总页数】11页(P696-706)【关键词】三维感应测井;各向异性地层;水平层状介质;直观解释【作者】肖加奇;张国艳;洪德成;杨善德【作者单位】中国石油长城钻探工程有限公司,北京 100176;中国石油长城钻探工程有限公司,北京 100176;吉林大学物理学院,长春130012;吉林大学物理学院,长春130012【正文语种】中文【中图分类】P6311 引言据估计,世界上大约有30%的油气储存于砂-泥岩薄交互层中[1],这种薄交互储层可等效为宏观的单轴各向异性地层[2](或称横向各向同性地层,简记为TI地层),探测和识别这类地层对于油气资源的开发有重要意义和现实需要.现有的轴向型感应测井仪器由于纵向分辨率还不够高,往往将这种薄交互油储层误认为是高含水饱和度层而漏掉[3].三维感应测井仪器是由三个彼此垂直的发射线圈和与之平行的三个接收线圈组成,能探测到地层的水平电导率和垂直电导率信息,可以从三维角度识别地层特性,对薄储层、复杂储层的探测具有先天的优势.快速正演模拟是研究测井响应特征和资料处理的必要计算工具.含井眼和纵向邻层的非均质地层中三维感应测井响应的求解是三维问题,没有正演解析解,只能用有限元和有限差分等数值计算方法求解[4-7],其计算工作量大,计算时间长.忽略井眼环境的水平层状地层模型可用来研究纵向邻层对测井响应的影响及经过井眼效应校正的测井资料的直观解释与非线性迭代反演.国内外已经发表了大量关于水平层状各向异性地层中三轴感应测井响应的快速算法研究,如数值模式匹配方法、传输线理论和传播矩阵法等[8-14],这些算法本质上都是计算层状各向异性介质中的电磁场,并且计算效率相当.本文针对水平和垂直线圈源(等价为磁偶极子源)激发的TE波和TM波,利用广义反射系数递推方法推导了纵向行波的递推公式,给出了水平层状各向异性介质中电磁场积分解析解的另一种表达形式,并利用快速汉克尔变换技术[15]实现其快速计算.三维感应测井响应与地层水平电导率、垂直电导率和井斜角同时有关,且非线性严重,加上邻层对不同分量测井曲线的影响又大不一样,其资料的直观解释很难取得好的效果[16].此外,除共轴分量外其余分量都与仪器方位角有关也增加了资料解释的困难.目前,三维感应测井资料处理除多频聚焦方法(MFF)外,主要是多参数非线性迭代反演方法[9,17],而测井资料的直观解释却鲜有报道.迭代反演方法可以获得更接近真值的地层参数,但其计算速度慢,并不能满足现场对测井资料解释的需要,属于测井资料后处理.资料直观解释计算速度快,有利于测井资料的现场处理,并可为进一步的非线性迭代反演提供好的初始值.本文数值模拟了三维感应测井仪器9个分量随测量深度和仪器方位角变化的响应曲线,并对模拟数据给出了基于组合量测井曲线的直观处理,得到了一些有意义的结果.2 三维感应测井响应快速计算2.1 TI介质中TE波和TM波的分解在TI介质的主轴坐标系O-XYZ中电导率可表示为对角张量diag(σh,σh,σv),垂直电导率σv方向与Z轴平行,水平电导率σh与XY平面平行.均匀TI介质中电磁场法向分量的Sommerfeld积分表达式为[12]:其中Jn()为n阶Bessel函数,kρ 为积分变量,kv,z=和分别为纵向波数和径向波数;=iωμσh=iωμσv,ω为发射线圈圆频率为各向异性系数,M为磁偶极子强度,下角标代表不同的发射方向;和分别为波数域电磁场法向分量.式(1)中包含三种模式波:由水平磁偶极子源(HMD)激发的TM波、TE波和由垂直磁偶极子源(VMD)激发的TE波.将(1)式代入麦克斯韦方程[18],有得到波数域中电磁场水平分量:2.2 层状介质中电磁波利用TE波和TM波分解技术,我们只需推导波数域中电磁场法向分量在地层中的分布,然后由式(3)和(4)即可求解电磁场的水平分量.图1为N个水平层状各向异性地层模型,层界面位置用dn(n=1,2,…,N-1)表示;地层n的磁导率μn为各向同性;电导率用张量σ^n=diag(σn,h,σn,h,σn,v)表示.图1 水平层状各向异性地层模型Fig.1 Model of horizontal layered anisotropy formation令源所在的位置为坐标原点,所在层为第m层,则任意层n的电磁场法向分量可表示为纵向行波与柱面波的叠加:F为纵向行波的传播项,上角标TM,h表示由水平磁偶极子源激发的TM波,其余类推;对不同模式波,传播项F的表达式略有不同,但可统一写为:(详细推导见附录A)其中,上角标“±”分别代表上行波和下行波;A为振幅,为边界处的广义反射系数.对于传播项将式(6)中的指数因子做替换kz→kh,z,振幅具体展开为:对于传播项,将式(6)中的指数因子做替换kz→kh,z,振幅具体展开为:对于传播项,将式(6)中的指数因子做替换kz →λkv,z,振幅具体展开为:由式(5)和(6)可以看出,如果知道每一层的广义反射系数R和纵向行波振幅A,即可求得电磁场法向分量.三种模式波振幅An有统一的递推关系,指数因子按式(7)—(9)方式替换可得到相应模式波的具体表达式.为了表述简明,我们略去递推公式中的上标和下标.不失一般性,任意层n中的传播项F可写为:上行波在m+1层中的振幅可表示为边界z=dm处的反射波和透射波振幅的迭加:其中Rm+1,m和Tm,m+1为m层对m+1层的狭义反射和透射系数,整理后同理,下行波在m-1层中的振幅为:重复上面的推导过程,我们得到其余相邻层振幅的递推关系:各向异性地层的广义反射系数可由相邻两层的狭义反射、透射系数表示[19],利用电磁场水平分量在边界处的连续条件,由式(3)和(4)得:由此,上述三种模式波在边界处狭义反射、透射系数的具体表达式可写为:类似于水平层状各向同性地层模型中广义反射系数的推导方式[19],可推导出水平层状各向异性地层中广义反射系数的递推公式并整理成统一的表达形式:为满足自然边界条件有 =RN-1,N 和R2,1.由此,将式(10)式代入(3)—(4),可得到任意层波数域电磁场的水平分量表达式:由式(5)和式(22)可以看出,水平磁偶极子源激发的磁场水平分量同时与地层水平电导率和垂直电导率有关,法向分量只与水平电导率有关;垂直磁偶极子源激发的磁场法向分量和水平分量都只与水平电导率有关,与垂直电导率无关.这些基本关系将有助于三维感应测井响应特征考察与资料处理.空间域中的电磁场积分表达式为:式(23)的右端为Sommerfeld积分表达式,其积分核是慢衰减的高震荡函数.由于水平层状地层模型的波数域电磁场分量是解析递推,计算速度快,所以空间域中电磁场分量的计算速度主要取决于Sommerfeld积分计算效率.陈桂波等[20]详细比较了快速汉克尔变换技术、高阶窗函数结合连分式展开技术以及直接积分等方法求解Sommerfeld积分的计算效率,发现快速汉克尔变换技术的计算效率最高,是高阶窗函数结合连分式展开技术的3倍,直接积分的27倍.本文也采用这种求解方法.2.3 三维感应测井视电导率响应分量三维感应测井仪器的线圈系结构(图2所示)由三个中心共点的、彼此垂直的发射线圈Tx,Ty,Tz和与其平行的三个接收线圈Rx,Ry,Rz(源距为L1)三个屏蔽线圈Bx,By,Bz(源距为L2)组成.屏蔽线圈可抵消发射线圈在接收线圈中产生的直耦信号,其绕线方向与接收线圈的相反.图2 三线圈系结构示意图Fig.2 The diagram of three-coil system当发射线圈系向周围发射正弦交流电时,三维感应测井仪器同时测量接收线圈系上的九个磁场分量,它们构成完整的二阶张量:为了考察倾斜井中三维感应测井响应,需引入三个坐标系:地层坐标系O-XYZ (介质主轴坐标系)、仪器坐标系 O-X′Y′Z′ 和线圈坐标系OX″Y″Z″[21].三个坐标系下的磁场张量满足如下旋转变换规则:其中α是仪器轴与地层法向分量之间的夹角,即井眼斜角;φ是仪器绕自身旋转的方位角,即仪器方位角.经补偿后,接收线圈测量的磁场强度可表示为:为便于测井响应与地层电参数之间的比较,通常将感应测井响应规格化为电导率量纲的测量,记为其中Im)为磁场的虚部,Kij为线圈系仪器系数.3 数值算例及资料处理3.1 算法检验快速汉克尔变换技术计算效率高,但其计算精度容易受仪器参数及地层参数变化的影响[20].下面我们考察一般条件下感应测井快速汉克尔变换技术的计算精度.感应测井仪器接收点距离发射源点一般为0.3~2.0m,发射源频率为10~100kHz,地层电导率变化范围0.001~10.0S/m.图3给出了均匀各向异性介质中利用快速汉克尔变换技术计算的磁场分量虚部与其代数解析解[23]的相对误差.图3(a、c、e)仪器发射频率f=25kHz,地层水平电导率σh=1.0s/m、σh/σv=2;(b、d、f)仪器发射点为坐标原点,接收点为(1,0,1).可以看出绝大部分相对误差都在0.1%以下,满足计算精度的要求.水平层状各向异性地层模型没有代数解析解,图4给出了三层模型中本文推导的解析解(曲线)与有限元法[7](点线)模拟的三维感应测井响应结果的比较.目标层厚为3m,地层电导率依次为:σ1,h=σ1,v=0.5S/m,σ2,h=1S/m,σ2,v=0.1S/m,σ3,h=0.5S/m和σ3,v=0.125S/m;井眼倾角α=60°,发射频率为f=25kHz,主接收线圈源距L1=39in,补偿线圈源距L2=27in,方位角φ=0°.在Pentium 2.9GHz计算机上,200个测量点计算所需总计算时间为0.55s.从图4可以看出,两种方法的结果完全一致,证明了本文的计算方法是正确的、高效的.图3 快速汉克尔变换技术的计算精度Fig.3 The calculation precision of fast Hankel transform technology3.2 三维感应测井资料处理方法由式(26)知,一般情况下倾斜井中三维感应测井响应的9个分量都不为零,且除轴向分量外,其余分量都是仪器方位角的函数.图5给出了在同时含有各向同性层和各向异性层的地层模型中随仪器旋进的测井响应数值模拟曲线.井眼斜角α=60°,地层模型层厚22.5m,仪器实际测量的行进距离(测量深度MD)是45m;假定沿测量深度的测量间隔是0.1m;每个测量间隔仪器方位角逆时针旋进5°;仪器参数同上.由图5a可以看出,共面分量视电导率响应σx″x″和σy″y″ 受邻层影响比共轴分量σz″z″ 的要大,在边界附近有“犄角”特性,甚至出现负响应,这主要是由层边界的电荷积累引起的.σx″x″ 和σy″y″ 的响应行为很接近,但受各向异性电阻率或纵向邻层的影响,其测井曲线并不重合;随着仪器方位角的旋进,σx″x″和σy″y″在厚层内测井曲线会出现波动.共轴分量视电导率响应σz″z″曲线比较光滑,其响应值与地层水平电导率的走势有一致的对应关系.由式(26)可知,σz″z″响应值与仪器方位角无关.图5b是交叉分量视电导率响应σx″y″ 和σy″x″,它们相互重叠并在零点附近波动.图5c和5d是交叉分量视电导率响应σx″z″、σy″z″ 和σz″x″、σz″y″,它们的响应幅度波动很大,特别是在地层边界附近.从上面的分析可以看出视电导率响应分量是地层各向异性电导率和几何角度的函数,单一的测井响应很难反映地层真实电导率的纵向分布和划分地层纵向边界.在资料处理上,现在普遍采取首先将测得的9条响应曲线变换成仪器坐标系(仪器方位角为零)中5条不为零的测井曲线再做资料解释.由于仪器方位角的变化范围是0~180°,现有的求解方法[21,24]都不能由单一测量点确定方位角所属的象限和逐点转化仪器坐标系中的测井响应分量.本文利用线圈坐标系中直接测量的响应曲线构造出三条与仪器方位角无关的响应曲线σ-、σ+和σapp.组合曲线的构造方法如下:图6给出了由图5中测井响应曲线构造的组合曲线σ-、σ+和σapp.由图6a可见,组合曲线σ-和σ+的形态显著地规则化了,所受邻层的影响显著减少了,与地层参数之间的关系变得简单和直观了:组合曲线σ-的峰值附近响应曲线明显的有规律的弯折是由屏蔽线圈造成的,峰值与地层边界位置相对应.σ-的测井响应与地层垂直电导率无关,峰值的大小与相邻两层水平电导率的反差成正比[25].组合曲线σ+在各向异性层中的响应值明显偏离零,在各向同性层中响应值趋近零,利用这一响应特征可直观识别各向异性地层.图6b中组合曲线σapp比共轴分量σz″z″的纵向分辨率明显提高;各向同性地层中的响应值与其真值很接近,各向异性层中的响应值处于水平电导率和垂直电导率之间,但比σz″z″更接近地层水平电导率,这是因为组合量σapp与电导率的非线性关系比σz″z″的更弱.综合利用这三条构造组合曲线,我们给出三维感应测井初步的直观解释流程:首先利用组合曲线σ-分层和确定边界位置;其次利用组合曲线σ+区分各向异性层和各向同性层;最后利用组合曲线σapp的标志点给出各向同性地层的视电导率及各向异性层中水平电导率视值.由于各向异性垂直电导率往往比较小,探测灵敏度低,它的视值解释方法与水平电导率的也不同,还需进一步研究.图6 构造组合曲线的数值模拟Fig.6 The numerical simulation of the combination curves图7 各向异性Oklahoma地层模型中的组合曲线Fig.7 The combination curves in TI Oklahoma modelOklahoma地层模型是检验电法测井资料处理方法的经典各向同性模型,为进一步检验本文的三维感应测井资料处理方法,我们将其中电导率较高的地层推广为各向异性层.图7给出了包含各向异性地层的Oklahoma模型三条组合曲线的处理结果.组合曲线σapp的纵向分辨率比σz″z″有明显改善,响应值与地层水平电导率更接近;结合组合曲线σ-可以确定出地层纵向边界位置,包括小反差的高阻层和层厚小于1.0m的薄层.靠近围岩的各向异性厚层中组合曲线σ+明显偏离零,而中间的高阻各向同性层中σ+响应值则接近零,为各向异性厚层的识别提供了判断依据.4 结论本文利用广义反射系数法推导了水平层状各向异性介质中电磁场的积分解析解,并利用快速汉克尔变换技术实现了其快速计算.三种模式波的递推公式可以写成统一的表达形式,推导过程简明.通过与已发表的计算结果相比较,验证了本文公式推导的正确性.倾斜井中数值模拟结果显示了线圈坐标系下的实测曲线响应行为复杂,很难由原始的测井曲线直观提取地层参数信息.本文利用实测曲线构造了响应行为更加规则的三条组合曲线,简化了测井响应与地层参数之间的关系,实现了三维感应测井响应资料的初步直观解释.附录A 发射源所在层纵向行波的表达式令发射源所在层为第m层,由(1)式,任意层n的电磁场法向分量可表示为:其中和U分别代表下行波和上行波,且由自然边界条件有DN=0,U1=0.以垂直磁偶极子激发的横电波波为例,将式(A4)改写为振幅与广义反射系数的表达式.发射源所在层的上行波可看作是下行波在下边界的广义反射波,因此,在边界z=dm-1处有:同理,在边界z=dm处有:联立式(A4)和(A5)求解有:其中将式(A6)—(A8)代入式(A4),经整理得:其中的具体表达式见(7a)—(7c).重复上面的推导过程,式(A2)可改写为:其中的具体表达式见(9a)—(9c).式(A3)可改写为:其中的具体表达式见(8a)—(8c).式(A10)—(A12)除指数因子外,表达形式相同,略去相应的上角标,可将统一写成表达式(6).参考文献(References)[1]党瑞荣,秦瑶,谢雁等.三分量感应测井系统研究.石油地球物理勘探,2006,41(4):484-488.Dang R R,Qin Y,Xie Y,et al.Study of 3-C inductionlog system.Oil Geophysical Prospecting (in Chinese),2006,41(4):484-488.[2]Klein J D,Martin P R,Allen D F.The petrophysics of electrically anisotropic reservoirs.The Log Analyst,1997,38(3):25-36.[3]Zhang Z Y,Yu L,Kriegsmauser B,et al.Simultaneous determination of relative angles and anisotropic resistivity using multicomponent induction logging data.SPWLA 42nd Annual Logging Symposium,Houston,2001. [4]王昌学,周灿灿,储昭坦等.电性各向异性地层频率域电磁响应模拟.地球物理学报,2006,49(6):1873-1883.Wang C X,Zhou C C,Chu Z T,etal.Modeling of electromagnetic responses in frequency domain to electrical anisotropic formations.Chinese J.Geophys.(in Chinese),2006,49(6):1873-1883.[5]Davydycheva S,Druskin V,Habashy T.An efficient finitedifference scheme for electromagnetic logging in 3D anisotropic in homogeneous media.Geophysics,2003,68(5):1525-2536.[6]孙向阳,聂在平,赵延文等.用矢量有限元方法模拟随钻测井仪在倾斜各向异性地层中的电磁响应.地球物理学报,2008,51(5):1600-1607.Sun X Y,Nie ZP,Zhao Y W,et al.The electromagnetic modeling of logging-while-drilling tool in tilted anisotropic formations using vector finite element method.Chinese J.Geophys.(in Chinese),2008,51(5):1600-1607.[7]Lv W G,Chu Z T,Zhao X Q,et al.Simulation of electromagnetic wave logging response in deviated wells based on vector finite element method.Chinese Physics Letters,2009,26(1):96-99.[8]汪宏年,陶宏根,姚敬金等.用模式匹配算法研究层状各向异性倾斜地层中多分量感应测井响应.地球物理学报,2008,51(5):1591-1599.Wang H N,Tao H G,Yao J J,et al.Study on the response of a multicomponent induction logging tool in deviated and layered anisotropic formations by using numerical mode matching method.Chinese J.Geophys.(in Chinese),2008,51(5):1591-1599.[9]Wang H N,Tao H G,Yao J J,et al.Fast multiparameter reconstruction of multicomponent induction well-logging datum in a deviated well in a horizontally stratified anisotropic formation.IEEE Transactions on Geoscience and Remote Sensing,2008,46(5):1525-1534.[10]Wei B J,Zhang G J,Liu Q H.Recursive algorithm and accurate computation of dyadic Green′s functions for stratified uniaxial anisotropic media.Sci.China Inf.Sci.,2008,51(1):63-80.[11]Zhong L L,Li J,Bhardwaj A,et putation of triaxial induction logging tools in layered anisotropic dipping formations.IEEETrans.Geosic.Remote.Sens.,2008,46(4),1148-1163.[12]范宜仁,杨震,邓少贵等.层状介质中大斜度井感应测井响应计算新方法.西南石油大学学报(自然科学版),2009,31(1):166-169.Fan Y R,Yang Z,Deng S G,et al.New computation method of induction logging responseof highly-deviated wells in stratified media.Journal of Southwest Petroleum University(Science &Technology Edition)(in Chinese),2009,31(1):166-169.[13]魏宝君,王甜甜,王颖.用磁流源并矢Green函数的递推矩阵方法计算层状各向异性地层中多分量感应测井响应.地球物理学报,2009,52(11):2920-2928.Wei B J,Wang T T,Wang puting the response of multi-component induction logging in layered anisotropic formation by the recursive matrix method for magneticcurrent-sourc e dyadic Green′s function.Chinese J.Geophys.(in Chinese),2009,52(11):2920-2928. [14]姚东华,汪宏年,杨守文等.用传播矩阵法研究层状正交各向异性地层中多分量感应测井响应.地球物理学报,2010,53(12):3026-3037.Yao D H,Wang H N,Yang S W,et al.Study on the responses of multi-component induction logging tool in layered orthorhombic anisotropy formations by using propagator matrix method.Chinese J.Geophys.(in Chinese),2010,53(12):3026-3037.[15]Andersion W puter program:numerical integration of related Hankel transforms of orders 0and 1by adaptive digital filtering.Geophysics,1979,44(7):1287-1305.[16]Wang H,Barber T, Morriss C,et al.Determining anisotropicformationresistivityatanyrelativedipusinga multiarray triaxial induction tool.∥SPE 103113presented at the 2006SPE Annual Technical Conference and Exhibition:San Antonio,Texas,U.S.A.,2006.[17]Barber T,Andersion B,Abubakar A,et al.Determining formationresistivity anisotropy in the presence of invasion.∥SPE 90256presented at the 2004SPE Annual Technical Conference and Exhibition:Houston,Texas,U.S.A.,2004.[18]Moran J H,Gianzero S.Effects of formation anisotropy on resistivity-logging measurements.Geophysics,1979,44(7):1266-1286.[19]Chen W C.Waves and Fields in Inhomogeneous Media.New York:Van Nostrand Reinhold,1990:48-65.[20]陈桂波,汪宏年,姚敬金等.水平层状各向异性介质中电磁场并矢Green函数的一种高效算法.物理学报,2009,58(3):1608-1618.Chen G B,Wang H N,Yao J J,et al.An efficient algorithm of the electromagnetic dyadic Green′s function in a horizontallayered anisotropic medium.Acta Physica Sinica(in Chinese),2009,58(3):1608-1618.[21]Zhang Z Y,Yu L M,Kriegsmauser B,et al.Simultaneous determination of relative angles and anisotropic resistivity using multicomponent induction logging data.∥SPWLA 42nd Annual Logging Symposium,Houston,2001,Paper Q.[22]Mallan R K, Torres-Verdin C.Effects of geometrical,environmental,and petrophysical parameters on multicomponent induction measurements acquired in high-angle wells.Petrophysics,2007,48(4):271-288.[23]Kennedy D,Peksen E,Zhdanov M.Foundations of tensor induction well-logging.Petrophysics,2001,44(6):588-610.[24]Wang H M,Barber T,Rosthal R,et al.Fast and rigorous inversion of triaxial induction logging data to determine formation resistivityanisotropy,bed boundary position,relative dip and azimuth angles.∥SEG International Exposition and Seventy-Third Annual Meeting,2003,22. [25]Minerbo G N,Omeragic D,Rosthal R A.Directional electromagnetic measurements insensitive to dip and anisotropy.Patent US6969994B2,2005.。
电法测井术语
电法测井3.1 电阻率测井 resistivity logging测量地层电阻率的测井方法。
3.2 岩石电阻率 resistivity of rock岩石的电学参数之一,是岩石阻抗电流通过其自身的特性。
3.3 岩石电导率 conductivity of rock岩石的电学参数之一,主要反映电流通过其自身的能力。
与岩石电阻率互为倒数关系。
3.4 钻井液电阻宰 mud resistivity钻井过程中使用的钻井液的电阻率。
3.5 泥饼电阻率 mud-cake resistivity由于井内压力与地层压力的压差,使钻井液中的固体颗粒附着在井壁上形成的泥饼的电阻率。
3.6 钻井液滤液电阻率 mud filtrate resistivity渗入渗透性地层的钻井液滤液电阻率。
3.7 冲洗带电阻率 flushed zone resistivity冲洗带地层的电阻率。
3.8 过渡带电阻率 transitional zone resistivity过渡带地层的电阻率。
3.9 侵入带电阻率 invaded zone resistivity侵入带地层的电阻率。
3.10 地层真电阻率 true formation resistivity原状地层的电阻率。
3.11 地层水电阻率 formation water resistivity原状地层孔隙中所含水的电阻率。
3.12 围岩电阻率 shoulder bed resistivity盖层及垫层的电阻率。
3.13 地层视电阻率 apparent formation resistivity受井筒、侵入带和围岩等测井环境的影响,地层真电阻率的测量结果。
3.14 水平电阻率 horizontal resistivity;parallel resistivity水平方向的电阻率。
3.15 垂直电阻率 perpendicular resistivity;vertical resistivity垂直方向的电阻率。
多分量感应测井理论和方法研究的开题报告
多分量感应测井理论和方法研究的开题报告
题目:多分量感应测井理论和方法研究
一、研究背景及意义
随着石油勘探开采工作的不断深入,传统的测井技术已经无法满足现代油气勘探的需求。
与传统测井技术相比,多分量感应测井技术具有更高的精度和分类能力,且对油气藏有更好的敏感性,因此在现代油气勘探中已经得到了广泛应用。
在未来的石油勘探中,多分量感应测井技术将继续发挥其重要的作用。
二、研究目标和内容
本文研究的目标是探讨多分量感应测井技术的理论基础和方法,以及其在油气勘探中的应用。
具体内容包括:
1.多分量感应测井技术的基本原理和方法;
2.多分量感应测井数据处理的方法和技术;
3.多分量感应测井技术在油气勘探中的应用;
4.多分量感应测井技术的优缺点及未来的发展方向。
三、研究方法
本文将采用文献综述法,对多分量感应测井技术的理论和应用进行梳理和分析。
首先对多分量感应测井技术的基本原理进行详细的介绍,包括多分量测量的方式、计算和数据处理方法等方面。
接着以实际应用中的数据为基础,对多分量感应测井技术在油气勘探中的应用进行深入研究和探讨。
最后,通过评价多分量感应测井技术的优缺点,分析其未来的发展方向。
四、预期成果
预期成果主要包括:
1.详细介绍多分量感应测井技术的基本原理和方法;
2.对多分量感应测井技术在油气勘探中的应用进行分析和总结;
3.评价多分量感应测井技术的优缺点,探讨其未来的发展方向。
三分量感应测井系统的理论研究
・开发设计・三分量感应测井系统的理论研究3闫敏杰 党瑞荣 袁阿明 谢 雁(光电油气测井与检测教育部重点实验室(西安石油大学) 陕西西安)摘 要:常规感应测井仪器只有Z方向的线圈,只能测量水平方向的平均电阻率。
三分量感应测井系统分别在X、Y、Z三个方向布置三组相互垂直的发射-接收线圈对,可直接测量地层的水平电阻率和垂直电阻率,通过测量两个交叉分量还可得出地层的倾角和方位,从而可以为储层评价提供更准确的信息。
文章设计了三分量感应测井系统的线圈系布局结构和排列方法以及线圈系的聚焦方式,主线圈系采用了发射聚焦方式,交叉线圈采用了接收聚焦方式。
通过线圈系匝数和距离的选取以及几何因子的计算,各向异性地层中的正、反演的研究及求解,提高了线圈系的水平探测深度和纵向分辨率,使三分量感应测井系统的探测性能达到最佳。
关键词:三分量;聚焦方式;几何因子;分辨率;探测性能中图法分类号:P631.8+11 文献标识码:B 文章编号:100429134(2009)01200092040 引 言感应测井是一种重要的测井方法,但长期以来,感应测井的基本理论是建立在均匀无限大地层假设基础上的[1],线圈系的轴线与井轴是平行或重合的,只能给出水平方向的一维电阻率,很容易低估储量,漏测储层[2]。
K lein验证并给出了常规测井仪器低估和漏测储层的实际事例[3~5]。
产生这种现象的原因是传统感应测井仪不能描述砂泥岩构成的薄交互储层所体现出的电各向异性的特性,仪器对这种薄交互储层给出的测量结果体现的是低电阻率特性,漏掉了高电阻率的油气储层。
三分量感应测井系统采用了轴线为X、Y、Z三个方向的线圈系,可给出水平电阻率和垂直电阻率,二者的比值可准确地描述地层的电各向异性的特性,对砂泥岩构成的薄交互储层,水平电阻率和垂直电阻率将有明显的差异,其中垂直电阻率体现的是高电阻率的油气层。
这就为储层的准确评价提供了重要的依据。
我国许多油田为低产、低渗、低阻油田,地质结构复杂,各向异性特性明显,而且随着钻井技术的进步及高效开采油气资源的背景需求,斜井、水平井的数目日益增多,对测井技术也提出了越来越高的要求,传统的感应测井技术遇到了新的挑战。
2001 感应测井几何因子理论及其应用研究
第 25 卷・第 6 期
仵 杰等: 感应测井几何因子理论及其应用研究
・419・
有趋肤效应校正, 深感应增加了三点或多点反仪器的响应, 随着几何因子理论和信号处 理技术的发展, 提出了利用几何因子进行测量信号处理 的多种方法。 传统双感应测井主要是围岩和趋肤效应校 正, 阵列感应测井是二维的环境影响和趋肤效应校正。
1 ( 1 - ik rT ) e ik rT + g 2 Do ll ( 2)
当它用于计算纵向和径向微分几何因子时, 除了对全空 间几何尺寸积分外还要对电导率积分。 田子立等人[ 3 ] 证 明在均匀地层中,M o ran 几何因子对全空间和电导率积 分等于均匀地层的精确解。
几何因子的应用
在实际应用中, 为了更好地描述感应测井测量信号 来源以及感应测井仪器的纵向分辨能力和径向探测特 性, 进一步由二维几何因子定义了一维几何因子, 它们分 别是纵向微分几何因子 g VD ( z ) =
) - Ρb ( z ) ]d z ′ [ Ρ( z ′
式中, Ρ- Ρb 表示背景电导率 Ρb 上的扰动; Ρba ( z ) 是地层 电导率为背景电导率 Ρb 时的视背景电导率; g Bo rn ( z , Ρb ) 是背景电导率 Ρb 的纵向Bo rn 几何因子。对于双线圈系, 二维 Bo rn 几何因子的表达式为
第 25 卷・第 6 期
仵 杰等: 感应测井几何因子理论及其应用研究
・417・
感应测井几何因子理论及其应用研究
仵 杰 庞巨丰
( 西安石油学院)
( 海军航空工程学院)
徐景硕
摘
要
仵杰, 庞巨丰, 徐景硕 . 感应测井几何因子理论及其应用研究 . 测井技术, 2001, 25 ( 6) : 417 ~ 422 首先讨论 Do ll、 M o ran、 G ianzero 和 Bo rn 等 4 种常用几何因子建立的物理条件及其应用范围; 然后研究这些几何因子在 实际应用中存在的问题。分析表明, Do ll 几何因子仅有理论意义, 一般情况下需用 G ianzero 或M o ran 几何因子分析仪器探测 特性和测量解释; G ianzero 和M o ran 几何因子的实部可用于描述实际测量信号的实部, 但虚部并未对应于信号的虚部; Bo rn 几何因子的实部和虚部近似地反映了非均匀地层对均匀背景地层的扰动, 所以, 背景电导率的选择是其应用的关键。 主题词: 几何因子理论 感应测井 局限性分析 应用
基于耦合标势与矢势的有限体积法模拟非均匀各向异性地层中多分量感应测井三维响应
基于耦合标势与矢势的有限体积法模拟非均匀各向异性地层中多分量感应测井三维响应张烨;汪宏年;陶宏根;杨守文【期刊名称】《地球物理学报》【年(卷),期】2012(055)006【摘要】为快速有效地研究、考察各向异性地层条件下多分量感应测井的响应特征,本文利用电场标势与矢势的有限体积法研制出三维频率域电磁场响应的数值模拟算法,克服由低频发射或高阻地层产生的低感应数问题,有效提高了三维电磁数值模拟算法的应用范围和计算效率.首先利用电场的标势与矢势将Maxwell方程转化为满足库仑规范条件的耦合势Helmholtz方程,以Yee氏交错非均匀网格中不同位置上的节点为中心建立四种控制体积单元,通过对控制体积单元中电磁场与电导率的积分平均实现耦合势方程和磁偶极子旋度的离散,并得到一个对角占优的大型稀疏复线性代数方程组,然后,通过不完全LU分解预处理和稳定双共轭梯度法快速求解离散方程.数值结果证明了该算法的有效性,并进一步考查了仪器偏心、倾斜井、垂直裂缝等复杂条件下多分量感应的响应特征.%In order to efficiently investigate the 3D responses of the multi-component induction logging (MCIL) tool in complex anisotropic formations, we apply a novel finite volume forward method to establish the 3D numerical simulation of the EM fields based on the coupled scalar-vector potentials. As the results, we overcome the low induction number problems (LINs) in both the low frequency domain and high resistivity formation, extend the application range of the 3D EM modeling and enhance its efficiency. We firstreformulate Maxwell s equation into Helmholtz equations in terms of coupled scalar-vector potentials with Coulomb gauge. Accordingrnto the different location of node in the Yee's non-uniform staggered grids, we set up four kinds of control volume cells. With volume integral averaging of both the EM fields and conductivity tensor on the control volume, we discretize the Helmholtz equations and the rotations of magnetic dipole sources successfully. After that, we obtain a large, sparse and complex linear system with a block diagonally dominant structure. A combination of bi-conjugate gradient stabilization (BICGSATB) with an incomplete LU-decomposition preconditioner is efficiently applied to iteratively solve the system. Numerical results validate the algorithm and we further investigate the characteristics of MCIL responses in the different cases such as tool eccentricity, dipping well and vertical fracture.【总页数】12页(P2141-2152)【作者】张烨;汪宏年;陶宏根;杨守文【作者单位】吉林大学物理学院,长春130012;吉林大学物理学院,长春130012;大庆钻探工程公司测井公司,大庆163412;吉林大学物理学院,长春130012【正文语种】中文【中图分类】P631【相关文献】1.层状非均质各向异性地层中多分量感应测井数值模拟与响应特征 [J], 汪建勋2.耦合势有限体积法模拟海洋移动导线源三维频率域电磁响应 [J], 周建美;李貅;戚志鹏3.非均质各向异性地层中方位随钻电磁测井响应三维有限体积法数值模拟算法∗[J], 王浩森;杨守文;白彦;陈涛;汪宏年4.用传播矩阵法研究层状正交各向异性地层中多分量感应测井响应 [J], 姚东华;汪宏年;杨守文;杨海亮5.倾斜各向异性地层中多分量电磁波测井响应三维时域有限差分(FDTD)算法 [J], 杨守文;汪宏年;陈桂波;姚东华因版权原因,仅展示原文概要,查看原文内容请购买。
关于感应测井纵向几何因子
关于感应测井纵向几何因子的报告,600字
感应测井纵向几何因子报告
感应测井是一种测量当前地下油水藏的动态和静态特性的一种方法。
它可以帮助地有关人员对油水藏体包裹进行识别和了解地质结构。
它是通过感应测井仪器在油水藏中产生电磁感应场,然后测量随时间变化的反射信号来实现的。
其中,最重要的参数就是纵向几何因子(longitudinal geometric factor)。
纵向几何因子是测井纵向方向上电磁功率的大小,它主要反映了油水藏的类型,其中有不同的孔隙和渗透率的类型。
它的大小受到油水藏的体积、井眼的位置和测井仪器倾斜角的影响。
它的变化会引起地下油水藏的反射信号变化,从而给出了更多有用的信息。
此外,在解决地质矿产勘探问题时,使用纵向几何因子还可以判断油水藏的密度和厚度。
它可以根据观察到的“折回”角度反映孔隙的密度,从而可以分析出油水藏的类型。
此外,它还可以用来判断油水藏的厚度,这在分析油水藏的地质结构方面非常重要。
由于纵向几何因子的重要性,在测量油水藏参数时应该特别注意纵向几何因子的变化,以及如何精确测量和计算它。
需要注意的是,纵向几何因子受多项因素的影响,包括油水藏体积、井眼位置和测井仪器倾斜角等,所以在使用它来测量和计算参数之前,必须先了解这些因素的影响。
总之,纵向几何因子是测井仪器用来测量反射信号的一项重要参数。
它可以用来充分了解油水藏的体积、井眼位置和测井仪器倾斜角等,从而更准确地测量和计算油水藏参数。
但是,受多种影响,在使用纵向几何因子前要进行充分的了解和研究,以准确测量反射信号。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
2.2
各向异性介质 Born 几何因子
本文中我们采用 Born 近似方法, 将二维 Doll
各向同性地层 Born 几何因子推广至三维空间各向 异性介质, 给出直角坐标系中各向异性介质多分 量感应测井三维 Born 几何因子. 根据 Born 近似方 法, 各向异性异常体存在时, 背景介质中任意点 r 处的磁场是源在背景介质中激发的一次场 Hb (r) 与各向异性异常体产生的散射场的叠加: H (r) ≈Hb (r) ∫ ˜ H (r, r′ ) · Eb (r′ ) d r′ δ σ ˆ (r ′ )G +
各向异性介质多分量感应测井三维 Born 几何因子 理论研究∗
王磊 1)2)† 范宜仁 1)2) 黄瑞 3) 韩玉娇 1)2) 巫振观 1)2) 邢东辉 1)2) 李炜 4)
1) (中国石油大学地球科学与技术学院, 青岛 2) (中国石油大学 CNPC 测井重点实验室, 青岛 3) (中国石油天然气勘探开发公司, 北京 266580) 266580)
100034) 300451)
4) (中海油田服务股份有限公司物探事业部, 天津
( 2015 年 5 月 23 日收到; 2015 年 8 月 3 日收到修改稿 )
几何因子理论被广泛运用于感应测井仪器设计、 响应特征分析及高分辨率曲线处理, 现有的感应测井几 何因子主要适用于二维各向同性介质, 难以满足各向异性介质多分量感应测井的研究需要. 本文基于 Born 近似方法将 Born 几何因子扩展至各向异性介质, 推导了各向异性介质多分量感应测井三维 Born 几何因子表 达式, 随后进一步考察了多分量感应测井对地层各向异性的敏感性和探测能力. 数值模拟结果表明, 各向异 性系数越大, 多分量感应测井几何因子空间分布越复杂, 其对各向异性的敏感性越高. 多分量感应测井仪器 在不同倾斜角度对各向异性均具有较高的敏感性. 直井条件下, xx 共面分量受各向异性影响严重, 井斜角在 40◦ —60◦ 时, xz /zx 交叉分量对各向异性的敏感性高, 水平井中 zz 同轴分量则对各向异性的探测能力最强. 各向异性介质三维 Born 几何因子弥补了现有几何因子理论的不足, 可为国内新型多分量感应测井仪器研发 及储层各向异性评价方法的发展提供有利条件.
0
x/in
100
-100 -100 0 y/ i n 100 -100
0
x/in
100
-100 -100 0 y/ i n 100 -100
0
x/in
100
(g) zx 100 z /in z /in 0 100 0
(h) zy 100 z /in 0
(i) zz
-100 -100 0 y/ i n 100 -100
斜井条件下, 假定仪器与地层 z 轴方向夹角为 θ, 当不考虑方位角变化时, 地层坐标系下单位磁偶 极子源产生的磁场张量 H 经两次坐标转换即可确 ˆ为 定仪器坐标下单位磁偶极子源产生的磁场 H cos θ 0 sin θ cos θ 0 − sin θ ˆ = H H 0 1 1 0 . (3) 0 0 − sin θ 0 cos θ sin θ 0 cos θ
网络出版时间:2015-11-02 15:27:23 网络出版地址:/kcms/detail/11.1958.O4.20151102.1527.036.html
物 理 学 报 Acta Phys. Sin.
Vol. 64, No. 23 (2015) 239301
∗ 国家自然科学基金 (批准号: 41474100)、 中央高校基本科研业务费 (批准号: 15CX06008A, 14CX02114A) 资助的课题. † 通信作者. E-mail: upcwanglei199133@
© 2015 中国物理学会 Chinese Physical Society 239301-1
因而可提供任意井斜角度下水平电阻率和垂直电 阻率及地层地质结构等信息 [25−28] . 目前对多分量感应测井的研究主要围绕仪器 宏观响应规律展开 [29−34] , 但关于不同空间分布对 多分量感应测井测量信号的相对贡献和敏感性等 方面研究极少 [35−38] . 感应测井中一般采用 Doll 几 何因子研究不同空间分布对测量信号的相对贡献, 但其未考虑趋肤效应的影响 [39] . Moran, Gianzero 进一步考虑了地层导电环涡流间的相互影响, 分别 提出了感应测井二维 Moran 几何因子和 Gianzero 几何因子 [40,41] . Zhang 等将感应测井响应作为参 数的泛函展开成 Taylor 级数, 给出了感应测井高阶 几何因子 [42,43] . Howard 等则将量子散射问题中的 Born 近似方法用于非均匀地层感应测井问题, 提
关键词: 多分量感应测井, 几何因子, 各向异性介质, Born 近似 PACS: 93.90.+y, 91.60.Pn, 93.85.–q DOI: 10.7498/aps.64.239301
1 引
言
电性各向异性广泛存在于砂泥岩薄互层、 页岩 储层等非常规油气藏 [1−4] , 其形成原因主要包括: 1) 方向性裂缝发育或小于测井仪器分辨率的薄互 层 [5,6] ; 2) 不同的沉积过程或者岩石颗粒的取向方 式不一致 [7] . 受各向异性的影响, 电测井仪器响应 变得更为复杂, 导致储层含油饱和度计算和测井 评价的困难 [8−16] . 针对各向异性地层探测和评价 的需要, 国外石油公司陆续推出了多分量感应测 井仪器, 如 Schlumberger 的 Rt-Scanner, Baker Atlas 的 3DEX 等 [17−20] . 与传统感应类测井仪器相 比 [21−24] , 多分量感应测井仪器测量 9 个磁场分量,
Vs
(4)
式中, Eb (r′ ) 为发射线圈在异常体内任意点 r′ 处产 生的电场. δ σ ˆ 为异常体电导率相对于背景介质的 ˜ H (r, r′ ) 代表 r′ 处的电流源在 变化, 张量格林函数 G r 处产生的磁场, (4) 式具体推导步骤见附录. (5) 式 为各向异性地层多分量感应测井三维 Born 几何因 子S 为
x 其中, Hy 的上标 x 和下标 y 分别代表发射线圈和接 收线圈的法线方向. 感应测井发射和接收线圈可等效为磁偶极子 源, 当忽略泥浆侵入等影响时, 地层模型简化为各 向异性介质, 在地层坐标系下 (x, y , z ), 电导率可 用张量形式 σ ˆ 表示, 即 σh 0 0 σ ˆ = 0 σh 0 , (2) 0 0 σv
↼xT֒ yT֒ zT↽
x x x Sx Sy Sz
(6)
y y y S = Sx Sy Sz , z z z Sx Sy Sz (5)
σbv r σbhห้องสมุดไป่ตู้z ᑀఀ̮᠏ x ↼xR֒ yR֒ zR↽
σv σh Vs पʹ
y
x 其中, Sy 上标 x 为发射线圈法线方向, 下标 y 为接 收线圈法线方向. 采用如 (3) 式所示的坐标变换, 可得到斜井条件下的多分量感应测井几何因子张 ˆ为 量S Sxx Sxy Sxz ˆ = S Syx Syy Syz Szx Szy Szz cos θ 0 − sin θ cos θ 0 sin θ = 0 1 0 S 0 1 0 , (7) sin θ 0 cos θ − sin θ 0 cos θ
式中, 水平方向和垂直方向电导率分别为 σh 和 σv , √ 各向异性系数 λ = σh /σv .
zϕ yϕ
xϕ
L
z σv θ σh y σh x
图1
多分量感应测井仪器示意图
Fig. 1. The configuration of multi-component induction logging tool.
图2 function.
Born 响应函数地层模型
Fig. 2. The configuration for modeling Born response
式中, S xz 表示仪器坐标系下 x 方向发射线圈 z 方 向接收线圈组合的 Born 几何因子.
由上式可知, Born 几何因子 S x (S y 和 S z 类似) 主要由三部分组成, 前两项为水平电导率的相对贡 献, 第三项则代表了垂直电导率的相对贡献. (6) 式 给出了地层坐标系中不同发射和接收线圈组合的 Born 几何因子张量 S 为
物 理 学 报 Acta Phys. Sin.
Vol. 64, No. 23 (2015) 239301
出了二维 Born 几何因子 [44] . 上述几何因子被广泛 运用于感应测井仪器设计、 响应特征分析及高分辨 率曲线处理, 但其只适用于同轴发射与同轴接收线 圈仪器. 为弥补二维几何因子的不足, Alumbaugh 等将几何因子推广至三维空间, 提出了适用于各个 分量的三维 Born 几何因子 [45−47] . 现有的几何因 子均以各向同性介质为基础, 尚不能满足斜井各向 异性地层仪器的探测特性分析与研究的需要. 因 此, 开展各向异性地层多分量感应测井几何因子理 论研究, 对进一步促进仪器设计与研发、 完善多分 量感应测井在非常规油气藏勘探开发中的应用具 有重要意义. 本文简述了多分量感应测井的基本原理, 分析 了各向同性介质多分量感应测井几何因子空间分 布及敏感性, 在各向同性介质 Born 几何因子的基 础上, 将其推广至各向异性介质, 给出了各向异性 介质多分量感应测井三维 Born 几何因子. 在此基 础上, 对比分析了不同井斜角度及各向异性系数条 件下, 空间分布对多分量感应测井的相对贡献和敏 感性. 本文的相关成果为国内多分量感应测井仪器 研发及测井解释评价提供了理论基础和技术支持.