地球物理学报  2016, Vol. 59 Issue (3): 1131-1142   PDF    
新型贴井壁式阵列方位侧向测井数值模拟与响应特性
高建申1, 孙建孟1, 于其蛟2, 崔利凯1, 姜艳娇1    
1. 中国石油大学(华东)地球科学与技术学院, 青岛 266580;
2. 中石化胜利石油工程有限公司测井公司, 山东 东营 257061
摘要: 在现有侧向测井仪器基础上,提出了一种新型贴井壁式阵列方位侧向测井电极系,该仪器能提供径向和周向方位电阻率测量.该电极系采用贴井壁测量方式,提供5种径向探测深度.利用三维有限元方法模拟了阵列方位侧向测井电极系的井眼影响特性、径向探测特性、纵向分层能力、方位分辨能力,并模拟其对井周地层、水平井非对称泥浆侵入和倾斜地层的响应.在导电泥浆中最大探测深度为1.23 m,纵向分辨率为0.3 m,可以识别出0.1 m薄层,方位分辨率为20°.贴井壁测量时,纵向分辨率不受泥浆和围岩电阻率的影响,能够准确测量井周方位电阻率,较不贴井壁测量具有很大优势,同时利用12条方位电阻率曲线能够反映出水平井泥浆非对称侵入特性,倾斜地层倾角和倾斜方向.
关键词: 阵列方位侧向测井     贴井壁     方位电阻率     有限元方法     倾斜地层    
Numerical simulation of the new sidewall array azimuthal laterolog and response characteristics
GAO Jian-Shen1, SUN Jian-Meng1, YU Qi-Jiao2, CUI Li-Kai1, JIANG Yan-jiao1    
1. School of Geosciences, China University of Petroleum, Qingdao 266580, China;
2. Well Logging Company of Shengli Petroleum Engineering Co., LTD, SINOPEC, Shandong Dongying 257061, China
Abstract: Based on existing lateral logging instruments,a new sidewall array azimuthal laterolog sonde is proposed,which allows measurement of radial and azimuthal resistivity with five radial investigation depths. Using the three-dimensional finite element method,we simulate the effects of borehole, radial investigation characteristics,longitudinal layering ability and azimuthal resolution. Moreover, the responses of the sonde in resistivity measurement around the borehole, in the horizontal well with asymmetric mud invasion and in dipping strata are also simulated. The results show that the biggest investigation depth of this device can reach 1.23 m for conductive mud. This sonde can also recognize 0.1 m thin bed without mud invasion, with 0.3 m longitudinal resolution and 20åzimuthal resolution. The sidewall way can accurately measure the resistivity around the borehole, basically eliminates the effect of the borehole and surrounding rock on the longitudinal resolution and has more advantages than the non-contact device.The twelve azimuthal resistivity curves can also reflect the characteristics of asymmetric mud invasion in the horizontal well and dipping strata.
Key words: Array azimuthal laterolog     Sidewall     Azimuthal resistivity     Finite element method     Dipping stratum    
1 引言

侧向测井在盐水泥浆和高阻地层条件下发挥重要作用.双侧向测井(Suau et al.,1972)提供深、浅电阻率测量,探测深度达到1.5 m以上,但纵向分辨率较低,为0.6~0.8 m(Chen et al.,1998).随着油田勘探开发的不断深入,对薄层、薄互层和复杂地层的识别与评价至关重要(Forsyth et al.,1993).为提高纵向分辨率,朱军等(2007)胡海涛等(2014)童茂松和宋建华(2014)对传统双侧向进行了改进,纵向分辨率提高到0.2~0.4 m.为实现井 周不同方位电阻率测量,ARI(Azimuthal Resistivity Imager)(Davies et al.,1994)在双侧向A2屏蔽电极上安装方位电极阵列,保持双侧向测井基础上,提供井周方位探测,纵向分辨率为0.2 m,并有文献报道ARI的现场应用(Al-Waheed and Audah,1994).但由于不 贴井壁,ARI严重受到仪器偏心和较大尺寸井眼影响(黄继贞等,2000),测量得到的视电阻率与地层真电阻率相差很大(杨韡,2002). HALS(High Resolution Azimuthal Laterolog Sonde)(Smits et al.,1995)将方位电极阵列放置在电极系中心位置,采用软件聚焦方式代替硬件聚焦,缩短电极系长度,提高了纵向分辨率,但也降低了径向探测深度,而且方位电极没有实现方位聚焦效果,未见大规模应用(黄继贞等,2000).在ARI和HALS基础上,黄继贞等(2001)提出了一种方位侧向测井仪,采用了贴井壁的测量方式,提高了方位探测能力,径向探测深度较小,主要用于井周不同方位非均质地层电阻率测量.

为得到更加丰富的泥浆侵入信息和消除Groningen效应(Woodhouse,1978)影响,HRLA(High-Resolution Laterolog Array tool)(Smits et al.,1998),HDLL(High-Definition Lateral Log)(Itskovich et al.,1998),Multi Laterolog(Zhou et al.,2008Li et al.,2013),HAL6505(High Resolution Array Laterolog Logging Tool)(贺飞等,2013)采用阵列化设计,增加电极数量,电流返回电极系自身,实现了多种不同径向深度探测模式,有效消除了Groningen效应,纵向分辨率提高到0.3 m,多篇文献报道了阵列侧向的应用情况(Hakvoort et al.,1998Griffiths et al.,1999Brambilla et al.,2013). 为了同时实现多种探测深度和方位测量,3D-LS(Yin et al.,2014)将方位电极阵列和屏蔽电极阵列结合,采用不贴井壁测量方式,提供径向阵列探测和方位探测,其纵向分辨率较低,测量方位地层真电阻率能力较差.

