水平井随钻电阻率与双侧向响应差异影响因素分析

doi:10.6038/pg20170519

地球物理学进展 ›› 2017, Vol. 32 ›› Issue (5): 1999-2008.doi: 10.6038/pg20170519

水平井随钻电阻率与双侧向响应差异影响因素分析

胡松¹, 郭洪波², 王昌学³, 孔强夫¹

1. 中国石化石油勘探开发研究院, 北京 100083;
2. 中国石油塔里木油田分公司勘探开发研究院, 新疆库尔勒 841000;
3. 中国石油勘探开发研究院, 北京 100083

收稿日期:2017-01-08 修回日期:2017-07-14 出版日期:2017-10-20 发布日期:2017-10-20
作者简介:胡松,男,1985年生,博士(后),主要从事水平井测井解释评价、电法数值模拟研究工作.(E-mail:husong.syky@sinopec.com)
基金资助:
国家科技重大专项“不同缝洞储集体测井精细评价”（2016ZX05014-002-001）、“致密复杂碎屑岩储层测井评价技术”（2016ZX05002-005-001）、“超深层碳酸盐岩缝洞性储层测井精细评价技术”（2017ZX05005005-005）和中国科学院先导A专项（XDA14010204）联合资助.

Influence factors analysis of response characteristics difference between dual laterolog resistivity and LWD resistivity in horizontal well

HU Song¹, GUO Hong-bo², WANG Chang-xue³, KONG Qiang-fu¹

1. Petroleum Exploration & Production Research Institute, SINOPEC, Beijing 100083, China;
2. Research Institute of Petroleum Exploration & Development, Tarim Oil Field Branch Company, Xinjiang Korla 841000, China;
3. Research Institute of Petroleum Exploration & Development, Petrochina, Bejing 100083, China

Received:2017-01-08 Revised:2017-07-14 Online:2017-10-20 Published:2017-10-20

摘要/Abstract

摘要：

在大斜度井和水平井中，不同电阻率测井系列测量的结果并非完全一致，测量结果的差异性分析对于正确认识和评价地层非常重要.依据双侧向测井和随钻电磁波测井原理，采用三维有限元素法模拟双侧向在水平井中的测井响应特征，采用积分方程法模拟随钻电磁波电阻率在水平井中的测井响应特征，通过对比分析值模拟结果，造成大斜度井和水平井随钻电阻率与双侧向电阻率测井响应差异的主要因素有：井斜、各向异性、地层界面以及泥浆侵入.随钻电阻率一般受泥浆侵入影响小，但是它在水平井中容易产生极化现象，有助于地层界面识别，同时受各向异性影响明显；而双侧向电阻率在水平地层中对地层界面不敏感，且受各向异性影响相对较弱，但会受泥浆侵入影响严重；正确认识这些差异，可以有效指导水平井的合理解释.

Abstract:

The measurement results of different resistivity logging series is not completely consistent in high angle and horizontal well due to the different measuring principle, the response characteristics difference analysis is very important for the correct understanding and evaluation of formation.Based on the logging principle of the dual laterolog and electromagnetic wave Logging While Drilling (LWD), adopted the 3-D finite element method to simulate the dual laterolog response characteristics in horizontal wells, and adopted the integral equation method to simulate the LWD resistivity response characteristics in horizontal wells. Through the comparison and analysis, the results showed that the response differences may be caused by the deviation, anisotropy, formation interface and mud intrusion. The LWD resistivity is little influenced by the mud invasion, but it is easy to produce polarization in horizontal wells which will help to formation interface recognition, at the same time it is severely affected by the anisotropic. The dual laterolog resistivity is not sensitive to the formation interface in horizontal well, and the affect by the anisotropic also relatively weak, but the affect by mud invasion is seriously. The correct understanding of these differences, can effectively guide the reasonable explanation of horizontal well.

中图分类号:

P631

胡松, 郭洪波, 王昌学, 孔强夫. 2017. 水平井随钻电阻率与双侧向响应差异影响因素分析. 地球物理学进展, 32(5): 1999-2008. doi: 10.6038/pg20170519.

HU Song, GUO Hong-bo, WANG Chang-xue, KONG Qiang-fu. 2017. Influence factors analysis of response characteristics difference between dual laterolog resistivity and LWD resistivity in horizontal well. Progress in Geophysics. 32(5): 1999-2008. doi: 10.6038/pg20170519.

