CN116556925A - 井孔椭圆度的计算方法、装置、计算设备及存储介质 - Google Patents
井孔椭圆度的计算方法、装置、计算设备及存储介质 Download PDFInfo
- Publication number
- CN116556925A CN116556925A CN202310507713.6A CN202310507713A CN116556925A CN 116556925 A CN116556925 A CN 116556925A CN 202310507713 A CN202310507713 A CN 202310507713A CN 116556925 A CN116556925 A CN 116556925A
- Authority
- CN
- China
- Prior art keywords
- depth point
- azimuth
- depth
- logging
- logging instrument
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 58
- 238000012545 processing Methods 0.000 claims abstract description 36
- 238000001914 filtration Methods 0.000 claims description 38
- 238000004891 communication Methods 0.000 claims description 16
- 238000004364 calculation method Methods 0.000 claims description 14
- 238000004458 analytical method Methods 0.000 claims description 12
- 238000001514 detection method Methods 0.000 claims description 7
- 239000002002 slurry Substances 0.000 claims description 7
- 238000005553 drilling Methods 0.000 abstract description 13
- 238000010586 diagram Methods 0.000 description 13
- 238000005259 measurement Methods 0.000 description 3
- 238000004590 computer program Methods 0.000 description 2
- 238000002592 echocardiography Methods 0.000 description 2
- 230000009897 systematic effect Effects 0.000 description 2
- 238000010276 construction Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000003628 erosive effect Effects 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000011022 operating instruction Methods 0.000 description 1
- 238000001028 reflection method Methods 0.000 description 1
- 230000035945 sensitivity Effects 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
- E21B47/00—Survey of boreholes or wells
-
- 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
- E21B47/00—Survey of boreholes or wells
- E21B47/08—Measuring diameters or related dimensions at the borehole
-
- 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
- E21B47/00—Survey of boreholes or wells
- E21B47/12—Means for transmitting measuring-signals or control signals from the well to the surface, or from the surface to the well, e.g. for logging while drilling
- E21B47/14—Means for transmitting measuring-signals or control signals from the well to the surface, or from the surface to the well, e.g. for logging while drilling using acoustic waves
-
- 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
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geology (AREA)
- Life Sciences & Earth Sciences (AREA)
- Mining & Mineral Resources (AREA)
- Environmental & Geological Engineering (AREA)
- Fluid Mechanics (AREA)
- Geophysics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Acoustics & Sound (AREA)
- Remote Sensing (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
本发明公开了一种井孔椭圆度的计算方法、装置、计算设备及存储介质,该方法包括:获取深度区间内各个深度点处各个方位的反射回波测井数据;根据各个深度点处各个方位的反射回波测井数据,计算各个深度点处各个方位上测井仪器中心与井壁的距离;根据各个深度点中的待处理深度点及其相邻的多个深度点处各个方位上测井仪器中心与井壁的距离以及各个方位的角度,进行椭圆拟合处理;根据拟合结果确定待处理深度点的井孔椭圆度。通过上述方式,采用了相邻多个深度点的数据进行椭圆拟合,仅仅利用随钻超声井径测井仪器即可以获取到井孔椭圆度,并且能够提高拟合井孔椭圆度的精准性。
Description
技术领域
本发明涉及测井技术领域,具体涉及一种井孔椭圆度的计算方法、装置、计算设备及存储介质。
背景技术
随钻超声井径测井利用了超声波反射法,可实时获取井眼直径。在钻井过程中,利用该技术可实时监测井筒的状况,为井眼的冲蚀和井眼的稳定性提供预警,从而为钻井安全提供必要测井技术参数,另外,精确的井径结果为随钻中子孔隙度的校正提供支持,对正确的随钻地层评价具有重要意义。
目前,随钻超声井径测井装置安装有三个按规律排列的换能器,测量时三个超声波换能器同时发射超声波,根据几何运算方法,就能确定井孔的直径、圆心以及其他相关参数。但是,实际的井眼可能不是圆形的,尤其在地应力的作用下,很多井眼都是椭圆井眼,现有的随钻超声井径测井仪器无法获取井孔椭圆度。
发明内容
鉴于上述问题,提出了本发明以便提供一种克服上述问题或者至少部分地解决上述问题的井孔椭圆度的计算方法、装置、计算设备及存储介质。
根据本发明的一个方面,提供了一种井孔椭圆度的计算方法,方法包括:
获取深度区间内各个深度点处各个方位的反射回波测井数据;
根据各个深度点处各个方位的反射回波测井数据,计算各个深度点处各个方位上测井仪器中心与井壁的距离;
根据各个深度点中的待处理深度点及其相邻的多个深度点处各个方位上测井仪器中心与井壁的距离以及各个方位的角度,进行椭圆拟合处理;
根据拟合结果确定待处理深度点的井孔椭圆度。
在一种可选的方式中,根据各个深度点处各个方位的反射回波测井数据,计算各个深度点处各个方位上测井仪器中心与井壁的距离进一步包括:
针对于任一深度点,对该深度点处各个方位的反射回波测井数据进行峰值检测,得到各个波形峰值的时间信息;
根据各个波形峰值的时间信息、泥浆声速以及测井仪器的半径,计算得到各个方位上测井仪器中心与井壁的距离。
在一种可选的方式中,方法进一步包括:
利用测井仪器对标准套管进行测量,获得标准套管的反射回波测井数据;根据标准套管的反射回波测井数据及内径,确定系统时间误差;
根据各个波形峰值的时间信息、泥浆声速以及测井仪器的半径,计算得到各个方位上测井仪器中心与井壁的距离进一步包括:
根据系统时间误差、各个波形峰值的时间信息、泥浆声速以及测井仪器的半径,计算得到各个方位上测井仪器中心与井壁的距离。
在一种可选的方式中,根据待处理深度点及其相邻的多个深度点处各个方位上测井仪器中心与井壁的距离以及各个方位的角度,进行椭圆拟合处理进一步包括:
针对于待处理深度点及其相邻的多个深度点中的任一深度点,根据该深度点处任一方位上测井仪器中心与井壁的距离以及该方位的角度,计算对应在预设坐标系中的坐标值;
根据各个坐标值进行椭圆拟合处理。
在一种可选的方式中,获取深度区间内各个深度点处各个方位的反射回波测井数据之后,方法还包括:
对各个深度点处各个方位的反射回波测井数据进行中值滤波处理;
根据各个深度点处各个方位的反射回波测井数据,计算各个深度点处各个方位上测井仪器中心与井壁的距离进一步包括:
根据各个深度点处的各个方位的滤波处理后的反射回波测井数据,计算各个深度点处各个方位上测井仪器中心与井壁的距离。
在一种可选的方式中,对各个深度点处各个方位的反射回波测井数据进行中值滤波处理进一步包括:
对于任一深度点,根据该深度点及其相邻的多个深度点处的任一方位的反射回波测井数据,对该深度点处的该方位的反射回波测井数据进行中值滤波处理。
在一种可选的方式中,对各个深度点处各个方位的反射回波测井数据进行中值滤波处理进一步包括:
对于任一深度点,获取各个方位的任一时刻的各个波形幅值中的波形幅度中值;将各个方位的该时刻的各个波形幅值分别减去该时刻的波形幅度中值。
根据本发明的另一方面,提供了一种井孔椭圆度的计算装置,装置包括:
获取模块,适于获取深度区间内各个深度点处各个方位的反射回波测井数据;
分析模块,适于根据各个深度点处各个方位的反射回波测井数据,计算各个深度点处各个方位上测井仪器中心与井壁的距离;
拟合模块,适于根据各个深度点中的待处理深度点及其相邻的多个深度点处各个方位上测井仪器中心与井壁的距离以及各个方位的角度,进行椭圆拟合处理;
处理模块,适于根据拟合结果确定待处理深度点的井孔椭圆度。
在一种可选的方式中,分析模块进一步适于:
针对于任一深度点,对该深度点处各个方位的反射回波测井数据进行峰值检测,得到各个波形峰值的时间信息;根据各个波形峰值的时间信息、泥浆声速以及测井仪器的半径,计算得到各个方位上测井仪器中心与井壁的距离。
在一种可选的方式中,分析模块进一步适于:
获取标准套管的反射回波测井数据;其中,标准套管的反射回波测井数据利用测井仪器对标准套管进行测量得到;根据标准套管的反射回波测井数据及内径,确定系统时间误差;
根据系统时间误差、各个波形峰值的时间信息、泥浆声速以及测井仪器的半径,计算得到各个方位上测井仪器中心与井壁的距离。
在一种可选的方式中,拟合模块进一步适于:
针对于待处理深度点及其相邻的多个深度点中的任一深度点,根据该深度点处任一方位上测井仪器与井壁的距离以及该方位的角度,计算对应在预设坐标系中的坐标值;根据各个坐标值进行椭圆拟合处理。
在一种可选的方式中,装置还包括:滤波模块,适于对各个深度点处各个方位的反射回波测井数据进行中值滤波处理;
分析模块进一步适于:根据各个深度点处的各个方位的滤波处理后的反射回波测井数据,计算各个深度点处各个方位上测井仪器中心与井壁的距离。
在一种可选的方式中,滤波模块进一步适于:对于任一深度点,根据该深度点及其相邻的多个深度点处的任一方位的反射回波测井数据,对该深度点处的该方位的反射回波测井数据进行中值滤波处理。
在一种可选的方式中,滤波模块进一步适于:对于任一深度点,获取各个方位的任一时刻的各个波形幅值中的波形幅度中值;将各个方位的该时刻的各个波形幅值分别减去该时刻的波形幅度中值。
根据本发明的又一方面,提供了一种计算设备,包括:处理器、存储器、通信接口和通信总线,所述处理器、所述存储器和所述通信接口通过所述通信总线完成相互间的通信;
所述存储器用于存放至少一可执行指令,所述可执行指令使所述处理器执行上述井孔椭圆度的计算方法对应的操作。
根据本发明的再一方面,提供了一种计算机存储介质,所述存储介质中存储有至少一可执行指令,所述可执行指令使处理器执行如上述井孔椭圆度的计算方法对应的操作。
根据本发明的井孔椭圆度的计算方法、装置、计算设备及存储介质,该方法包括:获取深度区间内各个深度点处各个方位的反射回波测井数据;根据各个深度点处各个方位的反射回波测井数据,计算各个深度点处各个方位上测井仪器中心与井壁的距离;根据各个深度点中的待处理深度点及其相邻的多个深度点处的各个方位上测井仪器中心与井壁的距离以及各个方位的角度,进行椭圆拟合处理;根据拟合结果确定待处理深度点的井孔椭圆度。通过上述方式,采用了相邻多个深度点的数据进行椭圆拟合,仅仅利用随钻超声井径测井仪器即可以获取到井孔椭圆度,并且能够提高拟合井孔椭圆度的精准性。
上述说明仅是本发明技术方案的概述,为了能够更清楚了解本发明的技术手段,而可依照说明书的内容予以实施,并且为了让本发明的上述和其它目的、特征和优点能够更明显易懂,以下特举本发明的具体实施方式。
附图说明
通过阅读下文优选实施方式的详细描述,各种其他的优点和益处对于本领域普通技术人员将变得清楚明了。附图仅用于示出优选实施方式的目的,而并不认为是对本发明的限制。而且在整个附图中,用相同的参考符号表示相同的部件。在附图中:
图1示出了本发明实施例提供的井孔椭圆度的计算方法的流程图;
图2示出了本发明实施例中多个深度点的波形分析结果的示意图;
图3示出了本发明另一实施例中基于深度中值滤波前后的波形的示意图;
图4示出了本发明另一实施例中基于方位中值滤波前后的波形的示意图;
图5示出了本申请另一实施例中标准套管与测井仪器的截面的示意图;
图6示出了本发明另一实施例中井孔椭圆度计算结果的示意图;
图7示出了本发明实施例提供的井孔椭圆度的计算装置的结构示意图;
图8示出了本发明实施例提供的计算设备的结构示意图。
具体实施方式
下面将参照附图更详细地描述本发明的示例性实施例。虽然附图中显示了本发明的示例性实施例,然而应当理解,可以以各种形式实现本发明而不应被这里阐述的实施例所限制。相反,提供这些实施例是为了能够更透彻地理解本发明,并且能够将本发明的范围完整的传达给本领域的技术人员。
图1示出了本发明实施例提供的井孔椭圆度的计算方法的流程图,如图1所示,该方法包括以下步骤:
步骤S110,获取深度区间内各个深度点处各个方位的反射回波测井数据。
使用测井仪器(随钻超声井径测井仪器)在深度区间内进行测井,在深度区间内每一个深度点处所采集的波形道数为N,即N个不同方位的反射回波测井数据,波形道数与实际的随钻超声井径测井仪器有关,方位与超声波换能器是一一对应的,一个方位的反射回波测井数据,即是通过相应超声波换能器发射超声波而采集的反射回波测井数据。随钻超声井径测井仪器一般包含3个或者4个超声波换能器,能够采集3道或4道不同方位的反射回波测井数据。
步骤S120,根据各个深度点处各个方位的反射回波测井数据,计算各个深度点处各个方位上测井仪器中心与井壁的距离。
对于每一个深度点,根据该深度点处采集的各个方位的反射回波测井数据,计算该深度点处各个方位上测井仪器中心与井壁之间的距离。
步骤S130,根据各个深度点中的待处理深度点及其相邻的多个深度点处各个方位上测井仪器中心与井壁的距离以及各个方位的角度,进行椭圆拟合处理。
当需要计算待处理深度点的椭圆度时,根据待处理深度点、该待处理深度点之上的相邻的K个深度点以及该待处理深度点之下的相邻的K个深度点(即总共2k+1个深度点)处的数据,进行椭圆拟合处理,每一个深度点的数据包括各个方位上测井仪器中心与井壁之间的距离以及各个方位的角度。
其中,方位的角度根据测井仪器本身的旋转角度以及超声波换能器在测井仪器上分布的角度而确定。
图2示出了本发明实施例中多个深度点的波形分析结果的示意图,其中,其中,深度点nDep-1表示待处理深度点nDep上方的第一个深度点,深度点nDep-2表示待处理深度点nDep上方的第二个深度点,深度点nDep+1表示待处理深度点nDep下方的第一个深度点,深度点nDep+2表示待处理深度点nDep下方的第二个深度点。如图2所示,钻铤在井孔中的不同位置。最外层的正方形为地层,中间的椭圆为井孔,最里面的圆形为钻铤示意图。钻铤上安装了3个超声波换能器,T1、T2、T3分别代表3个超声波换能器的方位,如图所示,在不同的深度点,仪器分布具有不同的方位角度。根据这5个深度点上所计算的3个测井仪器中心与井壁的距离以及3个超声波换能器的方位角度进行椭圆拟合。
步骤S140,根据拟合结果确定待处理深度点的井孔椭圆度。
井孔椭圆度的定义是:(最大直径-最小直径)/标准井眼直径,通过上述步骤拟合得到待处理深度点的椭圆信息,进一步计算待处理深度点的井孔椭圆度。
根据本申请实施例的井孔椭圆度的计算方法,获取深度区间内各个深度点处各个方位的反射回波测井数据;根据各个深度点处各个方位的反射回波测井数据,计算各个深度点处各个方位上测井仪器中心与井壁的距离;根据各个深度点中的待处理深度点及其相邻的多个深度点处各个方位上测井仪器中心与井壁的距离以及各个方位的角度,进行椭圆拟合处理;根据拟合结果确定待处理深度点的井孔椭圆度。通过上述方式,采用了相邻多个深度点的数据进行椭圆拟合,仅仅利用随钻超声井径测井仪器即可以获取到井孔椭圆度,并且能够提高拟合井孔椭圆度的精准性。
现有的随钻超声井径测井仪器所测量的超声反射回波信号中,经常会有发射脉冲余振的信号,该信号会干扰井眼反射回波信号的提取,进而导致井径计算错误,具体地,超声波换能器在激励过程中,会产生震荡,从而导致超声波信号中有余振信号,如果余振信号持续时间较长,或者超声波换能器距离井壁很近,又或者在重泥浆环境中超声反射回波信号衰减严重(在此情况下很容易检测到换能器余振信号),都会导致反射回波的峰值检测错误,从而无法获取真正的井径信息。
基于此,在一种可选的实施方式中,获得各个深度点处各个方位的反射回波测井数据之后,对各个深度点处各个方位的反射回波测井数据进行中值滤波处理,则在后续过程中,均以滤波后的反射回波测井数据计算用于拟合椭圆的相关参数。中值滤波处理是一种非线性数字滤波器技术,经常用于去除图像或者其它信号中的噪声。通过对反射回波信号进行中值滤波,能够消除超声波换能器余振的影响。
在一种可选的实施方式中,采用深度中值进行滤波处理,具体地:对于任一深度点,根据该深度点及其相邻的多个深度点处的任一方位的反射回波测井数据,对该深度点处的该方位的反射回波测井数据进行中值滤波处理。
根据该深度点及其相邻的多个深度点处的任一方位的反射回波测井数据,获取该深度点及其相邻的多个深度点的该方位上任一时刻的各个波形幅值,并确定各个波形幅值中的波形幅值中值,然后将该深度点处该方位、该时刻的波形幅值减去波形幅值中值,即得到该深度点处该方位的滤波处理后的反射回波测井数据。对于每个方位的反射回波测井数据,均按照同样的方法进行处理,得到该深度点的经过滤波处理的反射回波测井数据。
举例来说,为了对待处理深度点nDep的T1方位的反射回波测井数据进行中值滤波,分别获取深度点nDep-1、深度点nDep-2、待处理深度点nDep、深度点nDep+1以及深度点nDep+2的T1方位、S1时刻的波形幅值,并确定波形幅值中值,将待处理深度点nDep的T1方位的S1时刻的波形幅值减去该波形幅值中值,对于每一个时刻,均按照同样的方法进行处理,从而得到待处理深度点nDep的T1方位的滤波处理后的反射回波测井数据。
在另一种可选的实施方式中,采用方位中值进行滤波处理,具体地:对于任一深度点,获取各个方位的任一时刻的各个波形幅值中的波形幅度中值;将各个方位的该时刻的各个波形幅值分别减去该时刻的波形幅度中值。
举例来说,为了对待处理深度点nDep的三个方位的反射回波测井数据进行中值滤波,根据三个方位的反射回波测井数据,获取三个方位在S1时刻的波形幅值,并确定三个波形幅值中的波形幅度中值,然后,将三个方位对应在S1时刻的波形幅值减去S1时刻的波形幅度中值,对于每一个时刻均按照同样的方式进行处理,即得到待处理深度点nDep的三个方位的滤波处理后的反射回波测井数据。
图3示出了本发明另一实施例中基于深度中值滤波前后的波形的示意图,图中分别绘制了方位1、方位2、方位3的原始反射回波的波形,以及各个方位的经过中值滤波后的波形,以图3中虚线框内的波形为例,方位3原始波形中从起始时间开始的一段时间内存在余振信号波形,而经过深度中值滤波处理后,相应时段的波形曲线呈现为直线,即超声波换能器的余振信号已经大大消除了。
图4示出了本发明另一实施例中基于方位中值滤波前后的波形的示意图,图中分别绘制了方位1、方位2、方位3的原始反射回波的波形,以及各个方位的经过方位中值滤波后的波形,由图4可见,部分深度点的超声波换能器的余振信号已经削弱了,但是某些深度点效果一般,主要的原因是超声波换能器接收灵敏度不一致,余振的相位不可能完全一致。
在一种可选的实施方式中,根据各个深度点处各个方位的反射回波测井数据,计算各个深度点处各个方位上测井仪器中心与井壁的距离进一步包括:
针对于任一深度点,对该深度点处各个方位的反射回波测井数据进行峰值检测,得到各个波形峰值的时间信息;根据各个波形峰值的时间信息、泥浆声速以及测井仪器的半径,计算得到各个方位上测井仪器中心与井壁的距离。其中,峰值检测的结果可以取信号的包络或者直接取信号幅度峰值。
针对于任意一个深度点,根据N个方位的反射回波测井数据,分别获取其N个方位的反射回波波形的幅度最大的时间位置:TRi=max(Z),i=1,2……N,Z表示波形幅值;然后,分别确定N个方位上测井仪器中心与井壁的距离(用Di表示),具体计算公式如下:
其中,V表示泥浆声速,可以通过泥浆声速计进行测量,表示i方位上超声波换能器与井壁的距离,R表示测井仪器的半径。
在另一种可选的实施方式中,方法还包括:利用测井仪器对标准套管进行测量,获得标准套管的反射回波测井数据;根据标准套管的反射回波测井数据及内径,确定系统时间误差。则计算各个方位上测井仪器中心与井壁的距离的步骤进一步包括:根据系统时间误差、各个波形峰值的时间信息、泥浆声速以及测井仪器的半径,计算得到各个方位上测井仪器中心与井壁的距离。
信号最大幅值对应的时间或者包络信号对应的时间,均不是信号的起点,信号的起始点波形幅度小,很难测量到,因此,如果以波形最大幅值对应的时间或者包络信号对应的时间计算测井仪器中心与井壁的距离,则会引入一定的误差(系统时间误差)。基于此,为了屏蔽掉系统时间误差,通过在标准套管中进行测量,通过测量到的数据计算出系统时间误差。
图5示出了本申请另一实施例中标准套管与测井仪器的截面的示意图,大的圆形表示标准套管的截面,小的圆形表示测井仪器的截面,其中,标准套管为标准圆柱且其内径是已知的,以测井仪器的中心为原点建立直角坐标系,图中所示测井仪器具有3个超声波换能器,超声波换能器3与x轴的正轴的夹角为30度。
考虑到系统时间误差,i方位上测井仪器中心与标准套管的管壁的距离di的具体计算公式如下:
其中,V表示泥浆声速,ti表示对标准套管的i方位的反射回波测井数据进行分析所确定的波形峰值的时间,δt表示系统时间误差,R表示测井仪器的内径,右边第一项表示超声波换能器与管壁之间的距离。
根据上述三个方位的距离以及三个方位在坐标系中的角度,可以求出圆上三个点的坐标,具体为:
将上述三个坐标点代入到圆的求解方程,则可得到圆心及半径,具体计算公式如下:
其中,(x01,y01)表示圆心,r表示标准套管的内径。泥浆声速、标准套管对应的各个方位的波形峰值的时间、测井仪器的内径、标准套管的内径均是已知的,则可以基于上述方程进行求解,计算未知的系统时间误差。
则在该方式中,任一深度点处i方位上测井仪器中心与井壁之间的距离的具体计算公式如下:
其中,Di′表示i方位上测井仪器中心与井壁之间的距离,Ri表示波形峰值的时间,R表示测井仪器的半径,V表示泥浆声速,δt表示系统时间误差。该方式通过考虑系统时间误差,能够提高计算得到的测井仪器中心与井壁之间的距离的准确性。
在一种可选的实施方式中,根据该待处理深度点及其相邻的多个深度点处各个方位上测井仪器中心与井壁的距离以及各个方位的角度,进行椭圆拟合处理进一步包括:针对于待处理深度点及其相邻的多个深度点中的任一深度点,根据该深度点处任一方位上测井仪器与井壁的距离以及该方位的角度,计算对应在预设坐标系中的坐标值;根据各个坐标值进行椭圆拟合处理。由于在测井仪器在不同深度点的旋转角度不一致,因此,需要根据各个方位上测井仪器中心与井壁的距离与各个方位的角度转换为同一个坐标系下的坐标值,从而得到同一个椭圆上的多个点。
需要说明的是,待处理深度点及其相邻的多个深度点的数量与超声波换能器的数量之积不小于5,具体实施时,可将多个深度点的数量设置为较高,以便提高拟合的井眼椭圆的准确性。
在一种可选的实施方式中,采用最小二乘法进行椭圆拟合处理,具体拟合方法如下:
假设椭圆方程为:Ax2+Bxy+Cy2+Dx+Ey+1=0
设Pj(xj,yj)(j=1,2,……m)为椭圆轮廓上的测量点,m≥5,依据最小二乘原理,拟合的目标函数为:
欲使F最小,需要满足如下条件:
基于此,即可构建如下的方程:
通过解方程即可得到A、B、C、D以及E的值,进一步可计算得到椭圆的五个参数:位置参数(θ,x0,y0)以及形状参数(a,b),(x0,y0)表示椭圆圆心,θ表示偏心角,a表示椭圆长轴,b表示椭圆短轴。具体计算公式如下:
图6示出了本发明另一实施例井孔椭圆度计算结果的示意图,图中第一道61为椭圆长轴曲线,根据各个深度点处拟合得到的椭圆长轴的长度绘制得到,第二道62为椭圆短轴曲线,根据各个深度点处拟合得到的椭圆短轴的长度绘制得到,第三道63为椭圆度曲线,根据各个深度点处的椭圆度绘制得到,第四道64为多个深度点处的井孔二维截面图,采用本申请方法可以利用有限的数据得到井孔的椭圆度信息。
图7示出了本发明实施例提供的井孔椭圆度的计算装置的结构示意图,如图7所示,该装置包括:
获取模块71,适于获取深度区间内各个深度点处各个方位的反射回波测井数据;
分析模块72,适于根据各个深度点处各个方位的反射回波测井数据,计算各个深度点处的各个方位上测井仪器中心与井壁的距离;
拟合模块73,适于根据各个深度点中的待处理深度点及其相邻的多个深度点处的各个方位上测井仪器中心与井壁的距离以及各个方位的角度,进行椭圆拟合处理;
处理模块74,适于根据拟合结果确定待处理深度点的井孔椭圆度。
在一种可选的方式中,分析模块72进一步适于:
针对于任一深度点,对该深度点处各个方位的反射回波测井数据进行峰值检测,得到各个波形峰值的时间信息;根据各个波形峰值的时间信息、泥浆声速以及测井仪器的半径,计算得到各个方位上测井仪器中心与井壁的距离。
在一种可选的方式中,分析模块72进一步适于:
获取标准套管的反射回波测井数据;其中,标准套管的反射回波测井数据利用测井仪器对标准套管进行测量得到;根据标准套管的反射回波测井数据及内径,确定系统时间误差;
根据系统时间误差、各个波形峰值的时间信息、泥浆声速以及测井仪器的半径,计算得到各个方位上测井仪器中心与井壁的距离。
在一种可选的方式中,拟合模块73进一步适于:
针对于待处理深度点及其相邻的多个深度点中的任一深度点,根据该深度点处的任一方位上测井仪器中心与井壁的距离以及该方位的角度,计算对应在预设坐标系中的坐标值;根据各个坐标值进行椭圆拟合处理。
在一种可选的方式中,装置还包括:滤波模块,适于对各个深度点处各个方位的反射回波测井数据进行中值滤波处理;
分析模块72进一步适于:根据各个深度点处的各个方位的滤波处理后的反射回波测井数据,计算各个方位上测井仪器中心与井壁的距离。
在一种可选的方式中,滤波模块进一步适于:对于任一深度点,根据该深度点及其相邻的多个深度点处的任一方位的反射回波测井数据,对该深度点处的该方位的反射回波测井数据进行中值滤波处理。
在一种可选的方式中,滤波模块进一步适于:对于任一深度点,获取各个方位的任一时刻的各个波形幅值中的波形幅度中值;将各个方位的该时刻的各个波形幅值分别减去该时刻的波形幅度中值。
根据本申请实施例的井孔椭圆度的计算装置,获取深度区间内各个深度点处各个方位的反射回波测井数据;根据各个深度点处各个方位的反射回波测井数据,计算各个深度点处的各个方位上测井仪器中心与井壁的距离;根据各个深度点中的待处理深度点及其相邻的多个深度点处的各个方位上测井仪器中心与井壁的距离以及各个方位的角度,进行椭圆拟合处理;根据拟合结果确定待处理深度点的井孔椭圆度。通过上述方式,采用了相邻多个深度点的数据进行椭圆拟合,仅仅利用随钻超声井径测井仪器即可以获取到井孔椭圆度,并且能够提高拟合井孔椭圆度的精准性。
本发明实施例提供了一种非易失性计算机存储介质,所述计算机存储介质存储有至少一可执行指令,该计算机可执行指令可执行上述任意方法实施例中的井孔椭圆度的计算方法。
图8示出了本发明实施例提供的计算设备的结构示意图,本发明具体实施例并不对计算设备的具体实现做限定。
如图8所示,该计算设备可以包括:处理器(processor)802、通信接口(Communications Interface)804、存储器(memory)806、以及通信总线808。
其中:处理器802、通信接口804、以及存储器806通过通信总线808完成相互间的通信。通信接口804,用于与其它设备比如客户端或其它服务器等的网元通信。处理器802,用于执行程序810,具体可以执行上述用于计算设备的井孔椭圆度的计算方法实施例中的相关步骤。
具体地,程序810可以包括程序代码,该程序代码包括计算机操作指令。
处理器802可能是中央处理器CPU,或者是特定集成电路ASIC(ApplicationSpecific Integrated Circuit),或者是被配置成实施本发明实施例的一个或多个集成电路。计算设备包括的一个或多个处理器,可以是同一类型的处理器,如一个或多个CPU;也可以是不同类型的处理器,如一个或多个CPU以及一个或多个ASIC。
存储器806,用于存放程序810。存储器806可能包含高速RAM存储器,也可能还包括非易失性存储器(non-volatile memory),例如至少一个磁盘存储器。
在此提供的算法或显示不与任何特定计算机、虚拟系统或者其它设备固有相关。各种通用系统也可以与基于在此的示教一起使用。根据上面的描述,构造这类系统所要求的结构是显而易见的。此外,本发明实施例也不针对任何特定编程语言。应当明白,可以利用各种编程语言实现在此描述的本发明的内容,并且上面对特定语言所做的描述是为了披露本发明的最佳实施方式。
在此处所提供的说明书中,说明了大量具体细节。然而,能够理解,本发明的实施例可以在没有这些具体细节的情况下实践。在一些实例中,并未详细示出公知的方法、结构和技术,以便不模糊对本说明书的理解。
类似地,应当理解,为了精简本发明并帮助理解各个发明方面中的一个或多个,在上面对本发明的示例性实施例的描述中,本发明实施例的各个特征有时被一起分组到单个实施例、图、或者对其的描述中。然而,并不应将该公开的方法解释成反映如下意图:即所要求保护的本发明要求比在每个权利要求中所明确记载的特征更多的特征。更确切地说,如权利要求书所反映的那样,发明方面在于少于前面公开的单个实施例的所有特征。因此,遵循具体实施方式的权利要求书由此明确地并入该具体实施方式,其中每个权利要求本身都作为本发明的单独实施例。
本领域那些技术人员可以理解,可以对实施例中的设备中的模块进行自适应性地改变并且把它们设置在与该实施例不同的一个或多个设备中。可以把实施例中的模块或单元或组件组合成一个模块或单元或组件,以及此外可以把它们分成多个子模块或子单元或子组件。除了这样的特征和/或过程或者单元中的至少一些是相互排斥之外,可以采用任何组合对本说明书(包括伴随的权利要求、摘要和附图)中公开的所有特征以及如此公开的任何方法或者设备的所有过程或单元进行组合。除非另外明确陈述,本说明书(包括伴随的权利要求、摘要和附图)中公开的每个特征可以由提供相同、等同或相似目的的替代特征来代替。
此外,本领域的技术人员能够理解,尽管在此的一些实施例包括其它实施例中所包括的某些特征而不是其它特征,但是不同实施例的特征的组合意味着处于本发明的范围之内并且形成不同的实施例。例如,在权利要求书中,所要求保护的实施例的任意之一都可以以任意的组合方式来使用。
本发明的各个部件实施例可以以硬件实现,或者以在一个或者多个处理器上运行的软件模块实现,或者以它们的组合实现。本领域的技术人员应当理解,可以在实践中使用微处理器或者数字信号处理器(DSP)来实现根据本发明实施例的一些或者全部部件的一些或者全部功能。本发明还可以实现为用于执行这里所描述的方法的一部分或者全部的设备或者装置程序(例如,计算机程序和计算机程序产品)。这样的实现本发明的程序可以存储在计算机可读介质上,或者可以具有一个或者多个信号的形式。这样的信号可以从因特网网站上下载得到,或者在载体信号上提供,或者以任何其他形式提供。
应该注意的是上述实施例对本发明进行说明而不是对本发明进行限制,并且本领域技术人员在不脱离所附权利要求的范围的情况下可设计出替换实施例。在权利要求中,不应将位于括号之间的任何参考符号构造成对权利要求的限制。单词“包含”不排除存在未列在权利要求中的元件或步骤。位于元件之前的单词“一”或“一个”不排除存在多个这样的元件。本发明可以借助于包括有若干不同元件的硬件以及借助于适当编程的计算机来实现。在列举了若干装置的单元权利要求中,这些装置中的若干个可以是通过同一个硬件项来具体体现。单词第一、第二、以及第三等的使用不表示任何顺序。可将这些单词解释为名称。上述实施例中的步骤,除有特殊说明外,不应理解为对执行顺序的限定。
Claims (10)
1.一种井孔椭圆度的计算方法,其特征在于,所述方法包括:
获取深度区间内各个深度点处各个方位的反射回波测井数据;
根据各个深度点处各个方位的反射回波测井数据,计算所述各个深度点处各个方位上测井仪器中心与井壁的距离;
根据所述各个深度点中的待处理深度点及其相邻的多个深度点处各个方位上测井仪器中心与井壁的距离以及所述各个方位的角度,进行椭圆拟合处理;
根据拟合结果确定所述待处理深度点的井孔椭圆度。
2.根据权利要求1所述的方法,其特征在于,所述根据各个深度点处各个方位的反射回波测井数据,计算所述各个深度点处各个方位上测井仪器中心与井壁的距离进一步包括:
针对于任一深度点,对该深度点处各个方位的反射回波测井数据进行峰值检测,得到各个波形峰值的时间信息;
根据所述各个波形峰值的时间信息、泥浆声速以及所述测井仪器的半径,计算得到各个方位上测井仪器中心与井壁的距离。
3.根据权利要求2所述的方法,其特征在于,所述方法进一步包括:
利用所述测井仪器对标准套管进行测量,获得所述标准套管的反射回波测井数据;
根据所述标准套管的反射回波测井数据及内径,确定系统时间误差;
所述根据所述各个波形峰值的时间信息、泥浆声速以及所述测井仪器的半径,计算得到各个方位上测井仪器中心与井壁的距离进一步包括:
根据所述系统时间误差、所述各个波形峰值的时间信息、泥浆声速以及所述测井仪器的半径,计算得到各个方位上测井仪器中心与井壁的距离。
4.根据权利要求1-3中任一项所述的方法,其特征在于,所述根据所述待处理深度点及其相邻的多个深度点处各个方位上测井仪器中心与井壁的距离以及所述各个方位的角度,进行椭圆拟合处理进一步包括:
针对于待处理深度点及其相邻的多个深度点中的任一深度点,根据该深度点处任一方位上测井仪器中心与井壁的距离以及该方位的角度,计算对应在预设坐标系中的坐标值;
根据各个坐标值进行椭圆拟合处理。
5.根据权利要求1所述的方法,其特征在于,所述获取深度区间内各个深度点处各个方位的反射回波测井数据之后,所述方法还包括:
对各个深度点处各个方位的反射回波测井数据进行中值滤波处理;
所述根据各个深度点处各个方位的反射回波测井数据,计算所述各个深度点处各个方位上测井仪器中心与井壁的距离进一步包括:
根据各个深度点处的各个方位的滤波处理后的反射回波测井数据,计算各个深度点处各个方位上测井仪器中心与井壁的距离。
6.根据权利要求5所述的方法,其特征在于,所述对各个深度点处各个方位的反射回波测井数据进行中值滤波处理进一步包括:
对于任一深度点,根据该深度点及其相邻的多个深度点处的任一方位的反射回波测井数据,对该深度点处的该方位的反射回波测井数据进行中值滤波处理。
7.根据权利要求5所述的方法,其特征在于,所述对各个深度点处各个方位的反射回波测井数据进行中值滤波处理进一步包括:
对于任一深度点,获取各个方位的任一时刻的各个波形幅值中的波形幅度中值;将各个方位的该时刻的各个波形幅值分别减去该时刻的波形幅度中值。
8.一种井孔椭圆度的计算装置,其特征在于,所述装置包括:
获取模块,适于获取深度区间内各个深度点处各个方位的反射回波测井数据;
分析模块,适于根据各个深度点处各个方位的反射回波测井数据,计算所述各个深度点处各个方位上测井仪器中心与井壁的距离;
拟合模块,适于根据所述各个深度点中的待处理深度点及其相邻的多个深度点处各个方位上测井仪器中心与井壁的距离以及所述各个方位的角度,进行椭圆拟合处理;
处理模块,适于根据拟合结果确定所述待处理深度点的井孔椭圆度。
9.一种计算设备,包括:处理器、存储器、通信接口和通信总线,所述处理器、所述存储器和所述通信接口通过所述通信总线完成相互间的通信;
所述存储器用于存放至少一可执行指令,所述可执行指令使所述处理器执行如权利要求1-7中任一项所述的井孔椭圆度的计算方法对应的操作。
10.一种计算机存储介质,所述存储介质中存储有至少一可执行指令,所述可执行指令使处理器执行如权利要求1-7中任一项所述的井孔椭圆度的计算方法对应的操作。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310507713.6A CN116556925A (zh) | 2023-05-06 | 2023-05-06 | 井孔椭圆度的计算方法、装置、计算设备及存储介质 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310507713.6A CN116556925A (zh) | 2023-05-06 | 2023-05-06 | 井孔椭圆度的计算方法、装置、计算设备及存储介质 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN116556925A true CN116556925A (zh) | 2023-08-08 |
Family
ID=87497628
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310507713.6A Pending CN116556925A (zh) | 2023-05-06 | 2023-05-06 | 井孔椭圆度的计算方法、装置、计算设备及存储介质 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116556925A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117910883A (zh) * | 2024-03-18 | 2024-04-19 | 中海油田服务股份有限公司 | 测井数据质量检测方法、装置、计算设备及存储介质 |
-
2023
- 2023-05-06 CN CN202310507713.6A patent/CN116556925A/zh active Pending
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117910883A (zh) * | 2024-03-18 | 2024-04-19 | 中海油田服务股份有限公司 | 测井数据质量检测方法、装置、计算设备及存储介质 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US9322807B2 (en) | Ultrasonic signal time-frequency decomposition for borehole evaluation or pipeline inspection | |
AU2011382521B2 (en) | Acoustic transducer apparatus, systems, and methods | |
US9341731B2 (en) | Fracture detection using monopole source in acoustic well logging | |
RU2528279C1 (ru) | Наложение форм акустических сигналов с использованием группирования по азимутальным углам и/или отклонениям каротажного зонда | |
CN115992692B (zh) | 水泥环厚度测量方法、装置、电子设备及存储介质 | |
CN115144906B (zh) | 超声波测井的幅度偏心校正方法及装置 | |
CN116556925A (zh) | 井孔椭圆度的计算方法、装置、计算设备及存储介质 | |
CN115822564B (zh) | 套管厚度测量方法、装置、电子设备及计算机存储介质 | |
CN105074127A (zh) | 从声波波形中移除衰荡效应 | |
CN115324564B (zh) | 固井质量检测方法、装置、计算设备及存储介质 | |
CN106548456B (zh) | 校正超声图像轮廓中的偏心的方法 | |
CN115992691A (zh) | 基于超声兰姆波的固井质量检测方法及装置 | |
US9753170B2 (en) | Sourceless density determination apparatus, methods, and systems | |
CN106930758B (zh) | 一种随钻声波测井方法 | |
CN106767583B (zh) | 用于基桩检测声波透射法的纵向剖面等效桩径计算方法 | |
Mandal et al. | A new monocable circumferential acoustic scanner tool (CAST-M) for cased-hole and openhole applications | |
CN110824007B (zh) | 一种管桩裂缝检测方法及系统 | |
CN112904348B (zh) | 一种三维探测方法、装置、设备和存储介质 | |
CN116201526B (zh) | 微环隙检测方法、装置、计算设备及存储介质 | |
CN114139584A (zh) | 一种随钻声波测井套管波衰减反演方法和装置 | |
CN114135274A (zh) | 一种随钻方位固井质量扫描方法和装置 | |
CN116241239B (zh) | 基于远近单极的固井评价的方法、装置、设备及存储介质 | |
CN116660997B (zh) | 套管内外介质声阻抗反演方法、装置及电子设备 | |
CN111691876A (zh) | 一种利用声波测井对邻井成像的方法、装置及存储介质 | |
CN114779347A (zh) | 利用超声回波测井方法、装置、设备和存储介质 |
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 |