从2005年开始,Schlumberger,Baker Hughes和Halliburton公司在电磁波电阻率仪器的基础上陆续推出了自己的随钻方位电磁波电阻率仪器[1-3],并且随着油田勘探开发程度的提高,随钻方位电磁波电阻率仪器在地层评价和地质导向中的应用越来越广泛[4-6]。地层评价所用的幅度电阻率和相位电阻率是通过2个接收天线的电动势幅度比和相位差转换得到的,同时利用定向天线(倾斜或水平天线)测量得到的定向电动势来反演地层界面位置,因此天线测量精度会直接影响到地层评价和地质导向的应用效果。在目前的应用中,大部分只考虑了地层环境或地质因素对测量结果的影响[7],而实际测量过程中在某些地层条件下由于测量精度带来的误差可能要超过井眼环境因素对测量结果造成的影响。例如, 高电阻率地层的随钻电磁波电阻率测量值往往不再做环境影响校正,就是因为在高电阻率地层环境下的测量精度较低,导致视电阻率不能准确反映地层真实的电阻率。笔者采用数值模拟方法分析不同地层电阻率条件下仪器响应信号的特点,研究随钻方位电磁波电阻率仪器在不同地层环境下测量精度对电阻率以及界面距离反演结果的影响,以更好地评估测量资料和利用测量资料。
1 测量精度对电阻率影响模拟与传统的电磁波电阻率仪器相比,中石化胜利石油工程有限公司研发的随钻方位电磁波电阻率仪器(azimuthal multiple resistivity,AMR)增加了一个横向接收天线Rc[8](见图 1),轴向发射天线与横向接收天线彼此垂直,以仪器轴线为z轴,仪器坐标系如图 1所示,AMR可以测量zx或zy分量。
![]() |
图 1 AMR仪器结构示意 Fig.1 Structures of AMR tools |
AMR的工作频率为2 MHz和400 kHz,发射天线T1和T2、T3和T4分别关于轴向接收天线R1和R2对称,可以提供8条对称补偿的电磁波电阻率曲线;发射天线T5和T2关于横向接收天线Rc对称,可以提供2条对称补偿的定向电动势曲线,AMR利用横向天线测得的电动势反演界面距离来进行地质导向。
考虑仪器测量的磁场zx、zy和zz分量的响应特点,在层状地层模型(如图 2所示)中忽略井眼环境, 应用并矢格林函数求出任意方向磁偶极子的磁场z分量,推导出各地层的电磁场分布[9-11],进而计算出接收天线的电动势:
![]() |
(1) |
![]() |
(2) |
其中
![]() |
(3) |
![]() |
(4) |
式中:Ez为z方向的电场强度,V/m;Hz为z方向的磁场强度,A/m;μ为磁导率,H/m;ω为角频率,rad/m;φ为方位角,rad;Mh和Mv分别为水平磁矩和垂直磁矩,A·m2;Jn(·)为n阶Bessel函数; λ为积分变量;kzv和kzh分别为纵向波数和径向波数;σv和σh分别为垂直电导率和水平电导率,S/m。
随钻电磁波电阻率仪器是根据电磁波穿过不同物性地层(见图 2, 图中第m层地层特征参数包括电导率σm、磁导率μm和介电常数εm)时接收线圈感应电动势幅度比和相位的变化,反演得到地层的物性参数[12-15],转换得到电阻率与地层真实电阻率之间的相对误差:
![]() |
图 2 仪器及地层模型的空间关系 Fig.2 Spatial relationship between instrument and formation model |
![]() |
(5) |
![]() |
(6) |
![]() |
(7) |
式中:ΔA为幅度比,dB;Δφ为相位差,(°);UR1和UR2分别为2个接收线圈的感应电动势,V;φR1和φR2分别为2个接收线圈感应电动势的相位角,(°);δ为相对误差;Ra为测量的地层视电阻率,Ω·m;Rt为地层真电阻率,Ω·m。
随钻方位电磁波电阻率仪器幅度比、相位差与电导率的转换关系如图 3所示。由图 3可知,在低电导率地层(高电阻率地层),幅度比和相位差对电导率的分辨率很差,即仪器测得幅度比和相位差的微小误差都会引起电阻率发生巨大变化,因此仪器测量精度对转换电阻率影响较大。
![]() |
图 3 电磁波电阻率转换模板 Fig.3 Resistivity conversion template for electromagnetic wave |
通过转换关系可以模拟分析不同相位差和幅度比测量精度条件下的电阻率误差范围。接收天线相位差误差分别为±0.10°、±0.05°和±0.02°时转换电阻率的误差范围如图 4所示,模拟幅度比精度分别为±0.02 dB、±0.01 dB、±0.005 dB时的转换电阻率误差范围如图 5所示。
![]() |
图 4 相位差精度对电阻率影响模拟 Fig.4 Simulation of effect of phase difference accuracy on resistivity |
对电阻率的相对误差进行分析发现,接收天线幅度比和相位差的测量误差对地层电阻率的影响随着地层电阻率增大而增大。目前,国内外随钻电磁波电阻率仪器幅度比和相位差的测量精度大概为±0.005 dB和±0.02°,在该精度条件下,幅度比精度对电阻率误差的影响要远远大于相位差精度的影响,例如在电阻率100 Ω·m地层条件下,相位差精度为±0.02°时,根据电阻率转换关系,转换电阻率分别为100.88和97.31 Ω·m,误差在3%以内;幅度比精度为±0.005 dB时,转换电阻率分别为110.30和73.41 Ω·m,最大误差达到27%。而地层电阻率在1 000 Ω·m的情况下,幅度电阻率的误差已超过100%,不能有效测量地层的电阻率,因此即使仪器精度达到要求,在200 Ω·m的地层条件下幅度电阻率已经很难测出。
![]() |
图 5 幅度比测量精度对电阻率影响模拟 Fig.5 Simulation of effect of amplitude ratio accuracy on resistivity |
对比模拟结果可知,在目前天线测量精度条件下, 相位差电阻率精度要优于幅度比电阻率精度,因此相位差电率有效测量范围要大于幅度比电阻率的有效测量范围。在地层精细解释过程中,需要对测量结果进行环境影响因素校正,而当地层电阻率比较大时,电阻率的测量误差往往大于环境因素的影响,在该情况下再进行环境因素校正意义不大,在处理过程中需要特别注意。因此要想扩大随钻方位电阻率仪器的有效测量范围和提高测量精度,关键是要提高硬件电路对小信号的分辨精度。
2 界面距离反演结果影响模拟随钻方位电磁波电阻率仪器AMR并不能直接测量仪器到地层界面的距离,而是利用zx和zy方向的电动势对预先设定的地层模型进行反演来预测和判断地层边界,因此定向电动势的测量精度及地层模型的误差都会对反演结果带来影响。
2.1 定向电动势与界面距离响应关系模拟AMR的模型利用轴向发射天线T2和T5发射与接收天线Rc接收交叉耦合的电动势作为定向电动势信号,地层界面电阻率模型为1 Ω·m:50 Ω·m,仪器在高电阻率层中以接近水平姿态逐渐远离界面,其定向电动势信号幅度与界面距离关系如图 6所示。由图 6可看出,在同样地层条件下不考虑定向电动势信号的增益,工作频率越大,定向电动势信号的幅度越大[16-17]。定向信号幅度与界面距离呈近似的对数线性关系,且距离界面越近,对数线性的斜率越大,随着远离界面,对数线性斜率越小,说明随着远离界面,定向电动势信号精度对界面距离的影响越来越大。图 7为模拟层界面两侧地层电阻率对比度分别为1 Ω·m:2 Ω·m,1 Ω·m:5 Ω·m,1 Ω·m:10 Ω·m和1 Ω·m:50 Ω·m时,定向电动势响应与界面距离的关系。由图 7可看出,界面两侧电阻率对比度或电导率差越小,对数线性关系越明显;电阻率对比度越大,定向电动势信号幅度和有效界面探测距离也越大,在远离界面时线性关系斜率越小,测量精度对界面距离的影响也越大。
![]() |
图 6 不同频率下定向电动势信号幅度与界面距离的关系 Fig.6 Relationship between signal amplitude of directional electromotive force and interface distance under different frequencies |
![]() |
图 7 不同界面电阻率对比度下定向信号幅度与界面距离关系模拟 Fig.7 Simulation of relationship between directional signal amplitudes and interface distance with different resistivity contrast |
随钻方位电磁波仪是利用交叉电动势反演层界面距离的,而地质导向应用中影响界面反演距离的主要因素有交叉耦合电动势的测量精度、地层非均质性引起的模型噪声及地层电阻率的不确定性。地层非均质性引起的模型噪声很难定量模拟,本文主要讨论测量精度和电阻率不确定性对反演结果的影响。
将测井数据表示为向量形式,地层参数与仪器响应可以表示为[18-19]:
![]() |
(8) |
![]() |
(9) |
式中:R1和R2为界面两侧地层电阻率,Ω·m,可以通过邻井资料或随钻测井资料得到;h为仪器到地层界面的距离,m;y为随钻方位仪器响应;f为仪器响应函数;*表示转置;Φ为目标函数。
利用反演算法分析测量精度对界面距离的影响程度,反演采用的地层模型为2层地层模型,地层界面深度1 000 m,界面两侧地层电阻率分别为1 Ω·m和20 Ω·m,分别设定仪器横向接收天线小信号分辨率为20 nV和50 nV,利用AMR仪器Rc接收天线工作频率为400 kHz时的响应曲线进行反演,得到测量精度对预测界面距离的影响(见图 8(a),图中标出了真实界面距离10%的误差线)。由图 8(a)可看出,随着界面距离增大,接收天线测量精度对反演结果的影响也越来越大,即使是20 nV的小信号分辨率在界面距离大于2 m以后,反演误差仍然接近10%,在该模型中虽然界面距离在5 m时仍能反演出界面距离,但反演误差非常大。
![]() |
图 8 地层界面反演结果的影响因素模拟分析 Fig.8 Simulation and analysis of affecting factors on formation interface inversion results |
除了仪器测量精度的影响,地层模型中电阻率的不确定性也会给反演结果造成误差。在反演过程中,需建立地层电阻率模型,由于对地层认知程度不同,因此所采用的电阻率与真实地层电阻率之间存在一定偏差。在采用的反演模型中,低电阻率地层的电阻率保持1 Ω·m,高电阻率地层的电阻率分别设为10,20和30 Ω·m,反演结果如图 8(b)所示。由图 8(b)可知:在层界面距离较小时,高电阻率地层的电阻率对界面距离的反演结果影响较小,即无论高电阻率地层的电阻率是10 Ω·m,20 Ω·m还是30 Ω·m,反演界面距离的差别很小; 随着界面距离增大,反演得到的界面距离与真实界面距离(20 Ω·m时反演结果)的偏差越来越大;当高电阻率地层电阻率选取偏大时(30 Ω·m),反演得出的界面距离偏大,而高电阻率地层电阻率选取偏小时(10 Ω·m),反演得出的界面距离偏小。由于在远离界面时,仪器信号的测量精度以及地层电阻率选取的准确性对实际界面距离反演结果都有比较大的影响,因此反演过程中要充分考虑这2个因素。通过分析仪器测量精度及地层模型的不确定性对反演结果的影响程度,可以对反演结果的可靠性有一个合理的评估,从而有效降低地质导向的风险性。
3 结论1) 随钻方位电磁波仪器通过反演得到地层电阻率及界面距离,接收天线测量精度对电阻率及界面距离的确定有较大的影响,通过数值模拟可以定量确定不同地层环境条件下地层视电阻率及反演界面距离的误差范围。
2) 定向电动势幅度与界面距离呈近似对数线性关系,工作频率及界面两侧介质电阻率对比度越大,定向电动势幅度越大,随着界面距离增大,定向电动势测量精度对界面距离反演结果的影响程度增大。
3) 定向电动势信号测量精度、地层非均质性和地层电阻率的不确定性会对反演结果造成不利影响。因此准确进行地质导向的关键前提是如何评估定向电动势信号的测量精度和选取合适的地层模型。
[1] | LI Qiming, OMERAGIC D, CHOU L, et al.New directional electromagnetic tool for proactive geosteering and accurate formation evaluation while drilling[R].SPWLA-2005-UU, 2005. |
[2] | ZHANG Zhiyi, GONGUET C, RAJANI V, et al.Directional LWD resistivity tools and their business impacts[R].SPWLA-2008-FFFF, 2008. |
[3] | SEIFERT D, CHEMALI R, BITTAR M.The link between resistivity contrast and successful geosteering[R].SPWLA-2011-VVV, 2011. |
[4] |
张辛耘, 郭彦军, 王敬农.
随钻测井的昨天、今天和明天[J]. 测井技术, 2006, 30(6): 487–491.
ZHANG Xinyun, GUO Yanjun, WANG Jingnong. The logging while drilling:past, present and future[J]. Well Logging Technology, 2006, 30(6): 487–491. |
[5] |
杨震, 刘庆成, 岳步江, 等.
随钻电磁波测井中极化角的形成机理及其影响因素模拟分析[J]. 测井技术, 2010, 34(3): 10–14.
YANG Zhen, LIU Qingcheng, YUE Bujiang, et al. On mechanism of polarization angle of electromagnetic logging while drilling and its influence factors simulation[J]. Well Logging Technology, 2010, 34(3): 10–14. |
[6] |
李勇华, 杨锦舟, 杨震, 等.
随钻电阻率地层边界响应特征分析及应用[J]. 石油钻探技术, 2016, 44(6): 111–116.
LI Yonghua, YANG Jinzhou, YANG Zhen, et al. The analysis and application of formation interface response characteristic of the resistivity LWD tool[J]. Petroleum Drilling Techniques, 2016, 44(6): 111–116. |
[7] | NARDI G, MARTAKOV S, NIKITENKO M, et al.Evaluation of parameter uncertainty utilizing resolution analysis in reservoir navigation increases the degree of accuracy and confidence inwell-bore placement[R].SPWLA-2010-78038, 2010. |
[8] |
杨震, 杨锦舟, 韩来聚, 等.
随钻方位电磁波界面探测性能分析[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. DOI:10.7623/syxb201607012 |
[9] |
肖加奇, 张庚骥.
水平井和大斜度井中的感应测井响应计算[J]. 地球物理学报, 1995, 38(3): 396–404.
XIAO Jiaqi, ZHANG Gengji. Computation of induction logging response in horizontal and highly-deviated wells[J]. Acta Geophysica Sinica, 1995, 38(3): 396–404. |
[10] |
史晓锋.
水平井中随钻电阻率测量仪定位和预测地层界面的方法[J]. 测井技术, 2006, 30(2): 119–121.
SHI Xiaofeng. Methods for orientation of formation layerings and formation boundary prediction with LWD resistivity logs in HZ wells[J]. Well Logging Technology, 2006, 30(2): 119–121. |
[11] | ANDERSON B, BONNER S, LULING M G, et al. Response of 2-MHz LWD resistivity and wireline induction tools in dipping beds and laminated formation[J]. The Log Analyst, 1992, 33(5): 461–475. |
[12] | HOLDEN A J, THORSEN A K, GRAVEM T, et al.Application and interpretation of multiple advanced LWD measurements in horizontal wells[R].SPWLA-2006-WWW, 2006. |
[13] |
杨震, 杨锦舟, 韩来聚.
随钻方位电磁波电阻率成像模拟及应用[J]. 吉林大学学报(地球科学版), 2013, 43(6): 2035–2043.
YANG Zhen, YANG Jinzhou, HAN Laiju. Numerical simulation and application of azimuthal propagation resistivity imaging while drilling[J]. Journal of Jilin University(Earth Science Edition), 2013, 43(6): 2035–2043. |
[14] |
康俊佐, 邢光龙, 杨善德.
多传播电阻率测井的探测能力与资料处理方法[J]. 石油勘探与开发, 2006, 33(3): 345–350.
KANG Junzuo, XING Guanglong, YANG Shande. Detective ability and data processing method of the MPR logging[J]. Petroleum Exploration and Development, 2006, 33(3): 345–350. |
[15] |
孙向阳, 聂在平, 赵延文, 等.
用矢量有限元方法模拟随钻测井仪在倾斜各向异性地层中的电磁响应[J]. 地球物理学报, 2008, 51(5): 1600–1607.
SUN Xiangyang, NIE Zaiping, ZHAO Yanwen, et al. The electromagnetic modeling of logging-while-drilling tool in tilted anisotropic formations using vector finite element method[J]. Chinese Journal of Geophysics, 2008, 51(5): 1600–1607. |
[16] |
杨锦舟.
随钻方位电磁波仪器界面预测影响因素分析[J]. 测井技术, 2014, 38(1): 39–45.
YANG Jinzhou. Analysis on the affecting factors of prediction interface with azimuthal LWD electromagnetic tool[J]. Well Logging Technology, 2014, 38(1): 39–45. |
[17] |
杨锦舟, 魏宝君, 林楠.
倾斜线圈随钻电磁波电阻率测量仪器基本原理及其在地质导向中的应用[J]. 中国石油大学学报(自然科学版), 2009, 33(1): 44–49.
YANG Jinzhou, WEI Baojun, LIN Nan. Basic theory of electromagnetic wave resistivity measurement while drilling tool with tilted antenas and its application for geosteering[J]. Journal of China University of Petroleum(Edition of Natural Science), 2009, 33(1): 44–49. |
[18] |
魏宝君, 田坤, 张旭, 等.
定向电磁波传播随钻测量基本理论及其在地层界面预测中的应用[J]. 地球物理学报, 2010, 53(10): 2507–2515.
WEI Baojun, TIAN Kun, ZHANG Xu, et al. Physics of directional electromagnetic propagation measurement-while -drilling and its application for forecasting formation boundaries[J]. Chinese Journal of Geophysics, 2010, 53(10): 2507–2515. DOI:10.3969/j.issn.0001-5733.2010.10.024 |
[19] |
陈华, 范宜仁, 邓少贵, 等.
水平井中随钻电阻率实时确定地层界面方法[J]. 吉林大学学报(地球科学版), 2011, 41(5): 1623–1629.
CHEN Hua, FAN Yiren, DENG Shaogui, et al. Methods for real-time determination of formation boundary with LWD resistivity logs in horizontal wells[J]. Journal of Jilin University(Earth Science Edition), 2011, 41(5): 1623–1629. |