基于密度聚类的偶极子声波测井频散校正方法及相关装置制造方法及图纸

本发明专利技术公开了一种基于密度聚类的偶极子声波测井频散校正方法及相关装置，其中基于密度聚类的偶极子声波测井频散校正方法包括：对目标深度点的偶极子声波测井数据进行频散分析得到所述目标深度点对应的频散曲线；重复执行直至获取到目标层段的频散曲线组；利用密度聚类对所述频散曲线组进行处理，得到所述频散曲线组中的有效部分；对所述频散曲线组中的有效部分进行拟合，得到校正结果。本发明专利技术利用密度聚类对频散曲线进行聚类筛选，无需预先进行大量数据的训练，提升了有效数据的提取效率；通过对频散曲线中的有效部分进行拟合得到校正结果，在保证频散校正效率的同时提高了从实际资料中获取地层真实横波速度的准确率。际资料中获取地层真实横波速度的准确率。际资料中获取地层真实横波速度的准确率。

【技术实现步骤摘要】
基于密度聚类的偶极子声波测井频散校正方法及相关装置


[0001]本专利技术涉及地球物理声学测井领域，具体涉及一种基于密度聚类的偶极子声波测井频散校正方法及装置、计算设备及计算机存储介质。

技术介绍

[0002]偶极子声波测井采用偶极子激发，在井筒中激发挠曲波沿井内流体传播，偶极子波形具有明显的频散效应，时间域处理方法得到的时差结果是一段频率内的平均结果，无法反映真实的地层横波速度，故实际处理中需要对偶极子声波测井数据进行频散校正以得到地层真实的横波速度，有效的频散校正方法是偶极子声波测井有效应用的关键技术。
[0003]目前频散校正主要由获取模式波频散特性的频散分析方法和从频散特性中提取时差的方法构成。频散分析方法主要包括Prony方法、加权频谱相干法(WSS)、振幅相位估计法(APES)等，通过频散分析获得较高分辨率的频散曲线后，可用直方图统计方法和曲线拟合函数法对频散曲线进行频散校正。直方图统计方法是对从频率时差相关图中提取的频散曲线求时差的直方图，但数据会受噪声影响从而边缘使测量的结果有误差。曲线拟合函数法对频散曲线进行拟合，拟合后通过拟合曲线截止频率求取时差，但需要选择合适的拟合函数，对曲线质量要求较高。

技术实现思路

