一种地震品质因子估算方法、装置、设备及存储介质制造方法及图纸

本发明专利技术实施例公开了一种地震品质因子估算方法、装置、设备及存储介质。该方法包括：获取地震数据，并根据地震数据确定地震子波振幅对数谱以及震源子波主频；根据震源子波主频确定频带宽度，并根据震源子波主频和频带宽度确定地震子波振幅对数谱中的高频区域和低频区域；确定高频区域与低频区域的对数谱面积差，并根据对数谱面积差确定地震品质因子。本发明专利技术实施例的技术方案，解决了由于采集的地震数据中存在较多噪声，以及地震数据所对应的子波振幅谱不满足高斯分布假设而导致地震品质因子估算不准确的问题，提高了地震品质因子的估算精度以及估算结果的稳定性。

【技术实现步骤摘要】
一种地震品质因子估算方法、装置、设备及存储介质
本专利技术实施例涉及地震勘探
，尤其涉及一种地震品质因子估算方法、装置、设备及存储介质。
技术介绍
地震波在实际地层介质中传播过程中，几何扩散、投射、反射损失以及介质的非完全弹性等均会引起地震波的衰减。地震品质因子是描述地震波吸收衰减的重要参数，地震品质因子估算的准确性和稳定性对提高地震资料分辨率、地质体准确成像和解释、储层预测和油藏描述等具有重要意义。在当前工业界地震品质因子的估计方法主要有两种：1)谱比法(LogarithmSpectralRatio，LSR)通过在所选的频带范围内对接收子波(衰减后的地震波)和震源子波(衰减前的地震波)振幅谱比的对数进行线性拟合，得到拟合直线的斜率，根据斜率与品质因子的关系式对地震品质因子进行估算；2)质心频移法(CentroidFrequencyShift，CFS)根据地震波衰减前后质心频率的差异对地震品质因子进行估计。然而谱比法由于采用线性拟合方法，对噪声相对敏感，故其仅在理论上估计精度高，而对于包含高频和低频噪音的实际地震数据难以达到良好的估算精度；而质心频移法是以震源子波振幅谱满足高斯分布为基础的地震品质因子估算方法，其虽受到噪音干扰的程度小于谱比法，但由于实际地震数据中的子波振幅谱难以完全满足高斯分布假设，仍会引起地震品质因子估计的系统误差，降低了地震品质因子估计的准确性和稳定性。
技术实现思路
本专利技术提供一种地震品质因子估算方法、装置、设备及存储介质，以使地震品质因子的估算更加稳定和准确，降低了噪声对估值的干扰，减少了由于子波振幅谱不满足高斯分布假设而导致的系统误差。第一方面，本专利技术实施例提供了一种地震品质因子估算方法，包括：获取地震数据，并根据地震数据确定地震子波振幅对数谱以及震源子波主频；根据震源子波主频确定频带宽度，并根据震源子波主频和频带宽度确定地震子波振幅对数谱中的高频区域和低频区域；确定高频区域与低频区域的对数谱面积差，并根据对数谱面积差确定地震品质因子。第二方面，本专利技术实施例还提供了一种地震品质因子估算装置，包括：数据获取模块，用于获取地震数据，并根据地震数据确定地震子波振幅对数谱以及震源子波主频；区域确定模块，用于根据震源子波主频确定频带宽度，并根据震源子波主频和频带宽度确定地震子波振幅对数谱中的高频区域和低频区域；品质因子确定模块，用于确定高频区域与低频区域的对数谱面积差，并根据对数谱面积差确定地震品质因子。第三方面，本专利技术实施例还提供了一种计算机设备，所述计算机设备包括：一个或多个处理器；存储装置，用于存储一个或多个程序；当所述一个或多个程序被所述一个或多个处理器执行，使得所述一个或多个处理器实现如本专利技术任意实施例提供的地震品质因子估算方法。第四方面，本专利技术实施例还提供了一种包含计算机可执行指令的存储介质，所述计算机可执行指令在由计算机处理器执行时用于执行如本专利技术任意实施例提供的地震品质因子估算方法。本专利技术实施例通过获取地震数据，并根据地震数据确定地震子波振幅对数谱以及震源子波主频；根据震源子波主频确定频带宽度，并根据震源子波主频和频带宽度确定地震子波振幅对数谱中的高频区域和低频区域；确定高频区域与低频区域的对数谱面积差，并根据对数谱面积差确定地震品质因子。通过采用上述技术方案，根据震源子波主频在地震子波振幅对数谱中确定出频带宽度相等的高频区域和低频区域，使得地震数据中的有效信息可得到较好的保存并对地震数据中的高频和低频噪声进行一定程度的过滤；确定高频区域内的地震子波振幅对数谱面积以及低频区域内的地震子波振幅对数谱面积，进而根据上述两对数谱面积的差对地震品质因子进行估算，由于随机噪声的平均值为0，故在进行谱面积计算时可对其中的随机噪声也进行相加，对噪声起到了较好的抑制作用，降低了噪声对地震品质因子估算精度的影响；同时根据谱面积差对地震品质因子进行估算无需假设地震子波振幅谱满足高斯分布假设，减少了系统误差，解决了由于采集的地震数据中存在较多噪声，以及地震数据所对应的子波振幅谱不满足高斯分布假设而导致地震品质因子估算不准确的问题，提高了地震品质因子的估算精度以及估算结果的稳定性。附图说明图1是本专利技术实施例一中的一种地震品质因子估算方法的流程图；图2是本专利技术实施例二中的一种地震品质因子估算方法的流程图；图3是本专利技术实施例二中的一种地震子波振幅对数谱的示例图；图4是本专利技术实施例三中的一种合成衰减地震记录的示例图；图5是本专利技术实施例三中的三种方法地震品质因子估算结果的对比图；图6是本专利技术实施例三中的三种方法地震品质因子估计结果相对误差的对比图；图7是本专利技术实施例三中的一种含随机噪声的合成地震记录的示例图；图8是本专利技术实施例三中的三种方法的地震品质因子估计结果概率分布图；图9是本专利技术实施例三中的三种方法的地震品质因子估计结果均值示例图；图10是本专利技术实施例三中的三种方法的地震品质因子估计结果标准差示例图；图11是本专利技术实施例三中的中国某油田实际零偏VSP资料图；图12是本专利技术实施例三中的对数谱高频区与低频区面积差法对地震品质因子的估计结果示例图；图13是本专利技术实施例三中的实际地层速度的示例图；图14是本专利技术实施例四中的一种地震品质因子估算装置的结构示意图；图15是本专利技术实施例五中的一种计算机设备的结构示意图。具体实施方式下面结合附图和实施例对本专利技术作进一步的详细说明。可以理解的是，此处所描述的具体实施例仅仅用于解释本专利技术，而非对本专利技术的限定。另外还需要说明的是，为了便于描述，附图中仅示出了与本专利技术相关的部分而非全部结构。此外，在不冲突的情况下，本专利技术中的实施例及实施例中的特征可以互相组合。另外还需要说明的是，为了便于描述，附图中仅示出了与本专利技术相关的部分而非全部内容。在更加详细地讨论示例性实施例之前应当提到的是，一些示例性实施例被描述成作为流程图描绘的处理或方法。虽然流程图将各项操作(或步骤)描述成顺序的处理，但是其中的许多操作可以被并行地、并发地或者同时实施。此外，各项操作的顺序可以被重新安排。当其操作完成时所述处理可以被终止，但是还可以具有未包括在附图中的附加步骤。所述处理可以对应于方法、函数、规程、子例程、子程序等等。实施例一图1为本专利技术实施例一提供的一种地震品质因子估算方法的流程图，本实施例可适用于通过地震数据的子波振幅对数谱面积差对地震品质因子进行估算的情况，该方法可以由地震品质因子估算装置来执行，该地震品质因子估算装置可以由软件和/或硬件来实现，该地震品质因子估算装置可以配置在计算设备上，具体包括如下步骤：S101、获取地震数据，并根据地震数据确定地震子波振幅对数谱以及震源子波主频。在本实施例中，地震数据可理解为由检波器采集到的共反射点道集地震数据，可包括叠前地震数据和叠后地震数据。振幅谱本文档来自技高网...