参考文献

Anderson B, Bonner S, Luling M G, et al. 1992.Response of 2-MHz LWD resistivity and wireline induction tools in dipping beds and laminated formations[J]. The Log Analyst, 33(5):461-475.
Bilodeau B, De G, Wild T, et al. 2002.Integrating formation evaluation into earth modeling and 3-D petrophysics[C].//SPWLA 43rd Annual Logging Symposium. Oiso, Japan:Society of Petrophysicists and Well-Log Analysts.
Chen G B, Wang H N, Yao JJ, et al. 2009.Modeling of electromagnetic responses of 3-D electrical anomalous body in a layered anisotropic earth using integral equations[J]. Chinese J.Geophys. (in Chinese), 52(8):2174-2181, doi:10.3969/j.issn.0001-5733.2009.08.028.
Cozzolino K, da Conceição da Silva J. 2007. Synthetic focusing and simulation of dual laterolog tool in axisymmetric subsurface models[J]. Journal of Applied Geophysics, 61(2):102-110.
Davydycheva S, Druskin V, Habashy T. 2003. An efficient finite-difference scheme for electromagnetic logging in 3D anisotropic inhomogeneous media[J]. Geophysics, 68(5):1525-1536.
Gao J, Xie R H. 2000. 3D numerical forward modeling and fast correction of dual-laterolog for high angle deviated wells[J]. Petroleum Exploration and Development (in Chinese), 27(2):60-71.
Hagiwara T. 1994. Response of 2 MHz resistivity devices in thinly laminated formation(anisotropic resistivity and EM Log interpretation)[C]. SPE.
Hue Y K,Teixeira F L,Martin L S, et al. 2005. Three-Dimensional simulation of eccentric LWD tool response in boreholes through dipping formations[J]. IEEE Transactions on Geoscience and Remote Sensing, 43(2):257-268.
Li S J, Xiao C W, Wang H M, et al. 1996.Mathematical model of dual laterolog response to fracture and quantitativeinterpretation of fracture porosity[J]. ActaGeophysicaSinica(in Chinese), 39(6):845-852.
Liu D J, Ma Z H, Zhang Y Y, et al. 2013.Numerical simulation of the logging-while-drilling for fracture-vuggy reservoirs using adaptive FEM[J].Progress in Geophysics (in Chinese),28(2):928-936, doi:10.6038/pg20130246.
Nam M J, Pardo D, Torres-Verdín C. 2008.Simulation of DC dual-laterolog measurements in complex formations:A Fourier-series approach with nonorthogonal coordinates and self-adapting finite elements[J]. Geophysics, 74(1):E31-E43.
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[J]. Geophysics, 75(6):F143-F149.
Pardo D, Demkowicz L, Torres-Verdín C, et al. 2007. A self-adaptive goal-oriented hp-finite element method with electromagnetic applications. Part Ⅱ:Electrodynamics[J]. Comput. Methods Appl. Mech. Eng., 196(37-40):3585-3597.
Rodíguez-Rozas A N,Pardo D. 2016. A priori Fourier analysis for 2.5D finite elements simulations of logging-while-drilling (LWD) resistivity measurements[J]Procedia Computer Science, 80:782-791.
Shen J S. 2003. Modeling of 3-D electromagnetic responses in frequency domain by using staggered grid finite diferencemethod[J].Chinese J.Geophys.(in Chinese), 46(2):280-288, doi:10.3321/j.issn:0001-5733.2003.02.024.
Šolín P, ?erveny J, Dole?el I. 2008.Arbitrary-level hanging nodes and automatic adaptivity in the hp-FEM[J]. Math.Comput.Simul.,77(1):117-132.
Spitzer K. 1995.A 3-D finite-difference algorithm for DC resistivity modelling using conjugate gradient methods[J]. Geophysical Journal International, 123(3):903-914.
Sun X Y, Nie Z P, Zhao Y W, et al. 2008.The electromagnetic modeling of logging-while-drilling tool in tilted anisotropic formations using vector finite element method[J]. Chinese J.Geophys. (in Chinese), 51(5):1600-1607.
Tan M J, Gao J, Zou Y L, et al. 2012. Environment correction method of dual laterolog in directional well[J]. Chinese J.Geophys. (in Chinese), 55(4):1422-1432, doi:10.6038/j.issn.0001-5733.2012.04.038.
Tan M J, Zhang G J, Yun H Y, et al. 2007. 3-D numerical mode-matching (NMM) method for resistivity logging responses in nonsymmetricconditions[J]. Chinese J.Geophys. (in Chinese), 50(3):939-945, doi:10.3321/j.issn:0001-5733.2007.03.037.
Tang J T, Ren Z Y, Hua X R. 2007. Theoretical analysis of geo-electromagnetic modeling on Coulomb gauged potentials by adaptive finite element method[J]. Chinese J.Geophys. (in Chinese), 50(5):1584-1594, doi:10.3321/j.issn:0001-5733.2007.05.036.
Wang C X, Zhou C C, Chu Z T, et al. 2006. Modeling of electromagnetic responses in frequencydomain to electrical anisotropic formations[J]. Chinese J.Geophys. (in Chinese), 49(6):1873-1883, doi:10.3321/j.issn:0001-5733.2006.06.037.
Wang T, Fang Sheng. 2001. 3-D electromagnetic anisotropy modeling using finite differences[J]. Geophysics, 66(5):1386-1398.
Wang T, Hohmann G W. 1993. A finite-difference, time-domain solution for three-dimensional electromagnetic modeling[J]. Geophysics, 58(6):797-809.
Wang W, Wu X P. 2010. Rapid finite element resistivity forward modeling for 3D arbitrary anisotropic structures[J]. Progress in Geophysics (in Chinese),25(4):1365-1371, doi:10.3969/j.issn.1004-2903.2010.04.027.
Wei B J, Wang S S, Ou Y F, et al. 2011. Simulating the response of induction logging-while-drilling tools by vector eigenfunction expansion formulae for dyadic Green's functions[J]. Chinese J.Geophys. (in Chinese), 54(5):1391-1401, doi:10.3969/j.issn.0001-5733.2011.05.029.
Wei B J, Wang T T, Wang Y. 2009a. Computing the response of multi-component induction logging in layered anisotropic formation by the recursive matrix method for magnetic-current-source dyadic Green's function[J]. Chinese J.Geophys. (in Chinese), 52(11):2920-2928, doi:10.3969/j.issn.0001-5733.2009.11.028.
Wei B J, Wang Y, Wang T T. 2009b. The basic theory of electromagnetic wave resistivity instrument and its application to measurement of while-drilling[J]. Progress in Geophysics (in Chinese),24(2):774-781, doi:10.3969/j.issn.1004-2903.2009.02.057.
Wu X P, Wang T T. 2003. A 3-D finite-element resistivity forward modeling using conjugate gradient algorithm[J]. Chinese J.Geophys. (in Chinese), 46(3):428-432, doi:10.3321/j.issn:0001-5733.2003.03.023.
Xiao J Q, Zhang G J. 1996. Computation of dual laterolog responses in highly deviated and horizontal wells with 3-D finite element method[J]. Journal of the University of Petroleum,China (in Chinese), 20(1):24-28.
Yang Z, Fan Y R, Wen Y, et al. 2009.Three-dimensional numerical simulation of electromagnetic logging-while-drilling tool response in frequency domain[J]. Progress in Geophysics (in Chinese),24(5):1833-1838, doi:10.3969/j.issn.1004-2903.2009.05.038.
Zhang X, Wei B J, Liu K, et al. 2009. Computing electromagnetic wave resistivity MWD tool's response using accelerated iteration algorithm for integral equations[J]. Chinese J.Geophys. (in Chinese), 52(9):2394-2401, doi:10.3969/j.issn.0001-5733.2009.09.026.
Zhou Q. 2008. Log interpretation in high-deviation wells through user-friendly tool-response processing[C].//49th Annual Logging Symposium.Austin, Texas:Society of Petrophysicists and Well-Log Analysts.
陈桂波, 汪宏年, 姚敬金,等. 2009. 用积分方程法模拟各向异性地层中三维电性异常体的电磁响应[J]. 地球物理学报, 52(8):2174-2181, doi:10.3969/j.issn.0001-5733.2009.08.028.
楚泽涵,高杰, 黄隆基,等. 2007. 地球物理测井方法与原理[M]. 北京:石油工业出版社.
邓少贵. 2006. 大斜度井非均匀地层电测井响应研究[博士论文]. 青岛:中国石油大学(华东).
高杰, 谢然红. 2000. 大斜度井侧向测井三维正演数值模拟及曲线快速校正方法研究[J]. 石油勘探与开发, 27(2):60-71.
高杰, 刘传奇, 万金彬. 2012. 裂缝性储层双侧向测井响应临界角影响因素分析[J]. 测井技术, 36(5):456-459.
李善军, 肖承文, 汪涵明,等. 1996. 裂缝的双侧向测井响应的数学模型及裂缝孔隙度的定量解释[J]. 地球物理学报, 39(6):845-852.
李善军. 1998. 用三维有限元素法反演单一灰岩裂缝性储层的裂缝孔隙度和裂缝倾角[J]. 测井技术, 22(6):412-415.
刘得军, 马中华, 张颖颖,等. 2013. 缝洞型储层随钻测井自适应有限元法数值模拟[J]. 地球物理学进展, 28(2):928-936, doi:10.6038/pg20130246.
刘迪仁, 万文春, 夏培,等. 2012a. 不同井眼条件下的水平井双侧向测井响应[J]. 物探与化探, 36(3):422-425.
刘迪仁, 闫林辉, 夏培,等. 2012b. 水平井中砂泥岩互层的双侧向测井响应研究[J]. 西南石油大学学报(自然科学版), 34(3):78-82.
刘之的. 2006. 随钻测井响应反演方法及应用研究[博士论文]. 成都:西南石油大学.
沈金松. 2003. 用交错网格有限差分法计算三维频率域电磁响应[J]. 地球物理学报, 46(2):280-288, doi:10.3321/j.issn:0001-5733.2003.02.024.
孙向阳, 聂在平, 赵延文,等. 2008. 用矢量有限元方法模拟随钻测井仪在倾斜各向异性地层中的电磁响应[J]. 地球物理学报, 51(5):1600-1607.
谭茂金, 高杰, 邹友龙,等. 2012. 盐水泥浆条件下定向井双侧向测井环境校正方法研究[J]. 地球物理学报, 55(4):1422-1432, doi:10.6038/j.issn.0001-5733.2012.04.038.
谭茂金, 张庚骥, 运华云,等. 2007. 非轴对称条件下用三维模式匹配法计算电阻率测井响应[J]. 地球物理学报, 50(3):939-945, doi:10.3321/j.issn:0001-5733.2007.03.037.
谭永基, 于泳. 1997. 水平井和大斜度井的双侧向测井数值模拟[J]. 复旦学报(自然科学版), 36(6):633-638.
汤井田, 任政勇, 化希瑞. 2007. Coulomb规范下地电磁场的自适应有限元模拟的理论分析[J]. 地球物理学报, 50(5):1584-1594, doi:10.3321/j.issn:0001-5733.2007.05.036.
王昌学, 周灿灿, 储昭坦,等. 2006. 电性各向异性地层频率域电磁响应模拟[J]. 地球物理学报, 49(6):1873-1883, doi:10.3321/j.issn:0001-5733.2006.06.037.
汪涵明, 张庚骥, 李善军,等. 1995. 单一倾斜裂缝的双侧向测井响应[J]. 石油大学学报(自然科学版), 19(6):21-24.
王威, 吴小平. 2010. 电阻率任意各向异性三维有限元快速正演[J]. 地球物理学进展, 25(4):1365-1371, doi:10.3969/j.issn.1004-2903.2010.04.027.
魏宝君, 王莎莎, 欧永峰. 2011. 用并矢Green函数的矢量本征函数展开式模拟随钻感应测井仪器的响应[J]. 地球物理学报, 54(5):1391-1401, doi:10.3969/j.issn.0001-5733.2011.05.029.
魏宝君, 王甜甜, 王颖. 2009a. 用磁流源并矢Green函数的递推矩阵方法计算层状各向异性地层中多分量感应测井响应[J]. 地球物理学报, 52(11):2920-2928, doi:10.3969/j.issn.0001-5733.2009.11.028.
魏宝君, 王颖, 王甜甜. 2009b. 电磁波电阻率仪器的基本理论及其在随钻测量中的应用[J]. 地球物理学进展, 24(2):774-781, doi:10.3969/j.issn.1004-2903.2009.02.057.
吴小平, 汪彤彤. 2003. 利用共轭梯度算法的电阻率三维有限元正演[J]. 地球物理学报, 46(3):428-432, doi:10.3321/j.issn:0001-5733.2003.03.023.
肖加奇, 张庚骥. 1996. 水平井和大斜度井中双侧向测井响应的正演[J]. 石油大学学报(自然科学版), 20(1):24-28.
杨震, 范宜仁, 文艺,等. 2009. 三维频率域随钻电磁波测井数值模拟[J]. 地球物理学进展, 24(5):1833-1838, doi:10.3969/j.issn.1004-2903.2009.05.038.
雍世和, 张超谟. 2007. 测井数据处理与综合解释[M]. 东营:中国石油大学出版社.
张旭, 魏宝君, 刘坤,等. 2009. 利用积分方程的加速迭代算法计算随钻电磁波电阻率测量仪器的响应[J]. 地球物理学报, 52(9):2394-2401, doi:10.3969/j.issn.0001-5733.2009.09.026.

