一种二阶多重同步挤压多项式chirplet变换薄储层识别方法技术

本发明专利技术公开了一种二阶多重同步挤压多项式chirplet变换薄储层识别方法，包括以下步骤：S1、输入待分析的地震数据；S2、对多项式chirplet变换依次求取时间和频率方向的偏导，得到二阶瞬时频率估计算子；并将该估计算子进行多次迭代得到一种新的二阶多重瞬时频率估计函数；S3、利用二阶多重同步挤压多项式chirplet变换(MSSPCT)对地震数据进行时频变换，得到最终的时频结果；S4、对上述时频结果取模，得其高精度时频谱；S5、采用傅里叶变换获取过井道的地震道振幅谱，将振幅最大值设置为低频段，将振幅相对较低处设置为高频段；S6、从MSSPCT得到的时频谱中分别抽取高频段和低频段对应的频谱，通过两个共频率剖面中地震信号的衰减情况来表征该地震资料中的薄储层展布特征。特征。特征。

【技术实现步骤摘要】
一种二阶多重同步挤压多项式chirplet变换薄储层识别方法


[0001]本专利技术属于地震信号处理与地震储层识别领域，具体涉及一种二 阶多重同步挤压多项式chirplet变换薄储层识别方法。

技术介绍

[0002]时频分析是挖掘地震信号瞬时谱信息的重要方法，可用于储层的 精细识别。根据发展历程，时频分析方法可分为传统时频分析方法和 高精度时频分析方法两大类。传统时频分析方法包括短时傅里叶变换 (STFT)、小波变换(WT)、S变换(ST)及广义S变换(GST)等， 已被广泛应用于储层预测、烃类检测等领域；相较于传统方法，高精 度时频分析方法显著提高了时频分辨率，主要包括经验模态分解 (EMD)、同步挤压方法(SST)、同步提取方法(SET)等，在对实 际地震资料进行分频处理时，对信号刻画的精度更高，可以更为有效 地预测和识别薄储层。
[0003]除了前面提到的时频分析方法以外，chirplet变换(CT)是另一种 参数化的时频分析方法。CT对于线性调频信号能够获得很好的分析 效果。然而，当待分析信号瞬时频率为非线性函数时，例如强调频信本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种二阶多重同步挤压多项式chirplet变换薄储层识别方法，其特征在于，包括以下步骤：S1、输入待分析信号x(t)；S2、在该信号模型下，对多项式chirplet变换先进行时间方向的偏导，再对其沿频率方向求偏导，从而得到二阶瞬时频率的初步估计算子；根据上述瞬时频率估计算子，并以此类推，得到多次迭代后一种新的二阶多重瞬时频率估计函数；S3、根据同步挤压原理，在时频域上构建一种以信号瞬时频率为中心的同步挤压算子，用于挤压原时频谱瞬时频率曲线处的时频系数，得到二阶多重同步挤压多项式chirplet变换值；S4、对二阶多重同步挤压多项式chirplet变换的时频结果取模，得到每一点处的时频能量，从而得其高精度时频谱；S5、采用傅立叶变换获取过井道的地震道振幅谱，将振幅最大值处设置为低频段，将振幅相对较小值处设置为高频段；S6、从二阶多重同步挤压多项式chirplet变换得到的时频谱中分别抽取低频段和高频段对应的频谱，得到两个共频率剖面；通过对比两个共频率剖面中地震信号的衰减情况识别该地震资料中的薄储层展布特征。2.根据权利要求1所述的一种二阶多重同步挤压多项式chirplet变换薄储层识别方法，其特征在于，所述步骤S2中对多项式chirplet变换先进行时间方向求偏导，得到:其中，表示求关于时间变量t的偏导数，为在窗口g(t)和信号s(u)下的PCT，t为时间变量，ω为频率变量，A(t)和分别表示信号的瞬时振幅和瞬时相位函数，A
(1)
(t)为一阶导数，A
(2)
(t)为二阶导数，为在窗口tg(t)和信号s(u)下的PCT；为了得到...

【专利技术属性】
技术研发人员：陈辉，丁柏蓉，胡英，陈旭平，李军，
申请(专利权)人：成都理工大学，
类型：发明
国别省市：