另外,多名学者利用有限元数值模拟方法(李大潜等,1980Nam et al.,2009潘克家等,2013)对现有的双侧向、方位侧向和阵列侧向测井仪器进行了数值模拟研究,分析了多种地层环境对侧向测井响应的影响(Faivre,1993史謌,2004朱军和冯琳伟,2007邓少贵等,20062010a2010b2012邓少贵和李智强,2009; Nam et al.,2010谭茂金等,2012).

因此,集多种径向探测深度、较高纵向分辨率和良好的方位探测特性于一体的侧向测井仪器有待研究发展.针对这一现状,本文提出了一种贴井壁式阵列方位侧向测井电极系AALS(Array Azimuthal Laterolog Sonde).利用有限元数值模拟方法分析了电极系的探测特性,尤其是方位探测在井周电阻率测量,水平井泥浆侵入和倾斜地层中的应用,为侧向测井仪器发展提供支持.

2 电极结构及测量原理

阵列方位侧向测井电极系结构如图 1所示.图 1a为电极系整体结构示意图,电极系中心位置为主 电极A0,表面均匀分布着12个方位电极A0i(i=1,…,12), 为加强方位聚焦效果,在每个方位电极内镶嵌方位监督电极Mazi(i=1,…,12).每个方位电极对应圆心角度为20°,绝缘部分对应角度为10°.在主电极A0两侧对称分布着6对屏蔽电极A1、A1′、A2、A2′、A3、A3′、A4、A4′、A5、A5′、A6和A6′.为了增加聚焦效果,在主电极A0两侧增加对称分布的4对监督电极,从M1、M1′到M4、M4′.为了减少井眼影响,提高方位的探测能力,将方位电极、方位监督电极及监督电极M1、M1′和M2、M2′安装在12 个极板上,利用推靠臂贴井壁测量.图 1b为单个极板电极分布图,极板中间位置为方位监督电极,向外为方位电极,监督电极M1、M1′和M2、M2′对称分布在两端.

图 1 阵列方位侧向测井电极系结构 (a)电极阵列;(b)极板结构. Fig. 1 Structure of array azimuthal laterolog device (a) Electrode array; (b) Pad structure.

具体的工作原理如下:

模式Ⅰ:主电极A0发射主电流I0(实际为12个方位电极发射电流,下同),屏蔽电极A1、A1′发射与主电流同相位屏蔽电流,电流回流到屏蔽电极A2、A2′-A6、A6′,聚焦条件:VM1(M1′)=VM2(M2′),电流返回电极A2、A2′-A6、A6′等电位,探测深度最浅.

模式Ⅱ:主电极A0发射主电流I0,屏蔽电极A1、A1′、A2、A2′发射与I0同相位屏蔽电流,电流回流到屏蔽电极A3、A3′-A6、A6′,聚焦条件:VM1(M1′)=VM2(M2′)VM3(M3′)=VM4(M4′),电流返回电极A3、A3′-A6、A6′等电位.电流屏蔽作用较模式Ⅰ加强,探测深度增加.

模式Ⅲ:主电极A0发射主电流I0,屏蔽电极A1、A1′-A3、A3′发射与I0同相位屏蔽电流,电流回 流到屏蔽电极A4、A4′-A6、A6′,聚焦条件: VM1(M1′)=VM2(M2′)VM3(M3′)=VM4(M4′)VA2(A2′)=VA3(A3′),电流返 回电极A4、A4′-A6、A6′等电位.探测深度进一步加深.模式Ⅳ:主电极A0发射主电流I0,屏蔽电极A1、A1′-A4、A4′发射与I0同相位屏蔽电流,电流回 流到A5、A5′、A6、A6′,聚焦条件: VM1(M1′)=VM2(M2′)VM3(M3′)=VM4(M4′)VA2(A2′)=VA3(A3′)=VA4(A4′),电流返回电极A5、A5′、A6、A6′等电位.

模式Ⅴ:主电极A0发射主电流I0,屏蔽电极A1、A1′-A5、A5′发射与I0同相位屏蔽电流,电流回 流到A6、A6′,聚焦条件: VM1(M1′)=VM2(M2′)VM3(M3′)=VM4(M4′)VA2(A2′)=VA3(A3′)=VA4(A4′)=VA5(A5′).探测深度最大.

方位探测模式:在实现上述五个工作模式过程中,同时实现方位探测模式.主电极A0上12个方位电极各自独立发射电流Ii(i=1,…,12),通过调节电流Ii大小,在满足各模式聚焦条件同时还满足方位监督电极等电位,且与监督电极M1、M1′电位相同.各方位电流在轴向上受到屏蔽电流的聚焦作用,在周向上受到相邻方位电流的聚焦作用,实现方位探测.

3 数值模拟方法与地层模型

