CN115144906A - 超声波测井的幅度偏心校正方法及装置 - Google Patents
超声波测井的幅度偏心校正方法及装置 Download PDFInfo
- Publication number
- CN115144906A CN115144906A CN202211086231.XA CN202211086231A CN115144906A CN 115144906 A CN115144906 A CN 115144906A CN 202211086231 A CN202211086231 A CN 202211086231A CN 115144906 A CN115144906 A CN 115144906A
- Authority
- CN
- China
- Prior art keywords
- amplitude
- depth
- correction
- reflection echo
- preset
- 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.)
- Granted
Links
- 238000012937 correction Methods 0.000 title claims abstract description 122
- 238000000034 method Methods 0.000 title claims abstract description 44
- 238000003384 imaging method Methods 0.000 claims abstract description 34
- 238000004891 communication Methods 0.000 claims description 15
- 238000004364 calculation method Methods 0.000 claims description 10
- 239000002002 slurry Substances 0.000 claims description 10
- 238000000605 extraction Methods 0.000 claims description 6
- 238000010586 diagram Methods 0.000 description 11
- 229910000831 Steel Inorganic materials 0.000 description 3
- 239000010959 steel Substances 0.000 description 3
- 230000015572 biosynthetic process Effects 0.000 description 2
- 238000004590 computer program Methods 0.000 description 2
- 238000002592 echocardiography Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 239000012530 fluid Substances 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 230000001902 propagating effect Effects 0.000 description 2
- 239000000523 sample Substances 0.000 description 2
- 230000002238 attenuated effect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000011022 operating instruction Methods 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/40—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
- G01V1/44—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging using generators and receivers in the same well
- G01V1/48—Processing data
- G01V1/50—Analysing data
-
- 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
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Geology (AREA)
- Geophysics (AREA)
- Environmental & Geological Engineering (AREA)
- Mining & Mineral Resources (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- General Physics & Mathematics (AREA)
- Fluid Mechanics (AREA)
- Remote Sensing (AREA)
- Geochemistry & Mineralogy (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
本发明公开了一种超声波测井的幅度偏心校正方法及装置,方法包括:获取测井仪器在预设深度区间产生的测井数据;根据泥浆脉冲反射回波数据和到时数据,提取得到一次反射回波的第一波峰值的信号幅度和第一到时数据,以及二次反射回波的第二波峰值的信号幅度和第二到时数据,计算得到泥浆衰减因子;针对预设深度区域中任一深度点,获取该深度点上下预设窗长范围的井眼内壁反射回波幅度,计算得到该深度点的深度校正因子;根据深度校正因子对该深度点的井眼内壁反射回波幅度进行幅度校正;对预设深度区域各个深度点的井眼内壁反射回波幅度均完成幅度校正,得到幅度校正后的超声波成像曲线,从而解决测井仪器偏心导致的成像曲线图像明暗不均匀问题。
Description
技术领域
本发明实施例涉及地球物理勘探和石油测井领域,具体涉及一种超声波测井的幅度偏心校正方法及装置。
背景技术
超声波成像测井广泛应用于油气田勘探开发,利用波探头发射超声脉冲信号,超声波信号在井眼流体中传播,遇到井壁进行反射,反射回波信号被超声波探头接收。通过测量反射回波信号的幅度和到时可以评价井眼内壁的地质构造信息。测井仪器在测量过程中通常安装扶正器帮助居中,但在运动过程中不可避免会偏离井孔轴线。此时虽然井孔是圆形的,但由于探测器旋转轴不在井孔轴线,声波到达井壁并反射的路径不等,导致超声波到时成像出现明暗条纹;另外由于传播路径有差异,超声波信号衰减不同,导致超声波幅度成像不均匀。
如图1中所示,左侧为测井仪器居中测量时的超声波传播路径截面图,右侧为测井仪器偏心时超声波传播路径截面图。若测井仪器位于井眼中心,从换能器辐射的超声波信号入射到井壁上,会沿着原来的路径返回到换能器并被接收;若测井仪器偏心,则从换能器辐射的超声波信号入射到井壁上,超声波反射回来的路径与入射路径不一致,导致反射回来的超声波信号幅度降低,到时增加,测量得到的幅度与到时成像不能真实地反映地层的真实构造,对解释评价带来影响。测井仪器偏心时,扫描一周测量的超声反射回波信号,得到如图2所示的展示效果。其中,粗圆圈表示超声反射回波的到时,基于图2可见,圆圈的中心没有与井眼中心(圆点)重合,此时测井仪器偏心,测量得到的反射波幅度和到时不能真实反映地层的构造。
另外由于泥浆的衰减,不同距离下声波能量的损失也不相同,最终在幅度和旅时图像会出现两条竖直亮暗带,测井仪器偏心时测量到的幅度成像或到时成像可能掩盖一些有用的地质构造信息。
因此,急需一种超声波测井的幅度偏心校正方法,以解决由于测井仪器偏心导致的幅度与到时成像不能真实地反映地层的真实构造、反射波幅度成像图像的明暗不均匀等问题。
发明内容
鉴于上述问题,提出了本发明实施例以便提供一种克服上述问题或者至少部分地解决上述问题的超声波测井的幅度偏心校正方法及装置。
根据本发明实施例的一个方面,提供了一种超声波测井的幅度偏心校正方法,方法包括:
获取步骤,获取测井仪器在预设深度区间产生的测井数据;测井数据包括泥浆脉冲反射回波数据、井眼内壁反射回波幅度以及到时数据;
提取计算步骤,根据泥浆脉冲反射回波数据,提取得到一次反射回波的第一波峰值的信号幅度和第一到时数据,以及二次反射回波的第二波峰值的信号幅度和第二到时数据,计算得到泥浆衰减因子;
校正步骤,针对预设深度区域中任一深度点,获取该深度点上下预设窗长范围的井眼内壁反射回波幅度,计算得到该深度点的深度校正因子;根据深度校正因子对该深度点的井眼内壁反射回波幅度进行幅度校正;重复执行校正步骤,直至预设深度区域各个深度点的井眼内壁反射回波幅度完成幅度校正,得到幅度校正后的超声波成像曲线。
根据本发明实施例的另一方面,提供了一种超声波测井的幅度偏心校正装置,其包括:
获取模块,适于获取测井仪器在预设深度区间产生的测井数据;测井数据包括泥浆脉冲反射回波数据、井眼内壁反射回波幅度以及到时数据;
提取计算模块,适于根据泥浆脉冲反射回波数据,提取得到一次反射回波的第一波峰值的信号幅度和第一到时数据,以及二次反射回波的第二波峰值的信号幅度和第二到时数据,计算得到泥浆衰减因子;
校正模块,适于针对预设深度区域中任一深度点,获取该深度点上下预设窗长范围的井眼内壁反射回波幅度,计算得到该深度点的深度校正因子;根据深度校正因子对该深度点的井眼内壁反射回波幅度进行幅度校正;重复执行校正模块,直至预设深度区域各个深度点的井眼内壁反射回波幅度完成幅度校正,得到幅度校正后的超声波成像曲线。
根据本发明实施例的又一方面,提供了一种计算设备,包括:处理器、存储器、通信接口和通信总线,所述处理器、所述存储器和所述通信接口通过所述通信总线完成相互间的通信;
所述存储器用于存放至少一可执行指令,所述可执行指令使所述处理器执行上述超声波测井的幅度偏心校正方法对应的操作。
根据本发明实施例的再一方面,提供了一种计算机存储介质,所述存储介质中存储有至少一可执行指令,所述可执行指令使处理器执行如上述超声波测井的幅度偏心校正方法对应的操作。
根据本发明实施例的提供的超声波测井的幅度偏心校正方法及装置,通过获取测井仪器在预设深度区间的测井数据,计算泥浆衰减因子,对井眼内壁反射回波幅度进行偏心校正,得到整体偏心校正后的预设深度区间的超声波幅度成像曲线,从而解决测井仪器偏心导致的成像曲线图像明暗不均匀问题。
上述说明仅是本发明实施例技术方案的概述,为了能够更清楚了解本发明实施例的技术手段,而可依照说明书的内容予以实施,并且为了让本发明实施例的上述和其它目的、特征和优点能够更明显易懂,以下特举本发明实施例的具体实施方式。
附图说明
通过阅读下文优选实施方式的详细描述,各种其他的优点和益处对于本领域普通技术人员将变得清楚明了。附图仅用于示出优选实施方式的目的,而并不认为是对本发明实施例的限制。而且在整个附图中,用相同的参考符号表示相同的部件。在附图中:
图1示出了测井仪器居中以及偏心情况下超声波传播路径示意图;
图2示出了测井仪器偏心情况下超声反射回波到时的效果示意图;
图3示出了根据本发明一个实施例的超声波测井的幅度偏心校正方法的流程图;
图4示出了测井仪器设置在井眼中的示意图;
图5示出了超声波包络曲线示意图;
图6示出了原始反射波幅度与校正后幅度对比示意图;
图7示出了反射波幅度偏心校正后成果示意图;
图8示出了根据本发明一个实施例的超声波测井的幅度偏心校正装置的结构示意图;
图9示出了根据本发明一个实施例的一种计算设备的结构示意图。
具体实施方式
下面将参照附图更详细地描述本发明的示例性实施例。虽然附图中显示了本发明的示例性实施例,然而应当理解,可以以各种形式实现本发明而不应被这里阐述的实施例所限制。相反,提供这些实施例是为了能够更透彻地理解本发明,并且能够将本发明的范围完整的传达给本领域的技术人员。
图3示出了根据本发明一个实施例的超声波测井的幅度偏心校正方法的流程图,如图3所示,该方法包括如下步骤:
步骤S301,获取测井仪器在预设深度区间产生的测井数据。
本实施例中测井仪器可以采用如超声波成像测井仪,在测井时如图4所示,超声成像测井仪402位于充满泥浆的井眼401中,井眼401外为地层。超声成像测井仪402设置有两个超声波换能器,其中一个超声波换能器为超声换能器405,设置与井壁平行,不断向井壁发射超声波信号,通过测量井壁表面的超声波信号的幅度或者到时成像,用于评价井眼内壁的地质构造类型,如裂缝、层界面、孔洞等。另一个超声波换能器为泥浆声速计404,设置于井壁侧,向正上方发射超声波,在泥浆声速计404的正上方距离如L位置处固定设置反射块403,反射块403可以采用如钢板,超声波入射到反射块403后反射回来被泥浆声速计404接收。以上为举例说明,具体实施时可以根据实施情况设置,此处不做限定。
预设深度区间可以根据不同测井的实际情况设置,将测井仪器放置如井眼的预设深度区间,获取到该预设深度区间内产生的测井数据,测井数据包括如泥浆脉冲反射回波数据、井眼内壁反射回波幅度以及到时数据。如利用测井仪器中的泥浆声速计、超声换能器测量获取泥浆脉冲反射回波数据、井眼内壁反射回波幅度以及到时数据等。
步骤S302,根据泥浆脉冲反射回波数据,提取得到一次反射回波的第一波峰值的信号幅度和第一到时数据,以及二次反射回波的第二波峰值的信号幅度和第二到时数据,计算得到泥浆衰减因子。
利用泥浆声速计测量得到的泥浆脉冲反射回波数据,对泥浆脉冲反射回波数据进行如Hibert变换,可以得到对应的包络曲线,如图5所示,粗线条为包络曲线,细线条为泥浆脉冲信号。横轴为时间轴,单位为微秒,纵轴为幅度轴,单位为毫伏。根据包络曲线可以提取一次及二次超声脉冲反射回波信号,如图5中箭头所指向的第一次反射回峰波和第二次反射回峰波,对应的提取得到一次反射回波的第一波峰值的信号幅度和第一到时数据,以及二次反射回波的第二波峰值的信号幅度和第二到时数据。第一到时数据即图5中第一次反射回峰波对应的到时时间,第二到时数据即图5中第二次反射回峰波对应的到时时间。
计算确定泥浆脉冲信号入射到反射块的传播时间时,可以根据第二到时数据T2和第一到时数据T1的差值,计算确定传播时间T,如传播时间T=(T2-T1)/2。测井仪器与反射块的距离L即泥浆声速计距离反射块的距离L,距离L为设置测井仪器时可以获知。根据传播时间T以及距离L可以计算得到超声波在泥浆中的传播速度V,传播速度V=L/T。
超声波发射超声波信号入射到障碍物,如井壁后会再反射回来,发射的超声波信号幅度根据其在泥浆中的传播距离会产生衰减,反射后接收的超声波信号幅度与发射的超声波信号幅度、传播距离等存在关联关系。具体的,如A0为超声波发射时的信号幅度,P0为超声波发射后在泥浆中传播x距离后到达井壁时的信号幅度,P0=A0*e-ax。此处,a为衰减因子。Pr为P0从井壁反射时的信号幅度。超声波在井壁发生反射时在界面处存在能量损失,P0与Pr之间的关系为R=Pr/P0,R为声压反射系数,得到Pr=A0*R *e-ax。A1为接收到的超声波信号幅度(也即第一波峰值的信号幅度),为Pr又传播x距离后接收到的信号幅度,A1=Pr*e-ax=A0*R *e-2ax。即最终超声换能器接收所接收到的超声波信号不仅有传播路径产生的损失,还有井壁反射所产生的损失。声压反射系数R与反射界面相关,每个井眼中泥浆密度ρ,可以根据钻井液确定。根据传播速度V以及泥浆密度ρ,得到泥浆的声阻抗Z1,Z1=ρV。反射块的声阻抗Z2可以根据反射块的材质确定,如反射块采用钢板,反射块的声阻抗Z2可以采用钢的声阻抗,如47.6Mrayls。根据泥浆的声阻抗Z1和反射块的声阻抗Z2,可以计算得到声压反射系数R。R=(Z2-Z1)/(Z2+Z1)。
在本实施例中,由于A0在实际应用中较难准确获取,而超声波入射到反射块后会继续发生二次反射,再继续被超声换能器接收。因此,可以利用一次和二次反射回波,即第一波峰值的信号幅度、第二波峰值的信号幅度来确定泥浆衰减因子。第一波峰值的信号幅度A1和第二波峰值的信号幅度A2的关系为A2=A1*R *e-2ax。此处,考虑到距离在实际应用时存在测量不精准等问题,距离与传播时间T成正比,传播时间T可以通过测井仪器准确获知。衰减是基于泥浆这一介质引起,可以采用泥浆衰减因子μ,基于上述传播时间、泥浆衰减因子对A2的计算公式进行变形,得到A2=A1*R *e-T/μ。
根据传播时间T、第一波峰值的信号幅度A1、第二波峰值的信号幅度A2、传播速度V可以确定泥浆衰减因子。具体的,第一波峰值的信号幅度、第二波峰值的信号幅度可以从包络曲线中提取得到,如第一波峰值的信号幅度为A1,第二波峰值的信号幅度为A2。根据传播时间T、第一波峰值的信号幅度A1、第二波峰值的信号幅度A2以及声压反射系数R可以计算得到泥浆衰减因子μ。μ=-T/ln(A2/(A1*R))。
步骤S303,针对预设深度区域中任一深度点,获取该深度点上下预设窗长范围的井眼内壁反射回波幅度,计算得到该深度点的深度校正因子;根据深度校正因子对该深度点的井眼内壁反射回波幅度进行幅度校正。
针对预设深度区域中任一深度点,根据该深度点上下预设窗长范围的井眼内壁反射回波幅度,计算得到该深度点的深度校正因子。其中,不同深度点具有不同的深度校正因子,每个深度点的深度校正因子根据各自深度点上下预设窗长范围的井眼内壁反射回波幅度计算得到。具体的,先获取该深度点上下预设窗长范围内多个深度点在不同方位(如井眼360度的各个不同方位)的多个采集点的到时数据,如深度点为3887m,预设窗长范围为4m,获取深度点为3883m-3891m在不同方位的多个采集点的到时数据。根据多个深度点在不同方位的多个采集点的到时数据以及泥浆衰减因子计算得到该深度点的深度校正因子Correction。
其中,公式(1)中M为进行幅度校正的深度点,Time(θ,M)为方位θ,深度点M的反射波的到时数据,μ为泥浆衰减因子。N为测井仪器周向不同方位对应的N个采集点数,ndep为该深度点在进行幅度校正时预设窗长范围内的ndep个采集点数。其中,基于以上公式计算得到的是深度点M的深度校正因子。当深度点更改,需要基于以上公式重新计算对应更改后的深度点的深度校正因子。
计算得到深度校正因子Correction后,根据该深度点的深度校正因子Correction、泥浆衰减因子μ以及该深度点在预设方位的到时数据,对该深度点在预设方位的井眼内壁反射回波幅度进行幅度校正,得到该深度点在预设方位校正后的井眼内壁反射回波幅度,如利用以下公式进行幅度校正:
Ampc(θ)=Amp(θ)*e(Time(θ)/μ)/Correction(2)
其中,公式(2)中Ampc(θ)为方位θ在进行幅度校正后幅度值,Amp(θ)为方位θ基于测井仪器获取的反射波的幅度值。Time(θ)为方位θ的反射波的到时数据。
如图6所示,为深度点3886m的井眼内壁反射回波幅度与校正后幅度的示意图,其中,实线为使用测井仪器进行测井的原始幅度,虚线为校正后幅度。横轴代表不同方位的角度,单位为度;纵轴代表各个幅度值,单位为毫伏。图6中对比了深度点3886m处的井眼内壁反射回波幅度校正前后的结果,从图6中可见,校正前不同角度测量的幅度值变化程度较大,而校正后不同角度幅度值比较稳定。
步骤S304,判断预设深度区域内各个深度点的井眼内壁反射回波幅度是否均进行幅度校正。
对一个深度点完成幅度校正后,继续针对下一深度点进行幅度校正,如对3887m深度点进行幅度校正后,在继续对3888m深度点进行幅度校正,不同深度点使用不同深度校正因子进行幅度校正。判断预设深度区域内的各个深度点的井眼内壁反射回波幅度是否均进行幅度校正,若否,则重复执行步骤S303,直至完成预设深度区域各个深度点的井眼内壁反射回波幅度进行幅度校正;若是,则执行步骤S305,基于完成幅度校正的各个深度点的井眼内壁反射回波幅度还原超声波成像曲线。
具体的,如在对深度点进行幅度校正时,按照预设深度区域的各个深度点大小从小到大排序,从预设深度区域的最小深度开始进行幅度校正,在对深度点完成幅度校正,判断是否预设深度区域的最大深度,若否,则获取该深度点的下一个深度点继续进行幅度校正,直至完成预设深度区域各个深度点的井眼内壁反射回波幅度进行幅度校正;或者从预设深度区域的最大深度开始进行幅度校正,在对深度点完成幅度校正,判断是否预设深度区域的最小深度,若否,则获取该深度点的前一个深度点继续进行幅度校正,直至完成预设深度区域各个深度点的井眼内壁反射回波幅度进行幅度校正。以上为举例说明,具体根据实施情况设置,此处不做限定。
步骤S305,根据校正后的预设深度区域各个深度点的井眼内壁反射回波幅度得到幅度校正后的超声波成像曲线。
根据校正后的预设深度区域内各个深度点的井眼内壁反射回波幅度对应的生成幅度校正后的超声波成像曲线,如图7所示。图7中第一列为预设深度区域内的各个深度点,如从3883m-3891m,第二列为泥浆声速计测量的超声波信号曲线轨迹,第三列为测井仪器的原始数据得到的原始反射波幅度成像曲线,第四列为反射波到时成像曲线,第五列为基于校正后反射幅度成像曲线。由第一道超声波信号的波形曲线轨迹可见,波形中有两个波包,依次为超声波在反射块处的一次和二次反射回波信号。由第三列和第五列的幅度成像曲线对比可见,由于测井仪器偏心导致第三列原始反射波幅度成像时强时弱的问题,在第五列校正后反射幅度成像曲线均有明显的改善,幅度图像均匀,更有利用获取地层的地质构造信息。
根据本发明实施例提供的超声波测井的幅度偏心校正方法,通过获取测井仪器在预设深度区间的测井数据,计算泥浆衰减因子,对井眼内壁反射回波幅度进行偏心校正,得到整体偏心校正后的预设深度区间的超声波幅度成像曲线,从而解决测井仪器偏心导致的成像曲线图像明暗不均匀问题。
图8示出了本发明实施例提供的超声波测井的幅度偏心校正装置的结构示意图。如图8所示,该装置包括:
获取模块810,适于获取测井仪器在预设深度区间产生的测井数据;测井数据包括泥浆脉冲反射回波数据、井眼内壁反射回波幅度以及到时数据;
提取计算模块820,适于根据泥浆脉冲反射回波数据,提取得到一次反射回波的第一波峰值的信号幅度和第一到时数据,以及二次反射回波的第二波峰值的信号幅度和第二到时数据,计算得到泥浆衰减因子;
校正模块830,适于针对预设深度区域中任一深度点,获取该深度点上下预设窗长范围的井眼内壁反射回波幅度,计算得到该深度点的深度校正因子;根据深度校正因子对该深度点的井眼内壁反射回波幅度进行幅度校正;重复执行校正模块830,直至预设深度区域各个深度点的井眼内壁反射回波幅度完成幅度校正,得到幅度校正后的超声波成像曲线。
可选地,提取计算模块820进一步适于:
计算确定泥浆脉冲信号入射到反射块的传播时间;传播时间根据第一到时数据和第二到时数据确定;
根据传播时间以及测井仪器与反射块的距离,计算得到超声波在泥浆中的传播速度;
根据传播时间、第一波峰值的信号幅度、第二波峰值的信号幅度、传播速度计算确定泥浆衰减因子。
可选地,提取计算模块820进一步适于:
根据第二到时数据和第一到时数据的差值,计算确定传播时间。
可选地,提取计算模块820进一步适于:
根据传播速度以及泥浆密度,得到泥浆的声阻抗;
根据泥浆的声阻抗和反射块的声阻抗,计算得到声压反射系数;
根据传播时间、第一波峰值的信号幅度、第二波峰值的信号幅度以及声压反射系数计算得到泥浆衰减因子。
可选地,校正模块830进一步适于:
针对预设深度区域中任一深度点,获取该深度点上下预设窗长范围内多个深度点在不同方位的多个采集点的到时数据;
根据多个深度点在不同方位的多个采集点的到时数据以及泥浆衰减因子计算得到该深度点的深度校正因子;
根据该深度点的深度校正因子、泥浆衰减因子以及该深度点在预设方位的到时数据,对该深度点在预设方位的井眼内壁反射回波幅度进行幅度校正,得到该深度点在预设方位校正后的井眼内壁反射回波幅度。
可选地,校正模块830进一步适于:
判断预设深度区域内各个深度点的井眼内壁反射回波幅度是否均进行幅度校正;
若否,重复校正模块,直至完成预设深度区域各个深度点的井眼内壁反射回波幅度进行幅度校正;
若是,根据校正后的预设深度区域各个深度点的井眼内壁反射回波幅度得到幅度校正后的超声波成像曲线。
可选地,不同深度点具有不同的深度校正因子。
以上各模块的描述参照方法实施例中对应的描述,在此不再赘述。
本发明实施例还提供了一种非易失性计算机存储介质,计算机存储介质存储有至少一可执行指令,可执行指令可执行上述任意方法实施例中的超声波测井的幅度偏心校正方法。
图9示出了根据本发明实施例的一种计算设备的结构示意图,本发明实施例的具体实施例并不对计算设备的具体实现做限定。
如图9所示,该计算设备可以包括:处理器902、通信接口904、存储器906、以及通信总线908。
处理器902、通信接口904、以及存储器906通过通信总线908完成相互间的通信。
通信接口904,用于与其它设备比如客户端或其它服务器等的网元通信。
处理器902,用于执行程序910,具体可以执行上述超声波测井的幅度偏心校正方法实施例中的相关步骤。
具体地,程序910可以包括程序代码,该程序代码包括计算机操作指令。
处理器902可能是中央处理器CPU,或者是特定集成电路ASIC(ApplicationSpecific Integrated Circuit),或者是被配置成实施本发明实施例的一个或多个集成电路。计算设备包括的一个或多个处理器,可以是同一类型的处理器,如一个或多个CPU;也可以是不同类型的处理器,如一个或多个CPU以及一个或多个ASIC。
存储器906,用于存放程序910。存储器906可能包含高速RAM存储器,也可能还包括非易失性存储器(non-volatile memory),例如至少一个磁盘存储器。
程序910具体可以用于使得处理器902执行上述任意方法实施例中的超声波测井的幅度偏心校正方法。程序910中各步骤的具体实现可以参见上述超声波测井的幅度偏心校正实施例中的相应步骤和单元中对应的描述,在此不赘述。所属领域的技术人员可以清楚地了解到,为描述的方便和简洁,上述描述的设备和模块的具体工作过程,可以参考前述方法实施例中的对应过程描述,在此不再赘述。
在此提供的算法或显示不与任何特定计算机、虚拟系统或者其它设备固有相关。各种通用系统也可以与基于在此的示教一起使用。根据上面的描述,构造这类系统所要求的结构是显而易见的。此外,本发明实施例也不针对任何特定编程语言。应当明白,可以利用各种编程语言实现在此描述的本发明实施例的内容,并且上面对特定语言所做的描述是为了披露本发明实施例的较佳实施方式。
在此处所提供的说明书中,说明了大量具体细节。然而,能够理解,本发明的实施例可以在没有这些具体细节的情况下实践。在一些实例中,并未详细示出公知的方法、结构和技术,以便不模糊对本说明书的理解。
类似地,应当理解,为了精简本发明实施例并帮助理解各个发明方面中的一个或多个,在上面对本发明的示例性实施例的描述中,本发明实施例的各个特征有时被一起分组到单个实施例、图、或者对其的描述中。然而,并不应将该公开的方法解释成反映如下意图:即所要求保护的本发明实施例要求比在每个权利要求中所明确记载的特征更多的特征。更确切地说,如下面的权利要求书所反映的那样,发明方面在于少于前面公开的单个实施例的所有特征。因此,遵循具体实施方式的权利要求书由此明确地并入该具体实施方式,其特征在于每个权利要求本身都作为本发明的单独实施例。
本领域那些技术人员可以理解,可以对实施例中的设备中的模块进行自适应性地改变并且把它们设置在与该实施例不同的一个或多个设备中。可以把实施例中的模块或单元或组件组合成一个模块或单元或组件,以及此外可以把它们分成多个子模块或子单元或子组件。除了这样的特征和/或过程或者单元中的至少一些是相互排斥之外,可以采用任何组合对本说明书(包括伴随的权利要求、摘要和附图)中公开的所有特征以及如此公开的任何方法或者设备的所有过程或单元进行组合。除非另外明确陈述,本说明书(包括伴随的权利要求、摘要和附图)中公开的每个特征可以由提供相同、等同或相似目的的替代特征来代替。
此外,本领域的技术人员能够理解,尽管在此的一些实施例包括其它实施例中所包括的某些特征而不是其它特征,但是不同实施例的特征的组合意味着处于本发明的范围之内并且形成不同的实施例。例如,在下面的权利要求书中,所要求保护的实施例的任意之一都可以以任意的组合方式来使用。
本发明的各个部件实施例可以以硬件实现,或者以在一个或者多个处理器上运行的软件模块实现,或者以它们的组合实现。本领域的技术人员应当理解,可以在实践中使用微处理器或者数字信号处理器(DSP)来实现根据本发明实施例的一些或者全部部件的一些或者全部功能。本发明实施例还可以实现为用于执行这里所描述的方法的一部分或者全部的设备或者系统程序(例如,计算机程序和计算机程序产品)。这样的实现本发明实施例的程序可以存储在计算机可读介质上,或者可以具有一个或者多个信号的形式。这样的信号可以从因特网网站上下载得到,或者在载体信号上提供,或者以任何其他形式提供。
应该注意的是上述实施例对本发明实施例进行说明而不是对本发明进行限制,并且本领域技术人员在不脱离所附权利要求的范围的情况下可设计出替换实施例。在权利要求中,不应将位于括号之间的任何参考符号构造成对权利要求的限制。单词“包含”不排除存在未列在权利要求中的元件或步骤。位于元件之前的单词“一”或“一个”不排除存在多个这样的元件。本发明实施例可以借助于包括有若干不同元件的硬件以及借助于适当编程的计算机来实现。在列举了若干系统的单元权利要求中,这些系统中的若干个可以是通过同一个硬件项来具体体现。单词第一、第二、以及第三等的使用不表示任何顺序。可将这些单词解释为名称。上述实施例中的步骤,除有特殊说明外,不应理解为对执行顺序的限定。
Claims (10)
1.一种超声波测井的幅度偏心校正方法,其特征在于,方法包括:
获取步骤,获取测井仪器在预设深度区间产生的测井数据;所述测井数据包括泥浆脉冲反射回波数据、井眼内壁反射回波幅度以及到时数据;
提取计算步骤,根据泥浆脉冲反射回波数据,提取得到一次反射回波的第一波峰值的信号幅度和第一到时数据,以及二次反射回波的第二波峰值的信号幅度和第二到时数据,计算得到泥浆衰减因子;
校正步骤,针对所述预设深度区域中任一深度点,获取该深度点上下预设窗长范围的井眼内壁反射回波幅度,计算得到该深度点的深度校正因子;根据所述深度校正因子对该深度点的井眼内壁反射回波幅度进行幅度校正;重复执行所述校正步骤,直至所述预设深度区域各个深度点的井眼内壁反射回波幅度完成幅度校正,得到幅度校正后的超声波成像曲线。
2.根据权利要求1所述的方法,其特征在于,所述提取计算步骤进一步包括:
计算确定泥浆脉冲信号入射到反射块的传播时间;所述传播时间根据第一到时数据和第二到时数据确定;
根据所述传播时间以及测井仪器与反射块的距离,计算得到超声波在泥浆中的传播速度;
根据所述传播时间、第一波峰值的信号幅度、第二波峰值的信号幅度、传播速度计算确定泥浆衰减因子。
3.根据权利要求2所述的方法,其特征在于,所述计算确定泥浆脉冲信号入射到反射块的传播时间进一步包括:
根据第二到时数据和第一到时数据的差值,计算确定传播时间。
4.根据权利要求2所述的方法,其特征在于,所述根据所述传播时间、第一波峰值的信号幅度、第二波峰值的信号幅度、传播速度计算确定泥浆衰减因子进一步包括:
根据传播速度以及泥浆密度,得到泥浆的声阻抗;
根据泥浆的声阻抗和反射块的声阻抗,计算得到声压反射系数;
根据所述传播时间、第一波峰值的信号幅度、第二波峰值的信号幅度以及声压反射系数计算得到泥浆衰减因子。
5.根据权利要求1所述的方法,其特征在于,所述校正步骤进一步包括:
针对所述预设深度区域中任一深度点,获取该深度点上下预设窗长范围内多个深度点在不同方位的多个采集点的到时数据;
根据所述多个深度点在不同方位的多个采集点的到时数据以及泥浆衰减因子计算得到该深度点的深度校正因子;
根据该深度点的深度校正因子、泥浆衰减因子以及该深度点在预设方位的到时数据,对该深度点在预设方位的井眼内壁反射回波幅度进行幅度校正,得到该深度点在预设方位校正后的井眼内壁反射回波幅度。
6.根据权利要求1所述的方法,其特征在于,所述校正步骤进一步包括:
判断所述预设深度区域内各个深度点的井眼内壁反射回波幅度是否均进行幅度校正;
若否,重复所述校正步骤,直至完成所述预设深度区域各个深度点的井眼内壁反射回波幅度进行幅度校正;
若是,根据校正后的所述预设深度区域各个深度点的井眼内壁反射回波幅度得到幅度校正后的超声波成像曲线。
7.根据权利要求1-6中任一项所述的方法,其特征在于,不同深度点具有不同的深度校正因子。
8.一种超声波测井的幅度偏心校正装置,其特征在于,装置包括:
获取模块,适于获取测井仪器在预设深度区间产生的测井数据;所述测井数据包括泥浆脉冲反射回波数据、井眼内壁反射回波幅度以及到时数据;
提取计算模块,适于根据泥浆脉冲反射回波数据,提取得到一次反射回波的第一波峰值的信号幅度和第一到时数据,以及二次反射回波的第二波峰值的信号幅度和第二到时数据,计算得到泥浆衰减因子;
校正模块,适于针对所述预设深度区域中任一深度点,获取该深度点上下预设窗长范围的井眼内壁反射回波幅度,计算得到该深度点的深度校正因子;根据所述深度校正因子对该深度点的井眼内壁反射回波幅度进行幅度校正;重复执行所述校正模块,直至所述预设深度区域各个深度点的井眼内壁反射回波幅度完成幅度校正,得到幅度校正后的超声波成像曲线。
9.一种计算设备,其特征在于,包括:处理器、存储器、通信接口和通信总线,所述处理器、所述存储器和所述通信接口通过所述通信总线完成相互间的通信;
所述存储器用于存放至少一可执行指令,所述可执行指令使所述处理器执行如权利要求1-7中任一项所述的超声波测井的幅度偏心校正方法对应的操作。
10.一种计算机存储介质,其特征在于,所述存储介质中存储有至少一可执行指令,所述可执行指令使处理器执行如权利要求1-7中任一项所述的超声波测井的幅度偏心校正方法对应的操作。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211086231.XA CN115144906B (zh) | 2022-09-06 | 2022-09-06 | 超声波测井的幅度偏心校正方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211086231.XA CN115144906B (zh) | 2022-09-06 | 2022-09-06 | 超声波测井的幅度偏心校正方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115144906A true CN115144906A (zh) | 2022-10-04 |
CN115144906B CN115144906B (zh) | 2022-11-22 |
Family
ID=83415905
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211086231.XA Active CN115144906B (zh) | 2022-09-06 | 2022-09-06 | 超声波测井的幅度偏心校正方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115144906B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115992689A (zh) * | 2023-03-23 | 2023-04-21 | 中海油田服务股份有限公司 | 一种随钻超声成像测井装置和测井方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4930109A (en) * | 1988-12-29 | 1990-05-29 | Atlantic Richfield Company | Method and apparatus of measuring ultrasonic time travel information obtained from logging operations in a well borehole |
CN102169190A (zh) * | 2011-01-06 | 2011-08-31 | 中国科学院地质与地球物理研究所 | 一种调制补充子空间的井约束叠前弹性参数反演方法 |
CN103225502A (zh) * | 2013-04-16 | 2013-07-31 | 中国石油天然气集团公司 | 一种测量井径的系统及其方法 |
CN103616718A (zh) * | 2013-11-15 | 2014-03-05 | 中国石油天然气集团公司 | 适用于任意斜度套管井的扭杆式井下震源 |
CN105888655A (zh) * | 2016-05-05 | 2016-08-24 | 中国海洋石油总公司 | 一种随钻声波隔声结构 |
CN113779812A (zh) * | 2021-09-28 | 2021-12-10 | 中海油田服务股份有限公司 | 一种大斜度井环境中利用随钻测井数据的薄夹层识别方法 |
-
2022
- 2022-09-06 CN CN202211086231.XA patent/CN115144906B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4930109A (en) * | 1988-12-29 | 1990-05-29 | Atlantic Richfield Company | Method and apparatus of measuring ultrasonic time travel information obtained from logging operations in a well borehole |
CN102169190A (zh) * | 2011-01-06 | 2011-08-31 | 中国科学院地质与地球物理研究所 | 一种调制补充子空间的井约束叠前弹性参数反演方法 |
CN103225502A (zh) * | 2013-04-16 | 2013-07-31 | 中国石油天然气集团公司 | 一种测量井径的系统及其方法 |
CN103616718A (zh) * | 2013-11-15 | 2014-03-05 | 中国石油天然气集团公司 | 适用于任意斜度套管井的扭杆式井下震源 |
CN105888655A (zh) * | 2016-05-05 | 2016-08-24 | 中国海洋石油总公司 | 一种随钻声波隔声结构 |
CN113779812A (zh) * | 2021-09-28 | 2021-12-10 | 中海油田服务股份有限公司 | 一种大斜度井环境中利用随钻测井数据的薄夹层识别方法 |
Non-Patent Citations (1)
Title |
---|
王向阳,等: "超声波套损检测技术在煤层气井套管检测中的应用", 《中国科技论文》 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115992689A (zh) * | 2023-03-23 | 2023-04-21 | 中海油田服务股份有限公司 | 一种随钻超声成像测井装置和测井方法 |
Also Published As
Publication number | Publication date |
---|---|
CN115144906B (zh) | 2022-11-22 |
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 | |
US10253615B2 (en) | Method and a system for ultrasonic inspection of well bores | |
US9103196B2 (en) | Pipelined pulse-echo scheme for an acoustic image tool for use downhole | |
CN115992692B (zh) | 水泥环厚度测量方法、装置、电子设备及存储介质 | |
CN115144906B (zh) | 超声波测井的幅度偏心校正方法及装置 | |
CN115822564B (zh) | 套管厚度测量方法、装置、电子设备及计算机存储介质 | |
US5146432A (en) | Method for making cement impedance measurements with characterized transducer | |
CN105134170A (zh) | 一种用于评价套管井二界面水泥胶结质量的方法 | |
CN115992691A (zh) | 基于超声兰姆波的固井质量检测方法及装置 | |
US10551518B2 (en) | Method to correct eccentricity in ultrasonic image profiles | |
AU2018301988B2 (en) | Method and device for quantifying medium viscoelasticity | |
CN108571317A (zh) | 一种随钻测井井径测量系统 | |
US10126454B2 (en) | Method and system for fracture detection using acoustic waves | |
CN114151066B (zh) | 一种超声Lamb波测井井壁声学界面逆时偏移成像方法 | |
AU2016396055B2 (en) | Casing thickness estimation by frequency correlation | |
CN116556925A (zh) | 井孔椭圆度的计算方法、装置、计算设备及存储介质 | |
CN114139584A (zh) | 一种随钻声波测井套管波衰减反演方法和装置 | |
CN116660997B (zh) | 套管内外介质声阻抗反演方法、装置及电子设备 | |
US20180372688A1 (en) | Method for ultrasonic testing of an object | |
CN114135274A (zh) | 一种随钻方位固井质量扫描方法和装置 | |
CN112904348A (zh) | 一种三维探测方法、装置、设备和存储介质 | |
CN116241239B (zh) | 基于远近单极的固井评价的方法、装置、设备及存储介质 | |
CN116220667B (zh) | 基于超声兰姆波测井的直达波与反射波分离方法及装置 | |
CN115951423B (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |