一种地层品质因子反演方法及装置制造方法及图纸

本申请公开了一种地层品质因子反演方法及装置，该方法包括：获取检波点的垂直地震剖面VSP资料；对每个检波点的VSP资料进行非零偏校正，确定每个检波点的零偏VSP资料初至时间；根据每个检波点的零偏VSP资料初至时间和深度，确定每相邻两个检波点之间的层速度；根据每个检波点的VSP资料进行傅里叶变换，确定波场振幅谱，根据所述波场振幅谱确定每个检波点的质心频率；根据每个检波点的深度和质心频率，确定每相邻两个检波点之间的平均吸收系数；根据每相邻两个检波点之间的平均吸收系数和层速度，确定每相邻两个检波点之间的品质因子Q值。从而可以计算地下介质层之间的品质因子。

A Method and Device for Inversion of Formation Quality Factor

【技术实现步骤摘要】
一种地层品质因子反演方法及装置
本申请涉及地震资料处理
，尤其涉及一种地层品质因子反演方法及装置。
技术介绍
本部分旨在为权利要求书中陈述的本申请实施例提供背景或上下文。此处的描述不因为包括在本部分中就承认是现有技术。地震波在粘弹性非均匀的地下介质中传播时，由于受到地层的吸收，其不同频率成分的能量要受到不同程度的衰减损失。从地层的吸收程度上讲，吸收主要分为近地表吸收和潜水面以下固结岩层吸收两大类。而松散的近地表介质吸收了绝大部分的有效地震信号，相对而言地层对地震波高频成分的吸收越厉害，这是导致地震分辨率不高的主要原因之一，严重制约了目的层地震资料的品质，所以解决好近地表吸收补偿问题对石油地震勘探有着重要意义。近地表吸收补偿的关键是要获得比较准确的地层品质因子Q值，因此，如何准确求取地层的品质因子Q一直是物理学家研究的问题。目前主流的估算Q值的方法只能计算地层的平均品质因子，没有考虑到地下介质层间的品质因子的变化。
技术实现思路
本申请实施例提供一种地层品质因子反演方法及装置，用以计算地下介质层间的品质因子。第一方面，该地层品质因子反演方法包括：获取检波点的垂直地震剖面VSP资料；对每个检波点的VSP资料进行非零偏校正，确定每个检波点的零偏VSP资料初至时间；根据每个检波点的零偏VSP资料初至时间和深度，确定每相邻两个检波点之间的层速度；根据每个检波点的VSP资料进行傅里叶变换，确定波场振幅谱，根据所述波场振幅谱确定每个检波点的质心频率；根据每个检波点的深度和质心频率，确定每相邻两个检波点之间的平均吸收系数；根据每相邻两个检波点之间的平均吸收系数和层速度，确定每相邻两个检波点之间的品质因子Q值。第二方面，还提供一种地层品质因子反演装置，该装置包括：获取模块，用于获取检波点的垂直地震剖面VSP资料；校正模块，用于对每个检波点的VSP资料进行非零偏校正，确定每个检波点的零偏VSP资料初至时间；层速度确定模块，用于根据每个检波点的零偏VSP资料初至时间和深度，确定每相邻两个检波点之间的层速度；质心频率确定模块，用于根据每个检波点的VSP资料进行傅里叶变换，确定波场振幅谱，根据所述波场振幅谱确定质心频率；平均吸收系数确定模块，用于根据每个检波点的深度和质心频率，确定每相邻两个检波点之间的平均吸收系数；品质因子确定模块，用于根据每相邻两个检波点之间的平均吸收系数和层速度，确定每相邻两个检波点之间的品质因子Q值。第三方面，提供一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现上述第一方面的所述方法。第四方面，提供一种计算机可读存储介质，所述计算机可读存储介质存储有执行上述第一方面的所述方法的计算机程序。本申请实施例中，通过获取检波点的垂直地震剖面VSP资料；对每个检波点的VSP资料进行非零偏校正，确定每个检波点的零偏VSP资料初至时间；根据每个检波点的零偏VSP资料初至时间和深度，确定每相邻两个检波点之间的层速度；根据每个检波点的VSP资料进行傅里叶变换，确定波场振幅谱，根据所述波场振幅谱确定每个检波点的质心频率；根据每个检波点的深度和质心频率，确定每相邻两个检波点之间的平均吸收系数；根据每相邻两个检波点之间的平均吸收系数和层速度，确定每相邻两个检波点之间的品质因子Q值。从而可以计算地下介质层之间的品质因子。附图说明为了更清楚地说明本申请实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本申请的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。在附图中：图1为本申请实施例中提供的一种地层品质因子反演方法流程示意图；图2为本申请实施例中提供的一种三层模型的正演模型参数表；图3为本申请实施例中提供的一种Q值反演方法反演得到的Q随地层深度的变化示意图；图4为本申请实施例中提供的某VSP井资料的原始剖面示意图；图5为本申请实施例中提供的零偏VSP井资料预处理后提取出的下行直达波场剖面示意图；图6为本申请实施例中提供的某VSP井的层速度变化曲线示意图；图7为本申请实施例中提供的某VSP井主频随深度的变化曲线；图8为本申请实施例中提供的某VSP井的基于地质层位约束的Q值随深度变化曲线；图9为本申请实施例中提供的某VSP井未进行反Q滤波的时间偏移剖面示意图；图10为本申请实施例中提供的某VSP井进行了反Q滤波后的时间偏移剖面示意图；图11为本申请实施例中提供的反Q滤波补偿前后的叠前时间偏移剖面的振幅谱对比示意图；图12为本申请实施例提供的一种地层品质因子反演装置示意图。具体实施方式下面将结合本申请实施例中的附图，对本申请实施例中的技术方案进行清楚、完整地描述，所描述的实施例仅仅是本申请一部分实施例，而不是全部的实施例。基于本申请中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本申请保护的范围。下面参考本申请的若干代表性实施方式，详细阐释本申请的原理和精神。虽然本申请提供了如下述实施例或附图所示的方法操作步骤或装置结构，但基于常规或者无需创造性的劳动在所述方法或装置中可以包括更多或者更少的操作步骤或模块单元。在逻辑性上不存在必要因果关系的步骤或结构中，这些步骤的执行顺序或装置的模块结构不限于本申请实施例或附图所示的执行顺序或模块结构。所述的方法或模块结构的在实际中的装置或终端产品应用时，可以按照实施例或者附图所示的方法或模块结构进行顺序执行或者并行执行。垂直地震剖面(VerticalSeismicProfiling，VSP)是一种地震观测方法，与通常地面观测的地震剖面相对应。在地表附近的一些点上激发地震波，在沿井孔不同深度布置的一些多级多分量的检波点上进行观测。在垂直地震剖面中，因为检波器通过井置于地层内部，所以不仅能接收到自下而上传播的上行纵波和上行转换波，也能接收到自上而下传播的下行纵波及下行转换波，甚至能接收到横波。由于垂直地震剖面为地面激发、井中垂直方向接收的观测方式，其下行直达波受噪声的影响较小，较地面地震记录能更精确地反映地层的滤波作用，因此利用VSP直达下行波反演地层品质因子Q是最为有效的途径。在本实施方式中，执行所述地层品质因子反演的客体可以是具有逻辑运算功能的电子设备。所述电子设备可以是服务器和客户端。所述客户端可以为台式电脑、平板电脑、笔记本电脑、工作站等。当然，客户端并不限于上述具有一定实体的电子设备，其还可以为运行于上述电子设备中的软体。还可以是一种通过程序开发形成的程序软件，该程序软件可以运行于上述电子设备中。请参阅图1。图1示出了本申请实施例提供的一种地层品质因子反演方法的流程图。所述方法可以包括以下步骤。步骤101：获取检波点的垂直地震剖面VSP资料。步骤102：对每个检波点的VSP资料进行非零偏校正，确定每个检波点的零偏VSP资料初至时间。步骤103：根据每个检波点的零偏VSP资料初至时间和深度，确定每相邻两个检波点之间的层速度。步骤104：根据每个检波点的VSP资料进行傅里叶变换，确定波场振幅谱，根据所述波场振幅谱确定每个检波点的质心频率。步骤105：根据本文档来自技高网...

