一种确定叠加速度的方法、装置及系统制造方法及图纸

本申请实施例公开了一种确定叠加速度的方法、装置及系统，所述方法包括获取地震记录的地震反射数据，根据所述地震反射数据计算局部同相轴的斜率；根据所述局部同相轴的斜率计算得到零偏距时间与叠加速度的关系数据；计算确定的参考道与所述关系数据对应的地震道的局部相似系数；根据所述参考道与所述关系数据对应的地震道的局部相似系数确定所述关系数据的相似系数；将所述关系数据的相似系数作为所述关系数据的权重，根据加权聚类算法计算带有权重的关系数据点的中心点；根据所述中心点的速度值确定所述地震记录的叠加速度。利用本申请各个实施例，可以提高获得的叠加速度的精确度。

【技术实现步骤摘要】
一种确定叠加速度的方法、装置及系统
本专利技术涉及地震信号处理
，特别地，涉及一种确定叠加速度的方法、装置及系统。
技术介绍
速度估计是地震数据处理与成像中很重要的一个环节。一个准确的速度模型对于将界面反射时间转化为界面位置的成像至关重要。传统的速度分析方法是在速度谱上拾取速度，速度谱通过应用不同速度进行动校正，并计算它们对应的相干性。速度的拾取过程需要大量的人为工作而且伴随着人为误差，为了减少人为工作，大量的自动速度分析的方法被提出。目前，大部分自动速度分析方法在相干扫描速度谱上进行计算和拾取叠加速度，但是相干扫描自动拾取的方法由于速度谱存在多个局部极值，自动拾取速度精度不高。地震记录的同相轴的斜率因其包含了地层的速度信息，在自动估计叠加速度上具有较高的精确度。但地表与地下构造的复杂多变、数据激发和接收条件差等问题导致地震资料信噪比低，地震噪声的存在严重影响了直接通过同相轴斜率来估计叠加速度的精确度。因此，业内需要一种更加准确的估计叠加速度的方法。
技术实现思路
本申请实施例的目的在于提供一种确定叠加速度的方法、装置及系统，可以避免地震噪声的影响，提高获得的叠加速度的精确度。为达到上述目的，一方面，本申请实施例提供了一种确定叠加速度的方法，包括：获取地震记录的地震反射数据，根据所述地震反射数据计算局部同相轴的斜率；根据所述局部同相轴的斜率计算得到零偏距时间与叠加速度的关系数据；计算确定的参考道与所述关系数据对应的地震道的局部相似系数；根据所述参考道与所述关系数据对应的地震道的局部相似系数确定所述关系数据的相似系数；将所述关系数据的相似系数作为所述关系数据的权重，根据加权聚类算法计算带有权重的关系数据点的中心点；根据所述中心点的速度值确定所述地震记录的叠加速度。本申请实施例的确定叠加速度的方法，所述根据所述地震反射数据计算局部同相轴的斜率，包括：基于平面波分解算法计算局部同相轴斜率。本申请实施例的确定叠加速度的方法，所述将所述关系数据的相似系数作为所述关系数据的权重，根据加权聚类算法计算带有权重的关系数据点的中心点，包括：从所述关系数据点中随机选取K个点作为初始中心点；执行迭代步骤：将所述关系数据点根据上一次获得的中心点进行分类，确定所述关系数据点所属的类别基于下述公式根据本次确定的类别重新确定所述带有权重的关系数据点的中心点：其中，xp表示第p个关系数据点，wp表示xp的权重，Qi表示xp所属于的类别，μi表示Qi类的中心，iter表示迭代次数；当与的欧式距离小于预设阈值时，结束迭代；输出最后一次迭代的中心点确定所述最后一次迭代的中心点为所述关系数据点的中心点。本申请实施例的确定叠加速度的方法，所述根据所述中心点确定所述地震记录的叠加速度，包括：对所述中心点的速度值进行插值处理，确定所述地震记录的叠加速度。另一方面，本申请实施例还提供一种确定叠加速度的装置，包括：斜率计算模块，用于获取地震记录的地震反射数据，根据所述地震反射数据计算局部同相轴的斜率；速度计算模块，用于根据所述局部同相轴的斜率计算得到零偏距时间与叠加速度的关系数据；相似系数计算模块，用于计算确定的参考道与所述关系数据对应的地震道的局部相似系数，根据所述参考道与所述关系数据对应的地震道的局部相似系数确定所述关系数据的相似系数；叠加速度确定模块，用于将所述关系数据的相似系数作为所述关系数据的权重，根据加权聚类算法计算带有权重的关系数据点的中心点，根据所述中心点的速度值确定所述地震记录的叠加速度。本申请实施例的确定叠加速度的装置，所述斜率计算模块，包括：根据所述地震反射数据基于平面波分解算法计算局部同相轴斜率。本申请实施例的确定叠加速度的装置，所述叠加速度确定模块，包括：从所述关系数据点中随机选取K个点作为初始中心点；执行迭代步骤：将所述关系数据点根据上一次获得的中心点进行分类，确定所述关系数据点所属的类别基于下述公式根据本次确定的类别重新确定所述带有权重的关系数据点的中心点：其中，xp表示第p个关系数据点，wp表示xp的权重，Qi表示xp所属于的类别，μi表示Qi类的中心，iter表示迭代次数；当与的欧式距离小于预设阈值时，结束迭代；输出最后一次迭代的中心点确定所述最后一次迭代的中心点为所述关系数据点的中心点。本申请实施例的确定叠加速度的装置，对所述中心点的速度值进行插值处理，确定所述地震记录的叠加速度。本申请实施例的确定叠加速度的装置，包括处理器及用于存储处理器可执行指令的存储器，所述指令被所述处理器执行时实现包括以下步骤：获取地震记录的地震反射数据，根据所述地震反射数据计算局部同相轴的斜率；根据所述局部同相轴的斜率计算得到零偏距时间与叠加速度的关系数据；计算确定的参考道与所述关系数据对应的地震道的局部相似系数；根据所述参考道与所述关系数据对应的地震道的局部相似系数确定所述关系数据的相似系数；将所述关系数据的相似系数作为所述关系数据的权重，根据加权聚类算法计算带有权重的关系数据点的中心点；根据所述中心点的速度值确定所述地震记录的叠加速度。另一方面，本申请实施例还提供一种确定叠加速度的系统，包括上述任意一项所述的确定叠加速度的装置。本说明书一个或多个实施例提供的一种确定叠加速度的方法、装置及系统，可以通过局部同相轴斜率计算叠加速度与零偏距时间的关系数据，计算确定的参考道与所述关系数据对应的地震道的局部相似系数，根据所述参考道与所述关系数据对应的地震道的局部相似系数确定所述关系数据的相似系数。所述关系数据的相似系数越大，对应的地震道局部位置的信噪比越大，相应的叠加速度值越精确。将所述关系数据的相似系数作为所述关系数据的权重，计算带有权重的关系数据点的中心点，根据所述中心点的速度值确定所述地震记录的叠加速度。利用本实施例的方案，可以提高获得的叠加速度的精确度。附图说明为了更清楚地说明本申请实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本申请中记载的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。在附图中：图1为本说明书提供的一种确定叠加速度的方法实施例的流程示意图；图2为本说明书提供的一个实施例中实际数据的一个CMP道集的示意图(a)、以及其局部同相轴斜率示意图(b)；图3为本说明书提供的一个实施例中带有权重的数据点对应的叠加速度分布图(a)、以及利用加权聚类算法计算得到的叠加速度点和插值后的叠加速度曲线图(b)；图4为本说明书提供的一个实施例中通过本说明书实施例提供方法获得的叠加速度场示意图；图5为本说明书提供的一个实施例中通过常规方法确定的叠加速度进行动校正叠加的结果示意图(a)、以及通过本说明书实施例确定的叠加速度进行动校正叠加的结果示意图(b)；图6为本说明书提供的一个实施例中图5(a)所示矩形框中的放大剖面示意图(a)、以及图5(b)所示矩形框中的放大剖面示意图(b)；图7为本说明书提供的一种确定叠加速度的装置实施例的模块结构示意图。具体实施方式为了使本
的人员更好地理解本说明书中的技术方案，下面将结合本说明书一个或多个实施例中的附图，对本说明书一个或多个实施例中的技术方案进行清楚、本文档来自技高网...

【技术保护点】
一种确定叠加速度的方法，其特征在于，包括：获取地震记录的地震反射数据，根据所述地震反射数据计算局部同相轴的斜率；根据所述局部同相轴的斜率计算得到零偏距时间与叠加速度的关系数据；计算确定的参考道与所述关系数据对应的地震道的局部相似系数；根据所述参考道与所述关系数据对应的地震道的局部相似系数确定所述关系数据的相似系数；将所述关系数据的相似系数作为所述关系数据的权重，根据加权聚类算法计算带有权重的关系数据点的中心点；根据所述中心点的速度值确定所述地震记录的叠加速度。

【技术特征摘要】
1.一种确定叠加速度的方法，其特征在于，包括：获取地震记录的地震反射数据，根据所述地震反射数据计算局部同相轴的斜率；根据所述局部同相轴的斜率计算得到零偏距时间与叠加速度的关系数据；计算确定的参考道与所述关系数据对应的地震道的局部相似系数；根据所述参考道与所述关系数据对应的地震道的局部相似系数确定所述关系数据的相似系数；将所述关系数据的相似系数作为所述关系数据的权重，根据加权聚类算法计算带有权重的关系数据点的中心点；根据所述中心点的速度值确定所述地震记录的叠加速度。2.根据权利要求1所述的确定叠加速度的方法，其特征在于，所述根据所述地震反射数据计算局部同相轴的斜率，包括：基于平面波分解算法计算局部同相轴斜率。3.根据权利要求1所述的确定叠加速度的方法，其特征在于，所述将所述关系数据的相似系数作为所述关系数据的权重，根据加权聚类算法计算带有权重的关系数据点的中心点，包括：从所述关系数据点中随机选取K个点作为初始中心点；执行迭代步骤：将所述关系数据点根据上一次获得的中心点进行分类，确定所述关系数据点所属的类别基于下述公式根据本次确定的类别重新确定所述带有权重的关系数据点的中心点：其中，xp表示第p个关系数据点，wp表示xp的权重，Qi表示xp所属于的类别，μi表示Qi类的中心，iter表示迭代次数；当与的欧式距离小于预设阈值时，结束迭代；输出最后一次迭代的中心点确定所述最后一次迭代的中心点为所述关系数据点的中心点。4.根据权利要求1所述的确定叠加速度的方法，其特征在于，所述根据所述中心点确定所述地震记录的叠加速度，包括：对所述中心点的速度值进行插值处理，确定所述地震记录的叠加速度。5.一种确定叠加速度的装置，其特征在于，包括：斜率计算模块，用于获取地震记录的地震反射数据，根据所述地震反射数据计算局部同相轴的斜率；速度计算模块，用于根据所述局部同相轴的斜率计算得到零偏距时间与叠加速度的关系数据；相似系数计算模块，用于计算确定的参考道与所述关系数据对应的地震道的局部相似系数，根据所述参考...

【专利技术属性】
技术研发人员：刘国昌，李超，陈小宏，马继涛，刘兴业，
申请(专利权)人：中国石油大学北京，
类型：发明
国别省市：北京,11