井周界面电磁散射探测数值模拟

邓少贵, 张 盼, 王正楷, 袁习勇

(中国石油大学地球科学与技术学院,山东青岛 266580)

随钻电磁波测井是随钻地质导向与储层评价的重要方法[1-3]。传统随钻电磁波仪器发射天线和接收天线同轴设置,不具有方位探测能力,仪器工作频率为0.4～2 MHz,采用相位差和幅度比的方式刻度地层电阻率,获得地层界面信息,探测范围可达2～3 m[4]。随钻方位电磁波测井仪器,在多频、多源距基础上,采用同轴/倾斜/共面天线,实现对电磁场交叉分量的测量,仪器在高阻层且围岩为低阻层时,边界探测范围理论可达到5～6 m[5-9]。电磁波仪器界面探测主要受仪器测量精度、相对井斜角、地层电阻率以及界面两侧电阻率差异影响[10]。目前,电磁波测井基本靠发射和接收天线间的几何关系来分辨地层,即为了增加探测深度而增加仪器尺寸,但仪器过长会给实际应用带来不便[11]。井周存在地层界面时可以产生电磁散射,不同介质条件电磁散射和吸收情况不同,同种条件不同频率的电磁波也存在“色散”现象。散射波具有更好刻画地层细节的能力[12],地层边界处电磁波测井曲线“犄角”现象本质是电磁散射造成[13]。井中电磁波信号是一次场和散射场的叠加,但散射信号相对背景一次场弱很多(千分之几)。优化天线结构和参数,削弱地层背景一次场,提高异常体散射信噪比,可以增加随钻导向探边深度。笔者采用对称天线结构设计,提取井周异常体散射信号,根据不同发射频率、源距、电阻率对比度下信号响应特征,利用接收天线信号幅度差获得地层界面信息。

1 界面电磁散射响应计算方法

1.1 模型设定

采用单发双收对称天线结构(图1),发射和接收天线中心在同一轴线上,发射天线为T，两接收天线(R1、R2)到发射天线距离相等,线圈匝数均为1。在带井眼均匀介质中,发射天线产生电磁场具有对称性,两接收天线信号完全相同。当井周(含钻前)有界面存在时,电磁散射引起电磁场分布变化,散射信号到达两个接收天线时幅度和相位出现差异。仪器正前方(井眼)地层界面深度设为0,仪器到界面的垂直距离为D,发射天线到仪器前端距离为L,前探距离为d,d=D/cosα-L。除特殊说明外,以下算例均以此为模型。

图1 地层模型与线圈系结构Fig.1 Formation model and coil array configuration

1.2 计算方法

利用矢量电位法进行一维层状介质正演模拟算法的研究,将线圈等效为磁偶极子源(M),其电磁场空间分布可看成由水平磁偶极子源(Mx)和垂直磁偶极子源(Mz)单独产生的叠加。接收信号为背景介质一次矢量电位和界面电磁散射的二次矢量电位的叠加。引入矢量电位(F)描述磁偶极子源电磁场分布,F满足亥姆赫兹方程[14-15]:

(1)

式中,k、ω、μ和m分别为波数、角频率、磁导率和磁偶极子强度;r0和r分别为源点及测量点位置。

(2)

(3)

其中

ξ=(λ2-k2)1/2.

式中，λ为积分变量;ρ为源点到测量点的距离。当j层有源时,τj=iωμjmj;当j层无源时,τj=0。

(4)

(5)

水平分量分别为

(6)

(7)

式中,φ为仪器与地层夹角；Ai、Bi、Ci、Di为未知系数,可通过边界连续条件得到。根据式(2)～(7),计算得到分量信号:

(8)

(9)

由于磁场和矢量电位满足下式:

(10)

通过磁场求得感应电场为

E=iωμsH.

式中，s为线圈面积。

根据两个接收天线的感应电动势ER1、ER2,得到测量信号的幅度差和相位差,如

幅度差=|ER1|-|ER2|,

(11)

相位差=arg(ER1)-arg(ER2).

(12)

对比幅度差和相位差,可以识别地层界面。

1.3 界面响应

令发射和接收天线的源距为1 m,发射信号频率为400 kHz,介质Ⅰ和介质Ⅱ电阻率分别为100和10 Ω·m,仪器与地层界面垂直。两接收天线的信号幅度和相位如图2(a),当仪器离地层界面较远时,两接收天线的信号幅度和相位完全相同;当仪器靠近地层界面时,两接收天线中信号差异增大。如图2(b)所示,发射天线在界面附近时,接收天线信号差值达到峰值,其中信号幅度差约为7×10-7V,相位差约为4°;界面在仪器之前9 m或在仪器之后4.3 m,幅度差异降为1×10-9V。如采用实际仪器天线匝数(几十到上百匝),信号幅度会增大百倍以上。目前国外仪器的信号强度分辨能力为10 nV,该方法对界面探测范围可到达9 m,远高于现有仪器(线圈距243.84 cm)的探测深度5 m[16]。

图2 接收天线中信号及其差异Fig.2 Response and differences of two receives

2 影响因素模拟

2.1 发射频率

电磁波测井探测深度受到信号频率限制,也受介质电阻率本身性质影响。分别针对高阻地层和低阻地层的情况,探讨多频率下天线系统的探边特性：①高阻地层背景,介质Ⅰ、Ⅱ的电阻率分别为100和10 Ω·m,供电频率分别选取200、400、800 kHz和1 MHz;②低阻地层背景,介质Ⅰ、Ⅱ电阻率分别为10和1 Ω·m,发射频率分别选取40、80、100和200 kHz。假设源距为1 m,仪器与地层界面垂直。如图3所示,随着信号频率增大,接收天线信号差幅度明显增大,且信号变化率也明显加大,有利于信号检测,但信号衰减加快引起探测范围减小,不利于界面的远探测。

图3 不同频率下的信号Fig.3 Signal with different frequencies

供电频率选取需要考虑信号强度和信号变化速率两方面,也要考虑介质电阻率的影响。如在低阻背景条件下,发射频率高于100 kHz会造成边界探测距离明显减小;在高阻背景条件下,发射频率在800 kHz才会使边界探测距离明显减小;如使用相同的信号测量分辨率,低电阻率背景条件的探边距离要明显小于高阻背景条件的探边距离。所以根据地层电阻率不同需要选取不同频率信号,以达到好的探测效果。

2.2 线圈距

地层模型与2.1相同,发射频率分别为800和100 kHz,仪器源距分别选取0.5、1、1.2和2 m。如图4所示,随着源距增大,幅度差异信号峰值逐渐减小,但源距增大,界面处信号衰减变小,故离开界面位置较远时,接收天线信号差异较大,而使探测范围加大。考虑到信号强度、探测范围和线圈系长度等因素,源距从1、1.2 m增加到2 m,其探测深度增加效果并不明显。

图4 不同源距条件下的响应Fig.4 Response with different coil spacing

2.3 电阻率对比度

散射信号强弱受到界面两侧地层电阻率的影响。考虑两种条件:①高阻背景,背景介质Ⅰ电阻率为100 Ω·m,介质Ⅱ电阻率分别为10、20、30和50 Ω·m,天线源距为1 m,发射频率选取800 kHz;②低阻背景,背景介质Ⅰ电阻率为2、3、5和10 Ω·m,介质Ⅱ电阻率为1 Ω·m,发射频率为100 kHz。如图5所示,天线探测范围受到介质背景电阻率和电阻率对比度的影响,高阻地层信号衰减速度较慢,且天线探测范围较大;电阻对比度越大,界面散射效果越明显,差异信号幅度越大,探测范围加大。

图5 不同电阻率对比度下的信号Fig.5 Response with different resistivity contrasts

2.4 界面相对倾斜

