[发明专利]基于二阶同步挤压小波变换的地震资料时频分析方法

摘要

本发明公开一种基于二阶挤压小波变换的地震资料时频分析方法，首次在连续小波变换的基础上，引入二阶同步挤压算子，该变换通过对小波变换的系数在时间‑尺度域进行重排，从而压制了能量扩散，该方法可以得到一种具有更高时频分辨率的分析结果。相比于之前的同步挤压变换，改方法更加适合对地震信号这种带有较强频率调制的信号进行分析，进一步将其应用于实际地震数据的时频分析中，可以更加准确的界定地下碳氢储层的边界，更精细的刻画河道与断层等地质信息，从而有利于地震资料的进一步解释和井位确定。

主权项

1.基于二阶同步挤压小波变换的地震资料时频分析方法，其特征在于：二阶同步挤压小波变换包括以下步骤：1)对信号进行连续小波变换，得到小波变换谱W(a,b)，W(a,b)的时域和频域表达式分别为其中，ψ(t)为母小波函数，a和b分别表示尺度和平移因子，X(ω)和Ψ(ω)分别为x(t)和ψ(t)的傅里叶变换，*表示取共轭；基本母小波函数ψ(t)平方可积，无直流分量，并且满足以下容许性条件：2)计算每一个小波系数对应的参考频率其中3)将时间‑尺度域的小波系数重排至时间‑频率域，二阶同步挤压小波变换的形式为δ(·)表示单位冲击函数；通过公式(4),在时间‑尺度域所有和瞬时频率对应的小波系数沿着尺度方向进行纵向求和，然后将求和的结果放置在时间‑频率域w所对应的位置；基于以上二阶同步挤压小波变换的地震资料时频分析方法具体包括以下步骤：步骤1：在三维数据体中的某一地震剖面选取典型地震道，并对该地震数据进行二阶同步挤压小波变换，并找出异常区域所对应的频率成分；步骤2：对整个地震剖面逐道进行二阶同步挤压小波变换，提取出与异常区域相对应的频率切片；步骤3:对整个三维数据进行二阶同步挤压小波变换，得到对应的频率数据体，根据解释人员提供的层位信息，制作沿层切片，供地质人员解释。