基本曲波变换和压缩感知的不完备地震数据恢复.pdfVIP

中国地球物理201l ·547· (17)油气田与煤田地球物理勘探 基本曲波变换和压缩感知的不完备地震数据恢复 of seismicdatabasedoncurvelettransformand Recoveryincomplete compressedsensing 马坚伟1唐刚2汤文3 MaJianwei Wen TangGangTang 1．哈尔滨工业大学理学院应用数学研究所 哈尔滨150001；2．中国石油勘探开发研究院石油物探技术研 究所北京100083；3．清华大学航天航空学院地震波勘探开发研究所北京100084 1．研究背景 由于地形条件或者经济上的限制，地震勘探数据经常被进行不规则的(或者稀疏的)采样。这种不规 则采样与数据丢失的问题主要发生在三维地震勘探当中，在二维地震勘探当中也会偶尔出现。在陆地上， 这种情况可能起因于地震接收器的布置受到建筑、湖泊等复杂地形限制，也可能是起因于坏道或者严重污 染的道的出现。在海上勘探中，则可能由电缆的水平偏移等原因引起。由于许多现行的地震数据处理算法 都不能处理不规则采样的数据，并且从经济性和分辨率的提高等方面考虑，由较少的地震观测获得较多的 地震勘探资料具有重要的实际意义。 实际工作中处理不规则采样数据的方法是：①线形插值出邻近道：②忽略空缺道；③在三维地震数据 处理中不考虑实际中的不规则采样点，通过叠加将所有的道放在一个面元上，这些相对简单的处理结果是 常常会出现严重的假频，使大多数常规技术不能得到很好的运行。因而，寻找最优的地震数据规则化重建 方法，对获得的具有随机丢失地震道的数据进行恢复，使其成为密集而规则的排列，在地震数据处理、成 像中是至关重要的一步。 2．研究现状 地震数据规则化重建问题已经有了较深入研究，提出了许多方法，大致可分为以下五类：①数据投影 方法：包括偏移距延拓以及炮点延拓。②空间域插值方法：校正压制假频后利用采样定理的Sine函数插值、 最大相干倾角插值、功率谱差异倾斜相干插值等。③频率与插值方法：．肛域预测滤波插值、．肛域投影滤 波插值等。④基于特殊变换的方法：非均匀离散傅立叶变换，Radon变换，小波变换。⑤基于波动方程的 插值方法：NMO和反DMO交替进行数据规则化，AMO方法等。 3．本文内容 类方向小波变换，具有多尺度与多方向性的特点，它可以对地震数据进行稀疏表示【lJ。压缩感知【2】理论于 于经典Nyquist采样率的随机下采样中高精度地恢复出原信号。目前该理论已经广泛应用用医学成像、遥 感雷达、视频编码、科学计算、地震勘探等领域。对地震数据采集而言，压缩感知理论为我们提供了一个 “经济”的方案，即在数据采集端便以压缩的形式采集很少的数据，从而大大降低采集数据的成本和作业 时间，而将成本转移到数据恢复的计算中去，这在实际应用中很有意义。 根据压缩感知理论，先建立Curvelet稀疏约束的Ll范数最小化数学模型，然后用迭代阈值法或交替 方向法等快速凸优化方法求解该数学模型，并得到了良好的恢复效果13l。我们也对Lp范数非凸优化模型进 行的地震数据恢复的探讨。为了更有效的减少地震数据的获取，我们还比较了各种不同的测量矩阵(采样 格式)，提出了较为优化的伪随机采样方法14J。 参考文献 curvelet 18～133． transform[J】．IEEESignalProcess．Mag．。2010，27(2)：l 【l】J．Ma,G．Plonka，Thc from reconstruction uncertaintyprinciples：Exactsignal