采用三维有限元数值模拟方法计算阵列方位侧向测井响应.建立图 2所示的地层模型,地层划分为井眼、侵入带、原状地层和围岩四部分,在柱坐标系(rφz)中电位U(rφz)满足如下偏微分方程的边值问题:

电阻率R分布满足

边值条件:

在地层远端满足第一类边界条件:

在对称边界和绝缘边界满足第二类边界条件:

在地层交界面处满足连续边界条件:

在主电极、屏蔽电极表面满足等位面边界条件(李大潜等,1980):

在监督电极表面满足等位面边界条件:

其中,n表示边界法向方向,“ - ”、 “ + ”分别表示地层交界面两侧,i表示电极数量,Ci表示某一待定常数,ρi为对应电极半径,Γi为对应电极表面.

图 2 AALS数值模拟采用的地层模型 Fig. 2 Stratum model for AALS numerical simulation
4 探测特性及影响因素分析 4.1 井眼影响

选定Rm=0.1 Ωm,Rt变化范围从0.01~105 Ωm,井眼直径d从5~14 in(1 in=25.4 mm),图 3a图 3b分别给出了贴井壁条件下模式Ⅰ和模式Ⅴ的井眼影响图版,Ra1Ra5分别是模式Ⅰ、模式Ⅴ的视电阻率.整体上,井眼校正系数随着Rt增大先增大后减小,Rt在1~10 Ωm范围内,井眼校正系数达到最大值;随着井眼直径增大,井眼校正系数也增大.在整个Rtd 范围内,井眼校正系数控制在0.70~1.15,说明未经井眼校正处理的视电阻率范围在0.87Rt到1.43Rt.

图 3 井眼影响图版(1 in=25.4 mm) (a) 贴井壁,模式Ⅰ; (b) 贴井壁,模式Ⅴ; (c) 不贴井壁,模式Ⅰ; (d) 不贴井壁,模式Ⅴ; (e) 浅侧向; (f) 深侧向. Fig. 3 Borehole effect chart (a) Sidewall way, model Ⅰ; (b) Sidewall way, model Ⅴ; (c) Non-sidewall way, model Ⅰ; (d) Non-sidewall way, model Ⅴ; (e) LLS; (f) LLD.

保持电极系尺寸不变,图 3c图 3d分别给出了不贴井壁条件下模式Ⅰ和模式Ⅴ的井眼影响图版,与双侧向测井中浅侧向(视电阻率为RLLS)和深侧向(视电阻率为RLLD)类似,如图 3e图 3f.对比图 3a3d可知,采用贴井壁方式,主电流直接进入地层,视电阻率受井眼影响较小,而且当井眼尺寸较大时,浅探测模型Ⅰ中井眼影响曲线右端没有上升趋势.

4.2 径向探测特性

利用径向积分几何因子J表征径向探测深度,将J为0.5时对应的侵入半径ri作为测井仪器的探测深度(张庚骥,1996).模拟低侵地层,选择Rm=0.1 Ωm,Rxo=1 Ωm,Rt=10 Ωm,计算五种探测模式的径向积分几何因子,结果如图 4a所示.由图中看出,五种模式的探测深度依次增加,以J=0.5计算,最大探测深度为1.23 m.最大探测深度小于双侧向中深侧向的探测深度,这是因为径向探测深度与纵向分辨率相互制约,径向探测深度、纵向分辨率和电极系长度需要综合考量(Chen et al.,1998).

图 4 径向积分几何因子曲线 (d=8 in,Rm=0.1 Ωm,Rxo=1 Ωm) (a) 五种模式(Rt=10 Ωm); (b) 模式Ⅴ(Rt/Rxo=2,5,10,20,50). Fig. 4 Curves of radial integration geometrical factors(d=8 in, Rm=0.1 Ωm, Rxo=1 Ωm) (a) Model I to model Ⅴ with Rt=10 Ωm; (b) Model Ⅴ with Rt/Rxo=2,5, 10,20,50.

探测深度还受Rt/Rxo比值的影响,以模式Ⅴ为例,选择Rm=0.1 Ωm,Rxo=1 Ωm,图 4bRt/Rxo分别为2、5、10、20、50时对应的径向积分几何因子曲线.从图中看出,Rt/Rxo=2时对应的探测深度达到1.40 m,Rt/Rxo=5时,探测深度减小为1.21 m.

4.3 纵向探测特性

纵向分辨率主要指电极系对地层的分层能力.选定d=8 in,Rm=0.1 Ωm,围岩电阻率Rs=1 Ωm,为了与不贴井壁测量方式对比分析,选择目的层电阻率RtRs的比值分别为2、5、10、20、50、100、200、500、1000、2000、5000、10000等12个值.图 5给出了贴井壁测量方式模式Ⅰ(图 5a)和模式Ⅴ(图 5b)的分辨率曲线,从图中可以看出,目的层厚度H逐渐增大,围岩影响系数急剧减小,当围岩厚度为0.3 m时,Rt/Ra1Rt/Ra5的值(<1.3)约等于1,探测到目的层真电阻率.当围岩厚度大于0.3 m时,曲线近似为水平线,因此纵向分辨率为0.3 m.