【技术保护点】
1.一种地层品质因子反演方法，其特征在于，包括：获取检波点的垂直地震剖面VSP资料；对每个检波点的VSP资料进行非零偏校正，确定每个检波点的零偏VSP资料初至时间；根据每个检波点的零偏VSP资料初至时间和深度，确定每相邻两个检波点之间的层速度；根据每个检波点的VSP资料进行傅里叶变换，确定波场振幅谱，根据所述波场振幅谱确定每个检波点的质心频率；根据每个检波点的深度和质心频率，确定每相邻两个检波点之间的平均吸收系数；根据每相邻两个检波点之间的平均吸收系数和层速度，确定每相邻两个检波点之间的品质因子Q值。

【技术特征摘要】
1.一种地层品质因子反演方法，其特征在于，包括：获取检波点的垂直地震剖面VSP资料；对每个检波点的VSP资料进行非零偏校正，确定每个检波点的零偏VSP资料初至时间；根据每个检波点的零偏VSP资料初至时间和深度，确定每相邻两个检波点之间的层速度；根据每个检波点的VSP资料进行傅里叶变换，确定波场振幅谱，根据所述波场振幅谱确定每个检波点的质心频率；根据每个检波点的深度和质心频率，确定每相邻两个检波点之间的平均吸收系数；根据每相邻两个检波点之间的平均吸收系数和层速度，确定每相邻两个检波点之间的品质因子Q值。2.如权利要求1所述的方法，其特征在于，所述根据每个检波点的VSP深度和质心频率，确定每相邻两个检波点之间的平均吸收系数，包括：按照以下公式确定所述平均吸收系数：其中，αi表示第i个检波点与第i+1个检波点之间的平均吸收系数，σi2表示第i个检波点的方差，fi表示第i个检波点的质心频率，fi+1表示第i+1个检波点的质心频率，zi表示第i个检波点的深度，zi+1表示第i+1个检波点的深度，Δzi表示第i个检波点到第i+1个检波点之间的深度差，Δfi表示第i个检波点到第i+1个检波点之间的质心频率差。3.如权利要求1所述的方法，其特征在于，所述根据每相邻两个检波点之间的平均吸收系数和层速度，确定每相邻两个检波点之间的品质因子Q值，包括：按照以下公式确定所述Q值：其中，Qi表示第i个检波点到第i+1个检波点之间的品质因子，αi表示第i个检波点与第i+1个检波点之间的平均吸收系数，vi表示第i个检波点到第i+1个检波点之间的层速度。4.如权利要求1至3任一项所述的方法，其特征在于，所述根据每个检波点的VSP深度和质心频率，确定每两个检波点之间的平均吸收系数，包括：对所述零偏VSP资料进行几何扩散补偿以及波场分离，得到下行波场；根据所述下行波场中的每个检波点的VSP深度和质心频率，确定每相邻两个检波点之间的平均吸收系数；所述根据每相邻两个检波点之间的平均吸收系数和层速度，确定每相邻两个检波点之间的品质因子Q值，包括：根据所述下行波场中的每相邻两个检波点之间的平均吸收系数和层速度，确定每相邻两个检波点之间的Q值。5.一种地层品质因子反演装置，其特征在于，包括：获取模块，用于获取检波...

【专利技术属性】
技术研发人员：李斐，杨哲，赵玉莲，肖明图，刘伟明，谢俊法，
申请(专利权)人：中国石油天然气股份有限公司，
类型：发明
国别省市：北京,11