[0004]鉴于上述问题，提出了本专利技术以便提供一种克服上述频散校正准确率和效率较低问题的基于密度聚类的偶极子声波测井频散校正方法及装置、计算设备及计算机存储介质。
[0005]根据本专利技术的一个方面，提供了一种基于密度聚类的偶极子声波测井频散校正方法，包括：
[0006]步骤S1，对目标深度点的偶极子声波测井数据进行频散分析得到所述目标深度点对应的频散曲线；
[0007]步骤S2，重复执行所述步骤S1，直至获取到目标层段的频散曲线组；
[0008]步骤S3，利用密度聚类对所述频散曲线组进行处理，得到所述频散曲线组中的有效部分；
[0009]步骤S4，对所述频散曲线组中的有效部分进行拟合，得到校正结果。
[0010]在一种可选的方式中，所述频散曲线组的组数由所述目标层段的层段深度确定。
[0011]在一种可选的方式中，所述对目标深度点的偶极子声波测井数据进行频散分析得到所述目标深度点对应的频散曲线进一步包括：
[0012]利用加权频谱相干法对目标深度点的偶极子声波测井数据进行频散分析得到所述目标深度点对应的频散曲线。
[0013]在一种可选的方式中，所述利用密度聚类对所述频散曲线组进行处理，得到所述频散曲线组中的有效部分之前，所述方法还包括：
[0014]根据所述频散曲线组中的各个频散曲线的最大速度值、最小速度值、最大频率值、最小频率值，分别对所述各个频散曲线进行归一化处理。
[0015]在一种可选的方式中，所述根据所述频散曲线组中的各个频散曲线的最大速度值、最小速度值、最大频率值、最小频率值，分别对所述各个频散曲线进行归一化处理具体公式为：
[0016][0017]其中，x
min
为最小值；x
max
为最大值；x为原始数据；x
′
为归一化后的数据。
[0018]在一种可选的方式中，所述利用密度聚类对所述频散曲线组进行处理，得到所述频散曲线组中的有效部分进一步包括：
[0019]分别对所述频散曲线组中各个频散曲线设置半径和最小数目参数，利用DBSCAN聚类方法对所述各个频散曲线进行聚类分析，得到所述频散曲线组的有效部分。
[0020]在一种可选的方式中，所述对所述频散曲线组中的有效部分进行拟合，得到校正结果进一步包括：
[0021]对所述频散曲线组中各个频散曲线的有效频散点分别进行拟合得到多元线性回归模型，根据所述多元线性回归模型对所述各个频散曲线的有效频散点进行校正得到校正结果。
[0022]根据本专利技术的另一方面，提供了一种基于密度聚类的偶极子声波测井频散校正装置，包括：
[0023]频散分析模块，用于对目标深度点的偶极子声波测井数据进行频散分析得到所述目标深度点对应的频散曲线；重复执行，直至获取到目标层段的频散曲线组；
[0024]频散聚类模块，用于利用密度聚类对所述频散曲线组进行处理，得到所述频散曲线组中的有效部分；
[0025]频散校正模块，用于对所述频散曲线组中的有效部分进行拟合，得到校正结果。
[0026]根据本专利技术的又一方面，提供了一种计算设备，包括：处理器、存储器、通信接口和通信总线，所述处理器、所述存储器和所述通信接口通过所述通信总线完成相互间的通信；
[0027]所述存储器用于存放至少一可执行指令，所述可执行指令使所述处理器执行上述基于密度聚类的偶极子声波测井频散校正方法对应的操作。
[0028]根据本专利技术的再一方面，提供了一种计算机存储介质，所述存储介质中存储有至少一可执行指令，所述可执行指令使处理器执行如上述基于密度聚类的偶极子声波测井频散校正方法对应的操作。
[0029]本专利技术上述实施例提供的方案，对目标深度点的偶极子声波测井数据进行频散分析得到所述目标深度点对应的频散曲线；重复执行直至获取到目标层段的频散曲线组；利用密度聚类对所述频散曲线组进行处理，得到所述频散曲线组中的有效部分；对所述频散曲线组中的有效部分进行拟合，得到校正结果。本专利技术利用密度聚类对频散曲线进行聚类筛选，无需预先进行大量数据的训练，提升了有效数据的提取效率；通过对频散曲线中的有效部分进行拟合得到校正结果，在保证频散校正效率的同时提高了从实际资料中获取地层真实横波速度的准确率。
[0030]上述说明仅是本专利技术技术方案的概述，为了能够更清楚了解本专利技术的技术手段，
而可依照说明书的内容予以实施，并且为了让本专利技术的上述和其它目的、特征和优点能够更明显易懂，以下特举本专利技术的具体实施方式。
附图说明
[0031]通过阅读下文优选实施方式的详细描述，各种其他的优点和益处对于本领域普通技术人员将变得清楚明了。附图仅用于示出优选实施方式的目的，而并不认为是对本专利技术的限制。而且在整个附图中，用相同的参考符号表示相同的部件。在附图中：
[0032]图1示出了本专利技术实施例的基于密度聚类的偶极子声波测井频散校正方法的流程示意图；
[0033]图2示出了本专利技术另一个实施例的基于密度聚类的偶极子声波测井频散校正方法的流程示意图；
[0034]图3示出了本专利技术实施例的密度聚类示意图；
[0035]图4示出了本专利技术实施例的某一深度点的频散曲线示意图；
[0036]图5示出了本专利技术实施例的基于密度聚类的偶极子声波测井频散校正装置的结构示意图；
[0037]图6示出了本专利技术实施例的计算设备的结构示意图。
具体实施方式
[0038]下面将参照附图更详细地描述本专利技术的示例性实施例。虽然附图中显示了本专利技术的示例性实施例，然而应当理解，可以以各种形式实现本专利技术而不应被这里阐述的实施例所限制。相反，提供这些实施例是为了能够更本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于密度聚类的偶极子声波测井频散校正方法，其特征在于，包括：步骤S1，对目标深度点的偶极子声波测井数据进行频散分析得到所述目标深度点对应的频散曲线；步骤S2，重复执行所述步骤S1，直至获取到目标层段的频散曲线组；步骤S3，利用密度聚类对所述频散曲线组进行处理，得到所述频散曲线组中的有效部分；步骤S4，对所述频散曲线组中的有效部分进行拟合，得到校正结果。2.根据权利要求1所述的基于密度聚类的偶极子声波测井频散校正方法，其特征在于，所述频散曲线组的组数由所述目标层段的层段深度确定。3.根据权利要求1所述的基于密度聚类的偶极子声波测井频散校正方法，其特征在于，所述对目标深度点的偶极子声波测井数据进行频散分析得到所述目标深度点对应的频散曲线进一步包括：利用加权频谱相干法对目标深度点的偶极子声波测井数据进行频散分析得到所述目标深度点对应的频散曲线。4.根据权利要求1所述的基于密度聚类的偶极子声波测井频散校正方法，其特征在于，所述利用密度聚类对所述频散曲线组进行处理，得到所述频散曲线组中的有效部分之前，所述方法还包括：根据所述频散曲线组中的各个频散曲线的最大速度值、最小速度值、最大频率值、最小频率值，分别对所述各个频散曲线进行归一化处理。5.根据权利要求4所述的基于密度聚类的偶极子声波测井频散校正方法，其特征在于，所述根据所述频散曲线组中的各个频散曲线的最大速度值、最小速度值、最大频率值、最小频率值，分别对所述各个频散曲线进行归一化处理具体公式为：其中，x
min
为最小值；x
max
为最大值；x为原始数据；x
′
为归一化后的数据。6.根据权利要求1
‑
5中任一项...

【专利技术属性】
技术研发人员：张伟，钱玉萍，徐大年，张国华，崔维平，张传举，谢慧卓，丛平，侯振学，成家杰，张宏伟，范川，王俊华，
申请(专利权)人：中海油田服务股份有限公司，
类型：发明
国别省市：