随钻电磁波测井响应函数理论及探测特性

林发武, 田超国, 耿学杰, 司兆伟, 刘得芳, 庄东志

(中国石油冀东油田公司勘探开发研究院,河北唐山 064300)

随钻电磁波测井成为解决复杂油气储层开发的重要手段之一。如何掌握随钻电磁波测井探测特性及响应规律,获取其信号空间分布特征是关键。近十几年来部分学者对电磁波类测井探测机制进行了大量研究[1-7],在感应测井中几何因子理论能够有效指导仪器设计,定量分析响应特征等[8-11]。Doll几何因子将地层介质等效成由互不干扰的独立单元环组成,分析其对测井响应贡献的空间分布,以此来指导仪器设计及响应分析[12]。在此基础上,Howard等[13]将量子散射问题中的Born近似方法引入感应测井问题的求解,提出了二维Born几何因子;王磊等[14]针对多分量感应测井进行了三维Born几何因子的理论研究,分析了同轴分量及交叉分量在不同井斜角及地层各向异性下的响应规律。不同于感应测井直接探测电场信号,随钻电磁波类测井采用幅度比和相位差的刻度方式进行测量,传统的感应几何因子理论不再适用。邢光龙等[15-16]针对传播电阻率测井仪器提出了幅度比与相位差关于介电常数与电导率的响应函数,分析了响应信号空间分布规律,但未模拟地层各向异性及井斜角的影响,且未给出定量识别探测范围的方法。伪几何因子理论[17]可有效分析测井仪的径向探测深度,但对纵向及微观敏感性描述不足。目前国内外对于由幅度比、相位差刻度且对井周地层属性敏感的纵、横向微观及宏观几何贡献分布的精确描述研究较少。笔者从Born几何因子理论出发,从感应测井推广到随钻电磁波测井中,研究各向异性地层中电性异常体对接收信号的响应规律,提出一种能够刻画幅度比及相位差的响应函数方法,用于定量分析随钻电磁波类测井探测特性。

1 随钻电磁波响应函数

根据Born近似方法,当背景介质中存在异常体时,异常体内任意一点电场可以等效为源在背景介质中激发的一次场Eb(r)与异常体产生的散射场δE(r)的叠加:

E(r)=Eb(r)+δE(r).

(1)

其中背景场可通过解析方法获取[18],散射场可以表示为

(2)

式中,ω和μ0分别为角频率和地层磁导率;Vs为异常体边界;r和r′分别为接收点和异常点相对位置;Eb(r′)为发射线圈在异常体内某点处产生的电场;G(r,r′)为该异常点对接收信号的影响,可通过电偶极子的位函数推导得到。以x方向电偶极子源为例,其对应位函数为

(3)

其中

利用源的矢量位函数结合洛伦兹条件可以得到:

(4)

将式(3)带入式(4)可得

(5)

(6)

(7)

其中Gij代表i方向的电偶极子源在j方向产生的响应,因此异常体内某点产生的散射场可以表示为Eij=λ2EixGxj+gλ2EiyGyj+gizGzj,通过大地坐标系与仪器坐标系间变换关系,可以得仪器坐标系下接收线圈处散射场:

Es=Exxsin2φ+Ezzcos2φ+(Exz+Ezx)sinφcosφ.

(8)

随钻电磁波测井采用单发双收线圈结构(图1)，发射和接收线圈中心在同一轴线上，图中T表示发射线圈，R1和R2分别代表两接收线圈。利用两个电势信号V1和V2,计算两者之间的幅度比A和相位差Δθ来识别地层参数,可表示为

(9)

图1 随钻电磁波测井响应模型Fig.1 Modeling of electromagnetic logging while drilling response

对式(9)取微分可得随钻电磁波测井响应函数S:

(10)

整理可得

(11)

式中,dlnA和dΔθ分别为幅度比和相位差的响应函数,代表地层中任意地层单元对幅度比和相位差的贡献。

2 随钻电磁波测井探测特性分析

2.1 响应函数空间分布

