一种基于匹配追踪的Wigner高阶谱地震信号谱分解方法

摘要

一种基于匹配追踪的Wigner高阶谱地震信号谱分解方法，首先读入地震剖面，选定原子类型；然后选取一道地震数据；接着对信号进行复地震道分析以及对尺度因子进行全局搜索，确定原子的初始参数集；对参数集进行局部搜索，找到与信号最匹配的原子；计算最匹配原子的Wigner高阶谱的对角切片谱；计算该信号在最匹配原子方向上投影后的残差，并将其视为新的分解信号；对所有分解得到的原子的Wigner高阶谱的对角切片谱求和，作为该道地震数据的Wigner高阶谱时频谱，截取单频切片；对所有地震数据采用相同方法以得到谱分解的结果。本发明利用匹配追踪方法去除Wigner高阶谱的交叉项，可以获得时频聚集性更高的地震谱分解结果，为后续地震储层预测和流体识别提供更精确的信息。