图 5 纵向分辨率曲线(d=8 in,Rm=0.1 Ωm,Rs=0.1 Ωm) (a) 贴井壁,模式Ⅰ; (b) 贴井壁,模式Ⅴ; (c) 不贴井壁,模式Ⅰ; (d)不贴井壁,模式Ⅴ. Fig. 5 Vertical resolution curves (d=8 in, Rm=0.1 Ωm, Rs=0.1 Ωm) (a) Sidewall way, model Ⅰ; (b) Sidewall way, model Ⅴ; (c) Non-sidewall way, model Ⅰ; (d) Non-sidewall way, model Ⅴ.

另外,从图中还可以看出,Rt/Rs从2变化到10000,模式Ⅰ和模式Ⅴ的分辨率没有变化,皆为 0.3 m,说明采用贴井壁测量方式,纵向分层能力基 本不受围岩电阻率大小的影响,为了对比分析,做出采用不贴井壁测量方式对应的分辨率曲线,如图 5c图 5d所示.从图中可以看出,随着Rt/Rs增大,围岩影响系数对应的最低点向目的层厚度增大方向偏移,与双侧向测井围岩影响系数变化趋势一致(朱军等,2007).当Rt/Rs=10000时,目的层厚度为 1.0 m对应的 Rt/Ra 约为1,此时分辨率约为0.8 m. 经过对比可以得出,采用贴井壁方式,测量电流受泥浆和围岩分流作用很小,在分层能力上较不贴井壁测量有独特的优势.

选择d=8 in,Rm=0.1 Ωm,Rs=1 Ωm,地层厚度H依次为0.1、0.2、0.3、0.4、0.5、0.3、0.3、0.3、0.2、0.2 m,无泥浆侵入,第1到第5个地层Rt=10 Ωm,第6到第8个地层Rt依次为10、25、50 Ωm,第9、10地层Rt分别为5、10 Ωm.图 6为AALS在不同层厚地层中的响应曲线.从图 6中可以看出,AALS在地层厚度为0.1 m时也有明显响应,响应值约为0.2Rt;当地层厚度为0.3 m,即能探测到地层真电阻率.对后5个地层(模拟连续薄层)也保持良好的探测特性.

图 6 不同地层厚度下的AALS测井响应(无泥浆侵入) Fig. 6 Log responses of AALS for different stratum thickness (without mud invasion)
4.4 方位分辨率

方位分辨率考察方位电极的周向分辨能力,以角度表示.每个方位电极圆心角为20°,考察其中一个方位电极对方位地层的分辨能力,如图 7所示.选择Rm=0.1 Ωm,地层电阻率Rt为1 Ωm时,方位高阻层,方位电极对应地层电阻率RtazRt的比值分别为2、5、10、20、50;地层电阻率Rt为10 Ωm时,方位低阻层,Rtaz/Rt分别为0.1、0.05、0.01(Rtaz=Rm)、0.005、0.001,结果如图 8所示.

图 7 方位分辨率模型 Fig. 7 Azimuthal resolution model for AALS

图 8 方位分辨率曲线(d=8 in,Rm=0.1 Ωm) (a) Rt=1 Ωm; (b) Rt=10 Ωm. Fig. 8 Azimuthal resolution curves (d=8 in, Rm=0.1 Ωm)

图 8a中可以看出,φ从4°逐渐增加,模拟得到的方位电阻率Raaz缓慢增大,小于Rt,当φ从18°增大到20°时,Raaz发生急剧变化,以Rtaz/Rt=10为例,此时Raaz=9.5 Ωm 约等于Rtaz,当φ继续增大时,Raaz曲线增幅很小,趋于平缓.其原因是,随着φ从4°到40°,方位电极电压一直缓慢增加,而该方位电极电流在φ从18°变化到20°时,即方位电极边界正对方位地层边界时,方位电流大幅下降,然后一直缓慢下降.从图 8b可以得出,对于方位低阻地层,当Rtaz< Rm时,φ从18°增大到20°,Raaz才能发生较明显变化.模拟发现,随着φ增加,该方位电流值一直缓慢增加,电压值在φ=20°出现与Raaz相似的波动.从图 8中还可以得出,方位地层电阻率异常程度越高,与井周其他方位电阻率相差越大,电阻率曲线变化越明显,越容易识别方位地层.

4.5 井周电阻率测量

普通双侧向和阵列侧向电极系无方位电极,测量值是地层周向电阻率的加权平均响应(杨韡,2002).利用贴井壁方位电极实现井周电阻率测量,选择Rm=0.1 Ωm,各方位电极对应方位地层电阻率Rtaz大小不同,模拟各方位电极的电阻率响应,结果如图 9所示.图 9a模拟井周电阻率连续变化,图 9b模拟井周电阻率在某一方位发生较小幅度变化,图 9c模拟井周电阻率发生急剧减小,图 9d模拟井周电阻率发生急剧增大.从图 9中看出,测量得到方位电阻率Raaz曲线与Rtaz曲线变化趋势吻合,即使井周电阻率发生急剧变化(图 9c图 9d),RaazRtaz也 非常接近.图 9c中,5号方位电极对应Rtaz=10 Ωm,Raaz= 9.7 Ωm,误差为3%,8号方位电阻率对应 Rtaz=5 Ωm,Raaz=5.3 Ωm,误差为6%.因此,采用贴井壁方位电极可以很好地测量井周方位电阻率,而且对于井周电阻率急剧变化情况也有很好的测量效果,与现有文献结论一致(黄继贞等,2001杨韡,2002).