研究随钻电磁波测井的响应函数空间分布特性,发射线圈位于坐标原点,信号频率为2 MHz,发射线圈到两个接收线圈的距离分别为0.610和0.762 m,地层背景电阻率为10 Ω·m。图2给出了随钻电磁波测井幅度比和相位差的响应等值面,其中红色表示所有响应点值的90%;图3为沿x方向与z方向的积分平面。等值面及z方向积分平面表明,随钻电磁波测井响应函数在径向上具有旋转对称性,在井轴附近地层贡献较小,能够减弱井眼的影响;从x方向积分平面可以看出,发射线圈到接收线圈间区域贡献较小,而响应函数贡献主要集中在两接收线圈间,这表明其纵向分层能力与线圈距密切相关。对比幅度比和相位差积分平面表明,二者响应函数分布规律大致相同,幅度比敏感性空间分布较发散,相位差敏感性分布更为集中。

图2 幅度比及相位差响应函数等值面Fig.2 Response isosurface of amplitude ratio and phase shift

2.2 井斜角影响

图3 幅度比及相位差响应函数积分平面图Fig.3 Response function integral plane of amplitude ratio and phase shift

图4 不同井斜角条件下响应函数等值面Fig.4 Responses isosurface with different dipping angle

2.3 各向异性系数影响

图5 不同各向异性条件下响应函数等值面Fig.5 Response isosurface with different anisotropic coefficient

3 探测特性的影响因素

3.1 响应函数积分几何因子

对积分平面分别沿x轴与z轴进行积分,由于响应函数在径向上具有旋转对称性,x正半轴几何因子与径向几何因子相等,得到径向及纵向微分几何因子,如图6所示。径向微分几何因子曲线在井轴附近出现小段负敏感区域,幅度比信号尤为明显,负敏感区将导致井轴附近区对响应贡献率降低,井周地层对响应信号的贡献降低;纵向微分几何因子曲线在发射线圈与近接收线圈之间有小段曲线敏感度不变,在接收线圈附近敏感度高,相位差微分几何因子分布更为集中,能够有效区分薄层。对比幅度比和相位差敏感性曲线可以发现,幅度比信号径向探测深度更广,相位差具有更好的纵向分层能力。

图6 响应函数径向及纵向微分几何因子曲线Fig.6 Radial and vertical differential geometry factor of response function

为定量分析随钻电磁波测井径向探测深度及纵向分辨率的变化规律,引入径向积分几何因子Sr与纵向积分几何因子Sh,通过对微分几何因子进行积分,计算得到幅度比与相位差的径向与纵向积分几何因子曲线,定义表达式如下:

(12)

式中,Sx和Sz分别为响应函数在x>0和z>0空间内贡献的总和。

图7为随钻电磁波测井幅度比和相位差信号的径向积分几何因子与纵向积分几何因子变化曲线,图7中径向及纵向积分几何因子坐标原点分别为井轴及两接收线圈中点。类比伪几何因子[17],此处定义Sr=50%时对应的径向距离为径向探测深度,并且假定Sh=90%时对应的距离为纵向分辨率。观察图7可知,幅度比和相位差的径向探测深度分别为0.97和0.51 m,分辨率分别为0.76 和0.34 m,这也验证了幅度比信号具有更大的径向探测深度,而相位差信号能够有效分辨薄地层。

图7 随钻电磁波测井积分几何因子曲线Fig.7 Integrated geometric factor curve of electromagnetic logging while drilling

3.2 信号频率的影响

利用径向及纵向积分几何因子,在不同信号频率条件下,分析随钻电磁波测井径向探测深度和纵向分辨率的变化规律。仪器发射线圈到两个接收线圈的距离分别为0.610和0.762 m,背景地层电阻率Rb=10 Ω·m,信号频率分别选取200 kHz、400 kHz、800 kHz、1 MHz及2 MHz。图8(a)和(b)分别表示随钻电磁波测井幅度比和相位差信号的径向积分几何因子曲线。从图8中可以看出,随着信号频率的增加,幅度比和相位差信号曲线的探测范围都随之降低。对比图8(a)和图8(b)可以发现,幅度比信号的探测范围受频率的影响较大,且幅度比径向积分几何因子在小于0.5 m时贡献接近于零,说明幅度比信号受此范围内地层影响较小。图8(a)中幅度比信号在频率为200 kHz时,其径向探测深度近2 m,而当频率增加到2 MHz时,径向探测深度缩减到1 m以内,相对变化明显;图8(b)中相位差信号的径向探测深度明显小于幅度比信号,均小于1 m,且探测范围变化较小。在相同条件下,信号频率越低,信号的探测范围越大,且幅度比信号径向探测能力要强于相位差信号。

