随钻前探测电磁波测井响应影响因素分析_田园诗.pdf

测 井 技 术WELL LOGGING TECHNOLOGYVol.46 No.6 Dec 2022第46卷 第6期 2022年12月文章编号：10041338（2022）06067607随钻前探测电磁波测井响应影响因素分析田园诗1，2，朱军1，2，牒勇1，2，鲁晔1，2，刘刚1，2，关浩3（1.中国石油集团测井有限公司测井技术研究院，陕西 西安 710077；2.中国石油天然气集团有限公司测井重点实验室，陕西 西安 710077；3.中国石油集团测井有限公司长庆分公司，陕西 西安 723507）摘要：随着油气勘探开发领域由构造油气藏向缝洞、薄层、断层等非常规复杂油气藏转移，随钻电磁波边界探测技术在主动地质导向与精细储层刻画方面发挥着重要作用，能够用于探边、避水、预警、前探盐丘、预判高压地层与枯竭层等。随钻前探测电磁波测井仪器采用正交或水平天线，与传统电磁波仪器相比，大大提高了前探距离。阐述了一种随钻前探测电磁波测井仪器的典型天线系统以及相应的前探测响应信号计算方法，应用三维几何因子、镜像源等效法分析了不同磁场分量的前探测敏感性。计算了不同地层电阻率下随钻前探测电磁波测井的响应特征，分析了随钻前探测的影响因素，针对界面两侧地层电阻率对比度、地层电阻率背景值等环境影响因素以及仪器频率及源距等参数进行了分析。数值模拟结果表明，正交或倾斜天线系统的前探距离受天线距离、工作频率、地层电阻率对比度等多种因素影响，通过优化发射-接收天线的排列方式、频率、天线距离等参数，能够获得较大的前探距离。关键词：随钻电磁波测井；前探测；前探距离；各向异性；数值模拟中图分类号：P631.84 文献标识码：ADoi:10.16489/j.issn.10041338.2022.06.006Analysis of the Influence Factors on Look-Ahead Electromagnetic Wave Logging While DrillingTIAN Yuanshi 1,2,ZHU Jun 1,2,DIE Yong 1,2,LU Ye1,2,LIU Gang1,2,GUAN Hao3(1.Logging Technology Research Institute,China National Logging Corporation,Xian,Shaanxi,710077;2.Well Logging Key Laboratory,China National Petroleum Corporation,Xian,Shaanxi,710077;3.Changqing Branch,China National Logging Corporation,Xian,Shaanxi,723507)Abstract:With the transfer of oil exploration from structural reservoirs to very complex unconventional reservoir,such as fracture cave,thin-layer and fault,electromagnetic wave boundary detection technology while drilling plays an important role in geological guidance and reservoir characterization,which can be used in edge exploration,water avoidance,early warning,salt dome pre-exploration,predetermination of high pressure formation and depleted layer,etc.The look-ahead electromagnetic wave resistivity logging while drilling tool adopts orthogonal or horizontal antenna,which greatly improves the forward detection distance compared with the traditional electromagnetic wave logging tool.We describe the method and antenna system of typical look-ahead tool,and analyses look-ahead sensitivity of different magnetic field components through three-dimensional geometric factor and mirror source equivalent method in this paper.The response characteristics of look-ahead tool under different formation resistivities are calculated,and several influencing factors of measurement while drilling are analyzed.The environmental factors such as formation resistivity contrast on both sides of the interface,formation resistivity background value,instrument frequency,source distance and other parameters are analyzed.The numerical simulation results show that the look-ahead detection distance of orthogonal or inclined antenna system is affected by many factors such as antenna spacing,working frequency and formation resistivity contrast.By optimizing the parameters such as transmitting receiving antenna arrangement,frequency and antenna spacing,a large look-ahead detection distance can be obtainedKeywords:electromagnetic wave logging while drilling;look-ahead;look-ahead detection distance;anisotropy resistivity;numerical simulation基金项目：中国石油天然气集团有限公司项目“随钻远探测电磁波电阻率成像测井仪器研制”（2021ZG04）第一作者：田园诗，女，1990 年生，工程师，从事电法测井方法研究工作。E-mail：田园诗，等：随钻前探测电磁波测井响应影响因素分析第46卷 第6期677 0 引 言随钻前探测电磁波测井能够测量钻头前方的未钻地层信息，被业内公认为未来10年极具发展潜力的油气探勘开发新技术之一。目前国际上仅有斯伦贝谢公司于2019年推出 IriSphere随钻前探测电磁波测井仪器，实现最大30 m的前探距离，并在澳大利亚等地成功开展试验，检测出储层与盐层，识别出薄层，规避了钻井作业风险，优化了取心作业1-3。常规随钻电磁波测井仪器由轴向的发射和接收天线组成4-5，早期采用轴向的单发双收天线结构，如贝克休斯公司的双频电磁波测井仪（Dual Propaga-tion Resistivity Logging Tool，DPR）6、哈里伯顿公司的随钻电磁波电阻率测井仪（Electromagnetic Wave Resistivity Logging Tool，EWR）等。随着技术进步，常规随钻电磁波测井仪器逐渐发展成对称补偿的双发双收的天线结构，如斯伦贝谢公司的补偿双电阻率测井仪（Compensated Double Resistivity Logging Tool，CDR），这种结构能够实现井眼补偿功能，抵消探测器和电子线路引起的误差，改善仪器的精度，且测量点正好位于接收器的中点。后来为了获得多种探测深度，阵列结构的随钻电磁波仪器应运而生，如哈里伯顿公司的EWR-Phase4，采用四发双收六天线系结构，可得到8条曲线，实现深、中、浅、极浅4种不同探测深度的测量。分析常规随钻电磁波测井仪器的发展历程，可见其主要特点：最小天线单元均为轴向单发双收结构，测量分量仅有zz分量；尽管通过多源距、多频率阵列组合的天线能够获取不同探测深度，但由于其源距多在1 m以下，频率多为几百千赫兹到几兆赫兹，探测深度有限，不能实现精确实时的前探测7。随钻前探测电磁波测井仪器除了含有轴向线圈外，还配备倾斜或正交线圈，所测到的电压信号中包含方位信息的多个分量信号，并且具有较深的探测深度，因而，仪器对未钻地层界面方位和距离的探测能力有了极大提高8-11。应用随钻前探测电磁波测井仪器可以在井眼穿出储层之前及时采取纠正措施，使井眼轨迹保持在储层中，实现由“被动地质导向”到“主动地质导向”的升级。以斯伦贝谢公司的IriSphere为例，该仪器由1个发射短节、2个或2个以上的接收短节组成7；发射短节包括1个倾斜的发射天线，接收短节由3个相互正交的接收天线组成；根据现场作业需要可以将多个短节灵活组合，组合后其源距可达10 m以上；电磁波信号经过多频发射天线发出，在钻头前方地层中传播，到达多频接收天线，获取钻头前方地层电阻率信息，测量信号与其他模型响应对比，其偏移量能够指示前方界面的存在。随钻前探测测井技术为预测钻头前方未钻地层提供了手段，对于地质导向、风险预判具有重要意义。本文研究典型随钻前探测电磁波测井仪的天线系统特点，综合分析前探距离的有关定义，通过数值模拟探讨其影响因素，为随钻前探测电磁波测井仪器研制提供理论基础。1 随钻前探测电磁波测井方法1.1 随钻前探测电磁波测井仪的天线结构及测井方法随钻前探测电磁波测井仪器具有以下2种常见的天线组合：发射天线由倾斜天线或三轴正交天线组成，接收天线由三轴正交天线组成（见图1）。以仪器轴作为z轴，仪器工具面向上的方向作为x轴，法线为x、y、z这3个方向的三轴正交发射天线，记为Tx、Ty、Tz，法线为x、y、z这3个方向的三轴正交接收天线，记为Rx、Ry、Rz，倾斜发射天线记为T。zyxRyRzRxTyTzTxzyxTRyRzRx图1 随钻前探测电磁波测井仪2种天线排列示意图将仪器坐标系中接收线圈处的磁场记为H，其各个方向分量的矩阵形式可以表示为式（1）12。H=HHHHHHHHHxxxyxzyxyyyzzxzyzz（1）式中，Hmn表示m方向（x、y、z方向）发射线圈在n方向（x、y、z方向）接收线圈上产生的磁场分量，A/m。不同分量的磁场产生的感应电动势可表示为式（2），表征前探测响应的幅度比RP与相位差信号RA定义为式（3）和式（4）13-15。VI N A N A Hmnmn=-iTTTRR（2）2022年测 井 技 术678 RPVVVzzxxyy=-+arctan2（3）RAVVVzzxxyy=-+202lg（4）式中，Vmn为Hmn磁场分量产生的感应电动势，V；i为虚部算子；为角频率，rad*/s；为磁导率，H/m；IT为发射天线的电流，A；NT为发射天线匝数；AT为发射天线面积，m2；NR为接收天线匝数；AR为接收天线面积，m2。三轴正交发射-三轴正交接收天线组合能够直接测量到全部9个磁场分量，倾斜发射-三轴正交接收天线所测信号通常为2个磁场分量信号的叠加，通过信号分解能够得到全部9个磁场分量。其中Hxx、Hyy或者交叉耦合分量信号包含丰富的地质方位信息，更有利于地层界面探测16-18。1.2 Hxx、Hyy、Hzz分量的前探敏感性分析前探测响应信号RP主要由Hxx、Hyy、Hzz这3个磁场分量产生的感应电动势组成，基于三维几何因子理论计算分析这3个分量的三维空间响应12,19，图2为3个分量磁场信号90%贡献的等位面。设发射与接收天线之间源距为L（单位m），Hxx分量90%的贡献在z方向超过 3 L的空间范围，Hyy与Hxx在z方向上的贡献基本一致，而Hzz分量90%的贡献来自于z方向 2 L的空间范围内。对于直井中前探测信号，Hxx和Hyy分量比Hzz分量更具有前探敏感性。（c）Hzz三维几何因子图2 Hxx、Hyy、Hzz磁场分量的三维几何因子（a）Hxx三维几何因子（b）Hyy三维几何因子5 L4 L3 L2 L1 L0-3 L-2 L-1 L-4 L5 L4 L3 L2 L1 L0-3 L-2 L-1 L-4 L-5 L5 L5 L5 L4 L4 L4 L3 L3 L3 L2 L2 L2 L1 L1 L1 L000-3 L-3 L-3 L-2 L-2 L-2 L-1 L-1 L-1 L-4 L-4 L-4 L-5 L-5 L-5 LX/mX/mHxxHxxY/mZ/mZ/m5 L4 L3 L2 L1 L0-3 L-2 L-1 L-4 L5 L4 L3 L2 L1 L0-3 L-2 L-1 L-4 L-5 L5 L5 L5 L4 L4 L4 L3 L3 L3 L2 L2 L2 L1 L1 L1 L000-3 L-3 L-3 L-2 L-2 L-2 L-1 L-1 L-1 L-4 L-4 L-4 L-5 L-5 L-5 LX/mX/mHxxHxxY/mZ/mZ/m5 L4 L3 L2 L1 L0-3 L-2 L-1 L-4 L5 L4 L3 L2 L1 L0-3 L-2 L-1 L-4 L-5 L5 L5 L5 L4 L4 L4 L3 L3 L3 L2 L2 L2 L1 L1 L1 L000-3 L-3 L-3 L-2 L-2 L-2 L-1 L-1 L-1 L-4 L-4 L-4 L-5 L-5 L-5 LX/mX/mHxxHxxY/mZ/mZ/m1.3 前探测电磁波镜像法原理镜像法是一种求解电磁波在界面附近静电场问题的方法20。前探测电磁波测井主要适用于直井中探测前方界面是否存在，假设地层在仪器侧面方向上有规律地平面延伸。以仪器坐标系为参照，z方向是仪器的轴向并且指向钻进方向，x方向在纸面上向右，y方向垂直于纸面向外。*非法定计量单位，1 rad=（180/）（），下同对于图 3（a）所 示的发射-接收天线对，发射和接收都是z方向。仪器所在地层的电导率为2，邻层的电导率为1。如果仪器从一个高电阻率介质逼近一个低电阻率介质，并且考虑极限情况1趋近（理想低电阻率介质，如金属）、以及2趋近0（理想高电阻率介质，如空气），通过镜像法后，发射天线T的镜像源为T，可得到式（5）。对于图 3（b）所示的发射-接收天线对，发射和接收都是x方向，可得到式（6）。田园诗，等：随钻前探测电磁波测井响应影响因素分析第46卷 第6期679 HLDLzz=-[phone]（5）HLDLxx=-[phone]（6）式中，L为仪器发射与接收天线之间的源距，m；D为仪器前端与界面的位移，m。假设地层在仪器侧面方向上是均匀的，则 Hxx=Hyy，带入公式（2）可得-+=+-[phone]VVVD LD Lzzxxyy/（7）在这种情况下，仪器响应是 2 D/L 的函数。在理想的地层电导率下，仪器响应仅与D和L有关。当D远大于L时，前探信号极限趋近于2。2 趋近 01 趋近 RTTLD2 趋近 01 趋近 0RTLDD2 趋近 01 趋近RTTLD2 趋近 01 趋近 0RTLDD镜像法镜像法（a）Hzz的镜像源等效示意图图3 镜像法示意图（b）Hxx的镜像源等效示意图2 基于数值模拟的测井响应影响因素分析2.1 地层电阻率对前探距离的影响2.1.1 电阻率背景值对前探距离的影响随钻前探测电磁波测井技术主要用于测量钻头前方未钻地层的界面信息，建立直井单界面测井模型，令地层界面位于D=0处，地层1位于界面下方（D0）。地层1电阻率记为RT1，地层2电阻率记为RT2，地层电阻率对比度为地层2与地层1电阻率之间的比值，记为S。建立地层模型参数见表1。地层电阻率背景值即地层1电阻率分别为0.01、0.10、1.00 m和10.00 m，地层电阻率对比度均为100。表 1 地层模型参数模型地层电阻率/（m）地层电阻率对比度地层1地层210.011.0010020.1010.0010031.00100.00100410.001 000.00100数值模拟结果见图4。图 4（a）中，当对比度相同时，地层界面两侧的幅度比曲线均随着电阻率背景值的增大而增大。仪器在远离界面时幅度比基本平稳，电阻率越大，幅度比基值越大；从高电阻率地层靠近界面过程中，不同背景值的幅度比曲线先后出现衰减。仪器从高电阻率地层靠近低电阻率地层、或低电阻率地层靠近高电阻率地层时，幅度比曲线会由平稳基值发生偏移，这个偏移量能够指示前方界面的存1050-5-10-15-20-25-30-3520151050-5-10-15-20RA/dB1.00.90.80.70.60.50.40.30.20.1020RT1=0.01 mRT1=0.10 mRT1=1.00 mRT1=10.00 mRT1=0.01 mRT1=0.10 mRT1=1.00 mRT1=10.00 m15D/mRA/dBD/m1050图4 不同电阻率背景值下前探测响应曲线（a）不同电阻率背景值的前探测幅度比（b）归一化后不同电阻率背景值的幅度比2022年测 井 技 术680 在。受仪器本身测量精度的限制，5%以内的信号波动常被认为是误差，因此，令偏移量为5%时D作为前探距离。对高电阻率地层靠近界面时的响应做归一化处理，如图 4（b）所示。在电阻率对比度为100时，前探距离受背景值影响很大，背景值越大，前探距离越大。因而，后续计算均采用低电阻率、高电阻率这2种背景值进行对比。2.1.2 电阻率对比度对前探距离的影响建立的地层模型参数见表2，选取地层电阻率背景值为0.01 m和1.00 m这2种情况，通过数值模拟分别考察这2种地层电阻率对比度对前探距离的影响。将高电阻率地层中的前探测响应做归一化处理，如图5所示。图 5（a）中，在电阻率背景值较小时，随着电阻率对比度的增大，前探距离由3 m增大至9 m。图 5（b）中，在电阻率背景值较大时，前探距离几乎不随电阻率对比度的变化而变化，且前探距离约为13 m。1.00.90.80.70.60.50.40.30.20.1020S=10S=20S=50S=100S=1 000S=10S=20S=50S=100S=1 000151050RA/dB1.00.90.80.70.60.50.40.30.20.102015D/mRA/dBD/m1050图5 不同电阻率对比度下前探测响应曲线（b）归一化电阻率背景值为1.00 m的幅度比（a）归一化电阻率背景值为0.01 m的幅度比表 2 地层模型参数模型地层电阻率/（m）地层电阻率对比度地层1 地层2 10.010.101020.010.202030.010.505040.011.0010050.0110.001 00061.0010.001071.0020.002081.0050.005091.00100.00100101.001 000.001 0002.2 频率对前探距离的影响为考察不同频率对前探测测井响应的影响，选取界面两侧电阻率分别为0.01、1.00 m和1.00、100.00 m这2种模型，界面两侧电阻率对比度均为100。发射频率f为1 kHz1 MHz。将高电阻率地层中的前探测响应做归一化处理（见图6）。在2种地层电阻率模型下皆有以下特征：频率越高，幅度比在界面附近的衰减越快。在界面两侧地层电阻率分别为0.01、1.00 m时，前探距离约为16 m，频率越小，前探距离越大。在界面两侧地层电阻率为1.00、100.00 m时，前探距离为414 m，频率越小，前探距离越大。田园诗，等：随钻前探测电磁波测井响应影响因素分析第46卷 第6期681 1.00.90.80.70.60.50.40.30.20.1020F=1 kHzF=5 kHzF=10 kHzF=100 kHzF=500 kHzF=1 MHzF=1 kHzF=5 kHzF=10 kHzF=100 kHzF=500 kHzF=1 MHz151050RA/dB1.00.90.80.70.60.50.40.30.20.102015D/mRA/dBD/m1050图6 不同频率下前探测响应曲线（a）归一化地层电阻率分别为0.01、1.00 m的幅度比（b）归一化地层电阻率分别为1.00、100.00 m的幅度比2.3 源距对前探距离的影响为考察不同源距对前探测测井响应的影响，选取界面两侧地层电阻率分别为0.01、1.00 m和1.00、100.00 m这2种模型，发射-接收天线间的源距为110 m。将高电阻率地层中的前探测响应归一化（见图7）。在2种地层电阻率模型下皆有以下特征：源距越大，幅度比在界面附近的衰减量越大。在界面两侧地层电阻率为0.01、1.00 m时，前探距离约为35 m，源距越大，前探距离越大。在界面两侧地层电阻率为1.00、100.00 m时，前探距离为1115 m，源距越大，前探距离越大。RA/dB1.00.90.80.70.60.50.40.30.20.1020L=1 mL=3 mL=6 mL=10 m15D/m1050RA/dB1.00.90.80.70.60.50.40.30.20.1020L=1 mL=3 mL=6 mL=10 m15D/m1050图7 不同源距下前探测响应曲线（a）归一化地层电阻率为0.01、1.00 m的幅度（b）归一化地层电阻率1.00、100.00 m的幅度比3 结 论（1）前探测仪器的响应包含Hxx、Hyy、Hzz这3个磁场分量产生的感应电动势，其测量信号RA和RP对于仪器前方的界面非常敏感，在距离界面较远处就能够预判界面的存在。（2）前探测仪器从高电阻率地层靠近低电阻率地层、或低电阻率地层靠近高电阻率地层时，幅度比曲线会由平稳基值发生偏移，这个偏移量能够指示前方界面的存在。令偏移量为5%时D作为前探距离，在特定的电阻率对比度及仪器参数下，前探距离能够达到10 m以上。（3）前探距离受界面两侧地层电阻率对比度及地层电阻率背景值的影响，一般情况下，对比度越大，前探距离越大。前探距离受频率、源距等仪器参数影响，频率越低或源距越长，前探距离越长。（4）前探测电磁波镜像法在理想介质中大大简化了运算量，在其他介质中的计算精度与速度仍有待研究。2022年测 井 技 术682 参考文献：1 SEYDOUXJ,DENICHOU J M,AMIR I,et al.Real-time EM look-ahead:a maturing technology to decrease drilling risk in low inclination wells C/SPWLA 60th Annual Logging Symposium,Woodlands,Texas,2019.2 CONSTABLE M V,ANTONSEN F,STALHEIM S O,et al.Looking ahead of the bit while drilling:from vision to reality C/SPWLA 57th Annual Logging Symposium,Reykjavik,Iceland,2016.3 DERI C P,PETERNELL A,HEBERT J,et al.Reducing uncertainties in casing depth and formation evaluation using look-ahead logging while drilling technology C/SPE Annual Technical Conference and Exhibition,SanAntonio,Texas,2012.4 高杰，娄研，李可赛，等.随钻前视电磁波电阻率测井方法前期理论研究J.测井技术，2018，42（1）：20-24.5 LI K,GAO J,WU C,et al.Study on ultra-deep azimuthal electromagnrtic resistivity LWD tool by influence quantification on azimuthal depth of investigation and real signal J.Pure and Applied Geophysics,2018,175(12):4465-4482.6 CLARK B,LIILING M G,JUNDT J,et al.A dual depth resistivity measurement for FEWD C/SPWLA 29th Annual Logging Symposium,San Antonio,Texas,1988.7 牒勇，朱军，杨善森，等.随钻电磁波电阻率前探测特性分析J.测井技术，2021，45（3）：273-279.8 杨震，刘庆成，岳步江，等.随钻电磁波测井中极化角的形成机理及其影响因素模拟分析J.测井技术，2010，34（3）：210-214.9 BITTAR M,HVEDING F,CLEGG N,et al.Maximizing reservoir contact in the Oseberg field using a new azimuthal deep-reading technology C/SPE Annual Technical Conference and Exhibition,Denver,Colorado,2008.10 CHEMALI R,BITTAR M,HVEDING F,et al.Intergrating images from multiple depths of investigation and quantitative signal inversion in real time for accurate well placement C/International Petroleum Technology Conference,Kuala Lumpur,Malaysia,2008.11 李虎.复杂介质随钻方位电磁波测井数值模拟与应用基础研究D.青岛：中国石油大学（华东），2013.12 王磊.深探测多分量随钻电磁波测井理论与正反演研究D.青岛：中国石油大学（华东），2018.13 LI Q,OMERAGIC D,CHOU L,et al.New directional electromagnetic tool for proactive geosteering and accurate formation evaluation while drilling C/SPWLA 46th Annual Logging Symposium,New Orleans,Louisiana,2005.14 BITTAR M,CHEMALI R,PITCHER J,et al.Real-time proactive optimal well placement using geosignal and deep images C/Offshore Technology Conference,Houston,Texas,2010.15 OMERAGIC D,LI Q,CHOU L,et al.Deep directional electromagnetic measurements for optimal well placement C/SPE Annual Technical Conference and Exhibition,Dallas,Texas,2005.16 LI H,WANG H.Investigation of eccentricity effects and depth of investigation of azimuthal resistivity LWD tools using 3D fifinite difference method J.Journal of Petroleum Science and Engineering,2016,143:211-225.17 HAWKINS D,PHETTONGKAM N,NAKCHAMNAN N.Optimizing well placement in geosteering using an azimuthal resistivity tool in complex thin bed reservoirs in the gulf of Thailand C/SPWLA 56th Annual Logging Symposium,Long Beach,California,2015.18 ZHOU J,LI H,Rabinovich M,et al.Interpretation of azimuthal propagation resistivity measurements:modeling,inversion,application and discussion C/SPWLA 57th Annual Logging Symposium,Reykjavik,Iceland,2016.19 ZHU J,DIE Y,TIAN Y,et al.Deciphering the capabilities of look-ahead methods in LWD C/SPWLA 62nd Annual Logging Symposium Online,2021.20 ZHU J,DIE Y,TIAN Y,et al.Spatial sensitivity study essential to geosteering interpretation C/SPWLA 63rd Annual Logging Symposium,Stavanger,Norway,2022.（修改回稿日期：2022-09-27 编辑 叶春丽）欢迎订阅 欢迎赐稿