图 9 井周电阻率曲线(d=8 in,Rm=0.1 Ωm) (a) 连续变化; (b) 小幅度降低; (c) 急剧降低; (d)急剧增大. Fig. 9 Resistivity curves for different azimuths surrounding the borehole (d=8 in, Rm=0.1 Ωm) (a) Continuous change from 1 Ωm to 12 Ωm; (b) Slow decrease from about 100 Ωm to 50 Ωm; (c) Sharp drop from about 800 Ωm to around 10 Ωm; (d) Sharp increase from about 10 Ωm to above 100 Ωm.
4.6 非对称侵入响应

水平井泥浆侵入形态与垂直井泥浆侵入形态不 同,成“泪滴”型或“椭圆型”(周灿灿和王昌学,2006),如图 10a所示,随着侵入时间增加,井周泥浆侵入深度不均,因此水平井和垂直井的泥浆侵入电阻率响应不同.建立如图 10b所示的水平井泥浆侵入模型,井眼与侵入带剖面为不同心的圆状,对应中心点分别为O、O′,初始距离OO′=0.025 m,随着侵入半径ri增加,点O′向下偏移,与井眼中心的距离增大.3、4号方位电极正对最浅侵入一侧,9、10号方位 电极正对最深侵入一侧.选择d=8 in,Rm=0.1 Ωm,低侵地层:Rxo=1 Ωm,Rt=10 Ωm;高侵地层:Rxo=10 Ωm,Rt=1 Ωm,图 11为模式Ⅰ、模式Ⅴ的响应曲线.

图 10 水平井泥浆侵入 (a) 水平井中不同时间的泥浆侵入形态; (b) 水平井泥浆侵入模拟剖面图. Fig. 10 Schematic of mud invasion in horizontal well (a) Configuration of mud invasion in horizontal well; (b)Simulation profile for mud invasion in horizontal well.

图 11 AALS方位电阻率随侵入半径变化曲线(d=8 in,Rm=0.1 Ωm) 低侵地层(Rxo=1 Ωm, Rt=10 Ωm): (a)模式Ⅰ; (b) 模式Ⅴ. 高侵地层(Rxo=10 Ωm, Rt=1 Ωm): (c) 模式Ⅰ; (d)模式Ⅴ. Fig. 11 Azimuthal resistivity changes vs. invasion radius for AALS(d=8 in, Rm=0.1 Ωm) Low invasion (Rxo=1 Ωm, Rt=10 Ωm): (a)Model Ⅰ; (b)Model Ⅴ. High invasion (Rxo=10 Ωm, Rt=1 Ωm): (c)Model Ⅰ; (d)Model Ⅴ.

图 11a图 11b中可以看出,在低侵地层中,ri=0.05 m时,非对称侵入刚开始,测量电阻率值接近Rt,模式Ⅰ中3、4号电极测量电阻率值最大,9、10号电极测量电阻率值最小,最大值和最小值相差3.3 Ωm,模式Ⅴ中也出现类似效果,最大值和最小值相差2.5 Ωm.按照从3、4号电极到9、10号电极的顺序,其他方位电极测量电阻率值由大到小变化.随着侵入深度增加,测量电阻率值逐渐减小并接近Rxo,电阻率最大值和最小值之差逐渐减小.当ri=1 m时,模式Ⅰ中最大值与最小值之差为0.3 Ωm,模式Ⅴ中最大值与最小值之差为0.7 Ωm.

图 11c图 11d中可以看出,在高侵地层中,ri=0.05 m时,测量电阻率值接近Rxo,模式Ⅰ、模式Ⅴ中最大值与最小值分别相差0.5,0.8 Ωm,随着侵入深度增大,最大值与最小值之差变化不大,这与低侵地层中的响应不同,随着ri增加,测量电阻率值逐渐增大并接近Rxo.

另外,图 11Ra1Ra5分别为模式Ⅰ、模式Ⅴ中不考虑方位测量得到电阻率曲线.Ra1Ra5曲线位于12条方位电阻率曲线中间位置,反映出地层周向电阻率的平均大小,不能确定侵入最大和最小方向,提供的信息缺乏方位性.因此方位测量有利于侧向测 井在水平井中应用(Davies et al.,1994Smits et al.,1995).

4.7 地层倾角响应

当电极系穿过倾斜地层时,利用方位电极阵列可以反映地层倾角和倾斜方向,以下以模式Ⅰ为例.在图 2中,电极系轴向方向与倾斜地层法向方向的 夹角为地层倾角θ(邓少贵等,2010b),选择 井眼直径d=8 in,Rm=0.1 Ωm,无侵倾斜地层电阻率Rt=10 Ωm,围岩电阻率Rs=1 Ωm,倾斜地层厚度为1 m.

图 12是地层倾角分别为0°、10°、45°、80°、90°时,利用12个方位电极合成的电阻率曲线(无方位性),横坐标为电极系中心与倾斜地层中心的垂直距离.从图中可以看出,随着地层倾角增大,计算出的电阻率曲线中值减小,与双侧向和阵列侧向测井在倾斜地层中响应一致(Wang and Shen,2001; 邓少贵等,2010a).以半幅点计算,倾角为0°时,视地层厚度即为地层真实厚度1 m,倾角增大到80°,视地层厚度为4.5 m,与地层真实厚度1 m相差很大,即地层倾角越大,响应越不灵敏.