图8 不同发射频率条件下径向积分几何因子曲线Fig.8 Radial integral geometric factor curves under different frequencies

图9(a)和(b)分别表示随钻电磁波测井幅度比和相位差信号的纵向积分几何因子曲线,图中随着发射信号频率的增加,随钻电磁波测井信号的纵向分辨率增强。信号频率为200 kHz时,幅度比信号和相位差信号的纵向分辨率分别为1.56 和0.81 m,表明在相同信号频率条件下,相位差信号具有更好的纵向分层能力;而当信号增加到2 MHz时,纵向分辨能力均增强2～3倍,相位差纵向分辨率可达0.33 m。

图9 不同发射频率条件下纵向积分几何因子曲线Fig.9 Vertical integral geometric factor curves under different frequencies

3.3 线圈距的影响

考虑到实际随钻电磁波测井常采用阵列化设计,为分析不同线圈距对探测特性的影响,参照ARC312随钻电磁波仪器线圈布置进行模拟,令仪器发射频率为2 MHz,地层背景电阻率为10 Ω·m,发射线圈到两接收线圈中点距离分别为0.254、0.406、0.559、0.711和0.864 m,相邻两接收线圈间隔0.152 4 m。图10和图11表示了幅度比和相位差曲线在不同源距条件下径向及纵向积分几何因子的变化趋势。观察两图可发现,随着线圈距的增大,仪器的径向探测深度响应逐渐增加,而纵向分辨率逐渐减小。通过该方法可有效评估不同源距条件下探测特性,为随钻电磁波测井采用阵列化的线圈系设计提供定量分析,利用不同源距的线性组合,降低井眼、围岩等因素的影响,有效获取地层参数。

图11 不同线圈距下纵向积分几何因子曲线Fig.11 Response function integral curve with horizontal distance

3.4 地层电阻率的影响

发射线圈到两个接收线圈的距离分别为0.610和0.762 m,信号频率为2 MHz不变,选取不同背景地层电阻率Rb分别为1、5、10、50和100 Ω·m,分析Rb变化对随钻电磁波测井信号探测能力的影响。如图12所示,当电阻率由1 Ω·m增大到100 Ω·m时,幅度比信号径向探测深度由0.5 m增加到1.9 m,相位差信号径向探测深度由0.3 m增加到0.8 m。在电阻率较小时,幅度比径向积分几何因子在井轴附近存在小段负敏感区,能够有效减弱井周地层影响,随钻电阻率的增加影响逐渐消失;相位差径向几何因子会在1 m附近出现明显凸起,表明地层相对贡献在此处达到最大值,而后地层会产生负贡献,随着地层电阻率增加,峰值逐渐降低且距仪器距离增加。图13表明,随着电阻率的增加,仪器的纵向分辨能力逐渐降低。对比图12与图13可以发现,随着电阻率的增大,仪器探测范围加大,径向探测能力增强,纵向分层能力下降。这符合电阻率越大,趋肤深度越大,电磁波衰减慢,纵向分辨率受到影响更大。

图13 不同背景地层电阻率条件下纵向积分几何因子曲线Fig.13 Vertical integral geometric factor curves under different formation resistivity

4 结 论

(1)所提出的随钻电磁波测井响应函数可有效表征随钻电磁波幅度比和相位差信号的空间探测特性,数值模拟表明两者空间分布规律类似,均具旋转对称性,相位差信号相对幅度比信号而言更为集中,在相同条件下,幅度比信号具有更大的探测范围,而相位差信号纵向分辨率更强。

(2)各向异性系数及井斜角的加大均明显改变随钻电磁波测井响应函数的空间分布:在斜井及水平井条件下,随着各向异性系数的加大,响应函数在横向上呈现不规则扁平状,在接收和发射线圈附近逐渐出现负贡献,且随地层各向异性系数的增加而扩大;随着井斜角的加大,响应函数旋转对称性逐渐消失。

(3)在不同信号频率、地层电阻率及线圈距条件下,随钻电磁波测井幅度比与相位差信号探测范围均出现规律变化。随着信号频率的增加,地层电阻率加大及线圈距的增长,仪器的径向探测深度会随之加大,纵向分辨率降低,该方法对于仪器的探测特性研究及线圈系设计有着重要作用。