CN109765615A - 一种地层品质因子反演方法及装置 - Google Patents

Abstract

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

Description

一种地层品质因子反演方法及装置

技术领域

本申请涉及地震资料处理技术领域，尤其涉及一种地层品质因子反演方法及装置。

背景技术

本部分旨在为权利要求书中陈述的本申请实施例提供背景或上下文。此处的描述不因为包括在本部分中就承认是现有技术。

地震波在粘弹性非均匀的地下介质中传播时，由于受到地层的吸收，其不同频率成分的能量要受到不同程度的衰减损失。从地层的吸收程度上讲，吸收主要分为近地表吸收和潜水面以下固结岩层吸收两大类。而松散的近地表介质吸收了绝大部分的有效地震信号，相对而言地层对地震波高频成分的吸收越厉害，这是导致地震分辨率不高的主要原因之一，严重制约了目的层地震资料的品质，所以解决好近地表吸收补偿问题对石油地震勘探有着重要意义。

近地表吸收补偿的关键是要获得比较准确的地层品质因子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为本申请实施例提供的一种地层品质因子反演装置示意图。

具体实施方式

下面将结合本申请实施例中的附图，对本申请实施例中的技术方案进行清楚、完整地描述，所描述的实施例仅仅是本申请一部分实施例，而不是全部的实施例。基于本申请中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本申请保护的范围。

下面参考本申请的若干代表性实施方式，详细阐释本申请的原理和精神。

虽然本申请提供了如下述实施例或附图所示的方法操作步骤或装置结构，但基于常规或者无需创造性的劳动在所述方法或装置中可以包括更多或者更少的操作步骤或模块单元。在逻辑性上不存在必要因果关系的步骤或结构中，这些步骤的执行顺序或装置的模块结构不限于本申请实施例或附图所示的执行顺序或模块结构。所述的方法或模块结构的在实际中的装置或终端产品应用时，可以按照实施例或者附图所示的方法或模块结构进行顺序执行或者并行执行。

垂直地震剖面(Vertical Seismic Profiling，VSP)是一种地震观测方法，与通常地面观测的地震剖面相对应。在地表附近的一些点上激发地震波，在沿井孔不同深度布置的一些多级多分量的检波点上进行观测。在垂直地震剖面中，因为检波器通过井置于地层内部，所以不仅能接收到自下而上传播的上行纵波和上行转换波，也能接收到自上而下传播的下行纵波及下行转换波，甚至能接收到横波。由于垂直地震剖面为地面激发、井中垂直方向接收的观测方式，其下行直达波受噪声的影响较小，较地面地震记录能更精确地反映地层的滤波作用，因此利用VSP直达下行波反演地层品质因子Q是最为有效的途径。

在本实施方式中，执行所述地层品质因子反演的客体可以是具有逻辑运算功能的电子设备。所述电子设备可以是服务器和客户端。所述客户端可以为台式电脑、平板电脑、笔记本电脑、工作站等。当然，客户端并不限于上述具有一定实体的电子设备，其还可以为运行于上述电子设备中的软体。还可以是一种通过程序开发形成的程序软件，该程序软件可以运行于上述电子设备中。

请参阅图1。图1示出了本申请实施例提供的一种地层品质因子反演方法的流程图。所述方法可以包括以下步骤。

步骤101：获取检波点的垂直地震剖面VSP资料。

步骤102：对每个检波点的VSP资料进行非零偏校正，确定每个检波点的零偏VSP资料初至时间。

步骤103：根据每个检波点的零偏VSP资料初至时间和深度，确定每相邻两个检波点之间的层速度。

步骤104：根据每个检波点的VSP资料进行傅里叶变换，确定波场振幅谱，根据所述波场振幅谱确定每个检波点的质心频率。

步骤105：根据每个检波点的深度和质心频率，确定每相邻两个检波点之间的平均吸收系数。

步骤106：根据每相邻两个检波点之间的平均吸收系数和层速度，确定每相邻两个检波点之间的品质因子Q值。

在步骤102中，对地震波的初至拾取即计算地震波到检波点的时间，由于VSP井资料炮点位置跟井口位置有一段距离，需要通过非零偏校正从而得到每个检波点的零偏VSP资料初至时间。

在步骤103之前，需对VSP井资料进行预处理，为了保持数据的振幅相对关系不改变、频率和相位特征不损失，本申请实施例采用包括但不限于球面扩散补偿法对VSP井资料进行处理，旨在去除地震波由于几何扩散所引起的振幅损失，从而恢复由于地下介质因素引起的振幅相对变化关系。利用分频压制面波也可以消除地滚波强能量对有用数据的影响。

上述对VSP资料进行预处理的步骤，旨在不改变数据的振幅相对关系、数据的频率和相位特征。上述方法仅为示例，本申请实施例不对此做具体限定。

在步骤103中可以包括直达波初至时间拾取、再计算平均速度，进而可以求取每一个检波距段的层速度。其中，根据初至时间和深度，求取平均速度以及层速度。深度为检波器在井中的深度。

可选地，在步骤103之后，通过VSP资料的解释成果中提供的层速度数据表与步骤103中求得的层速度数据进行比对，这两组层速度数据在极浅层处往往有偏差，而井点层速度的准确度直接影响了井点Q值的精度，因此需要综合对比分析，最终得到一套相对准确的层速度数据。

考虑到地震波向下传播的过程中，零偏VSP资料的每个检波点接收到的下行波场，同时可以当做是下一个检波点的震源波场，本申请实施例提供的品质因子反演方法可以根据零偏VSP资料逐道反演任意两个检波点之间地层的时变Q值，这样的Q值计算方法具有更高的分辨率。因此在步骤105中，可以按照以下公式(1)确定所述平均吸收系数：

其中，α_i表示第i个检波点与第i+1个检波点之间的平均吸收系数，σ_i ²表示第i个检波点的方差，f_i表示第i个检波点的质心频率，f_i+1表示第i+1个检波点的质心频率，z_i表示第i个检波点的深度，z_i+1表示第i+1个检波点的深度，Δz_i表示第i个检波点到第i+1个检波点之间的深度差，Δf_i表示第i个检波点到第i+1个检波点之间的质心频率差。

进一步在步骤106中，可以按照以下公式(2)确定所述Q值：

其中，Q_i表示第i个检波点到第i+1个检波点之间的品质因子，α_i表示第i个检波点与第i+1个检波点之间的平均吸收系数，v_i表示第i个检波点到第i+1个检波点之间的层速度。

在一种可能的实施方式中，可以首先对零偏VSP资料进行几何扩散补偿以及波场分离，得到下行波场；根据所述下行波场中的每个检波点的VSP深度和质心频率，确定每相邻两个检波点之间的平均吸收系数。进一步地，根据所述下行波场中的每相邻两个检波点之间的平均吸收系数和层速度，确定每相邻两个检波点之间的Q值。

其中通过波场分离提取VSP井资料高品质的下行波场的方法包括但不限于中值滤波法，中值滤波法可以较好的消除地震记录中随机噪声，最大限度保护有效信息。由于上行波和下行波的视速度不同，下行波速度为正值，上行波的速度为负值，可以通过此差异来进行波场分离。当沿着深度方向进行中值滤波时，下行波的强振幅保持不变，而上行波强振幅将被滤掉。从而可以提取出VSP资料中品质较高的下行波场。

图2示出了本申请实施例提供的一个三层模型的正演模型参数表，其震源位于地表，检波器设置间距为10米，共60道。根据各检波点间Q值与深度关系，再结合实际地层情况，根据各检波点的层速度情况对检波点进行分层，便于后续数据的处理。

图3示出了本申请实施例提供的反演得到的Q值随地层深度的变化示意图。通过本申请实施例提供的品质因子反演方法，对每一口VSP井的下行波场利用离散化频移法，计算出每一个检波器位置处的波场主频极其方差，反演出该井点每一个检波距段的地层Q值。可以看出反演值与理论值较为吻合。

进一步，可以对工区所有VSP井资料与地面地震数据井震标定，利用VSP井合成记录或者VSP走廊叠加剖面，进行时间和深度的标定，使得时间域地震数据中包含深度信息，相当于为VSP井地质分层和地震数据解释层位建立桥梁，便于VSP井资料求取的Q值能够准确应用在地震数据上，与地震数据匹配。

进一步，可以结合工区地震解释层位，对每口VSP井计算的地层Q值进行地质层位划分，并利用地质层位约束，以及层速度趋势对地质层位分界面处存在异常极值的井口地层Q值进行适当删减平滑；以及根据Q值与深度关系，结合实际地层情况，参考层速度的分层情况对Q值曲线进行分层处理，使其与层速度变化趋势相符合；最终计算得到每个地质层位内的平均Q值，从而获得较为准确的地层Q值。

进一步，可以建立全工区近地表Q场。具体地，在计算得到每个VSP井点地质层位平均Q值基础上，将工区所有VSP井点解释后的地层Q值结果运用克里金方法进行插值内插，并做适当平滑，最终建立全工区近地表Q值场。

反Q滤波是将因Q因素而衰减损耗的地震波的能量补充回来，反Q滤波可看作是地震波的逆向传播或者是偏移过程，将地表记录到的数据反向延拓并剔除地球介质对入射地震波所产生的Q值滤波作用。进一步，可以利用Futterman地球衰减频率恒定Q模型，基于零偏VSP资料得到的近地表Q值场，应用基于波动方程的反Q滤波算法，对地面地震数据进行反Q滤波，最终补偿地震数据由于近地表吸收引起的高频信息损失。

本申请实施例还提供了一个具体的实施例来对上述地层品质因子反演的方法进行详细说明，然而，该具体的实施例仅是为了更好地说明本申请实施例，并不构成对本申请实施例的不当限定。

将该方法在柴达木盆地腹部某工区进行了实际资料处理应用。该区地表条件复杂，地貌类型主要有沙漠、盐碱地、沼泽、草场戈壁和少量山地，表层的复杂多变的地貌对地震资料的吸收衰减较为明显，研究目标区资料分辨率较低，高频损失明显，因此需要开展近地表吸收补偿处理，恢复地震资料的高频信息。

工区内共有VSP井13口，井施工方式为：检波点之间距离10m，首道到井口的距离为30m，激发炮距离井口59.5m，接收点深度为从1350m到2400米不等，偏移距与接收点深度之比远小于十分之一，进行零偏校正之后可当做零偏的VSP资料进行处理。

利用本申请实施例中的地层品质因子反演方法，对该工区进行了近地表吸收衰减补偿和分辨率提高处理，具体实施方式如下：

第一步，对获取的13口VSP资料逐一进行炮偏分析，通过非零偏校正得到每口井准确的零偏VSP井资料初至到达时间。

第二步，对每口VSP井资料进行预处理，该工区VSP井资料信噪比较低，图4是本申请实施例提供的某口VSP井资料的原始剖面示意图；选用基于地表一致性的分频减去法来压制面波，并对去噪前后以及去除的面波在频率域进行监控，保证去噪过程有效信号的振幅频率以及相位特征无损伤，最后通过中值滤波波场分离法提取出每口VSP井高品质的下行直达波场。图5示出了本申请实施例提供的零偏VSP井资料预处理后提取出的下行直达波场剖面示意图。

第三步，通过上述VSP井资料预处理计算得到每口井层速度数据，部分VSP井提供了层速度数据表，将两者进行比对。例如100米以上的极浅层速度差异比较大，VSP井解释成果提供的层速度均比处理计算的值要小，再往下的层速度越来越趋近，由于处理计算过程消除了非零偏影响，可以使得结果更为准确。最终综合两者得到每口VSP井处相对准确的层速度数据，图6示出了本申请实施例提供的某口VSP井的层速度随深度的变化曲线。

第四步，对每一口VSP井下行直达波场，计算出每一个检波器接收到的波场主频及其方差，图7示出了某口VSP井主频随深度的变化曲线。利用离散化频移法，计算出该井点每一检波距段(每相邻两个检波点)的层间Q值。进一步，利用VSP走廊叠加剖面对VSP井跟井旁地震数据进行井震时间深度标定，并根据地震解释层位对计算得到的井点地层Q值进行地质层位划分，并利用地质层位约束以及层速度趋势，对地质层位分界面处存在异常极值的井点地层Q值进行适当删减平滑，根据Q值与深度关系，结合实际地层情况，参考层速度的分层情况对Q值曲线进行分层处理，使其与层速度变化趋势相符合，最终计算得到每个地质层位内的平均Q值，从而获得较为准确的地层Q值。图8示出了某口VSP井计算得到的基于地质层位约束的Q值随深度变化曲线。

第五步，在每个VSP井点地质层位平均Q值基础上，将工区所有VSP井点地层Q值运用克里金方法进行插值，并做适当平滑，建立全工区近地表Q值场。

第六步，利用全工区近地表Q值场，对地震数据应用基于波动方程的反Q滤波算法，最终补偿地震数据由于近地表吸收引起的高频信息损失。图9示出了本申请实施例提供的某口井未进行反Q滤波的时间偏移剖面示意图，图10则示出经过反Q滤波的时间偏移剖面示意图，对比图9和图10可以看到，反Q滤波后整个剖面的频率明显提高，分辨率也提高，地层的断点位置更加清楚。

图11示出了地面地震叠前数据反Q滤波补偿前后的叠前时间偏移剖面的振幅谱分析。从图中可以看出，反Q滤波处理后，中、高频成分的能量得到补偿,地震信号的频带宽度得到拓宽，主频向高频端移动，分辨率明显提高。

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

基于同一技术构思，本申请实施例中还提供了一种地层品质因子反演装置，如下面的实施例所述。由于地层品质因子反演装置解决问题的原理与地层品质因子反演方法相似，因此地层品质因子反演装置的实施可以参见地层品质因子反演方法的实施，重复之处不再赘述。以下所使用的，术语“单元”或者“模块”可以实现预定功能的软件和/或硬件的组合。尽管以下实施例所描述的装置较佳地以软件来实现，但是硬件，或者软件和硬件的组合的实现也是可能并被构想的。图12示出了本申请实施例的地层品质因子反演装置的一种结构框图，如图12所示，该装置包括：获取模块1201、校正模块1202、层速度确定模块1203、质心频率确定模块1204、平均吸收系数确定模块1205和品质因子确定模块1206。下面对该结构进行说明。

获取模块1201，用于获取检波点的垂直地震剖面VSP资料。

校正模块1202，用于对每个检波点的VSP资料进行非零偏校正，确定每个检波点的零偏VSP资料初至时间。

层速度确定模块1203，用于根据每个检波点的零偏VSP资料初至时间和深度，确定每相邻两个检波点之间的层速度。

质心频率确定模块1204，用于根据每个检波点的VSP资料，通过傅里叶变换得到波场振幅谱，根据所述波场振幅谱确定质心频率。

平均吸收系数确定模块1205，用于根据每个检波点的深度和质心频率，确定每相邻两个检波点之间的平均吸收系数。

品质因子确定模块1206，用于根据每相邻两个检波点之间的平均吸收系数和层速度，确定每相邻两个检波点之间的品质因子Q值。

平均吸收系数确定模块1205具体用于按照上述公式(1)确定所述平均吸收系数。

品质因子确定模块1206具体用于按照上述公式(2)确定所述品质因子。

所述平均吸收系数确定模块1205用于对所述零偏VSP资料进行几何扩散补偿以及波场分离，得到下行波场；根据所述下行波场中的每个检波点的VSP深度和质心频率，确定每相邻两个检波点之间的平均吸收系数；进一步，所述品质因子确定模块1206具体用于：根据所述下行波场中的每相邻两个检波点之间的平均吸收系数和层速度，确定每相邻两个检波点之间的Q值。

在另外一个实施例中，还提供了一种软件，该软件用于执行上述实施例及优选实施方式中描述的技术方案。

在另外一个实施例中，还提供了一种存储介质，该存储介质中存储有上述软件，该存储介质包括但不限于：光盘、软盘、硬盘、可擦写存储器等。

本说明书中上述方法的各个实施例均采用递进的方式描述，各个实施例之间相同相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。相关之处参见方法实施例的部分说明即可。除非另有定义，本文所使用的所有的技术和科学术语与属于本申请的技术领域的技术人员通常理解的含义相同。本文中在本申请的说明书中所使用的术语只是为了描述具体的实施方式的目的，不是旨在于限制本申请。

需要说明的是，尽管在附图中以特定顺序描述了本申请方法的操作，但这并非要求或者暗示必须按照该特定顺序来执行这些操作，或是必须执行全部所示的操作才能实现期望的结果。附加地或备选地，可以省略某些步骤，将多个步骤合并为一个步骤执行，和/或将一个步骤分解为多个步骤执行。

虽然本申请提供了如实施例或流程图的方法操作步骤，但基于常规或者无创造性的手段可以包括更多或者更少的操作步骤。实施例中列举的步骤顺序仅仅为众多步骤执行顺序中的一种方式，不代表唯一的执行顺序。在实际中的装置或客户端产品执行时，可以按照实施例或者附图所示的方法顺序执行或者并行执行(例如并行处理器或者多线程处理的环境，甚至为分布式数据处理环境)。术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、产品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、产品或者设备所固有的要素。在没有更多限制的情况下，并不排除在包括所述要素的过程、方法、产品或者设备中还存在另外的相同或等同要素。

本领域内的技术人员应明白，本申请的实施例可提供为方法、系统、或计算机程序产品。因此，本申请可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。

本申请是参照根据本申请实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。

这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。

这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。

以上所述的具体实施例，对本申请的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上所述仅为本申请的具体实施例而已，并不用于限定本申请的保护范围，凡在本申请的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本申请的保护范围之内。

Claims (10)

1.一种地层品质因子反演方法，其特征在于，包括：

获取检波点的垂直地震剖面VSP资料；

对每个检波点的VSP资料进行非零偏校正，确定每个检波点的零偏VSP资料初至时间；

根据每个检波点的零偏VSP资料初至时间和深度，确定每相邻两个检波点之间的层速度；

根据每个检波点的VSP资料进行傅里叶变换，确定波场振幅谱，根据所述波场振幅谱确定每个检波点的质心频率；

根据每个检波点的深度和质心频率，确定每相邻两个检波点之间的平均吸收系数；

根据每相邻两个检波点之间的平均吸收系数和层速度，确定每相邻两个检波点之间的品质因子Q值。

2.如权利要求1所述的方法，其特征在于，所述根据每个检波点的VSP深度和质心频率，确定每相邻两个检波点之间的平均吸收系数，包括：

按照以下公式确定所述平均吸收系数：

其中，α_i表示第i个检波点与第i+1个检波点之间的平均吸收系数，σ_i ²表示第i个检波点的方差，f_i表示第i个检波点的质心频率，f_i+1表示第i+1个检波点的质心频率，z_i表示第i个检波点的深度，z_i+1表示第i+1个检波点的深度，Δz_i表示第i个检波点到第i+1个检波点之间的深度差，Δf_i表示第i个检波点到第i+1个检波点之间的质心频率差。

3.如权利要求1所述的方法，其特征在于，所述根据每相邻两个检波点之间的平均吸收系数和层速度，确定每相邻两个检波点之间的品质因子Q值，包括：

按照以下公式确定所述Q值：

其中，Q_i表示第i个检波点到第i+1个检波点之间的品质因子，α_i表示第i个检波点与第i+1个检波点之间的平均吸收系数，v_i表示第i个检波点到第i+1个检波点之间的层速度。

4.如权利要求1至3任一项所述的方法，其特征在于，所述根据每个检波点的VSP深度和质心频率，确定每两个检波点之间的平均吸收系数，包括：

对所述零偏VSP资料进行几何扩散补偿以及波场分离，得到下行波场；

根据所述下行波场中的每个检波点的VSP深度和质心频率，确定每相邻两个检波点之间的平均吸收系数；

所述根据每相邻两个检波点之间的平均吸收系数和层速度，确定每相邻两个检波点之间的品质因子Q值，包括：

根据所述下行波场中的每相邻两个检波点之间的平均吸收系数和层速度，确定每相邻两个检波点之间的Q值。

5.一种地层品质因子反演装置，其特征在于，包括：

获取模块，用于获取检波点的垂直地震剖面VSP资料；

校正模块，用于对每个检波点的VSP资料进行非零偏校正，确定每个检波点的零偏VSP资料初至时间；

层速度确定模块，用于根据每个检波点的零偏VSP资料初至时间和深度，确定每相邻两个检波点之间的层速度；

质心频率确定模块，用于根据每个检波点的VSP资料进行傅里叶变换，确定波场振幅谱，根据所述波场振幅谱确定质心频率；

平均吸收系数确定模块，用于根据每个检波点的深度和质心频率，确定每相邻两个检波点之间的平均吸收系数；

品质因子确定模块，用于根据每相邻两个检波点之间的平均吸收系数和层速度，确定每相邻两个检波点之间的品质因子Q值。

6.如权利要求5所述的装置，其特征在于，所述平均吸收系数确定模块具体用于按照以下公式确定所述平均吸收系数：

其中，α_i表示第i个检波点与第i+1个检波点之间的平均吸收系数，σ_i ²表示第i个检波点的方差，f_i表示第i个检波点的质心频率，f_i+1表示第i+1个检波点的质心频率，z_i表示第i个检波点的深度，z_i+1表示第i+1个检波点的深度，Δz_i表示第i个检波点到第i+1个检波点之间的深度差，Δf_i表示第i个检波点到第i+1个检波点之间的质心频率差。

7.如权利要求5所述的装置，其特征在于，所述品质因子确定模块具体用于按照以下公式确定所述品质因子：

其中，Q_i表示第i个检波点到第i+1个检波点之间的品质因子，α_i表示第i个检波点与第i+1个检波点之间的平均吸收系数，v_i表示第i个检波点到第i+1个检波点之间的层速度。

8.如权利要求5至7任一项所述的装置，其特征在于，所述平均吸收系数确定模块具体用于：

对所述零偏VSP资料进行几何扩散补偿以及波场分离，得到下行波场；

根据所述下行波场中的每个检波点的VSP深度和质心频率，确定每相邻两个检波点之间的平均吸收系数；

所述品质因子确定模块具体用于：

根据所述下行波场中的每相邻两个检波点之间的平均吸收系数和层速度，确定每相邻两个检波点之间的Q值。

9.一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述计算机程序时实现权利要求1至4任一所述方法。

10.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质存储有执行权利要求1至4任一所述方法的计算机程序。