图 12 倾斜地层12个方位电极合成电阻率响应(d=8 in,Rm=0.1 Ωm,H=1 m) Fig. 12 Synthetic resistivity responses using twelve azimuthal electrodes in dipping strata (d=8 in, Rm=0.1 Ωm, H=1 m)

虽然利用12个方位电极合成的电阻率曲线能够反映地层倾角,与双侧向、阵列侧向测井在倾斜地层中响应一致,但是由于不具有方位性,不能反映出地层倾斜方向(倾斜地层法线绕井轴旋转).

利用方位电极的方位性,根据方位电极提供的方位电阻率曲线变化规律可以判断地层倾斜方向.以倾角80°为例,模拟条件同上.图 13中曲线(标号为A)是模拟得出的12条方位电阻率曲线.从图中可以看出,7号方位电极电阻率首先发生变化,表示倾斜地层首先与7号方位电极接触.随后在7号电极两侧的方位电极6号、5号、4号和8号、9号、10号按先后顺序发生变化,而且与7号电极越远,变化幅度逐渐减小.随着倾斜地层移动,12条方位电阻率曲线形态与合成电阻率曲线(图 12)一致.在曲线另一侧,当电极系中心逐渐远离倾斜地层时,1号方位电极最后达到最高值,说明倾斜地层最后从1号方位电极位置与方位电极分离.据此可以判断地层的倾斜方向.同理,可以判断倾斜方位低阻层(Rt=1 Ωm,Rs=10 Ωm,其他条件同上),如图 13中曲线(标号为B).

图 13 80°倾斜地层中12个方位电极的电阻率响应 Fig. 13 Resistivity responses of twelve azimuthal electrodes in 80° dipping strata

选定1号、4号、7号、10号方位电极分别对应正东、正北、正西、正南方向,根据电极方位性及方位曲线形态确定地层倾斜方向,地层倾斜角分别为10°、45°、80°,倾斜方向分别为北偏东60°、南偏东30°和正西方向,利用最先、最后幅度发生变化的方位电阻率确定地层倾斜方向,结果如图 14所示.从图 14还可以看出,当倾角为10°时,曲线凸起值小于 中点值;当倾角为45°时,两者接近;当倾角为80°时,曲线凸起值大于中点,据此可以反映出地层倾角大小,结合倾角响应曲线(图 12),可以判断倾斜地层厚度.

图 14 倾斜地层中AALS方位响应 (a) θ=10°; (b) θ=45°; (c) θ=80°. Fig. 14 Azimuthal responses of AALS in dipping strata
5 结论

本文提出了一种贴井壁式阵列方位侧向测井电极系,介绍了电极系结构和工作原理,利用三维有限元数值模拟方法研究了电极系的探测特性,得到以下结论:

(1)与常规不贴井壁测量对比,贴井壁测量可以减小井眼影响.井眼直径从5~14 in,AALS的井眼影响系数控制在0.7~1.15范围内.

(2)AALS实现了5种探测深度,在低侵地层中最大探测深度为1.23 m.另外,侵入带电阻率与目的层电阻率的比值影响探测深度的大小,侵入带电阻率越小,探测深度越小.

(3)AALS在无侵中可以识别0.1 m薄层,当层厚为0.3 m时,即可准确测量地层电阻率.对比不贴井壁测量方式,采用贴井壁测量方式使得测量电流直接进入地层,大幅度降低泥浆和围岩对纵向分层能力的影响.目的层电阻率与围岩电阻率比值从2变化到10000,纵向分辨率基本不变.

(4)AALS利用12个方位电极,方位分辨率为20°.采用贴井壁的测量方式,可以准确测量出井周不同方位的电阻率,而且对井周电阻率异常变化的低阻和高阻地层也可以准确响应.

(5)利用方位电阻率曲线可以判断水平井泥浆非对称侵入.12个方位电阻率中,最大值与最小值之差,随着侵入半径增加,在低侵地层中逐渐减小,在高侵地层中基本保持不变.

(6)在倾斜地层中,利用12条方位电阻率曲线先后变化形态可以判断地层倾角大小和倾斜方向.倾斜地层“迈过”方位电极前后,对应方位电阻率曲线形态发生突变,倾角越大,突变越剧烈.

参考文献
[1] Al-Waheed H H, Audah T M. 1994. Applications of the azimuthal resistivity imager tool in Saudi Arabia.//SPE 28439, Annual Technical Conference and Exhibition. New Orleans, Louisiana:Society of Petroleum Engineers Inc., 811-818.
[2] Brambilla F, Corley B, Garcia A, et al. 2013. Improved laterolog resistivity measurements provide real-time true formation resistivity and invasion parameters.//OMC-2013-100, Offshore Mediterranean Conference and Exhibition. Ravenna, Italy:Offshore Mediterranean Conference, 1-11.
[3] Chen Y H, Chew W C, Zhang G J. 1998. A novel array laterolog method. The Log Analyst, 39(5):23-33.
[4] Davies D H, Faivre O, Gounot M T, et al. 1994. Azimuthal resistivity imaging:A new generation laterolog. SPE Formation Evaluation, 9(3):165-174.
[5] Deng S G, Tong Z Q, Fan Y R. 2006. Numerical simulation of dual laterolog response in tilted anisotropic formation. Acta Petrolei Sinica (in Chinese), 27(3):61-64.
[6] Deng S G, Li Z Q. 2009. Simulation of array laterolog response of fracture in fractured reservoir. Earth Science-Journal of China University of Geosciences(in Chinese), 34(5):841-847.
[7] Deng S G, Xu R W, Jiang J L. 2010a. Study on the array laterolog response of heterogeneous formation in deviated well. Well Logging Technology (in Chinese), 34(2):130-134.
[8] Deng S G, Li Z Q, Fan Y R, et al. 2010b. Numerical simulation of mud invasion and its array laterolog response in deviated wells. Chinese J. Geophys. (in Chinese), 53(4):994-1000, doi:10.3969/j.issn.0001-5733.2010.04.024.
[9] Deng S G, Mo X X, Lu C L, et al. 2012. Numerical simulation of the dual laterolog response to fractures and caves in fractured-cavernous formation. Petroleum Exploration and Development (in Chinese), 39(6):706-712.
[10] Faivre O. 1993. Fracture evaluation from quantitative azimuthal resistivities.// SPE 26434, Annual Technical Conference and Exhibition. Houston, Texas:Society of Petroleum Engineers Inc.,179-192.
[11] Forsyth D, Nawawi H, Ho T C. 1993. Review of techniques for the interpretation and evaluation of thin sand sequences.//SPE 25357, SPE Asia Pacific Oil and Gas conference. February, Singapore:Society of Petroleum Engineers Inc., 347-355.
[12] Griffiths R, Smits J W, Faivre O, et al. 1999. Better saturation from new array laterolog.// SPWLA 40th Annual Logging Symposium.Oslo, Norway:Society of Petrophysicists & Well Log Analysts, 1-14.
[13] Hakvoort R G, Fabris A, Frenkel M A, et al. 1998. Field measurements and inversion results of the high-definition lateral log.// SPWLA 39th Annual Logging Symposium.Keystone, Colorado:Society of Petrophysicists & Well Log Analysts, 1-12.
[14] He F, Ma X, Feng L W, et al. 2013. HAL 6505 array laterolog logging tool. Oil Forum (in Chinese),(2):59-62.
[15] Hu H T, Bai Q J, Xiao Z S, et al. 2014. Numerical simulation of high resolution laterolog tool based on computed focusing. Progress in Geophysics(in Chinese), 29(5):2178-2182, doi:10.6038/pg20140528.
[16] Huang J Z, Zhu J, Li M X, et al. 2000. A survey on azimuthal resistivity imaging abroad. Well Logging Technology(in Chinese), 24(3):170-175.
[17] Huang J Z, Zhu J, Li M X, et al. 2001. A new azimuthal resistivity imaging tool. China Patent:CN00226552.
[18] Itskovich G B, Mezzatesta A G, Strack K M, et al. 1998. High-definition lateral log-resistivity device:basic physics and resolution.// SPWLA 39th Annual Logging Symposium.Keystone, Colorado:Society of Petrophysicists & Well Log Analysts, 1-12.
[19] Li D Q, Zheng S M, Tan Y J, et al. 1980. Application of Finite Element Method in Electric Well Logging(in Chinese). Beijing:Petroleum Industry Press, 77-79.
[20] Li S J, Martin L E-S, Bittar S M, et al. 2013. Multi-array laterolog tools and methods with differential voltage measurements. US Patent:20130234718 A1.
[21] Nam M J, Pardo D, Torres-Verdín C. 2009. Simulation of DC dual-laterolog measurements in complex formations:A Fourier-series approach with nonorthogonal coordinates and self-adapting finite elements. Geophysics, 74(1):E31-E43.
[22] Nam M J, Pardo D, Torres-Verdín C. 2010. Assessment of Delaware and Groningen effects on dual-laterolog measurements with a self-adaptive HP finite-element method. Geophysics, 75(6):F143-F149.
[23] Pan K J, Wang W J, Tang J T, et al. 2013. Mathematical model and fast finite element modeling of high resolution array lateral logging. Chinese J. Geophys.(in Chinese), 56(9):3197-3211, doi:10.6038/cjg20130932.
[24] Shi G, He T, Wu Y Q, et al. 2004. A study on the dual laterolog response to fractures using the forward Numerical modeling. Chinese J. Geophys.(in Chinese), 47(2):359-363.
[25] Smits J W, Benimeli D, Dubourg I, et al. 1995. High resolution from a new laterolog with azimuthal imaging.//SPE 30584, Annual Technical Conference and Exhibition. Dallas, Texas:Society of Petroleum Engineers Inc., 563-576.
[26] Smits J W, Dubourg I, Luling M G, et al. 1998. Improved resistivity interpretation utilizing a new array laterolog tool and associated inversion processing.//SPE 49328, Annual Technical Conference and Exhibition. New Orleans, Louisiana:Society of Petroleum Engineers Inc., 831-843.
[27] Suau J, Grimaldi P, Poupon A, et al. 1972. The dual Laterolog-Rxo tool.//SPE 4018, Fall Meeting of the Society of Petroleum Engineers of AIME. San Antonio, Texas:Society of Petroleum Engineers Inc., 1-44.
[28] Tan M J, Gao J, Zou Y L, et al. 2012. Environment correction method of dual laterolog in directional well. Chinese J. Geophys. (in Chinese), 55(4):1422-1432, doi:10.6038/j.issn.0001-5733.2012.04.038.
[29] Tong M S, Song J H. 2014. Numeric simulation of 0.2 m vertical-resolution dual laterolog tool. Progress in Geophysics(in Chinese), 29(5):2251-2257, doi:10.6038/pg20140539.
[30] Wang H M, Shen L C. 2001. Dual laterolog responses in anisotropic crossbedding formation. Petrophysics, 42(6):624-632.
[31] Woodhouse R. 1978. The laterolog Groningen phantom can cost you money.//SPWLA 19th Annual Logging Symposium. El Paso, Texas:Society of Petrophysicists & Well Log Analysts, 1-12.
[32] Yang W. 2002. ARI log responses to 3D anisotropic formations. Well Logging Technology(in Chinese), 26(1):30-34.
[33] Yin C F, Ke S Z, Xu W, et al. 2014. 3D laterolog array sonde design and response simulation. Applied Geophysics, 11(2):223-234.
[34] Zhang G J. 1996. Electrical Logging(in Chinese). Dongying:Petroleum University Press.
[35] Zhou C C, Wang C X. 2006. Technology review on the log interpretation of horizontal well. Progress in Geophysics(in Chinese), 21(1):152-160.
[36] Zhou Z Q, Corley B H, Khokhar R, et al. 2008. Adaptive borehole corrections accounting for eccentricity for array laterologs.//SPE 114704, Annual Technical Conference and Exhibition. Denver, Colorado, USA:Society of Petroleum Engineers Inc., 1-18.
[37] Zhu J, Feng L W. 2007. Analysis on numeric simulation of high-resolution dual laterolog(DLL) response. Oil Geophysical Prospecting(in Chinese), 42(4):457-462.
[38] Zhu J, Feng L W, Li J H, et al. 2007. A new high resolution dual laterolog logging method. Well Logging Technology(in Chinese), 31(2):118-123.
[39] 邓少贵, 仝兆岐, 范宜仁. 2006. 各向异性倾斜地层双侧向测井响应数值模拟. 石油学报, 27(3):61-64.
[40] 邓少贵, 李智强. 2009. 裂缝性储层裂缝的阵列侧向测井响应数值模拟. 地球科学——中国地质大学学报, 34(5):841-847.
[41] 邓少贵, 徐悦伟, 蒋建亮. 2010a. 倾斜井非均匀地层的阵列侧向测井响应研究. 测井技术, 34(2):130-134.
[42] 邓少贵, 李智强, 范宜仁等. 2010b. 斜井泥浆侵入仿真及其阵列侧向测井响应数值模拟. 地球物理学报, 53(4):994-1000, doi:10.3969/j.issn.0001-5733.2010.04.024.
[43] 邓少贵, 莫宣学, 卢春利等. 2012. 缝-洞型地层缝洞的双侧向测井响应数值模拟. 石油勘探与开发, 39(6):706-712.
[44] 贺飞, 马骁, 冯琳伟等. 2013. HAL6505阵列侧向测井仪. 石油科技论坛,(2):59-62.
[45] 胡海涛, 白庆杰, 肖占山等. 2014. 基于合成聚焦的高分辨率侧向测井仪数值模拟研究.地球物理学进展, 29(5):2178-2182, doi:10.6038/pg20140528.
[46] 黄继贞, 朱军, 李妙霞等. 2000. 国外方位电阻率测井调查分析. 测井技术, 24(3):170-175.
[47] 黄继贞, 朱军, 李妙侠等. 2001. 一种新型方位侧向测井仪. 中国:CN00226552.
[48] 李大潜, 郑宋穆, 谭永基等. 1980. 有限元素在电法测井中的应用. 北京:石油工业出版社, 77-79.
[49] 潘克家, 王文娟, 汤井田等. 2013. 高分辨率阵列侧向测井的数学模型及有限元快速正演. 地球物理学报, 56(9):3197-3211, doi:10.6038/cjg20130932.
[50] 史謌, 何涛, 仵岳奇等. 2004. 用正演数值计算方法开展双侧向测井对裂缝的响应研究. 地球物理学报, 47(2):359-363.
[51] 谭茂金, 高杰, 邹友龙等. 2012. 盐水泥浆条件下定向井双侧向测井环境校正方法研究. 地球物理学报, 55(4):1422-1432, doi:10.6038/j.issn.0001-5733.2012.04.038.
[52] 童茂松, 宋建华. 2014. 0.2 m分辨率双侧向测井仪器数值模拟. 地球物理学进展, 29(5):2251-2257, doi:10.6038/pg20140539.
[53] 杨韡. 2002. 三维各向异性地层中方位电阻率测井的响应. 测井技术, 26(1):30-34.
[54] 周灿灿, 王昌学. 2006. 水平井测井解释技术综述. 地球物理学进展, 21(1):152-160.
[55] 朱军, 冯琳伟. 2007. 高分辨率双侧向测井响应数值模拟分析. 石油地球物理勘探, 42(4):457-462.
[56] 朱军, 冯琳伟, 李剑浩等. 2007. 一种新型的高分辨率双侧向测井方法. 测井技术, 31(2):118-123.