[1]	屈念念,李家斌,张西君,刘俊. 基于重磁资料研究贵州省深部构造特征[J]. 地球物理学进展, 2019, 34(5): 1785-1793.
[2]	安鹏,曹丹平,赵宝银,杨晓利,张明. 基于LSTM循环神经网络的储层物性参数预测方法研究[J]. 地球物理学进展, 2019, 34(5): 1849-1858.
[3]	陈稳,薛国强,陈卫营,周楠楠,卢云飞. SOTEM多分量激电响应特性分析[J]. 地球物理学进展, 2019, 34(5): 1859-1865.
[4]	豆辉,徐逸鹤. 黏性逆时偏移成像研究进展[J]. 地球物理学进展, 2019, 34(5): 1866-1877.
[5]	毛庆辉,王鹏,曾隽. 水力压裂微地震事件定位方法综述[J]. 地球物理学进展, 2019, 34(5): 1878-1886.
[6]	尚新民,王延光,崔庆辉,金昌昆,赵翠霞,滕厚华,赵胜天. 准南缘山前带精细近地表建模技术研究[J]. 地球物理学进展, 2019, 34(5): 1887-1892.
[7]	石战战,夏艳晴,王元君,周怀来,庞溯,池跃龙. 基于L₁-L₁范数稀疏表示的共偏移距道集地震数据重建方法[J]. 地球物理学进展, 2019, 34(5): 1893-1899.
[8]	王迪,陈敏,王健伟,盛志超,严曙梅,马美媛,李东哲,郑欣. 基于T₂和毛管压力分类转换的渗透率计算新方法[J]. 地球物理学进展, 2019, 34(5): 1900-1909.
[9]	王海波,张伟,张宏,肖虎,高云路,顾剑华,唐美斌,韩峰. 高精度可控震源在深反射地震采集中的应用[J]. 地球物理学进展, 2019, 34(5): 1910-1916.
[10]	王静, 张军华, 王延光, 杨勇, 冯德永, 武刚, 黄德峰. 特征值相干理论诠释及效果比较[J]. 地球物理学进展, 2019, 34(5): 1917-1923.
[11]	王维红,包培楠,陈国飞,林弘喆. 基于最小平方的横波速度拟合及应用[J]. 地球物理学进展, 2019, 34(5): 1924-1929.
[12]	王亚,王孟华,田建章,王成泉,彭玲丽,冯小英,周兴海,张浩. 一种适用于深潜山及内幕的有效储层地震预测方法——以廊固凹陷杨税务潜山为例[J]. 地球物理学进展, 2019, 34(5): 1930-1937.
[13]	吴曲波,潘自强,陈金勇,Mubarak,Hani. 利用瞬态瑞雷面波法探测浅层玄武岩三维分布——以沙特Sabkhah Ad Dumathah地区为例[J]. 地球物理学进展, 2019, 34(5): 1938-1944.
[14]	严小丽,康慧敏,王光杰,何国丽,李红领. AMT方法在鳌山卫花岗岩地区深部地热构造勘探中的应用[J]. 地球物理学进展, 2019, 34(5): 1945-1953.
[15]	尹奇峰,潘冬明,郭全仕,何情. 基于快速模拟退火算法的井中微地震事件定位反演[J]. 地球物理学进展, 2019, 34(5): 1954-1961.

水平井随钻电阻率与双侧向响应差异影响因素分析

Influence factors analysis of response characteristics difference between dual laterolog resistivity and LWD resistivity in horizontal well

PDF (PC)

可视化

摘要/Abstract

引用本文

使用本文

参考文献

相关文章 15

编辑推荐

多维度评价

本文评价