CN111058834B - 基于瞬变多分量感应测井的各向异性地层倾角确定方法 - Google Patents
基于瞬变多分量感应测井的各向异性地层倾角确定方法 Download PDFInfo
- Publication number
- CN111058834B CN111058834B CN201911237750.XA CN201911237750A CN111058834B CN 111058834 B CN111058834 B CN 111058834B CN 201911237750 A CN201911237750 A CN 201911237750A CN 111058834 B CN111058834 B CN 111058834B
- Authority
- CN
- China
- Prior art keywords
- stratum
- transient
- dip angle
- induction logging
- anisotropic
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 230000006698 induction Effects 0.000 title claims abstract description 35
- 230000001052 transient effect Effects 0.000 title claims abstract description 30
- 238000000034 method Methods 0.000 title claims abstract description 26
- 230000005284 excitation Effects 0.000 claims abstract description 12
- 230000008859 change Effects 0.000 claims abstract description 11
- 238000000605 extraction Methods 0.000 claims abstract description 7
- 238000005259 measurement Methods 0.000 claims abstract description 7
- 238000011426 transformation method Methods 0.000 claims abstract description 4
- 230000015572 biosynthetic process Effects 0.000 claims description 39
- 230000000007 visual effect Effects 0.000 claims description 9
- 229920000535 Tan II Polymers 0.000 claims description 3
- 230000005540 biological transmission Effects 0.000 claims description 3
- 230000035699 permeability Effects 0.000 claims description 3
- 239000011159 matrix material Substances 0.000 claims description 2
- 230000036962 time dependent Effects 0.000 claims description 2
- 238000011156 evaluation Methods 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 8
- 230000008569 process Effects 0.000 description 4
- 230000007704 transition Effects 0.000 description 4
- 238000004422 calculation algorithm Methods 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000006978 adaptation Effects 0.000 description 1
- 238000007792 addition Methods 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000009792 diffusion process Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Images
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH DRILLING; MINING
- E21B—EARTH DRILLING, e.g. DEEP DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B47/00—Survey of boreholes or wells
- E21B47/02—Determining slope or direction
- E21B47/026—Determining slope or direction of penetrated ground layers
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
- G01V3/18—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation specially adapted for well-logging
- G01V3/30—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation specially adapted for well-logging operating with electromagnetic waves
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
- G01V3/38—Processing data, e.g. for analysis, for interpretation, for correction
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/30—Assessment of water resources
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Environmental & Geological Engineering (AREA)
- Geophysics (AREA)
- Remote Sensing (AREA)
- General Physics & Mathematics (AREA)
- Mining & Mineral Resources (AREA)
- Fluid Mechanics (AREA)
- Geochemistry & Mineralogy (AREA)
- Electromagnetism (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开了一种基于瞬变多分量感应测井的各向异性地层倾角确定方法,所述方法包括:步骤1:建立横向各向同性地层模型,设计多分量的发射接收天线组合模式;步骤2:选择下阶跃信号源作为瞬变多分量电磁波测井的激励源,测量激励源关断后地层中的纯二次场;步骤3:采用逆拉普拉斯变换法将频率域多分量感应测井张量转换到时间域,得到瞬变多分量感应测井的感应电动势;步骤4:构建瞬变多分量感应测井测量信号与地层倾角之间的代数关系;步骤5:根据时间域中各向异性地层倾角的变化特征,确定地层倾角提取方法。本发明能够解决目前各向异性地层倾角反演的繁琐耗时和多解性问题,可为各向异性地层地层评价和地质导向快速准确地提供倾角信息。
Description
技术领域
本发明涉及油气勘探开发领域,具体涉及一种基于瞬变多分量感应测井的各向异性地层倾角确定方法。
背景技术
电性各向异性一直是困扰砂泥岩薄互层、页岩储层等非常规油气藏勘探与开发的重要难题,受各向异性的影响,地层电阻率的空间分布与方向有关,具体表现为地层层理方向电阻率远小于垂向电阻率。传统感应测井只能测量直井中的水平方向电阻率,利用阿尔奇公式评价的含水饱和度过高,容易导致有效储层的漏失。因此各大石油公司陆续推出了多分量感应测井仪器,与传统感应测井不同,多分量感应测井曲线无法提供电阻率、边界等信息的直观显示,一般需要反演实现地层参数的提取,具体为根据邻井或其他测井资料提供的先验信息建立初始地层模型,通过不断对比多分量感应测井实测曲线与正演结果,实现反演模型与实际地质情况的逐步逼近。反演过程繁琐,计算量大,且经常存在多解性问题,采用多频聚焦的方法处理zz、xx和yy分量提取地层各向异性,借助多分量感应测井同轴、共面和交叉分量的组合可同时实现各向异性地层倾角和各向异性系数的线性表达,不需要通过反演。但地层参数提取方法有很大的应用局限,各向异性系数较小、电导率较大、仪器工作频率较高的情况都将导致算法的应用效果变差。瞬变电磁法,是通过发射线圈激励电流,在地层中产生一次场,某一时刻激励源突然关断,通过测量由此产生的涡流二次场的扩散和衰减过程,提取地质体的电性参数,在近地表地球物理勘探。与当前的直流源和时谐源电测井方法相比,瞬变电磁波测井具有探测距离远;只测量纯二次场信息,不受一次场的干扰;仪器尺寸较小等优点,具有广阔的应用前景,基于此,本发明提供一种瞬变多分量感应测井的各向异性地层倾角确定方法,通过研究同轴、共面、交叉分量在时域内不同的衰减特性,建立瞬变多分量电磁波测井地层倾角的代数提取方法,可以避免复杂的反演过程,对各向异性地层倾角的确定有重要意义。
发明内容
针对现有技术需要通过复杂反演来提取地层倾角的问题,本发明提出了一种基于瞬变多分量感应测井的各向异性地层倾角确定方法,设计合理,应用简单,具有良好的前景。
基于瞬变多分量感应测井的各向异性地层倾角确定方法,所述方法包括以下步骤:
步骤1:建立各向异性地层模型,地层坐标体系为(x,y,z),设计多分量的发射接收天线组合模式;
步骤2:选择下阶跃信号源作为瞬变多分量电磁波测井的激励源,测量激励源关断后地层中的纯二次场;
步骤3:采用逆拉普拉斯变换法将频率域多分量感应测井磁场张量转换到时间域,进一步得到瞬变多分量感应测井的感应电动势;
步骤4:构建瞬变多分量感应测井测量信号与地层倾角之间的代数关系;
步骤5:根据时间域中各向异性地层视倾角的变化特征,确定地层真实倾角提取方法。
进一步地,所述步骤1中,建立横向各向同性的地层模型,采用多分量发射接收模式,发射线圈是三个相互正交的线圈,接收线圈采用同样相互正交的布置方式,一次发射能够同时测得xx、yy、zz、xy、xz、yx、yz、zx和zy 9个分量。
进一步地,所述步骤3具体包括以下步骤:
式中,n为线圈匝数,S为线圈面积,μ为地层磁导率。
进一步地,所述步骤4具体包括以下步骤:
步骤4.1:确定测量信号与地层倾角之间的表达式;
式中,α为地层真实倾角,Ψ为地层倾角误差,V'表示仪器坐标系(x’,y’,z’)下的感应电动势,仪器坐标系下的感应电动势V'可由地层坐标系下的感应电动势V通过旋转矩阵得到;
式中αapp为视地层倾角。
进一步地,在步骤4中,地层坐标系下的感应电动势V与仪器坐标系下的感应电动势V’存在如下关系:
由此可以得到:
进一步地,在步骤5中,视地层倾角αapp随时间的变化大致分为三个阶段,早期,过渡期和晚期,其中,早期和晚期的视地层倾角αapp等于地层真实倾角α;过渡期的视地层倾角αapp偏离真实值较大。
有益效果:本发明通过将瞬变多分量感应测井多分量测量信号的简单代数组合即可提取各向异性地层倾角,设计合理,原理简单,计算速度快,能够解决目前各向异性地层倾角反演的繁琐耗时和多解性问题,可为各向异性地层地层评价和地质导向快速准确地提供倾角信息具有广阔的应用前景。
附图说明
图1为本发明基于瞬变多分量感应测井的各向异性地层倾角确定方法的天线结构示意图;
图2为本发明采用的下阶跃脉冲激励源示意图;
图3为本发明具体实施例的视地层倾角示意图;
图4为本发明具体实施例的地层真实倾角为30°时视地层倾角随垂向电阻率的变化示意图;
图5为本发明具体实施例的地层真实倾角为60°时视地层倾角随垂向电阻率的变化示意图;
图6为本发明具体实施例的地层真实倾角为30°时视地层倾角随横向电阻率的变化示意图;
图7为本发明具体实施例的地层真实倾角为60°时视地层倾角随横向电阻率的变化示意图;
图8为本发明具体实施例的地层真实倾角为60°时视地层倾角随源距的变化示意图。
图中,1、下阶跃电流,2、电流、3、上阶跃电流。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅用以解释本发明,并不用于限定本发明,即所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。
下面结合附图以及具体实施方式对本发明作进一步详细说明:
步骤1:建立各向异性地层模型,设计多分量的发射接收天线组合模式;
图1为本发明提供的基于瞬变多分量感应测井的各向异性地层倾角确定方法的天线结构示意图,采用三发三收的多分量线圈系结构,发射线圈T1、T2和T3两两正交布置,接受线圈R1、R2和R3布置方式与发射线圈相同,一次发射能够同时测得介质的xx、yy、zz、xy、xz、yx、yz、zx和zy 9个分量。
图2所示为本发明采用的下阶跃脉冲激励源示意图,下阶跃电流1与上阶跃电流3的关系,向地层中发射电流2并激励足够长的时间,之后关断电流以测量地层中的纯二次场,
步骤3:采用逆拉普拉斯变换法,将频率域多分量感应测井张量转换到时间域,得到瞬变多分量感应测井的感应电动势;
式中,n为线圈匝数,S为线圈面积,μ为磁导率,进而推导出时域的感应电动势Vzz(t)、Vzx(t)、Vxz(t)和Vxx(t)的表达式为:
步骤4:构建瞬变多分量电磁波测井探测信号与视地层倾角之间的代数关系。
步骤4.1:仪器坐标系下的感应电动势V'可由地层坐标系下的感应电动势V存在如下关系:
V'zz(t)-V'xx(t)=(cos2α-sin2α)[Vzz(t)-Vxx(t)]+2sinαcosα[Vzx(t)+Vxz(t)]
V'xz(t)+V'zx(t)=(cos2α-sin2α)[Vzx(t)+Vxz(t)]+2sinαcosα[Vxx(t)-Vzz(t)]
V'zz(t)+V'xx(t)=Vzz(t)+Vxx(t)
由此可以得到:
步骤5:根据时间域中各向异性地层视倾角的变化特征,确定地层真实倾角提取方法。
不同的地层倾角情况下,构建一个垂向电阻率Rv=8Ω·m,水平电阻率Rh=2Ω·m,源距为40in的均匀各向异性地层模型,得到的视地层倾角值如图3所示。视地层倾角值在早期反应真实的地层倾角,在过渡期时偏离真实值,在晚期时重新接近真实值;地层倾角越大,过渡期来的越晚且计算所耗时间更短,在水平井中,过渡期非常小并且视地层倾角在整个时段几乎都反应了地层倾角的真实值,说明该算法在大角度/水平井中更为适用。
不同地层电阻率条件下,在一个地层倾角为30°,源距为40in的均匀各向异性地层模型中,水平电阻率不变,垂向电阻率变化得到的视地层倾角示意图如图4所示,垂向电阻率不变,水平电阻率变化得到的视地层倾角示意图如图6所示;在一个地层倾角为60°,源距为40in的均匀各向异性地层模型中,水平电阻率不变,垂向电阻率变化得到的视地层倾角示意图如图5所示,垂向电阻率不变,水平电阻率变化得到的视地层倾角示意图如图7所示。结果表明垂向电阻率对视地层倾角值影响较小而水平电阻率影响较大,随水平电阻率增加,过渡期出现的时间变早。
不同源距情况下,构建一个垂向电阻率Rv=8Ω·m,水平电阻率Rh=2Ω·m,地层倾角为60°的均匀各向异性地层模型,得到的视地层倾角值如图8所示。可以看出源距越长,过渡期出现得越早,因此想要尽早确定倾角需要更长的源距,而如果在晚期确定倾角则可以利用短源距实现。
当然,上述说明并非是对本发明的限制,本发明也并不仅限于上述举例,本技术领域的技术人员在本发明的实质范围内所做出的变化、改型、添加或替换,也应属于本发明的保护范围。
Claims (2)
1.基于瞬变多分量感应测井的各向异性地层倾角确定方法,其特征在于:所述方法包括以下步骤:
步骤1:建立各向异性地层模型,地层坐标体系为(x,y,z),设计多分量的发射接收天线组合模式;
步骤2:选择下阶跃信号源作为瞬变多分量电磁波测井的激励源,测量激励源关断后地层中的纯二次场;
步骤3:采用逆拉普拉斯变换法将频率域多分量感应测井磁场张量转换到时间域,进一步得到瞬变多分量感应测井的感应电动势;
步骤4:构建瞬变多分量感应测井测量信号与地层倾角之间的代数关系;
步骤5:根据时间域中各向异性地层视倾角的变化特征,确定地层真实倾角提取方法;
所述步骤1中,建立横向各向同性的地层模型,采用多分量发射接收模式,发射线圈是三个相互正交的线圈,接收线圈采用同样相互正交的布置方式,一次发射能够同时测得xx、yy、zz、xy、xz、yx、yz、zx和zy 9个分量;
所述步骤3具体包括以下步骤:
式中,n为线圈匝数,S为线圈面积,μ为地层磁导率;
所述步骤4具体包括以下步骤:
步骤4.1:确定测量信号与地层倾角之间的表达式;
式中,α为地层真实倾角,Ψ为地层倾角误差,V′表示仪器坐标系(x’,y’,z’)下的感应电动势,仪器坐标系下的感应电动势V′由地层坐标系下的感应电动势V通过旋转矩阵得到;
式中αapp为视地层倾角;
在步骤4中,V与V’存在如下关系:
V′zz(t)-V′xx(t)=(cos2α-sin2α)[Vzz(t)-Vxx(t)]+2sinαcosα[Vzx(t)+Vxz(t)]
V′xz(t)+V′zx(t)=(cos2α-sin2α)[Vzx(t)+Vxz(t)]+2sinαcosα[Vxx(t)-Vzz(t)] (6)
V′zz(t)+V′xx(t)=Vzz(t)+Vxx(t)
由此得到:
式中,V所代表的为地层坐标系下的感应电动势,V’所代表的为仪器坐标系下的感应电动势。
2.根据权利要求1所述的基于瞬变多分量感应测井的各向异性地层倾角确定方法,其特征在于,在步骤5中,视地层倾角αapp随时间的变化分为三个阶段,早期,过渡期和晚期,其中,早期和晚期的视地层倾角αapp等于地层真实倾角α;过渡期的视地层倾角αapp偏离真实值较大。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911237750.XA CN111058834B (zh) | 2019-12-06 | 2019-12-06 | 基于瞬变多分量感应测井的各向异性地层倾角确定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911237750.XA CN111058834B (zh) | 2019-12-06 | 2019-12-06 | 基于瞬变多分量感应测井的各向异性地层倾角确定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111058834A CN111058834A (zh) | 2020-04-24 |
CN111058834B true CN111058834B (zh) | 2023-05-02 |
Family
ID=70300144
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911237750.XA Active CN111058834B (zh) | 2019-12-06 | 2019-12-06 | 基于瞬变多分量感应测井的各向异性地层倾角确定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111058834B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111794738B (zh) * | 2020-07-07 | 2023-07-25 | 中国石油天然气集团有限公司 | 一种基于阵列感应与各向异性联合反演地层倾角的方法 |
CN112814668B (zh) * | 2021-01-07 | 2023-06-20 | 河南理工大学 | 一种时间域电磁测井的地层倾角估算方法 |
CN113504573A (zh) * | 2021-07-06 | 2021-10-15 | 北京航空航天大学 | 一种基于聚焦磁场的井周电阻率测量装置 |
CN114089428B (zh) * | 2021-10-29 | 2022-08-12 | 中国矿业大学 | 一种时间域电磁张量测井方法 |
CN117328862B (zh) * | 2023-12-01 | 2024-02-13 | 齐鲁工业大学(山东省科学院) | 应用于随钻方位测井仪的全息测井方法和发射接收电路 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5115198A (en) * | 1989-09-14 | 1992-05-19 | Halliburton Logging Services, Inc. | Pulsed electromagnetic dipmeter method and apparatus employing coils with finite spacing |
EP2841692A1 (en) * | 2013-07-18 | 2015-03-04 | Halliburton Energy Services, Inc. | Detecting boundary locations of multiple subsurface layers |
CN107748393A (zh) * | 2017-11-10 | 2018-03-02 | 西南石油大学 | 一种基于数值模拟的地层倾角对电阻率影响的校正方法 |
CN109209354A (zh) * | 2018-10-15 | 2019-01-15 | 中国石油大学(华东) | 一种时间域瞬变电磁波测井边界远探测方法 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8112227B2 (en) * | 2004-06-15 | 2012-02-07 | Baker Hughes Incorporated | Processing of multi-component induction measurements in a biaxially anisotropic formation |
US7392137B2 (en) * | 2004-06-15 | 2008-06-24 | Baker Hughes Incorporated | Determination of formation anistrophy, dip and azimuth |
GB2451974B (en) * | 2006-04-26 | 2010-10-20 | Baker Hughes Inc | Method and apparatus for correcting underestimation of formation anistropy ratio |
BRPI0817711B1 (pt) * | 2007-11-30 | 2018-12-18 | Exxonmobil Upstream Res Co | método para detectar planos de falha ou de estratificação cruzada, e, método para produzir hidrocarbonetos a partir de uma região da subsuperfície |
-
2019
- 2019-12-06 CN CN201911237750.XA patent/CN111058834B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5115198A (en) * | 1989-09-14 | 1992-05-19 | Halliburton Logging Services, Inc. | Pulsed electromagnetic dipmeter method and apparatus employing coils with finite spacing |
EP2841692A1 (en) * | 2013-07-18 | 2015-03-04 | Halliburton Energy Services, Inc. | Detecting boundary locations of multiple subsurface layers |
CN107748393A (zh) * | 2017-11-10 | 2018-03-02 | 西南石油大学 | 一种基于数值模拟的地层倾角对电阻率影响的校正方法 |
CN109209354A (zh) * | 2018-10-15 | 2019-01-15 | 中国石油大学(华东) | 一种时间域瞬变电磁波测井边界远探测方法 |
Non-Patent Citations (1)
Title |
---|
沈金松 ; 郭乃川 ; .各向异性层状介质中视电阻率与磁场响应研究.地球物理学报.(第05期),全文. * |
Also Published As
Publication number | Publication date |
---|---|
CN111058834A (zh) | 2020-04-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111058834B (zh) | 基于瞬变多分量感应测井的各向异性地层倾角确定方法 | |
Li et al. | Investigation of eccentricity effects and depth of investigation of azimuthal resistivity LWD tools using 3D finite difference method | |
Dannowski et al. | Estimation of water content and porosity using combined radar and geoelectrical measurements | |
WO2020078003A1 (zh) | 一种时间域瞬变电磁波测井边界远探测方法 | |
Lee et al. | Numerical modeling of eccentered LWD borehole sensors in dipping and fully anisotropic Earth formations | |
US20070256832A1 (en) | Method of analyzing a subterranean formation and method of producing a mineral hydrocarbon fluid from the formation | |
CN104166168B (zh) | 一种井中和地面联合激发电磁数据采集方法 | |
RU2326414C1 (ru) | Способ для использования прибора многокомпонентного индукционного каротажа при управлении параметрами бурения и при интерпретации результатов измерений удельного электрического сопротивления в горизонтальных скважинах | |
Giannino et al. | Electromagnetic methods in geophysics: applications in GeoRadar, FDEM, TDEM, and AEM | |
CN112814668B (zh) | 一种时间域电磁测井的地层倾角估算方法 | |
CN107191181B (zh) | 一种基于电磁散射的井周界面探测方法 | |
CN104360401A (zh) | 一种瞬变电磁b场确定地下目标体地质信息方法 | |
Lin et al. | First evidence of the detection of an underground nuclear magnetic resonance signal in a tunnel | |
NO20181100A1 (en) | Enhanced two dimensional inversion | |
Lin et al. | Rapid and high-resolution detection of urban underground space using transient electromagnetic method | |
CN108169802A (zh) | 一种粗糙介质模型的时域电磁数据慢扩散成像方法 | |
CN107939385B (zh) | 定量计算极化值及应用的方法 | |
Onegova et al. | 3D simulation of transient electromagnetic field for geosteering horizontal wells | |
CN106610509B (zh) | 一种瞬变电磁资料的时域处理方法 | |
Dang et al. | Multi-coil array for long-distance cross-well electromagnetic detection | |
Liu et al. | A study on directional resistivity logging-while-drilling based on self-adaptive hp-FEM | |
Eltsov et al. | Beyond steel casing: detecting zonal isolation in the borehole environment | |
Dutta et al. | Modeling tools for drilling, reservoir navigation, and formation evaluation | |
Chen et al. | Characteristic analysis for dual-induction logging response in highly deviated wells and horizontal wells | |
Yuan et al. | Deep-detection of formation boundary using transient multicomponent electromagnetic logging measurements |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |