非平稳信号时频分析与地震频谱成像研究

论文摘要

为了提高油气储层预测和流体识别精度,非平稳地震信号的时频域去噪方法和高精度时频分析研究日趋活跃。本文以非平稳地震信号的分析和处理为研究课题,重点研究了非平稳地震信号的去噪技术、基于分数阶傅立叶变换的高精度时频分析方法、基于稀疏正则约束的稀疏时频分析方法等。主要研究内容分为以下几个方面:（1）研究了非平稳信号时频分析方法和数据预处理方法。时频分析方法和数据预处理是分析非平稳信号的重要技术。由于非平稳信号源的复杂多变性,时频分析和数据预处理问题是非常有挑战性和活跃的研究话题。第一部分调研了各类经典的时频分析和数据预处理方法。（2）提出一种基于交叠组稀疏的广义全变分去噪方法。该方法重点讨论了基于交叠组稀疏的改进广义全变分方法。广义全变分模型是全变分模型的推广并被证明是能够有效去除全变分模型阶梯效应的方法,然而,该模型独立处理每个像素,忽略了图像的结构相似先验。因此,广义全变分模型对于大幅度噪声并不鲁棒。本部分的研究动机是通过利用图像的结构相似性进一步提高广义全变分模型的去噪性能。通过引入交叠组稀疏到广义全变分模型,挖掘图像一阶和二阶领域差分梯度信息的结构相似性,从而达到提高广义全变分模型对重噪声污染的鲁棒性。为求解提出模型,使用了加速重启的交替乘子迭代法（Alternating Direction Method of Multipliers,ADMM）,将复杂的多约束问题转为若干子问题。为避免空域大型矩阵的运算,将差分算子视为卷积形式,然后利用快速傅立叶变换和卷积定理求解提出模型。针对不同噪声的地震信号进行去噪实验发现如下结论:1)提出模型特别对平滑区域的重噪声有较好的去噪效果;2)带重启的加速交替乘子迭代法算法能求解提出模型,比ADMM效率高;3)组合块的大小需要合理选择,从而获得最佳的去噪性能。（3）提出一种基于分数阶傅立叶变换和贪婪策略的高精度时频分析方法。Wigner-Ville分布是一种在地震信号处理领域具有高分辨率的重要时频分析技术,然而它被大量交叉项干扰。为了在去除Wigner-Ville分布交叉项的同时保持其良好的聚集性,本文提出一种基于分数阶傅立叶变换和贪婪策略的多向自适应模糊域窗函数。首先从Cohen类与Gabor变换的关系出发,并充分利用贪婪策略和分数阶傅立叶变换的旋转特性获得模糊域自适应多方向窗,把最优分数解Gabor变换的单方向、一维的最优窗函数推广为二维的、多向模糊域窗函数。通过这种方法,在处理多成分信号时,多向窗函数能准确匹配Wigner-Ville分布的交叉项。利用贪婪策略,提出方法能将最优方向以及其他子方向充分考虑,从而避免了最优分数阶Gabor变换的局部聚集现象。（4）提出一种基于一阶原始对偶优化的稀疏时频分析方法。时频分析广泛应用于多种工程领域。但是传统时频分析方法存在分辨率低或者交叉项干扰。为解决上述问题,提出一种基于L1范数约束的稀疏时频分析方法,从而满足信号局部频谱的稀疏先验。首先阐述了稀疏频谱与短时测量的关系,提出了局部时间频谱反演模型。然后,利用一阶原始对偶方法求解提出模型。通过这种方法,使得重构的频谱变得稀疏。一方面,提出算法的聚集性在L1范数约束下变得稀疏,另一方面,由于提出方法基于短时傅立叶变换和凸优化技术,可以避免交叉项干扰。为了反映算法效果,文章分别对理论信号和实际地震信号进行实验并对比了其他先进的时频分析方法。结果显示,提出方法相比于其他对比算法能获得更加准确的时频谱分布。（5）提出一种基于Lp伪范数和交替乘子迭代法的稀疏时频分析方法。该方法中,视短时截断数据为稀疏表示中的观测信号并设计了一种字典矩阵,建立起短时测量和稀疏频谱之间的关系。基于这种关系和Lp伪范数描述的稀疏约束,稀疏时频表示模型得以建立,然后利用交替乘子迭代法求解提出模型。通过若干合成信号、一道实际地震信号和一组含有天然气的地震剖面进行实验。这些实例都显示提出方法相比于一些先进的时频分析方法能够获得更高分辨率的时频图。因此提出方法对于地震勘探具有重要意义。（6）提出一种基于匹配追踪算法的稀疏时频分析方法。基于匹配追踪方法的稀疏时频表示有效避免了短时傅立叶变换的低分辨率问题,保持了局部信号的稀疏频谱先验。实验表明,提出方法能获得高分辨率时频谱,相比一些先进的时频分析方法更具有竞争优势,这对地震信号谱分解具有重要意义。

