[发明专利]一种地震属性体提取的方法及系统

说明书

本发明提供了一种地震属性体提取的方法与系统，其中，该方法包括：根据地震信号的子波特征，选取预设的小波函数作为基阶滤波函数；变化小波函数的带宽形成具有相同中心频率和不同滤波带宽的滤波函数集合，记为超小波基；将地震信号的每一道分别与超小波基进行卷积运算得到相应基阶小波主频的超小波变换信号；改变基阶滤波函数的主频，获得不同主频对应的超小波变换信号；根据不同主频对应的超小波变换信号，获得地震信号的属性体能量分布图。本发明通过选取可调的小波函数构建超小波基，能够获取高分辨率的能量图，进而提高提取地震属性体的准确率，为后续地震断裂位置、断裂的轮廓等方面的识别提供了可靠的数据支撑。

技术领域

本发明涉及地震探测技术领域，特别是涉及一种地震属性体提取的方法及系统。

背景技术

地震数据通常是以时间和空间同时采样的信号，其中蕴含探测区地下各种特征信息。基于时频分析技术的地震信号分析方法是获得地震属性体特征的有效方法之一，如提取勘探地震剖面中的断裂甚至裂缝信息、二维时间截面中的河道信息等。

地震信号一般属于非线性非平稳信号，传统的信号时频分析方法有短时傅立叶变换、小波变换等，这些时频方法具有一定的时间(空间)分辨和频率分辨，在地震数据属性体分析中得到了广泛应用。而受限于海森堡测不准原理，上述方法不能同时在时间和频率上均能够取得最优分辨，一定程度上影响了对属性体特征的识别。

发明内容

为了克服现有技术的不足，本发明的目的是提供一种地震属性体提取的方法及系统，以提高属性体在不同频率下的特征分辨能力。

为实现上述目的，本发明提供了如下方案：

一种地震属性体提取的方法，包括：

步骤1：获取待处理的地震信号；

步骤2：根据所述地震信号的子波特征，选取预设的小波函数作为基阶滤波函数；

步骤3：变化所述小波函数的带宽形成具有相同中心频率和不同滤波带宽的滤波函数集合，记为超小波基；

步骤4：将所述地震信号的每一道分别与所述超小波基进行卷积运算得到相应基阶小波主频的超小波变换信号；

步骤5：改变基阶滤波函数的主频，返回步骤2，获得不同主频对应的超小波变换信号；

步骤6：根据所述不同主频对应的超小波变换信号，获得所述地震信号的属性体能量分布图。

优选的，所述基阶滤波函数为：

其中，ψ_f,c(t)表示基阶滤波函数，f表示中心频率，c为周期数，B_c为有效带宽，t表示时间。

优选的，所述步骤3：变化所述小波函数的带宽形成具有相同中心频率和不同滤波带宽的滤波函数集合，记为超小波基，包括：

按照等比序列或等差序列设定所述小波函数的带宽形成具有相同中心频率和不同滤波带宽的滤波函数集合，记为超小波基。

优选的，所述步骤4：将所述地震信号的每一道分别与所述超小波基进行卷积运算得到相应基阶小波主频的超小波变换信号，包括：

采用公式：

将所述地震信号的每一道x_i分别与所述超小波基进行卷积运算得到相应基阶小波主频的超小波变换信号；其中，*表示复卷积算子，R[SL_f,o]表示超小波变换信号，o表示阶次，c_i表示基阶小波的带宽参数。

本发明还提供了一种地震属性体提取的系统，包括：