【技术保护点】
1.一种地震品质因子估算方法，其特征在于，包括：/n获取地震数据，并根据所述地震数据确定地震子波振幅对数谱以及震源子波主频；/n根据所述震源子波主频确定频带宽度，并根据所述震源子波主频和所述频带宽度确定所述地震子波振幅对数谱中的高频区域和低频区域；/n确定所述高频区域与所述低频区域的对数谱面积差，并根据所述对数谱面积差确定地震品质因子。/n

【技术特征摘要】
20210108 CN 20211002600471.一种地震品质因子估算方法，其特征在于，包括：
获取地震数据，并根据所述地震数据确定地震子波振幅对数谱以及震源子波主频；
根据所述震源子波主频确定频带宽度，并根据所述震源子波主频和所述频带宽度确定所述地震子波振幅对数谱中的高频区域和低频区域；
确定所述高频区域与所述低频区域的对数谱面积差，并根据所述对数谱面积差确定地震品质因子。


2.根据权利要求1所述的方法，其特征在于，所述根据所述叠前地震数据确定地震子波振幅对数谱，包括：
将所述地震数据进行子波拾取，确定地震子波；
将所述地震子波进行傅里叶变换，确定地震子波振幅谱；
将所述地震子波振幅谱取对数，确定地震子波振幅对数谱。


3.根据权利要求1所述的方法，其特征在于，所述根据所述震源子波主频确定频带宽度，包括：
根据所述地震数据的信噪比、所述震源子波主频以及预设映射关系确定频带比例；
将所述震源子波主频与所述频带比例的乘积确定为频带宽度。


4.根据权利要求1所述的方法，其特征在于，所述根据所述震源子波主频和所述频带宽度确定所述地震子波振幅对数谱中的高频区域和低频区域，包括：
将所述震源子波主频与所述频带宽度的和确定为高频边界，并将所述震源子波主频与所述频带宽度的差确定为低频边界；
将所述震源子波主频与所述高频边界间的区域确定为所述地震子波振幅对数谱中的高频区域；
将所述震源子波主频与所述低频边界间的区域确定为所述地震子波振幅对数谱中的低频区域。


5.根据权利要求1所述的方法，其特征在于，所述地震子波振幅对数谱包括震源子波振幅对数谱和接收子波振幅对数谱，所述确定所述高频区域与所述低频区域的对数谱面积差，包括：
将所述震源子波振幅对数谱和所述接收子波振幅对数谱分别在所述高频区域积分，确定震源子波高频谱面积和接收子波高频谱面积；
将所述震源子波振幅对数谱和所述接收子波振幅对数谱分别在所述低频区域积分，确定震源子波低频谱面积和接收子波低频谱面积；

【专利技术属性】
技术研发人员：杨登锋，刘军，彭光荣，吴静，白海军，刘道理，李振升，杜晓东，
申请(专利权)人：中海石油中国有限公司深圳分公司，
类型：发明
国别省市：广东;44