[发明专利]一种时频域等效Q场的获取方法及装置

说明书

本发明提供了一种时频域等效Q场的获取方法及装置，包括：对获取的叠前地震数据进行分频，获得分频地震数据；对分频地震数据沿时间/深度方向统计所处频段的振幅能量；根据振幅能量获取所处频段沿时间/深度方向的等效Q值；根据获取的等效Q值重构单道Q值曲线。本申请结合了基于频率信息和基于振幅信息进行Q场建模方法的优点，实现了能够人为控制有效频带范围和加权系数并使求取的Q场具有高精度和适应性较强的效果。

技术领域

本申请属于石油地震勘探中的地震数据处理技术领域，具体地讲，涉及一种时频域等效Q场的获取方法及装置。

背景技术

由于大地吸收衰减效应的影响，地震信号在传播过程中会产生高频成分衰减、子波相位畸变等现象，这种与地下介质有关的吸收效应通常用介质品质因子Q标识。为了恢复振幅衰减和相位畸变，通常采用反Q滤波的方法模拟地震波传播的反过程，并用考虑吸收补偿因素的QPSDM偏移算法(含Q吸收补偿的叠前深度偏移算法——Q compensationprestack depth migration)进行成像，两种方法均需要首先求取空间Q场。

现今，Q场的求取方法主要可分为两大类：基于数据频率信息的方法和基于振幅信息的方法。基于数据频率信息的方法不受振幅保真度、信噪比低等问题的干扰，包括：Q同步迭代建模方法和反射波质心频率偏移建模方法，但是这种方法一般假设Q值变化与不同频率分量无关，利用数据主频变化估算Q值；基于振幅信息的方法由于不受地震数据类型的限制可以应用于三维数据的Q建模，包括：反射波衰减旅行时的Q层析建模方法，但是这种方法由于受到射线方法的局限性会对Q层析造成误差累计，进而影响Q层析的精度。

发明内容

本申请提供了一种时频域等效Q场的获取方法及装置，以至少解决现有技术中在进行Q场求取的过程中受到射线类方法的局限性导致对Q层析造成误差累计进而影响Q层析的精度的问题。

根据本申请的一个方面，提供了一种时频域等效Q场的获取方法，包括：

对获取的叠前地震数据进行分频，获得分频地震数据；

对分频地震数据沿时间/深度方向统计所处频段的振幅能量；

根据振幅能量获取所处频段沿时间/深度方向的等效Q值；

根据获取的等效Q值重构单道Q值曲线。

在一实施例中，对获取的叠前地震数据进行分频，获得分频地震数据，包括：

构建余弦滤波函数，并利用余弦滤波函数对叠前地震数据进行分频获得分频地震数据及分频重构系数。

在一实施例中，根据振幅能量获取所处频段沿时间/深度方向的等效Q值，包括：

采用高阶自然指数曲线对振幅能量进行拟合，获得拟合后的振幅能量；

根据拟合后的振幅能量获得所处频段沿时间/深度方向的等效Q值。

在一实施例中，根据获取的等效Q值重构单道Q值曲线，包括：

利用分频重构系数和等效Q值获得Q值曲线；

根据Q值曲线获得全区等效Q场；

对全区等效Q场进行平滑后输出。

根据本申请的另一个方面，还提供了一种时频域等效Q场的确定装置，包括：

分频处理单元，用于对获取的叠前地震数据进行分频，获得分频地震数据；

振幅能量统计单元，用于对分频地震数据沿时间/深度方向统计所处频段的振幅能量；

等效Q值获取单元，用于根据振幅能量获取所处频段沿时间/深度方向的等效Q值；

重构单元，用于根据获取的等效Q值重构单道Q值曲线。