一种基于稀疏窗S变换的分频-重构波阻抗反演方法*

冉喜阳 周怀来 张益明 李雷豪 杨吉鑫

(1.成都理工大学地球物理学院 四川成都 610059； 2.成都理工大学油气藏地质及开发工程国家重点实验室 四川成都 610059；3.四川省地震局 四川成都 610041； 4.中海油研究总院有限责任公司 北京 100028)

在地震勘探中，不同频率的地震反射波特征反映了不同的地层岩性、物性及含流体性质[1],故相比常规地震资料，分频技术的应用能够从地震数据中获取更多的有效信息，有利于对目标体进行更全面的分析解释，进而提高储层预测精度,降低解释多解性。地震分频技术的基础是时频分析方法，目前常见的时频分析方法包括短时傅里叶变换、Gabor变换、Wigner Ville时频分布、小波变换、S变换等[2]。

为了从地震数据中获得更多的有效信息，研究人员对分频技术做了很多的研究工作[3-12]。前人利用各种时频分析方法对测井数据和地震数据进行分频处理，发现不同频率尺度地层特征在此基础上进行波阻抗反演获得了较高分辨率的反演结果，但在应用过程中也存在一些问题。王振卿 等[5]对塔中西部台内数据用广义S变换进行谱分解，使用带通滤波处理得到15～25 Hz优势频段，对其进行波阻抗反演，虽然提高了识别精度，应用效果良好，但是在滤波过程中也滤掉了大量的有效波，忽视了部分细节。

另外，在常规的时频分析方法中，由于短时Fourier变换的窗函数固定，不能适应非平稳信号中不同频率成分特性的分析；S变换通过改进时间窗函数的尺度因子，具备多尺度分析能力，但S变换中基本小波是固定的，使其在应用中受到了限制；……