CN116398126B - 基于套管开窗水平井小直径随钻电阻率的数据处理方法 - Google Patents
基于套管开窗水平井小直径随钻电阻率的数据处理方法 Download PDFInfo
- Publication number
- CN116398126B CN116398126B CN202310664210.XA CN202310664210A CN116398126B CN 116398126 B CN116398126 B CN 116398126B CN 202310664210 A CN202310664210 A CN 202310664210A CN 116398126 B CN116398126 B CN 116398126B
- Authority
- CN
- China
- Prior art keywords
- data
- resistivity
- sub
- parameters
- measured
- 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
- 238000005553 drilling Methods 0.000 title claims abstract description 66
- 238000003672 processing method Methods 0.000 title claims abstract description 15
- 238000005259 measurement Methods 0.000 claims abstract description 119
- 238000012937 correction Methods 0.000 claims abstract description 70
- 238000000034 method Methods 0.000 claims abstract description 47
- 238000004364 calculation method Methods 0.000 claims description 33
- 230000035772 mutation Effects 0.000 claims description 15
- 238000009825 accumulation Methods 0.000 claims description 9
- 230000002159 abnormal effect Effects 0.000 claims description 6
- 230000005540 biological transmission Effects 0.000 claims description 5
- 239000011435 rock Substances 0.000 abstract description 11
- 239000002689 soil Substances 0.000 abstract description 2
- 230000015572 biosynthetic process Effects 0.000 description 9
- 238000005755 formation reaction Methods 0.000 description 9
- 238000010586 diagram Methods 0.000 description 5
- 230000000694 effects Effects 0.000 description 5
- 230000005856 abnormality Effects 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 239000012530 fluid Substances 0.000 description 3
- 238000003801 milling Methods 0.000 description 3
- XEEYBQQBJWHFJM-UHFFFAOYSA-N Iron Chemical compound [Fe] XEEYBQQBJWHFJM-UHFFFAOYSA-N 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 239000004568 cement Substances 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 230000009545 invasion Effects 0.000 description 2
- 238000000691 measurement method Methods 0.000 description 2
- 230000011218 segmentation Effects 0.000 description 2
- 230000003139 buffering effect Effects 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 229910052742 iron Inorganic materials 0.000 description 1
- 238000000053 physical method Methods 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B49/00—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
- E21B49/08—Obtaining fluid samples or testing fluids, in boreholes or wells
- E21B49/087—Well testing, e.g. testing for reservoir productivity or formation parameters
-
- 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)
- Engineering & Computer Science (AREA)
- Geology (AREA)
- Mining & Mineral Resources (AREA)
- Physics & Mathematics (AREA)
- Environmental & Geological Engineering (AREA)
- Fluid Mechanics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明涉及土层或岩石钻进技术领域,具体为一种基于套管开窗水平井小直径随钻电阻率的数据处理方法,包括:通过先获取随钻测井过程产生的实测数据,其中,所述实测数据包括测量参数与初始电阻率参数,然后根据实测数据的测量参数与电阻率参数对应的分布曲线,构建邻域扩展比,并且基于邻域扩展比对应的窗口半径与纠正尺度,计算电阻率修正值,最后通过所述电阻率修正值对所述电阻率参数进行修正,确认最终电阻率参数,能够根据测量参数带来的影响在不同时间段内与视电阻率分布关系体现的增减变化,对视电阻率进行精准修正,以获得更精准视电阻率,进而降低钻井的风险,提高钻井效率。
Description
技术领域
本发明涉及土层或岩石钻进技术领域,具体为一种基于套管开窗水平井小直径随钻电阻率的数据处理方法。
背景技术
套管开窗,是指在已下套管的井段内进行侧钻时,必须先在套壁上开一个“窗口”,才能把井眼侧钻到套管外的地层里去,这个工艺过程称为套管开窗。开窗前,须在井内打一段水泥塞,造成人工井底,然后在水泥塞上固定一个导斜器。开窗所用的钻具为:钻杆下面接少量钻铤,再接铣锥。开始铣进时,加很小钻压或不加钻压,靠铣锥转动,在套管壁上铣出窗口,当窗口足够长时,可适当扩窗和打捞碎铁。最后换牙轮钻头及其他造斜工具继续造斜侧钻。
套管开窗水平井技术是在定向井、水平井等技术基础上发展的钻井技术,得益于提高油田采收率、节约套管使用费、节约开发成本等优点被推广应用。在钻井过程中,会采用随钻测井技术对钻井过程的数据进行测量与监控。随钻测井是在测井领域的常用技术之一,能够在水平井钻井过程中同时测量相关参数,并将测量结果实时传输至地面的测量技术。
在随钻测井过程中,由于地层环境复杂多变的,层厚、围岩、各向异性、钻井液侵入、井眼环境、仪器与地层的相对夹角等因素都会对地层真实电阻率信息的测量结果产生较大的影响,以导致随钻测井的测量精度较低,进而导致钻井风险增大。
发明内容
鉴于以上内容,有必要提供一种基于套管开窗水平井小直径随钻电阻率的数据处理方法,可以提高随钻测井的测量精度,以获得更精准视电阻率,进而降低钻井的风险,提高钻井效率。
本申请第一方面提供一种基于套管开窗水平井小直径随钻电阻率的数据处理方法,应用于套管开窗水平井的测井领域,所述方法包括:获取随钻测井过程产生的实测数据,其中,所述实测数据包括测量参数与初始电阻率参数;根据实测数据的测量参数与电阻率参数对应的分布曲线,构建邻域扩展比;基于邻域扩展比对应的窗口半径与纠正尺度,计算电阻率修正值;通过所述电阻率修正值对所述电阻率参数进行修正,确认最终电阻率参数。
在其中一种实施例中,所述获取随钻测井过程产生的实测数据,其中,所述实测数据包括测量参数与初始电阻率参数,具体包括:通过脉冲传输方式,接受随钻测井过程产生的测量参数与初始电阻率参数,其中,所述测量参数包括井深、井斜角、方位、测井仪器的倾角、层厚度。
在其中一种实施例中,所述根据实测数据的测量参数与电阻率参数对应的分布曲线,构建邻域扩展比,具体包括:根据不同的实测数据的测量参数与初始电阻率参数同时序下的数值分布情况,形成分布曲线;通过所述分布曲线确认所述实测数据对应的影响值,以构建邻域扩展比。
在其中一种实施例中,所述测量参数包括若干个子测量参数,对应的,所述通过所述分布曲线确认所述实测数据对应的影响值,以构建邻域扩展比,具体包括:将所述实测数据分成若干份相同长度的子区间数据,以获取若干个子区间数据的子测量参数的实测影响值;根据所述若干个子区间数据的子测量参数的实测值,计算每份子区间数据的若干测量参数的平均影响值;将所述若干个子区间数据的子测量参数的实测影响值与平均影响值做差,获取若干个数据偏移度;根据所述若干个数据偏移度与预先构建的标准计算公式,计算邻域扩展比。
在其中一种实施例中,所述根据所述若干个数据偏移度与预先构建的标准计算公式,计算邻域扩展比,具体包括:
;
;
;
其中,子测量参数,m为子区间数据,j是第i种子测量参数的第m个子区间数据
内第j个接收时刻,是指第m个子区间数据第i种子测量参数的实测值,是第m个子区间数
据的第i种子测量参数若干次测量的平均影响值,是每份子区间数据的第i种子测量参数
测量值的分布方差;是第j个接收时刻第i种子测量参数第m个子区间数据的数据偏
移度;
是第m个子区间的区间显著度,是偏移度序列、的自相关
系数,其中,,t是子区间数据的长度,是第m个子区
间数据第i种子测量参数第t个接收时刻的数据偏移度;
是第i种子测量参数第m个子区间的邻域扩展比是第i种子测量参数第m个子区
间数据的邻域扩展比,、分别是第j个接收时刻的幅度比电阻率差值、相位差
电阻率差值,、分别是第m个子区间数据的幅度比电阻率区间均值、相位差电
阻率区间均值。
在其中一种实施例中,所述基于邻域扩展比对应的窗口半径与纠正尺度,计算电阻率修正值,具体包括:根据所有子测量参数在若干个子区间数据的邻域扩展比,确认每个数据接收时刻的窗口半径;通过所述每个数据接收时刻的窗口半径,计算每个数据接受时刻的纠正尺度;根据所述窗口半径和纠正尺度,获取接收时刻的电阻率修正值。
在其中一种实施例中,所述根据所有子测量参数在若干个子区间数据的邻域扩展比,确认每个数据接收时刻的窗口半径,具体包括:将所有子测量参数在若干个子区间数据的邻域扩展比输入第一计算公式,计算若干个子区间数据的扩展累加比;根据第二计算公式与所述若干个子区间数据的扩展累加比,计算每个数据接收时刻的窗口半径;对应的:
所述第一计算公式为:
;
所述第二计算公式为:
;
其中,所述是第m个子测量参数的扩展累加比,n是子测量参数的种类数量,是
第i种子测量参数所有子区间数据对应的邻域扩展比均值;r是窗口半径的初始值,是第j
个数据接收时刻的窗口半径,是第j个数据接收时刻的时间,是第d个累加突变点的时
间,是累加突变点的数量。
在其中一种实施例中,所述通过所述每个数据接收时刻的窗口半径,计算每个数据接受时刻的纠正尺度,具体包括:将所述每个数据接收时刻的窗口半径输入第三计算公式,计算每个数据接受时刻的拟合权重;根据第四计算公式与所述每个数据接受时刻的拟合权重,计算每个数据接受时刻的纠正尺度;对应的:
所述第三计算公式为:
;
所述第四计算公式为:
;
其中,所述是第j个接收时刻的拟合权重,、分别是第k个、第j个接收时刻的
时间,是调参因子;是第j个接收时刻的纠正尺度,是第j个接收时刻的窗口内接收时
刻的数量,是第k个接收时刻的异常得分值。
在其中一种实施例中,所述通过所述电阻率修正值对所述初始电阻率参数进行修正,确认最终电阻率参数,具体包括:
通过预设规则,根据对应的所述电阻率修正值分别对不同层次的初始电阻率参数进行修正,确认不同层次的最终电阻率参数。
本申请实施例通过先获取随钻测井过程产生的实测数据,其中,所述实测数据包括测量参数与初始电阻率参数,然后根据实测数据的测量参数与电阻率参数对应的分布曲线,构建邻域扩展比,并且基于邻域扩展比对应的窗口半径与纠正尺度,计算电阻率修正值,最后通过所述电阻率修正值对所述电阻率参数进行修正,确认最终电阻率参数。通过实测数据的测量参数与电阻率参数对应的分布曲线构建的邻域扩展比,进一步获得对应的窗口半径与纠正尺度计算电阻率修正值,能够根据测量参数带来的影响在不同时间段内与视电阻率分布关系体现的增减变化,对视电阻率进行精准修正,以获得更精准视电阻率,进而降低钻井的风险,提高钻井效率。
附图说明
图1是本申请实施例的基于套管开窗水平井小直径随钻电阻率的数据处理方法的流程示意图。
图2是本申请实施例的基于套管开窗水平井小直径随钻电阻率的数据处理方法的第一子流程示意图。
图3是本申请实施例的基于套管开窗水平井小直径随钻电阻率的数据处理方法的第二子流程示意图。
图4是本申请实施例的基于套管开窗水平井小直径随钻电阻率的数据处理方法的第三子流程示意图。
图5是本申请实施例的基于套管开窗水平井小直径随钻电阻率的数据处理方法的第四子流程示意图。
图6是本申请实施例的基于套管开窗水平井小直径随钻电阻率的数据处理方法的第五子流程示意图。
图7是本申请实施例的优选实施例的一个分布曲线示意图。
图8是本申请实施例的优选实施例的另一个分布曲线示意图。
具体实施方式
在本申请实施例的描述中,“示例性”、“或者”、“例如”等词用于表示作例子、例证或说明。本申请实施例中被描述为“示例性”或者“例如”的任何实施例或设计方案不应被解释为比其它实施例或设计方案更优选或更具优势。确切而言,使用“示例性”、“或者”、“例如”等词旨在以具体方式呈现相关概念。
除非另有定义,本文所使用的所有的技术和科学术语与属于本申请中的技术领域的技术人员通常理解的含义相同。本申请的说明书中所使用的术语只是为了描述具体的实施例的目的,不是旨在于限制本申请。应理解,本申请中除非另有说明,“/”表示或的意思。例如,A/B可以表示A或B。本申请中的“和/或”仅仅是一种描述关联对象的关联关系,表示可以存在三种关系。例如,A和/或B,可以表示:单独存在A,同时存在A和B,单独存在B三种情况。“至少一个”是指一个或者多个。“多个”是指两个或多于两个。例如,a、b或c中的至少一个,可以表示:a,b,c,a和b,a和c,b和c,a、b和c七种情况。
另外需要说明的是,本申请的说明书和权利要求书及附图中的术语“第一”、“第二”是用于部别类似的对象,而不是用于描述特定的顺序或先后次序。本申请实施例中公开的方法或流程图所示出的方法,包括用于实现方法的一个或多个步骤,在不脱离权利要求的范围的情况下,多个步骤的执行顺序可以彼此互换,其中某些步骤也可以被删除。
本申请实施例首先提出一种基于套管开窗水平井小直径随钻电阻率的数据处理方法。请参阅图1,基于套管开窗水平井小直径随钻电阻率的数据处理方法,应用于套管开窗水平井的测井领域,所述方法包括:
S101、获取随钻测井过程产生的实测数据,其中,所述实测数据包括测量参数与初始电阻率参数。
其中,所述随钻测井是在测井领域的常用技术之一,是指能够在水平井钻井过程中同时测量相关参数,并将测量结果实时传输至地面的测量技术。在套管开窗水平井的随钻测井过程中,将测井仪器直接安装在钻头后方一定间隔距离的位置,间隔距离的大小与地层和钻头相关,优选的,本案间隔距离大小可以取经验值5米。优选的,所述测井仪器是指电阻率传感器,例如EWR-Phase4仪器、CPR仪器等。
其中,所述实测数据是指随钻测井过程中能够获取的有关钻井过程的数据,所述实测数据包括测量参数与初始电阻率参数,所述测量参数是指与该井相关的一些角度、深度、位置以及厚度等物理测量因素,所述初始电阻率参数是指根据测井仪器所反馈的当下的有关地层的电阻率,由于地层环境复杂多变的,在随钻测井过程中,层厚、围岩、各向异性、钻井液侵入、井眼环境、仪器与地层的相对夹角等因素都会对地层真实电阻率信息的测量结果产生较大的影响。
在钻井的同时,利用测井仪器采集实测数据,所述实测数据包括井深、井斜角、方位、测井仪器的倾角、层厚度,并利用脉冲传输方式将实测数据传输至地表。
在本申请的一种实施例中,所述获取随钻测井过程产生的实测数据,其中,所述实测数据包括测量参数与初始电阻率参数,具体包括:
通过脉冲传输方式,接受随钻测井过程产生的测量参数与初始电阻率参数,其中,所述测量参数包括井深、井斜角、方位、测井仪器的倾角、层厚度。
S102、根据实测数据的测量参数与初始电阻率参数对应的分布曲线,构建邻域扩展比。
其中,不同的地层环境对应的电阻率参数不同,测量参数包括多种,可以包括井深、井斜角、方位、测井仪器的倾角、层厚度,且每种测量参数的改变,对电阻率参数的测量结果都造成不同的影响,且不同地层深度,测量参数对电阻率参数的测量的影响大小也不同。
需要说明的是,地面上地表接收端将实测数据中的幅度比信号转化成幅度比电阻
率、将得到的相位差信号转化成相位差电阻率,通过、来反映不同地层视电阻
率的情况。原理是通过发射电极或发射线圈向地层发射电磁波,再由二个接收天线接收来
自地层的电磁波的相位差和幅度比,测量的相位差和幅度比与地层的电阻率和介电常数之
间存在函数关系,这样就可以得到地层的电阻率和介电常数。由相位差得到的电阻率称为
相位差电阻率,由幅度比得出的电阻率为幅度比电阻率,介电常数亦同理。视电阻率,是用
来反映岩石和矿石导电性变化的参数。在地下岩石电性分布不均匀(有两种或两种以上导
电性不同的岩石或矿石)或地表起伏不平的情况下,若仍按测定均匀水平大地电阻率的方
法和计算公式求得的电阻率称之为视电阻率,以符号ρs表示,单位和电阻率相同,为Ω·m。
在地下存在多种岩石的情况下用电阻率法测得的电阻率,不是某一种岩石的真电阻率。它除受各种岩石电阻率的综合影响外,还与岩、矿石的分布状态(包括—些构造因素)、电极排列等具体情况有关,所以称它为视电阻率。
另外,随着测量井垂直深度的不同,垂直电阻率与水平电阻率也随着改变,因
为垂直电阻率与水平电阻率也随着改变,对不同地层视电阻率的影响也不同。
因此,建立不同测量参数与电阻率参数对应的分布曲线,则能够将不同测量参数对不同地层视电阻率的影响进行体现。则,可选定一个测量参数,并保持其他测量参数,随着垂直深度的增大,建立与初始电阻率参数对应的分布曲线,可表征不同地层视电阻率的影响。
在一个优选实施例中,选定测量参数为井斜角,并保持其余的测量参数不变,随着
井垂直深度变化,将实测数据中的井斜角与幅度比电阻率、相位差电阻率形成分布
曲线,参考图7与图8。不难看出在层界面处,、均出现极角化,在目的层处,、基
本不变,即井斜角对目的层视电阻率的拟合影响较小。
除此之外,在小角度倾斜井以及直井中,幅度比电阻率、相位差电阻率主要
反映水平电阻率,此时、两者之间的差别较小,且大小基本相等;在大斜度井甚至
水平井中,幅度比电阻率、相位差电阻率会逐渐趋于垂直电阻率的大小,大于,甚至会超过。不难看出在层界面处,、均出现极角化;在目的层处,、基本不变,即井斜角对目的层视电阻率的拟合影响较小。
因此,分别对实测数据的不同的测量参数与初始电阻率参数形成不同的分布曲线,以此来构建邻域扩展比,用于表征不同的测量参数对视电阻率拟合结果的影响程度。
S103、基于邻域扩展比对应的窗口半径与纠正尺度,计算电阻率修正值。
其中,在获取所述实测数据时,测井仪器利用脉冲传输方式将实测数据传输至地面上地表接收端,则地表接收端获取所述实测数据包括多个接收时刻,所述窗口半径则与所述多个接收时刻进行对应,窗口半径是反应当前接收时刻所接收的实测数据的误差大小,则,接收时刻的窗口半径越小,当前接收时刻所对应的接收数据误差概率越大。将所有接收时刻的窗口半径作为特征值,利用LOF算法计算每个接收时刻的异常得分,以此来构建纠正尺度,所述纠正尺度用于表征每个接收时刻接收的实测数据造成误差量的大小,换句话说,纠正尺度反映了每个接收时刻接收的实测数据造成误差量的大小。
在获取到所述窗口半径与纠正尺度之后,将所述窗口半径与纠正尺度作为预设算法的输入参数,以获取最终的电阻率修正值,所述预设算法可以是加权回归Lowess算法,将窗口半径作为局部加权回归Lowess算法的窗口长度,将纠正尺度作为局部加权回归Lowess算法中的权值函数,利用局部加权回归Lowess算法获取TVD值与电阻率的拟合曲线,并利用所获拟合曲线对传输至地表的实测数据中不同层的电阻率进行修正。
S104、通过所述电阻率修正值对所述电阻率参数进行修正,确认最终电阻率参数。
其中,在获取到所述电阻率修正值之后,针对电阻率参数进行修正进行对应修正,以获取最终的电阻率参数。
本申请实施例通过先获取随钻测井过程产生的实测数据,其中,所述实测数据包括测量参数与初始电阻率参数,然后根据实测数据的测量参数与电阻率参数对应的分布曲线,构建邻域扩展比,并且基于邻域扩展比对应的窗口半径与纠正尺度,计算电阻率修正值,最后通过所述电阻率修正值对所述电阻率参数进行修正,确认最终电阻率参数。通过实测数据的测量参数与电阻率参数对应的分布曲线构建的邻域扩展比,进一步获得对应的窗口半径与纠正尺度计算电阻率修正值,能够根据测量参数带来的影响在不同时间段内与视电阻率分布关系体现的增减变化,对视电阻率进行精准修正,以获得更精准视电阻率,进而降低钻井的风险,提高钻井效率。
在本申请的一种实施例中,并参阅图2,所述步骤S102:所述根据实测数据的测量参数与初始电阻率参数对应的分布曲线,构建邻域扩展比,具体包括:
S201、根据不同的实测数据的测量参数与初始电阻率参数同时序下的数值分布情况,形成分布曲线。
其中,所述不同的实测数据的测量参数与初始电阻率参数同时序下的数值分布情况,是指随着某个参数的改变,不同的实测数据的测量参数与初始电阻率参数在同时序下的数值表征。则,可选定一个测量参数,并保持其他测量参数,随着垂直深度的增大,建立与初始电阻率参数对应的分布曲线,可表征不同地层视电阻率的影响。
S202、通过所述分布曲线确认所述实测数据对应的影响值,以构建邻域扩展比。
其中,所述通过所述分布曲线确认所述实测数据对应的影响值,是指不同的实测数据的测量参数对所述初始电阻率参数的影响程度序列值,用于表征具体的测量参数对所述视电阻率的影响程度。
在本申请的一种实施例中,并参阅图3,所述测量参数包括若干个子测量参数,对应的,S202:所述通过所述分布曲线确认所述实测数据对应的影响值,以构建邻域扩展比,具体包括:
S301、将所述实测数据分成若干份相同长度的子区间数据,以获取若干个子区间数据的子测量参数的实测影响值。
其中,将所述实测数据分成等长度数据的子区间数据,且所述子区间数据为若干份,且所述测量参数也包括若干份,则分别对应获取若干份的子区间数据中的若干份子测量参数的实测影响值。
S302、根据所述若干个子区间数据的子测量参数的实测值,计算每份子区间数据的若干测量参数的平均影响值。
其中,在获取若干份的子区间数据中的若干份子测量参数的实测影响值后,对应获取每个子区间数据的若干测量参数的平均影响值。举个例子,其中一个子区间数据对应的若干测量参数分别为井深、井斜角、方位三个,则对所述井深、井斜角、方位三个测量参数的实测影响值进行球平均数,求得该子区间数据的平均影响值。
S303、将所述若干个子区间数据的子测量参数的实测影响值与平均影响值做差,获取若干个数据偏移度。
其中,所述数据偏移度用于表征测量值偏离整体分布的程度。举个例子,假设测井过程中地表接收相邻两次实测数据的时间间隔为t,将N次实测数据分为M个区间长度为t的子区间数据,例如第一子区间为[1,t],即包含第1次至第t次接收的实测数据。对于每个子区间数据,计算子区间内每种子测量参数的均值,其次将实测数据的测量值减去对应的均值,得到每种子测量参数的差分序列。例如对于第一子区间数据中的井斜角,计算t个井斜角的均值,将t个测量值减去均值得到对应的井斜角差分序列。
S304、根据所述若干个数据偏移度与预先构建的标准计算公式,计算邻域扩展比。
其中,所述标准计算公式为:
;
其中,子测量参数,m为子区间数据,j是第i种子测量参数的第m个子区间数据
内第j个接收时刻,是指第m个子区间数据第i种子测量参数的实测值,是第m个子区间数
据的第i种子测量参数若干次测量的平均影响值,是每份子区间数据的第i种子测量参数
测量值的分布方差;是第j个接收时刻第i种子测量参数第m个子区间数据的数据偏
移度;
是第m个子区间的区间显著度,是偏移度序列、的自相关
系数,其中,,t是子区间数据的长度,是第m个子区
间数据第i种子测量参数第t个接收时刻的数据偏移度;
是第i种子测量参数第m个子区间的邻域扩展比是第i种子测量参数第m个子区
间数据的邻域扩展比,、分别是第j个接收时刻的幅度比电阻率差值、相位差
电阻率差值,、分别是第m个子区间数据的幅度比电阻率区间均值、相位差电
阻率区间均值。
需要解析的是,的值越大,第j个接收时刻的测量值偏离整体分布的程度越
大。求自相关系数的过程为公知技术,具体过程不再赘述。的值越大,第m个子区间数据在
所有区间中的显著程度越大,对数据拟合结果的影响越大。、分别是第m个子
区间的幅度比电阻率区间均值、相位差电阻率区间均值,所述幅度比电阻率区间均值是将
第m个子区间的幅度比电阻率相加求均值得到。
需要说明的是,邻域扩展比反映了不同因素对视电阻率拟合结果的影响程度。第j
个接收时刻的测量值偏离整体分布的程度越大,的值越大,的值越大,第m个
子区间在所有区间中越突出,偏移度序列在相邻两个时间段内的变化越大,自相关程度降
低,的值越小,的值越大;第j个接收时刻的幅度比电阻率差值、相位差电
阻率差值越大,的值越大、的值越大,第j个接收时刻越有可能靠近极角化区
域,的值越大,第m个子区间数据的幅度比电阻率区间均值与相位差电阻
率区间均值的大小越接近,对目的层的影响越小,的值越小,即的值越
大,第m个子区间数据对视电阻率拟合结果的影响越大。邻域扩展比考虑了每种测量因素在
不同时间段内与视电阻率分布关系相关关系的增减变化,其有益效果在于能够利用子测量
参数的极角化区域精准评估每个时间时刻实测数据对视电阻率拟合结果的影响程度。
在本申请的一种实施例中,并参阅图4,所述步骤S103:所述基于邻域扩展比对应的窗口半径与纠正尺度,计算电阻率修正值,具体包括:
S401、根据所有子测量参数在若干个子区间数据的邻域扩展比,确认每个数据接收时刻的窗口半径。
S402、通过所述每个数据接收时刻的窗口半径,计算每个数据接受时刻的纠正尺度。
S403、根据所述窗口半径和纠正尺度,获取接收时刻的电阻率修正值。
其中,将所述实测数据分为若干份子区间数据,并获取所有子测量参数在若干个子区间数据的邻域扩展比,再获取若干份个子区间数据的对应的接收时刻的窗口半径,能够将数据进行精准与具体化,进一步提高拟合的精准程度。
具体的,参照附图5,所述S401、根据所有子测量参数在若干个子区间数据的邻域扩展比,确认每个数据接收时刻的窗口半径,具体包括:
S501、将所有子测量参数在若干个子区间数据的邻域扩展比输入第一计算公式,计算若干个子区间数据的扩展累加比;
S501、根据第二计算公式与所述若干个子区间数据的扩展累加比,计算每个数据接收时刻的窗口半径;对应的:
所述第一计算公式为:
;
所述第二计算公式为:
;
其中,所述是第m个子测量参数的扩展累加比,n是子测量参数的种类数量,是
第i种子测量参数所有子区间数据对应的邻域扩展比均值;r是窗口半径的初始值,是第j
个数据接收时刻的窗口半径,是第j个数据接收时刻的时间,是第d个累加突变点的时
间,是累加突变点的数量。
需要解析的是,所述累加突变点是利用BG序列分割算法获取的n种因素偏移度序
列自相关系数组成的序列突变点集合的交集中的突变点。例如第i种测量参数偏移度序列
自相关系数组成的单突变序列为,第n种测量参数偏移度序列自相关系数组成的单突变
序列为,分别利用BG序列分割算法获取n个单突变序列中的突变点,上述突变点的交集即
为累加突变点。的值越大,接收时刻越靠近累加突变点的时间,接收时刻的窗口半径越
小。
至此,得到所有实测数据接收时刻的窗口半径。
能够理解的是,根据接收时刻的窗口半径获取接收时刻相邻时刻出现测量误差数据的概率,如果测井过程中测井仪器出现误差或钻井液浓度等因素发生变化导致实测数据出现异常,那么接收上述异常实测数据的接收时刻的窗口半径相对偏小。除此之外,测井过程中的地层环境复杂,造成实测数据异常的因素不会瞬时消失,通常情况下会经过一定的缓冲周期恢复至正常,即较小的窗口半径的接收时刻具有较大的局部异常密度,则在进行拟合时,应该优先降低拟合点附近实测数据的数据误差,并使得拟合点附近实测数据具有较大的权重,这样才能尽可能小的减少测量误差对拟合曲线的影响,提高电阻率的预测精度。
具体的,参照附图6,所述S402、所述通过所述每个数据接收时刻的窗口半径,计算每个数据接受时刻的纠正尺度,具体包括:
S601、将所述每个数据接收时刻的窗口半径输入第三计算公式,计算每个数据接受时刻的拟合权重;
S602、根据第四计算公式与所述每个数据接受时刻的拟合权重,计算每个数据接受时刻的纠正尺度;对应的:
所述第三计算公式为:
;
所述第四计算公式为:
;
其中,所述是第j个接收时刻的拟合权重,、分别是第k个、第j个接收时刻的
时间,是调参因子;是第j个接收时刻的纠正尺度,是第j个接收时刻的窗口内接收时
刻的数量,是第k个接收时刻的异常得分值。
需要解析的是,的作用是防止分母为0,的大小取经验值0.01。的值越大,
第k个时刻对第j个接收时刻电阻率拟合结果的影响越大。的值越大,第j个接收时刻的窗
口内接收时刻实测数据的异常程度越高。
能够理解的是,纠正尺度反映了每个接收时刻接收的实测数据造成误差量的大
小。在对每个接收时刻的实测数据进行拟合时,第k个、第j个接收时刻之间的时间差越小,
第k个时刻对第j个接收时刻电阻率拟合结果的影响越大,的值越小,的值越大,
窗口内相邻时刻的异常得分越大,的值越大,即的值越大,相邻时刻的接受数据靠
近突变点,在视电阻率的分布曲线中越靠近极角化区域,接受的修正尺度应该越大,使得测
量误差的降到最小。纠正尺度考虑了每个接收时刻拟合时相邻时刻受测量误差异常程度,
其有益效果在于能够利用各子测量参数与视电阻率分布曲线中各底层界面的极角化区域,
提高每个接收时刻电阻率的拟合精度。
在本申请的一种实施例中,所述步骤S104:所述通过所述电阻率修正值对所述初始电阻率参数进行修正,确认最终电阻率参数,具体包括:
通过预设规则,根据对应的所述电阻率修正值分别对不同层次的初始电阻率参数进行修正,确认不同层次的最终电阻率参数。
其中,在获取到所述电阻率修正值之后,可根据预设规则对不同层次的初始电阻率参数进行修正,修正不同层次因不同测量参数造成的电阻率误差。所述预设规则可为加权回归Lowess算法,将窗口半径作为局部加权回归Lowess算法的窗口长度,将纠正尺度作为局部加权回归Lowess算法中的权值函数,则可并利用所获拟合曲线对传输至地表的实测数据中不同层的电阻率进行修正。
需要进一步说明的是,根据上述方法获取到最终电阻率参数后,即说明能够进行提供更加有效的地质导向信息,使得钻头尽可能在高阻含油层中钻进,进而降低钻井的风险,提高钻井效率。
本申请实施例通过先获取随钻测井过程产生的实测数据,其中,所述实测数据包括测量参数与初始电阻率参数,然后根据实测数据的测量参数与电阻率参数对应的分布曲线,构建邻域扩展比,并且基于邻域扩展比对应的窗口半径与纠正尺度,计算电阻率修正值,最后通过所述电阻率修正值对所述电阻率参数进行修正,确认最终电阻率参数。通过实测数据的测量参数与电阻率参数对应的分布曲线构建的邻域扩展比,进一步获得对应的窗口半径与纠正尺度计算电阻率修正值,能够根据测量参数带来的子测量参数在不同时间段内与视电阻率分布关系体现的增减变化,对视电阻率进行精准修正,以获得更精准视电阻率,进而降低钻井的风险,提高钻井效率。
对于本领域技术人员而言,显然本申请不限于上述示范性实施例的细节,而且在不背离本申请的精神或基本特征的情况下,能够以其他的具体形式实现本申请。因此,无论从哪一点来看,均应将本申请上述的实施例看作是示范性的,而且是非限制性的,本申请的范围由所附权利要求而不是上述说明限定,因此旨在将落在权利要求的等同要件的含义和范围内的所有变化涵括在本申请内。
Claims (2)
1.一种基于套管开窗水平井小直径随钻电阻率的数据处理方法,应用于套管开窗水平井的测井领域,其特征在于,所述方法包括:
获取随钻测井过程产生的实测数据,其中,所述实测数据包括测量参数与初始电阻率参数;
根据实测数据的测量参数与初始电阻率参数对应的分布曲线,构建邻域扩展比;
基于邻域扩展比对应的窗口半径与纠正尺度,计算电阻率修正值;
通过所述电阻率修正值对所述电阻率参数进行修正,确认最终电阻率参数;
所述获取随钻测井过程产生的实测数据,其中,所述实测数据包括测量参数与初始电阻率参数,具体包括:
通过脉冲传输方式,接收随钻测井过程产生的测量参数与初始电阻率参数,其中,所述测量参数包括井深、井斜角、方位、测井仪器的倾角、层厚度;
所述根据实测数据的测量参数与初始电阻率参数对应的分布曲线,构建邻域扩展比,具体包括:
根据不同的实测数据的测量参数与初始电阻率参数同时序下的数值分布情况,形成分布曲线;
通过所述分布曲线确认所述实测数据对应的影响值,以构建邻域扩展比;
所述测量参数包括若干个子测量参数,对应的,所述通过所述分布曲线确认所述实测数据对应的影响值,以构建邻域扩展比,具体包括:
将所述实测数据分成若干份相同长度的子区间数据,以获取若干个子区间数据的子测量参数的实测影响值;
根据所述若干个子区间数据的子测量参数的实测值,计算每份子区间数据的若干测量参数的平均影响值;
将所述若干个子区间数据的子测量参数的实测影响值与平均影响值做差,获取若干个数据偏移度;
根据所述若干个数据偏移度与预先构建的标准计算公式,计算邻域扩展比;
所述根据所述若干个数据偏移度与预先构建的标准计算公式,计算邻域扩展比,具体包括:
;
其中,为子测量参数,m为子区间数据,j是第i种子测量参数的第m个子区间数据内第j个接收时刻,/>是指第m个子区间数据第i种子测量参数的实测值,/>是第m个子区间数据第i种子测量参数若干次测量的平均影响值,/>是每份子区间数据的第i种子测量参数测量值的分布方差;/>是第j个接收时刻第i种子测量参数第m个子区间数据的数据偏移度;
是第m个子区间的区间显著度,/>是偏移度序列/>、/>的自相关系数,其中,/>,t是子区间数据的长度,/>是第m个子区间数据第i种子测量参数第t个接收时刻的数据偏移度;
是第i种子测量参数第m个子区间的邻域扩展比是第i种子测量参数第m个子区间数据的邻域扩展比,/>、/>分别是第j个接收时刻的幅度比电阻率差值、相位差电阻率差值,/>、/>分别是第m个子区间数据的幅度比电阻率区间均值、相位差电阻率区间均值;
所述基于邻域扩展比对应的窗口半径与纠正尺度,计算电阻率修正值,具体包括:
根据所有子测量参数在若干个子区间数据的邻域扩展比,确认每个数据接收时刻的窗口半径;
通过所述每个数据接收时刻的窗口半径,计算每个数据接受时刻的纠正尺度;
根据所述窗口半径和纠正尺度,获取接收时刻的电阻率修正值;
所述根据所有子测量参数在若干个子区间数据的邻域扩展比,确认每个数据接收时刻的窗口半径,具体包括:
将所有子测量参数在若干个子区间数据的邻域扩展比输入第一计算公式,计算若干个子区间数据的扩展累加比;
根据第二计算公式与所述若干个子区间数据的扩展累加比,计算每个数据接收时刻的窗口半径;对应的:
所述第一计算公式为:
;
所述第二计算公式为:
;
其中,所述是第m个子测量参数的扩展累加比,n是子测量参数的种类数量,/>是第i种子测量参数所有子区间数据对应的邻域扩展比均值;r是窗口半径的初始值,/>是第j个数据接收时刻的窗口半径,/>是第j个数据接收时刻的时间,/>是第d个累加突变点的时间,/>是累加突变点的数量;
所述通过所述每个数据接收时刻的窗口半径,计算每个数据接受时刻的纠正尺度,具体包括:
将所述每个数据接收时刻的窗口半径输入第三计算公式,计算每个数据接受时刻的拟合权重;
根据第四计算公式与所述每个数据接受时刻的拟合权重,计算每个数据接受时刻的纠正尺度;对应的:
所述第三计算公式为:
;
所述第四计算公式为:
;
其中,所述是第j个接收时刻的拟合权重,/>、/>分别是第k个、第j个接收时刻的时间,是调参因子;/>是第j个接收时刻的纠正尺度,/>是第j个接收时刻的窗口内接收时刻的数量,/>是第k个接收时刻的异常得分值。
2.根据权利要求1所述的基于套管开窗水平井小直径随钻电阻率的数据处理方法,其特征在于,所述通过所述电阻率修正值对所述初始电阻率参数进行修正,确认最终电阻率参数,具体包括:
通过预设规则,根据对应的所述电阻率修正值分别对不同层次的初始电阻率参数进行修正,确认不同层次的最终电阻率参数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310664210.XA CN116398126B (zh) | 2023-06-07 | 2023-06-07 | 基于套管开窗水平井小直径随钻电阻率的数据处理方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310664210.XA CN116398126B (zh) | 2023-06-07 | 2023-06-07 | 基于套管开窗水平井小直径随钻电阻率的数据处理方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN116398126A CN116398126A (zh) | 2023-07-07 |
CN116398126B true CN116398126B (zh) | 2023-08-25 |
Family
ID=87009035
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310664210.XA Active CN116398126B (zh) | 2023-06-07 | 2023-06-07 | 基于套管开窗水平井小直径随钻电阻率的数据处理方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116398126B (zh) |
Citations (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5886526A (en) * | 1996-06-19 | 1999-03-23 | Schlumberger Technology Corporation | Apparatus and method for determining properties of anisotropic earth formations |
CN101139925A (zh) * | 2006-09-08 | 2008-03-12 | 西南石油大学 | 一种随钻测试储层参数特性并实时调整钻井措施的方法 |
CN102678106A (zh) * | 2012-05-02 | 2012-09-19 | 中国电子科技集团公司第二十二研究所 | 随钻电磁波电阻率测井仪器的数据处理方法 |
CN103410504A (zh) * | 2013-07-22 | 2013-11-27 | 中国石油天然气股份有限公司 | 一种水平井/大斜度井确定地层真电阻率的方法及装置 |
CN105484740A (zh) * | 2015-12-04 | 2016-04-13 | 中国石油天然气集团公司 | 一种用于探测地层复电阻率的多频探测装置 |
WO2016057122A1 (en) * | 2014-10-10 | 2016-04-14 | Halliburton Energy Services, Inc. | Formation resistivity measurement apparatus, systems, and methods |
CN105545301A (zh) * | 2015-12-24 | 2016-05-04 | 西安石油大学 | 一种不同泥浆体系下复杂储层流体性质识别方法 |
CN107045154A (zh) * | 2017-02-08 | 2017-08-15 | 中国海洋石油总公司 | 一种水平井环境中的识别地层产状的方法和装置 |
CN107762490A (zh) * | 2016-08-17 | 2018-03-06 | 中国石油化工股份有限公司 | 一种水平井基于双侧向测井真电阻率反演方法 |
CN108005646A (zh) * | 2017-12-14 | 2018-05-08 | 中国石油大学(华东) | 基于随钻电磁波测井资料的地层各向异性电阻率提取方法 |
CN108019206A (zh) * | 2017-11-16 | 2018-05-11 | 中国石油集团长城钻探工程有限公司 | 一种高介电常数下随钻电磁波电阻率仪器量程扩展方法 |
CN108073765A (zh) * | 2017-12-11 | 2018-05-25 | 中海石油(中国)有限公司湛江分公司 | 一种水平井常规随钻测井地层界面识别与边界距反演方法 |
CN108547611A (zh) * | 2018-03-01 | 2018-09-18 | 杭州迅美科技有限公司 | 水平井复杂环境随钻电磁波电阻率测井快速仿真方法 |
CN109138975A (zh) * | 2018-07-10 | 2019-01-04 | 中国海洋石油集团有限公司 | 一种基于时移测井数据的求解相渗特征曲线的新方法 |
CN110005398A (zh) * | 2019-04-04 | 2019-07-12 | 中国石油大学(北京) | 随钻电磁波电阻率测井仪器设计参数确定方法及装置 |
CN110685600A (zh) * | 2018-06-20 | 2020-01-14 | 中国石油化工股份有限公司 | 一种用于地质导向的钻头调整预测方法 |
CN113779812A (zh) * | 2021-09-28 | 2021-12-10 | 中海油田服务股份有限公司 | 一种大斜度井环境中利用随钻测井数据的薄夹层识别方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7301338B2 (en) * | 2001-08-13 | 2007-11-27 | Baker Hughes Incorporated | Automatic adjustment of NMR pulse sequence to optimize SNR based on real time analysis |
US8035392B2 (en) * | 2008-10-17 | 2011-10-11 | Baker Hughes Incorporated | Method and apparatus for while-drilling transient resistivity measurements |
-
2023
- 2023-06-07 CN CN202310664210.XA patent/CN116398126B/zh active Active
Patent Citations (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5886526A (en) * | 1996-06-19 | 1999-03-23 | Schlumberger Technology Corporation | Apparatus and method for determining properties of anisotropic earth formations |
CN101139925A (zh) * | 2006-09-08 | 2008-03-12 | 西南石油大学 | 一种随钻测试储层参数特性并实时调整钻井措施的方法 |
CN102678106A (zh) * | 2012-05-02 | 2012-09-19 | 中国电子科技集团公司第二十二研究所 | 随钻电磁波电阻率测井仪器的数据处理方法 |
CN103410504A (zh) * | 2013-07-22 | 2013-11-27 | 中国石油天然气股份有限公司 | 一种水平井/大斜度井确定地层真电阻率的方法及装置 |
WO2016057122A1 (en) * | 2014-10-10 | 2016-04-14 | Halliburton Energy Services, Inc. | Formation resistivity measurement apparatus, systems, and methods |
CN105484740A (zh) * | 2015-12-04 | 2016-04-13 | 中国石油天然气集团公司 | 一种用于探测地层复电阻率的多频探测装置 |
CN105545301A (zh) * | 2015-12-24 | 2016-05-04 | 西安石油大学 | 一种不同泥浆体系下复杂储层流体性质识别方法 |
CN107762490A (zh) * | 2016-08-17 | 2018-03-06 | 中国石油化工股份有限公司 | 一种水平井基于双侧向测井真电阻率反演方法 |
CN107045154A (zh) * | 2017-02-08 | 2017-08-15 | 中国海洋石油总公司 | 一种水平井环境中的识别地层产状的方法和装置 |
CN108019206A (zh) * | 2017-11-16 | 2018-05-11 | 中国石油集团长城钻探工程有限公司 | 一种高介电常数下随钻电磁波电阻率仪器量程扩展方法 |
CN108073765A (zh) * | 2017-12-11 | 2018-05-25 | 中海石油(中国)有限公司湛江分公司 | 一种水平井常规随钻测井地层界面识别与边界距反演方法 |
CN108005646A (zh) * | 2017-12-14 | 2018-05-08 | 中国石油大学(华东) | 基于随钻电磁波测井资料的地层各向异性电阻率提取方法 |
CN108547611A (zh) * | 2018-03-01 | 2018-09-18 | 杭州迅美科技有限公司 | 水平井复杂环境随钻电磁波电阻率测井快速仿真方法 |
CN110685600A (zh) * | 2018-06-20 | 2020-01-14 | 中国石油化工股份有限公司 | 一种用于地质导向的钻头调整预测方法 |
CN109138975A (zh) * | 2018-07-10 | 2019-01-04 | 中国海洋石油集团有限公司 | 一种基于时移测井数据的求解相渗特征曲线的新方法 |
CN110005398A (zh) * | 2019-04-04 | 2019-07-12 | 中国石油大学(北京) | 随钻电磁波电阻率测井仪器设计参数确定方法及装置 |
CN113779812A (zh) * | 2021-09-28 | 2021-12-10 | 中海油田服务股份有限公司 | 一种大斜度井环境中利用随钻测井数据的薄夹层识别方法 |
Non-Patent Citations (1)
Title |
---|
随钻电磁波电阻率测井曲线分离关系研究;岳喜洲;马明学;李国玉;;石油管材与仪器;第02卷(第02期);第53页-第56页+第59页 * |
Also Published As
Publication number | Publication date |
---|---|
CN116398126A (zh) | 2023-07-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
EP2320251B1 (en) | Multilevel workflow method for extracting resistivity anisotropy data from 3D induction measurements | |
US6353321B1 (en) | Uncompensated electromagnetic wave resistivity tool for bed boundary detection and invasion profiling | |
US6594584B1 (en) | Method for calculating a distance between a well logging instrument and a formation boundary by inversion processing measurements from the logging instrument | |
US10605072B2 (en) | Well ranging apparatus, systems, and methods | |
EP2836861B1 (en) | Resistivity logging system and method employing ratio signal set for inversion | |
US10358911B2 (en) | Tilted antenna logging systems and methods yielding robust measurement signals | |
EP2941534B1 (en) | Deep formation evaluation systems and methods | |
RU2582477C1 (ru) | Электромагнитный способ получения азимутального угла падения | |
US8749243B2 (en) | Real time determination of casing location and distance with tilted antenna measurement | |
AU2013394401B2 (en) | Detecting boundary locations of multiple subsurface layers | |
US20120068712A1 (en) | System and method for using a resistivity tool with wired drill pipe and one or more wells | |
Kuckes et al. | New electromagnetic surveying/ranging method for drilling parallel horizontal twin wells | |
US20060017443A1 (en) | Deep reading propagation resistivity tool for determination of distance to a bed boundary with a transition zone | |
US7663968B2 (en) | Method of processing geological data | |
US20180347353A1 (en) | Improved mci logging for processing downhole measurements | |
CN116411950A (zh) | 一种随钻前探测电磁波天线***及测量方法 | |
CN116398126B (zh) | 基于套管开窗水平井小直径随钻电阻率的数据处理方法 | |
US20180329105A1 (en) | Utilizing Diverse Excitation Sources in Electromagnetic Ranging | |
CN110374579A (zh) | 随钻方位电磁波测井深浅电阻率预测控制方法 | |
US11885925B2 (en) | System and methods for evaluating a formation using pixelated solutions of formation data | |
WO2022266614A9 (en) | Resistivity determination from one transmitter and one receiver antennas | |
WO2017019031A1 (en) | Assessment of formation true dip, true azimuth, and data quality with multicomponent induction and directional logging | |
JP2002296347A (ja) | トンネル掘進機の地中探査レーダの探査距離測定方法及び測定装置 |
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 | ||
PE01 | Entry into force of the registration of the contract for pledge of patent right |
Denomination of invention: Data processing method for small diameter resistivity while drilling in horizontal wells based on casing window opening Effective date of registration: 20231115 Granted publication date: 20230825 Pledgee: China Construction Bank Corporation Dongying nonggao District sub branch Pledgor: Shandong Wan Yang Petroleum Technology Co.,Ltd. Registration number: Y2023980065946 |
|
PE01 | Entry into force of the registration of the contract for pledge of patent right |