当地层界面在天线旁边时,随着仪器轴与地层界面夹角减小,界面电磁散射轴向分量降低,会使与发射天线同轴设置的接收天线信号减弱。设置倾斜(45°)接收天线,结合同轴天线,分析不同倾斜界面电磁散射信号变化。介质Ⅰ的电阻率为100 Ω·m,介质Ⅱ电阻率为10 Ω·m;天线源距取1 m,发射频率取800 kHz。结果显示,当仪器与地层垂直时,接收天线差异信号幅度最大,随着夹角的减小,信号幅度和探边距离(地层界面到发射天线的距离)都显著降低,如图6(a);倾斜天线信号差异幅度和探测范围变化较小,如图6(b),随着仪器与界面夹角减小,信号差异先增大再减小。同轴天线和倾斜天线的响应差异反映了地层界面散射信号分量变化,与地层界面位置、角度、介质电阻率对比度有关。

3 界面相对位置的确定

对比同轴和倾斜天线的信号,根据轴向和径向分量之间关系,确定仪器与地层界面夹角。采用分量组合天线结构,如图7所示,R1和R2为同轴天线,R3和R4为倾斜天线,源距分别为0.8和1 m,发射频率为800和100 kHz。分别针对高阻地层和低阻地层情况进行分析:①高阻地层背景,介质Ⅰ、Ⅱ的电阻率分别为100和10 Ω·m;②低阻地层背景,介质Ⅰ、Ⅱ电阻率分别为10和1 Ω·m,仪器均位于介质Ⅰ中。

图6 不同夹角条件下的信号Fig.6 Response with different angles

图7 仪器线圈系结构Fig.7 Configuration of transmitters and receivers

根据接收天线R1与R2、R3与R4之间的信号差异确定界面相对位置和倾斜角度。如图8所示,横坐标为同轴天线信号差异强度,纵坐标为倾斜天线信号差异强度,竖实线表示相对倾斜角度,横虚线表示界面位置。当地层界面与天线轴线夹角较大时,模拟地层界面在仪器前方,交会图对相对倾斜角度的分辨能力相对较弱,但对界面位置确定有利。当地层界面与天线轴线夹角较小时,模拟地层界面在仪器侧方,交会图对夹角的分辨能力相对较强,但对界面位置响应相对降低。根据图1所示模型中仪器和地层界面位置关系,可预测界面在仪器前所在距离(准向前看)。

图8 界面位置识别图版Fig.8 Identification chart for location of interface

4 结 论

(1)采用天线对称布局方式,利用接收天线信号差,可以较好消除背景一次场,反映界面电磁散射特点,在有限的源距条件下达到更优的探测效果。

(2)探边距离受电阻率对比度、发射频率、天线倾角等因素影响,增大电阻率对比度、降低发射频率以及天线倾斜都会导致探边距离加大。

(3)利用交会图法识别仪器与界面的距离及夹角,识别效果简洁,可以准确地判断界面的位置,有效指导仪器准确着陆及地质导向。

(4)散射信号相对一次场信号弱很多,进一步提升散射信号强度和信噪比,是电磁波测井向前/远看的关键。

致谢衷心感谢哈里伯顿公司李善军博士的悉心指导。

[1] 宋殿光. 随钻电磁波测井仪的数值模拟及地质导向应用研究[D].北京:中国地质大学(北京),2014.

SONG Dianguang. The study on numerical simulation and geosteering application of electromagnetic logging while drilling tools[D]. Beijing: China University of Geosciences(Beijing),2014.

[2] EDWARDS M, WATCHARANANTAKUL R, FARID S, et al. Advanced high definition multi-layer geosteering technology to overcome horizontal drilling challenges in complex reservoirs: proceedings of Abu Dhabi International Petroleum Exhibition & Conference, Abu Dhabi, 7-10 November，2016[C].Dallas, TX:Society of Petroleum Engineers,2016.

[3] FARHAD K, MOHAMMAD R R, ALI S. Investigating induction log response in the presence of natural fractures[J]. Journal of Petroleum Science and Engineering, 2016,145:357-369

[4] RABINOVICH M, LE F, LOFTS J, et al. The vagaries and myths of look-around deep resistivity measurements while drilling[J]. Petrophysics, 2012,53(2):86-101.

[5] 王磊,范宜仁,黄瑞,等.各向异性介质多分量感应测井三维Born几何因子理论研究[J].物理学报,2015,64(23):438-448.

WANG Lei, FAN Yiren, HUANG Rui, et al. Three dimensional born geometrical factor of multi-component induction logging in anisotropic media[J]. Acta Physica Sinica, 2015,64(23):438-448.

[6] 魏宝君,王成园,党峰,等.用柱状成层各向异性介质的并矢Green函数模拟多分量感应测井仪器的响应[J].中国石油大学学报(自然科学版),2016,40(2):59-69.

WEI Baojun, WANG Chengyuan, DANG Feng, et al. Simulating responses of multi-component induction logging tools by dyadic Greens functions in cylindrically stratified anisotropic media[J]. Journal of China University of Petroleum(Edition of Natural Science), 2016,40(2):59-69.

[7] ZHOU J, LI H. Interpretation of azimuthal propagation resistivity measurements: modeling, inversion, application and discussion: proceedings of SPWLA 57th Annual Logging Symposium, Reykjavik, Iceland, 25-29 June, 2016 [C]. Houston, TX:Society of Petrophysicists and Well-Log Analysts, 2016.

[8] 宋殿光,段宝良,魏宝君,等.倾斜线圈随钻电磁波电阻率测量仪器的响应模拟及应用[J]. 中国石油大学学报(自然科学版),2014,38(2):67-74.

SONG Dianguang, DUAN Baoliang, WEI Baojun, et al. Basic theory of electromagnetic wave resistivity measurement while drilling tool with tilted antennas and its application for geosteering[J]. Journal of China University of Petroleum(Edition of Natural Science),2014,38(2): 67-74.

[9] CHRIS B, JOHB H, PETER E, et al. Navigating and imaging in complex geology with azimuthal propagation resistivity while drilling[R]. SPE 102637, 2006.

[10] 杨震,杨锦舟,杨涛.随钻方位电磁波仪器补偿测量方法研究[J]. 中国石油大学学报(自然科学版),2015,39(3):62-69.

YANG Zhen, YANG Jinzhou, YANG Tao. Research on azimuthal electromagnetic tool while drilling measuring method ofcompensation[J]. Journal of China University of Petroleum(Edition of Natural Science),2015,39(3):62-69.

[11] 杨锦舟,魏宝君,林楠.倾斜线圈随钻电磁波电阻率测量仪器基本原理及其在地质导向中的应用[J]. 中国石油大学学报(自然科学版),2009,33(1):44-49.

YANG Jinzhou, WEI Baojun, LIN Nan. Response simulation and application of electromagnetic wave resistivity measurement-while-drilling tool with tilted coil[J]. Journal of China University of Petroleum(Edition of Natural Science),2009,33(1):44-49.

[12] 罗明璋.脉冲电磁法理论研究与硬件实现[D]. 武汉:长江大学, 2012.

LUO Mingzhang. Research on pulsed electromagnetic method and the hardware implementation[D]. Wuhan: Yangtze University, 2012.

[13] 杨震,刘庆成,岳步江,等.随钻电磁波测井中极化角的形成机理及其影响因素模拟分析[J].测井技术,2010,34(3):210-214.

YANG Zhen, LIU Qingcheng, YUE Bujiang. On mechanism of polarization angle of electromagnetic logging while drilling and its influence factors simulation[J]. Well Logging Technology, 2010,34(3):210-214.

[14] 李虎.复杂介质随钻方位电磁波测井数值模拟与应用基础研究[D].青岛:中国石油大学,2013.

LI Hu. Simulation and applied basic research of directional electromagnetic LWD tool in complex media[D]. Qingdao:China University of Petroleum, 2013.

[15] 范宜仁,巫振观,王磊,等.大斜度井各向异性地层双感应测井响应特征研究[J]. 测井技术,2016,40(3):262-269,275.

FAN Yiren, WU Zhenguan, WANG Lei, et al. Response characteristics of dual induction logging in anisotropic formation and high deviated wells[J]. 2016,40(3):262-269,275.

[16] 杨震,杨锦舟,韩来聚,等.随钻方位电磁波界面探测性能分析[J]. 石油学报,2016,37(7):930-938.

YANG Zhen, YANG Jinzhou, HAN Laiju, et al. Interface detection performance analysis of azimuthal electromagnetic while drilling[J]. Acta Petrolei Sinica,2016,37(7):930-938.