基于快速Curvelet变化和压缩感知技术地震数据插值方法.pdfVIP

202 SPG/SEG 深圳 2011 国际地球物理会议论文集 基于快速 Curvelet 变化和压缩感知技术的地 震数据插值方法 邵 彧 詹 毅 李合群 （中国石油东方地球物理公司物探技术研究中心） 摘要 本文阐述了一种基于快速Curvelet变换和压缩感知技术的地震数据插值方法，此方法利用 Curvelet变换多尺度、多方向、局部性，及可处理非稳态信号等特性，把地震数据按频率、角度和空 间坐标在Curvelet域进行分解，并根据压缩感知理论在变换空间对数据进行插值重建。基于稀疏变换 的缺失地震数据插值的基本原理可认为是在反问题正则化理论框架下，在L 范数约束的求L 范数最小 2 1 值的过程，输出结果为均匀的完整数据。此方法采用变角度变频率阈值估计和自适应阈值滤波来提升 插值处理的效果，降低残差。最后将此插值方法应用于理论模型和实际地震数据的插值试算，进而证 明其有效性。 关键词 Curvelet 变化 压缩感知 稀疏反演 插值 震反射波信号的物理性质，通过数据之间的相 1 前言 关作用，实现对海上拖缆数据，OBS数据的插 野外采集时由于环境复杂和检波器故障等 值和外推。该类方法的特点是处理过程完全基 原因，会使采集到的数据中出现废道。室内处 于采集资料。（2 ）基于预测误差滤波器(PEFs) 理环节中，通常会将出现的废道剔除，进而影 的插值方法：这类方法主要利用预测误差滤波 响了覆盖次数的均匀性。目前常规叠前数据的 器在倾角谱的处理，实现缺道的插值恢复。其 插值方法，应用最广泛的主要还是基于 NMO 中具有代表性的方法是Antoine Guitton 和Jon 的插值方法。此类方法的基本原理是：首先对 Claerbout在 2009 年提出的在金字塔域的插值 实际采集数据做 NMO 校正拉平，然后采用水 方法[4] 。（3 ）基于求解某目标函数最小化：这 平叠加方法进行插值恢复，最后做反 NMO 变 类方法属于图像处理的方法，具有代表性的分 换。这种方法的不足之处是：正反变换并不是 别是Liu 和Sacchi 于 2005 年提出的Minimum 完全可逆的，同时由于存在动校拉伸的问题， Weighted Norm Interpolation 方法 [6] ，以及 容易造成原始信号的改变。 Herrmann 和Hennenfent 2008 年提出的基于 为了解决以上问题，近几年全球很多科研 Curvelet变换的办法[5] 。前一种方法不需要知道 机构都在不断研究新的插值方法，这些较新的 地下介质速度，它首先对采集数据进行傅里叶 方法大概可以分为以下三类：（1）基于地震波 变换，然后通过加权求取采集数据的L p 范数并 反射信号物理性质的插值方法：这类方法主要 使其最小化，经过若干次迭代最终得到插值结 基于地震波干涉原理，其中具有代表性的是 果。后一种方法则基于压缩感知理论，此方法 Wang 、Luo和Schuster于 2009 年提出的干涉插 通过将采集信号投影到某个空间（比如说 值方法[7] 。此类方法的基本原理是：利用了地 Fourie