论文目录

• 摘要
• abstract
• 缩略词表
• 主要符号表
• 第一章 绪论
•   1.1 研究背景及意义
•   1.2 国内外研究现状和发展态势
•     1.2.1 流体识别方法概述
•     1.2.2 非平稳地震信号预处理技术
•     1.2.3 地震信号时频分析常见方法
•   1.3 本文主要研究内容及技术路线
•   1.4 论文结构安排
• 第二章 信号时频分析与数据预处理基础
•   2.1 离散傅立叶变换及其快速实现
•   2.2 STFT与海森伯格测不准原理
•   2.3 Cohen类时频分布理论简介
•   2.4 稀疏表示理论
•   2.5 数据的稀疏正则化预处理
•   2.6 本章小结
• 第三章 基于交叠组稀疏广义全变分的地震数据时频重建
•   3.1 概述
•   3.2 预备知识
•     3.2.1 二阶TGV模型
•     3.2.2 二维矩阵交叠组稀疏邻近算子
•   3.3 模型及其求解方法
•     3.3.1 基于交叠组稀疏的改进TGV模型
•     3.3.2 模型的ADMM求解
•     3.3.3 模型的加速重启ADMM求解
•   3.4 实验
•     3.4.1 实验平台与评价指标
•     3.4.2 合成地震信号去噪性能对比测试
•     3.4.3 实际资料去噪性能测试
•   3.5 本章小结
• 第四章 基于FrFT和贪婪策略的多向窗时频分析
•   4.1 概述
•   4.2 最优分数域Gabor算法与其局部聚集效应
•   4.3 模糊域多向窗函数设计
•     4.3.1 Cohen分布与Gabor变换的关系
•     4.3.2 模糊域多向窗设计
•   4.4 数值实验
•     4.4.1 模型一：单调频信号
•     4.4.2 模型二：多个线性调频信号组合
•     4.4.3 模型三：非线性调频信号
•   4.5 实际资料处理及应用
•   4.6 本章小结
• 第五章 基于稀疏时频分析的地震信号谱分解
•   5.1 概述
•   5.2 预备知识
• 0和P₁建模'>    5.2.1 稀疏重构的P₀和P₁建模
•     5.2.2 Chambolle-Pock一阶原始对偶方法
•   5.3 提出方法
•   5.4 数值实验
•     5.4.1 评价指标
•     5.4.2 理论信号测试
•   5.5 实际资料处理及应用
•     5.5.1 单道地震信号测试
•     5.5.2 地震信号谱分解中的应用
•   5.6 本章小结
• 第六章 基于Lp伪范数的地震信号稀疏时频分析
•   6.1 概述
•   6.2 提出方法
•     6.2.1 基于Lp伪范数约束的稀疏频谱重构模型
•     6.2.2 基于ADMM和 Lp收缩算子的求解方法
•   6.3 数值实验
•     6.3.1 评价指标
•     6.3.2 模型一:单成分信号
•     6.3.3 模型二:抛物线频率调制信号
•     6.3.4 模型三:多成分信号
•   6.4 地震信号谱分解中的应用
•     6.4.1 单道地震信号测试
•     6.4.2 地震信号谱分解
•   6.5 本章小结
• 第七章 基于匹配追踪的地震信号稀疏时频分析
•   7.1 提出方法
•   7.2 实际资料处理及应用
•     7.2.1 单道地震信号测试
•     7.2.2 地震信号谱分解
•   7.3 本章小结
• 第八章 总结与展望
•   8.1 全文工作总结
•   8.2 主要贡献与创新点
•   8.3 工作展望
• 致谢
• 参考文献
• 攻读博士学位期间的成果

文章来源

类型: 博士论文

作者: 陈颖频

导师: 彭真明

关键词: 非平稳信号,分数阶傅立叶变换,时频分析,稀疏正则约束,稀疏频谱成像

来源: 电子科技大学

年度: 2019

分类: 基础科学,工程科技Ⅰ辑

专业: 地质学,地质学,地球物理学,石油天然气工业,矿业工程

单位: 电子科技大学

基金: 国家自然科学基金项目“复杂地震信号分数域频谱成像理论及应用研究”(项目编号:41274127),项目研究年限为2013年1月至2016年12月

分类号: P631.44;P618.13

总页数: 146

文件大小: 10297K

下载量: 363