CN116660997B - 套管内外介质声阻抗反演方法、装置及电子设备 - Google Patents

套管内外介质声阻抗反演方法、装置及电子设备 Download PDF

Info

Publication number
CN116660997B
CN116660997B CN202310965604.9A CN202310965604A CN116660997B CN 116660997 B CN116660997 B CN 116660997B CN 202310965604 A CN202310965604 A CN 202310965604A CN 116660997 B CN116660997 B CN 116660997B
Authority
CN
China
Prior art keywords
waveform
acoustic impedance
ultrasonic echo
casing
sleeve
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
Application number
CN202310965604.9A
Other languages
English (en)
Other versions
CN116660997A (zh
Inventor
陶爱华
王文梁
张勇
程林波
李�亨
刘汇鑫
李�学
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Oilfield Services Ltd
Original Assignee
China Oilfield Services Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by China Oilfield Services Ltd filed Critical China Oilfield Services Ltd
Priority to CN202310965604.9A priority Critical patent/CN116660997B/zh
Publication of CN116660997A publication Critical patent/CN116660997A/zh
Application granted granted Critical
Publication of CN116660997B publication Critical patent/CN116660997B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/40Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging

Landscapes

  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Acoustics & Sound (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • Remote Sensing (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Length Measuring Devices Characterised By Use Of Acoustic Means (AREA)
  • Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)

Abstract

本发明公开了一种套管内外介质声阻抗反演方法、装置及电子设备,涉及超声脉冲回波固井质量评价技术领域。其中,所述方法包括:获取探头在井筒预设空间位置处的实测超声回波波形和理论仿真波形,对所述实测超声回波波形和所述理论仿真波形进行调整,得到目标仿真波形,根据所述实测超声回波波形和所述目标仿真波形反演所述套管的套管厚度,根据所述套管厚度、所述实测超声回波波形的包络特征和所述目标仿真波形的包络特征同时反演所述套管内泥浆声阻抗和所述套管外水泥声阻抗,可以提高反演结果准确性。

Description

套管内外介质声阻抗反演方法、装置及电子设备
技术领域
本发明涉及超声波固井质量评价领域,具体涉及一种套管内外介质声阻抗反演方法、装置及电子设备。
背景技术
油气井进行固井后需要对固井水泥质量进行检测,超声扫描成像测井利用脉冲-回波检测方法测量套后环空介质的声阻抗,是评价固井水泥胶结质量的重要手段。现有计算套管内外介质声阻抗的方法主要有以下两种:
(1)利用套管共振波归一化幅度的对数与水泥声阻抗之间的线性关系建立套管后水泥声阻抗的计算公式。为求得水泥声阻抗,在测量井段内需要找到自由套管作为参考点。如果没有自由套管,则无法计算出准确的水泥声阻抗。另一方面,在斜井或水平井中,由于泥浆沉淀使得泥浆周向不均匀,此方法无法准确计算圆周上各个方位的水泥声阻抗。
(2)通过共振波中心频率计算套管壁厚,用雷克子波作为仿真声源,建立套管内外介质多层流体的一维平面声波模型,计算仿真波形共振波包络的截距和斜率特征,并与实测波形的截距和斜率特征联立建立方程组,利用二分法迭代依次计算出泥浆声阻抗和水泥声阻抗。该方法忽略了套管的固体特性,无法模拟出探头有限尺寸及套管的声场三维特征;采用的声源无法模拟超声探头在高温高压条件下的频谱变化。因此由理论仿真波形引入的误差将导致反演结果的不准确。
发明内容
鉴于上述问题,提出了本发明以便提供一种克服上述问题或者至少部分地解决上述问题的套管内外介质声阻抗反演方法、装置及电子设备。
根据本发明的一个方面,提供了一种套管内外介质声阻抗反演方法,包括:
获取探头在井筒预设空间位置处的实测超声回波波形和理论仿真波形;
对所述实测超声回波波形和所述理论仿真波形进行调整,得到目标仿真波形;
根据所述实测超声回波波形和所述目标仿真波形反演所述套管的套管厚度;
根据所述套管厚度、所述实测超声回波波形的包络特征和所述目标仿真波形的包络特征同时反演所述套管内泥浆声阻抗和所述套管外水泥声阻抗。
根据本发明的另一个方面,提供了一种套管内外介质声阻抗反演装置,包括:
波形获取模块,用于获取探头在井筒预设空间位置处的实测超声回波波形和理论仿真波形;
波形调整模块,用于对所述实测超声回波波形和所述理论仿真波形进行调整,得到目标仿真波形;
厚度反演模块,用于根据所述实测超声回波波形和所述目标仿真波形反演所述套管的套管厚度;
声阻抗反演模块,用于根据所述套管厚度、所述实测超声回波波形的包络特征和所述目标仿真波形的包络特征同时反演所述套管内泥浆声阻抗和所述套管外水泥声阻抗。
根据本发明的另一方面,提供了一种电子设备,包括:处理器、存储器、通信接口和通信总线,所述处理器、所述存储器和所述通信接口通过所述通信总线完成相互间的通信;
所述存储器用于存放至少一可执行指令,所述可执行指令使所述处理器执行本发明所述的套管内外介质声阻抗反演方法对应的操作。
根据本发明的另一方面,提供了一种计算机存储介质,所述存储介质中存储有至少一可执行指令,所述可执行指令使处理器执行本发明所述的套管内外介质声阻抗反演方法对应的操作。
根据本发明公开的套管内外介质声阻抗反演方法、装置及电子设备,通过获取探头在井筒预设空间位置处的实测超声回波波形和理论仿真波形,对所述实测超声回波波形和所述理论仿真波形进行调整,得到目标仿真波形,根据所述实测超声回波波形和所述目标仿真波形反演所述套管的套管厚度,根据所述套管厚度、所述实测超声回波波形的包络特征和所述目标仿真波形的包络特征同时反演所述套管内泥浆声阻抗和所述套管外水泥声阻抗,可以提高反演结果准确性。
上述说明仅是本发明技术方案的概述,为了能够更清楚了解本发明的技术手段,而可依照说明书的内容予以实施,并且为了让本发明的上述和其它目的、特征和优点能够更明显易懂,以下特举本发明的具体实施方式。
附图说明
通过阅读下文优选实施方式的详细描述,各种其他的优点和益处对于本领域普通技术人员将变得清楚明了。附图仅用于示出优选实施方式的目的,而并不认为是对本发明的限制。而且在整个附图中,用相同的参考符号表示相同的部件。在附图中:
图1示出了根据本发明实施例一提供的一种套管内外介质声阻抗反演方法的流程示意图;
图2示出了根据本发明实施例二提供的一种套管内外介质声阻抗反演装置的结构示意图;
图3示出了根据本发明实施例四提供的一种电子设备的结构示意图。
具体实施方式
下面将参照附图更详细地描述本发明的示例性实施例。虽然附图中显示了本发明的示例性实施例,然而应当理解,可以以各种形式实现本发明而不应被这里阐述的实施例所限制。相反,提供这些实施例是为了能够更透彻地理解本发明,并且能够将本发明的范围完整的传达给本领域的技术人员。
实施例一
图1示出了根据本发明实施例一提供的一种套管内外介质声阻抗反演方法的流程示意图。本实施例的执行主体为本发明实施例提供的套管内外介质声阻抗反演装置,该装置可以通过软件或硬件实现,该装置可以集成在超声脉冲反射法测井仪中。如图1所示,该方法包括:
步骤S11,获取探头在井筒预设空间位置处套管中的实测超声回波波形和理论仿真波形。
其中,预设空间位置为套管井井眼内某一具体位置,例如井眼内600米的位置,在该位置采用探头扫描一周可以得到多个方位处套管的超声脉冲回波即实测超声回波波形。
其中,理论仿真波形可以通过高精度的井下理论仿真实现。
步骤S12,对实测超声回波波形和理论仿真波形进行调整,得到目标仿真波形。
由于探头自身的频率响应以及井下高温高压环境的影响,即使使用高精度正演技术仍然难以计算出与实际波形完全一样的波形。为此本实施例提出了频谱校正技术,校正探头自身频率响应和井下环境的影响。具体的,基于实测超声回波波形与理论仿真波形提取的声源频谱的差异,校正理论仿真波形的频谱特性,这样可以提升理论仿真波形的准确度。
步骤S13,根据实测超声回波波形和目标仿真波形反演套管的套管厚度。
具体的,可以根据实测超声回波波形和目标仿真波形的群延时曲线和半高宽反演套管的套管厚度。其中,群延时曲线的峰值频率主要由反射波的厚度和材料声速决定,而半高宽主要由泥浆声阻抗决定。具体实现参见下述可选实施例中的相关描述。
步骤S14,根据套管厚度、实测超声回波波形的包络特征和目标仿真波形的包络特征同时反演套管内泥浆声阻抗和套管外水泥声阻抗。
其中,该包络特征可以为部分包络特征,也可以为完整包络特征。具体的,根据套管厚度、实测超声回波波形的部分包络特征和目标仿真波形的部分包络特征可以同时得到套管内泥浆声阻抗和套管外水泥声阻抗,或者,根据套管厚度、实测超声回波波形的完整包络特征和目标仿真波形的完整包络特征同时得到套管内泥浆声阻抗和套管外水泥声阻抗。
由此可见,本实施例通过获取探头在井筒预设空间位置处套管中的实测超声回波波形和理论仿真波形,对实测超声回波波形和理论仿真波形进行调整,得到目标仿真波形,根据实测超声回波波形和目标仿真波形反演套管的套管厚度,根据套管厚度、实测超声回波波形的包络特征和目标仿真波形的包络特征同时反演套管内泥浆声阻抗和套管外水泥声阻抗,可以提高反演结果准确性。
在一个可选实施例中,步骤S11具体包括:
步骤S111,根据套管参数、探头参数、水泥参数和泥浆参数构建三维声波传播模型。
其中,套管参数至少包括套管外径、套管厚度、套管最小厚度、套管密度、纵波速度、横波速度;探头参数至少包括探头的扫描半径、探头晶片的宽与高、探头的固有频率、变迹参数、探头在套管中的位置、探头偏心距、探头方位角;泥浆参数至少包括密度、纵波速度;水泥参数至少包括密度、纵波速度、横波速度。根据上述参数构建三维声波传播模型。
步骤S112,通过所述三维声波传播模型确定所述套管的理论仿真波形。
其中,三维声波传播模型采用高精度正演技术模拟生成套管的理论仿真波形。另外,可以将理论仿真波形形成波形库或频谱库存储下来,在后面反演计算时可以直接从库中提取,这样可提高迭代反演的速度。
在一个可选实施例中,步骤S12具体包括:
步骤S121,基于所述实测超声回波波形和所述理论仿真波形的声源频谱差异,对所述理论仿真波形的频谱特性进行校正。
对实测超声回波波形和所述理论仿真波形分别提取表面回波窗和有用全波窗进行滤波,将滤波的波形数据中共振波后面的二次反射波去掉,保留有用的接收信号,施加快速傅里叶变换,确定有效的频率范围,并对该范围内的频谱进行趋势拟合,剔除凹陷区域频值点,进行精确拟合,并通过施加反傅里叶变换得到声源信号,将重建的声源信号与理论仿真声源信号进行对比,求取频谱校正系数,待反演时基于实测超声回波波形和所述理论仿真波形提取的声源频谱之间的校正系数,校正理论仿真波形的频谱特性,该步骤可提升理论仿真波形的准确度,校正理论仿真波形的频谱特性,可提升理论仿真波形的准确度。
步骤S122,对所述实测超声回波波形和校正后的理论仿真波形进行波形重整,以规则化波形的形态。
具体的,对实测超声回波波形和校正后的理论仿真波形应用自适应汉宁窗长度的波形重整技术,分别提取表面回波窗和有用全波窗,计算傅里叶变换,并归一化频谱,通过群延时计算共振波中心频率,在频域设计高斯幅度谱叠加在归一化频谱中,使共振波中心频率的泄漏峰位于频谱的中心,最后施加反傅里叶变换获得时域波形,消除频谱中的异常数据点,使得时域波形的包络更加光滑,提高反演结果的稳定性。
在一个可选实施例中,步骤S13具体包括:
步骤S131,根据所述实测超声回波波形和所述目标仿真波形分别确定各自的群延时曲线和半高宽。
其中,群延时曲线的峰值频率主要由套管的厚度和材料声速决定,而半高宽主要由泥浆声阻抗决定。
步骤S132,根据所述实测超声回波波形的群延时曲线和半高宽与所述目标仿真波形的群延时曲线和半高宽建立迭代方程。
由于套管内的泥浆与套管外介质(通常为水泥)不同时,其声阻抗群延时曲线形状也不同,可以用峰值高度和半高宽度/>两个参数来进行表征。峰值高度是峰的最高点至两侧最低点平均值的距离,单位为/>s;半高宽度是从最高点向两侧下降至峰值高度/>时的频率点之差,单位为/>。假定理论仿真波形为/>,其中,Zm为泥浆声阻抗,d为反射板厚度,则仿真理论波形的群延时曲线表达为Gd(Zm,d),群延时曲线的峰值频率可表达为/>,半高宽表达为
泥浆的声速可以通过回波时间和探头与反射板之间的距离计算得到,反射板的声学参数预先已知。对于实际测量得到的反射板回波信号即实测超声回波波形,可计算其群延时的峰值频率和半高宽,分别记为F0和Hw0,则参数反演可转化为求解如下方程组:
结合牛顿迭代公式:,其中,
其中,,/>,/>为雅可比矩阵。采用求解方程组的迭代公式得:
其中,由于雅可比矩阵没有解析的表达式,所以实际计算中采用差分方法数值计算导数。由于声阻抗对群延时频率的影响微乎其微,所以雅可比矩阵可写成,迭代公式可以将厚度和泥浆声阻抗分离出来。又/>,所以,迭代公式为:
通过上述公式可以独立反演套管厚度和水泥声阻抗。
步骤S133,根据所述迭代方程对所述套管厚度和泥浆声阻抗进行反演。
在一个可选实施例中,步骤S14具体包括:
步骤S141,分别对所述目标仿真波形和所述实测超声回波波形的共振波包络进行线性回归直线拟合,根据各自拟合直线的截距和斜率建立迭代方程。
具体的,假设目标仿真波形为,其中,套管厚度d可以由厚度反演模块23预先反演出来,因此目标仿真波形简化为/>,其中,/>为水泥声阻抗。目标仿真波形的共振波包络拟合直线的截距为/>,斜率为/>,对于实测超声回波波形,计算其共振波包络取对数后拟合直线的截距和斜率,分别记为/>和/>,则参数反演转化为如下方程组:
采用求解方程组的迭代公式后得:
或,基于所述套管厚度、所述目标仿真波形和所述实测超声回波波形的全包络建立迭代方程。
具体的,假定目标理论仿真波形为,其中,套管厚度d可以由厚度反演模块23先反演出来,包络特征记为/>,其中T为采样三时间间隔,。对于测量的套管回波即实测超声回波波形的包络特征记为/>,采用最小二乘优化方法建立迭代方程组,其中,最小二乘误差为:
其中,雅可比矩阵为/>,采用求解方程组的迭代公式为:
步骤S142,基于所述迭代方程同时反演所述套管内泥浆声阻抗和所述套管外水泥声阻抗。
上述实施例在构建迭代反演算法时,最大限度的利用了回波信号的完整(或部分)包络所含信息,避免了对包络直线拟合可能漏失的信息。
在一个可选实施例中,上述方法还包括:
步骤S15,对反演结果中套管内泥浆声阻抗和套管外水泥声阻抗的无效值进行优化。
具体的,可以基于深度方向和时间采样方向的数据优化算法,对反演计算结果Zm、Zc、d进行数据优化处理。根据待优化点周围的9点和25点去除无效值后采取均值替代,减少瑕疵,平滑图像。
实施例二
图2示出了根据本发明实施例二提供的一种套管内外介质声阻抗反演装置的结构示意图。该装置包括:波形获取模块21、波形调整模块22、厚度反演模块23和声阻抗反演模块24;其中,
波形获取模块21用于获取探头在井筒预设空间位置处套管中的实测超声回波波形和理论仿真波形;
波形调整模块22用于对所述实测超声回波波形和所述理论仿真波形进行调整,得到目标仿真波形;
厚度反演模块23用于根据所述实测超声回波波形和所述目标仿真波形反演所述套管的套管厚度;
声阻抗反演模块24用于根据所述套管厚度、所述实测超声回波波形的包络特征和所述目标仿真波形的包络特征同时反演所述套管内泥浆声阻抗和所述套管外水泥声阻抗。
进一步的,所述波形获取模块21具体用于:根据套管参数、探头参数、水泥参数和泥浆参数构建三维声波传播模型;通过所述三维声波传播模型确定所述套管的理论仿真波形。
进一步的,所述波形调整模块22具体用于:基于所述实测超声回波波形和所述理论仿真波形的声源频谱差异,对所述理论仿真波形的频谱特征进行校正;对所述实测超声回波波形和校正后的理论仿真波形进行波形重整,以规则化波形的形态。
进一步的,所述声阻抗反演模块24具体用于:分别对所述目标仿真波形和所述实测超声回波波形的共振波包络进行线性回归直线拟合,根据各自的拟合直线的截距和斜率建立迭代方程;或基于所述套管厚度、所述目标仿真波形和和所述实测超声回波波形的全包络建立迭代方程;基于所述迭代方程同时反演所述套管内泥浆声阻抗和所述套管外水泥声阻抗。
进一步的,所述声阻抗反演模块24具体用于:分别对所述目标仿真波形和所述实测超声回波波形的共振波全包络进行采样,得到各自的采样包络特征;基于所述套管厚度、所述目标仿真波形的采样包络特征、所述实测超声回波波形的采样包络特征建立迭代方程;基于所述迭代方程同时反演所述套管内泥浆声阻抗和所述套管外水泥声阻抗。
进一步的,所述厚度反演模块23具体用于:根据所述实测超声回波波形和所述目标仿真波形分别确定各自的群延时曲线和半高宽;根据所述实测超声回波波形的群延时曲线和半高宽以及所述目标仿真波形的群延时曲线和半高宽建立迭代方程;根据所述迭代方程对所述套管厚度和泥浆声阻抗进行反演。
进一步的,所述装置还包括:优化模块25;
所述优化模块25用于对反演结果中套管内泥浆声阻抗和套管外水泥声阻抗的无效值进行优化。
本实施例所述的套管内外介质声阻抗反演装置用于执行上述实施例一所述的套管内外介质声阻抗反演方法,其工作原理与技术效果类似,这里不再赘述。
实施例三
本发明实施例三提供了一种非易失性计算机存储介质,所述计算机存储介质存储有至少一可执行指令,该计算机可执行指令可执行上述任意方法实施例中的套管内外介质声阻抗反演方法。
实施例四
图3示出了根据本发明实施例四提供的一种电子设备的结构示意图。本发明具体实施例并不对电子设备的具体实现做限定。
如图3所示,该电子设备可以包括:处理器302、通信接口304、存储器306、以及通信总线308。
其中:处理器302、通信接口304、以及存储器306通过通信总线308完成相互间的通信。通信接口304,用于与其它设备比如客户端或其它服务器等的网元通信。处理器302,用于执行程序310,具体可以执行上述方法实施例中的相关步骤。
具体地,程序310可以包括程序代码,该程序代码包括计算机操作指令。
处理器302可能是中央处理器CPU,或者是特定集成电路ASIC,或者是被配置成实施本发明实施例的一个或多个集成电路。电子设备包括的一个或多个处理器,可以是同一类型的处理器,如一个或多个CPU;也可以是不同类型的处理器,如一个或多个CPU以及一个或多个ASIC。
存储器306,用于存放程序310。存储器306可能包含高速RAM存储器,也可能还包括非易失性存储器,例如至少一个磁盘存储器。
程序310具体可以用于使得处理器302执行上述任意方法实施例中的套管内外介质声阻抗反演方法。
在此提供的算法或显示不与任何特定计算机、虚拟系统或者其它设备固有相关。各种通用系统也可以与基于在此的示教一起使用。根据上面的描述,构造这类系统所要求的结构是显而易见的。此外,本发明实施例也不针对任何特定编程语言。应当明白,可以利用各种编程语言实现在此描述的本发明的内容,并且上面对特定语言所做的描述是为了披露本发明的最佳实施方式。
在此处所提供的说明书中,说明了大量具体细节。然而,能够理解,本发明的实施例可以在没有这些具体细节的情况下实践。在一些实例中,并未详细示出公知的方法、结构和技术,以便不模糊对本说明书的理解。
类似地,应当理解,为了精简本发明并帮助理解各个发明方面中的一个或多个,在上面对本发明的示例性实施例的描述中,本发明实施例的各个特征有时被一起分组到单个实施例、图、或者对其的描述中。
本领域那些技术人员可以理解,可以对实施例中的设备中的模块进行自适应性地改变并且把它们设置在与该实施例不同的一个或多个设备中。可以把实施例中的模块或单元或组件组合成一个模块或单元或组件,以及此外可以把它们分成多个子模块或子单元或子组件。除了这样的特征和/或过程或者单元中的至少一些是相互排斥之外,可以采用任何组合对本说明书中公开的所有特征以及如此公开的任何方法或者设备的所有过程或单元进行组合。除非另外明确陈述,本说明书中公开的每个特征可以由提供相同、等同或相似目的的替代特征来代替。
此外,本领域的技术人员能够理解,尽管在此的一些实施例包括其它实施例中所包括的某些特征而不是其它特征,但是不同实施例的特征的组合意味着处于本发明的范围之内并且形成不同的实施例。 M 。
本发明的各个部件实施例可以以硬件实现,或者以在一个或者多个处理器上运行的软件模块实现,或者以它们的组合实现。本领域的技术人员应当理解,可以在实践中使用微处理器或者数字信号处理器(DSP)来实现根据本发明实施例的一些或者全部部件的一些或者全部功能。本发明还可以实现为用于执行这里所描述的方法的一部分或者全部的设备或者装置程序(例如,计算机程序和计算机程序产品)。这样的实现本发明的程序可以存储在计算机可读介质上,或者可以具有一个或者多个信号的形式。这样的信号可以从因特网网站上下载得到,或者在载体信号上提供,或者以任何其他形式提供。
应该注意的是上述实施例对本发明进行说明而不是对本发明进行限制,并且本领域技术人员在不脱离本发明保护范围的情况下可设计出替换实施例。本发明可以借助于包括有若干不同元件的硬件以及借助于适当编程的计算机来实现。在列举了若干装置中,这些装置中的若干个可以是通过同一个硬件项来具体体现。单词第一、第二、以及第三等的使用不表示任何顺序。可将这些单词解释为名称。上述实施例中的步骤,除有特殊说明外,不应理解为对执行顺序的限定。

Claims (6)

1.一种套管内外介质声阻抗反演方法,其特征在于,包括:
获取探头在井筒预设空间位置处的实测超声回波波形和理论仿真波形;
对所述实测超声回波波形和所述理论仿真波形进行调整,得到目标仿真波形;
根据所述实测超声回波波形和所述目标仿真波形反演所述套管的套管厚度;
根据所述套管厚度、所述实测超声回波波形的包络特征和所述目标仿真波形的包络特征同时反演所述套管内泥浆声阻抗和所述套管外水泥声阻抗;
其中,对所述实测超声回波波形和所述理论仿真波形进行调整,包括:
基于所述实测超声回波波形和所述理论仿真波形的声源频谱差异,对所述理论仿真波形的频谱特征进行校正;
对所述实测超声回波波形和校正后的理论仿真波形进行波形重整,以规则化波形的形态;
其中,根据所述套管厚度、所述实测超声回波波形的包络特征和所述目标仿真波形的包络特征同时反演所述套管内泥浆声阻抗和所述套管外水泥声阻抗,包括:
分别对所述目标仿真波形和所述实测超声回波波形的共振波包络进行线性回归直线拟合,根据各自的拟合直线的截距和斜率建立迭代方程;或基于所述套管厚度、所述目标仿真波形和所述实测超声回波波形的全包络建立迭代方程;基于所述迭代方程同时反演所述套管内泥浆声阻抗和所述套管外水泥声阻抗;或
分别对所述目标仿真波形和所述实测超声回波波形的共振波全包络进行采样,得到各自的采样包络特征;基于所述套管厚度、所述目标仿真波形的采样包络特征、所述实测超声回波波形的采样包络特征建立迭代方程;基于所述迭代方程同时反演所述套管内泥浆声阻抗和所述套管外水泥声阻抗;
其中,根据所述实测超声回波波形和所述目标仿真波形反演所述套管的套管厚度,包括:
根据所述实测超声回波波形和所述目标仿真波形分别确定各自的群延时曲线和半高宽;根据所述实测超声回波波形的群延时曲线和半高宽以及所述目标仿真波形的群延时曲线和半高宽建立迭代方程;根据所述迭代方程对所述套管厚度和泥浆声阻抗进行反演。
2.根据权利要求1所述的方法,其特征在于,所述获取探头在井筒预设空间位置处的理论仿真波形,包括:
根据套管参数、探头参数、水泥参数和泥浆参数构建三维声波传播模型;
通过所述三维声波传播模型确定所述套管的理论仿真波形。
3.根据权利要求1或2所述的方法,其特征在于,所述方法还包括:
对反演结果中套管内泥浆声阻抗和套管外水泥声阻抗的无效值进行优化。
4.一种套管内外介质声阻抗反演装置,其特征在于,包括:
波形获取模块,用于获取探头在井筒预设空间位置处套管中的实测超声回波波形和理论仿真波形;
波形调整模块,用于对所述实测超声回波波形和所述理论仿真波形进行调整,得到目标仿真波形;
厚度反演模块,用于根据所述实测超声回波波形和所述目标仿真波形反演所述套管的套管厚度;
声阻抗反演模块,用于根据所述套管厚度、所述实测超声回波波形的包络特征和所述目标仿真波形的包络特征同时反演所述套管内泥浆声阻抗和所述套管外水泥声阻抗;
其中,所述波形调整模块具体用于:基于所述实测超声回波波形和所述理论仿真波形的声源频谱差异,对所述理论仿真波形的频谱特征进行校正;对所述实测超声回波波形和校正后的理论仿真波形进行波形重整,以规则化波形的形态;
其中,所述声阻抗反演模块具体用于:分别对所述目标仿真波形和所述实测超声回波波形的共振波包络进行线性回归直线拟合,根据各自的拟合直线的截距和斜率建立迭代方程;或基于所述套管厚度、所述目标仿真波形和所述实测超声回波波形的全包络建立迭代方程;基于所述迭代方程同时反演所述套管内泥浆声阻抗和所述套管外水泥声阻抗;或,分别对所述目标仿真波形和所述实测超声回波波形的共振波全包络进行采样,得到各自的采样包络特征;基于所述套管厚度、所述目标仿真波形的采样包络特征、所述实测超声回波波形的采样包络特征建立迭代方程;基于所述迭代方程同时反演所述套管内泥浆声阻抗和所述套管外水泥声阻抗;
其中,所述厚度反演模块具体用于:根据所述实测超声回波波形和所述目标仿真波形分别确定各自的群延时曲线和半高宽;根据所述实测超声回波波形的群延时曲线和半高宽以及所述目标仿真波形的群延时曲线和半高宽建立迭代方程;根据所述迭代方程对所述套管厚度和泥浆声阻抗进行反演。
5.一种电子设备,包括:处理器、存储器、通信接口和通信总线,所述处理器、所述存储器和所述通信接口通过所述通信总线完成相互间的通信;
所述存储器用于存放至少一可执行指令,所述可执行指令使所述处理器执行如权利要求1-3中任一项所述的套管内外介质声阻抗反演方法对应的操作。
6.一种计算机存储介质,所述存储介质中存储有至少一可执行指令,所述可执行指令使处理器执行如权利要求1-3中任一项所述的套管内外介质声阻抗反演方法对应的操作。
CN202310965604.9A 2023-08-02 2023-08-02 套管内外介质声阻抗反演方法、装置及电子设备 Active CN116660997B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310965604.9A CN116660997B (zh) 2023-08-02 2023-08-02 套管内外介质声阻抗反演方法、装置及电子设备

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310965604.9A CN116660997B (zh) 2023-08-02 2023-08-02 套管内外介质声阻抗反演方法、装置及电子设备

Publications (2)

Publication Number Publication Date
CN116660997A CN116660997A (zh) 2023-08-29
CN116660997B true CN116660997B (zh) 2023-09-29

Family

ID=87714037

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202310965604.9A Active CN116660997B (zh) 2023-08-02 2023-08-02 套管内外介质声阻抗反演方法、装置及电子设备

Country Status (1)

Country Link
CN (1) CN116660997B (zh)

Citations (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2569476A1 (fr) * 1984-08-24 1986-02-28 Schlumberger Prospection Procede et dispositif pour evaluer la qualite du ciment entourant le tubage d'un puits
CN1038489A (zh) * 1988-06-09 1990-01-03 中国科学院声学研究所 声波全面评价水泥胶结固井质量的方法
US7753117B2 (en) * 2008-04-04 2010-07-13 Schlumberger Technology Corporation Tool and method for evaluating fluid dynamic properties of a cement annulus surrounding a casing
CN102128028A (zh) * 2010-01-12 2011-07-20 同济大学 用于评估水泥与地层界面胶结质量的超声信号反演方法
CN106383173A (zh) * 2016-09-13 2017-02-08 中国海洋石油总公司 一种水泥声阻抗计算方法和装置
CN106595538A (zh) * 2017-01-13 2017-04-26 长江大学 一种检测套管外水泥与套管间裂隙宽度的系统及方法
EP3179277A1 (en) * 2015-12-11 2017-06-14 Services Pétroliers Schlumberger Resonance-based inversion of acoustic impedance of annulus behind casing
CN110806601A (zh) * 2019-11-08 2020-02-18 中国海洋石油集团有限公司 套内外物质阻抗计算、关系确定、固井评价方法和设备
CN111007565A (zh) * 2019-12-24 2020-04-14 清华大学 三维频率域全声波成像方法及装置
CN111335888A (zh) * 2020-03-28 2020-06-26 西安石油大学 一种确定地质层中井筒性质的方法
CN114740532A (zh) * 2022-03-18 2022-07-12 电子科技大学 一种水泥声阻抗、泥浆声阻抗的计算方法
CN114941521A (zh) * 2022-06-16 2022-08-26 中海油田服务股份有限公司 一种超声波测井方法、装置及设备
CN115324564A (zh) * 2022-10-11 2022-11-11 中海油田服务股份有限公司 固井质量检测方法、装置、计算设备及存储介质
CN115822564A (zh) * 2023-02-20 2023-03-21 中海油田服务股份有限公司 套管厚度测量方法、装置、电子设备及计算机存储介质
CN115992692A (zh) * 2023-03-23 2023-04-21 中海油田服务股份有限公司 水泥环厚度测量方法、装置、电子设备及存储介质
CN115992691A (zh) * 2023-02-20 2023-04-21 中海油田服务股份有限公司 基于超声兰姆波的固井质量检测方法及装置

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9784875B2 (en) * 2014-01-31 2017-10-10 Schlumberger Technology Corporation Method to estimate cement acoustic wave speeds from data acquired by a cased hole ultrasonic cement evaluation tool
US10012749B2 (en) * 2014-10-20 2018-07-03 Schlumberger Technology Corporation Fast model based inversion of acoustic impedance of annulus behind casing
US11359488B2 (en) * 2019-03-12 2022-06-14 Baker Hughes Oilfield Operations Llc Self-calibrated method of determining borehole fluid acoustic properties

Patent Citations (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2569476A1 (fr) * 1984-08-24 1986-02-28 Schlumberger Prospection Procede et dispositif pour evaluer la qualite du ciment entourant le tubage d'un puits
CN1038489A (zh) * 1988-06-09 1990-01-03 中国科学院声学研究所 声波全面评价水泥胶结固井质量的方法
US7753117B2 (en) * 2008-04-04 2010-07-13 Schlumberger Technology Corporation Tool and method for evaluating fluid dynamic properties of a cement annulus surrounding a casing
CN102128028A (zh) * 2010-01-12 2011-07-20 同济大学 用于评估水泥与地层界面胶结质量的超声信号反演方法
EP3179277A1 (en) * 2015-12-11 2017-06-14 Services Pétroliers Schlumberger Resonance-based inversion of acoustic impedance of annulus behind casing
CN106383173A (zh) * 2016-09-13 2017-02-08 中国海洋石油总公司 一种水泥声阻抗计算方法和装置
CN106595538A (zh) * 2017-01-13 2017-04-26 长江大学 一种检测套管外水泥与套管间裂隙宽度的系统及方法
CN110806601A (zh) * 2019-11-08 2020-02-18 中国海洋石油集团有限公司 套内外物质阻抗计算、关系确定、固井评价方法和设备
CN111007565A (zh) * 2019-12-24 2020-04-14 清华大学 三维频率域全声波成像方法及装置
CN111335888A (zh) * 2020-03-28 2020-06-26 西安石油大学 一种确定地质层中井筒性质的方法
CN114740532A (zh) * 2022-03-18 2022-07-12 电子科技大学 一种水泥声阻抗、泥浆声阻抗的计算方法
CN114941521A (zh) * 2022-06-16 2022-08-26 中海油田服务股份有限公司 一种超声波测井方法、装置及设备
CN115324564A (zh) * 2022-10-11 2022-11-11 中海油田服务股份有限公司 固井质量检测方法、装置、计算设备及存储介质
CN115822564A (zh) * 2023-02-20 2023-03-21 中海油田服务股份有限公司 套管厚度测量方法、装置、电子设备及计算机存储介质
CN115992691A (zh) * 2023-02-20 2023-04-21 中海油田服务股份有限公司 基于超声兰姆波的固井质量检测方法及装置
CN115992692A (zh) * 2023-03-23 2023-04-21 中海油田服务股份有限公司 水泥环厚度测量方法、装置、电子设备及存储介质

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
声阻抗技术在泡沫水泥固井质量评价中的应用;李兆阳;;勘探地球物理进展(第05期);全文 *
多功能超声成像测井仪在套损及固井质量评价中的应用;孙志峰;王春艳;陈洪海;刘西恩;;石油天然气学报(第06期);全文 *
套管中模式波的响应特征;陈雪莲;唐晓明;刘临政;;应用声学(第01期);全文 *
微环厚度对垂直入射超声脉冲反射波影响分析;孙志峰;王春艳;罗瑜林;陶爱华;陈洪海;;测井技术(第03期);全文 *
有效测量套管井泄漏Lamb波衰减的实验研究;李疾翎;陶爱华;刘临政;程林波;王明辉;;应用声学(第01期);全文 *

Also Published As

Publication number Publication date
CN116660997A (zh) 2023-08-29

Similar Documents

Publication Publication Date Title
US11300547B2 (en) Ultrasound matrix inspection
US9322807B2 (en) Ultrasonic signal time-frequency decomposition for borehole evaluation or pipeline inspection
CN115144906B (zh) 超声波测井的幅度偏心校正方法及装置
CN110806601B (zh) 套内外物质阻抗计算、关系确定、固井评价方法和设备
CN115992692B (zh) 水泥环厚度测量方法、装置、电子设备及存储介质
US5146432A (en) Method for making cement impedance measurements with characterized transducer
CN115822564B (zh) 套管厚度测量方法、装置、电子设备及计算机存储介质
TWI686180B (zh) 運動資訊獲取方法及裝置
CN115324564B (zh) 固井质量检测方法、装置、计算设备及存储介质
CN115992691A (zh) 基于超声兰姆波的固井质量检测方法及装置
CA3070622A1 (en) Method and device for quantifying viscoelasticity of a medium
CN108508093A (zh) 一种工件缺陷高度的检测方法及系统
Mandal et al. A new monocable circumferential acoustic scanner tool (CAST-M) for cased-hole and openhole applications
CN116660997B (zh) 套管内外介质声阻抗反演方法、装置及电子设备
CN116556925A (zh) 井孔椭圆度的计算方法、装置、计算设备及存储介质
WO2023123952A1 (zh) 利用das井中地震数据反演地层波阻抗的方法及装置
CN112904348B (zh) 一种三维探测方法、装置、设备和存储介质
CN116220667B (zh) 基于超声兰姆波测井的直达波与反射波分离方法及装置
WO2017015954A1 (zh) 一种地震信号处理方法、装置和系统
US10578589B2 (en) System and method for ultrasound inspection with time reversal
Gunarathne et al. Development of a synthetic A-scan technique for ultrasonic testing of pipelines
CN116241239B (zh) 基于远近单极的固井评价的方法、装置、设备及存储介质
CN116201526B (zh) 微环隙检测方法、装置、计算设备及存储介质
CN113837163B (zh) 一种基于三维探地雷达的隧道监测方法、系统以及存储介质
Liang et al. Maximum non-Gaussianity parameters estimation of ultrasonic echoes and its application in ultrasonic non-destructive evaluation

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