基于峭度最大化的地震盲反褶积方法.pdfVIP

第51卷第 1期 石 油 物 探 Vol_51，No．1 2012年 1月 GE0PHYSICALPR0SPEC INGHClRPETR0LE1M Jan．，2012 文章编号：1000—1441(2012)01—0030—07 基于峭度最大化的地震盲反褶积方法 蔡连芳，田学民 (中国石油大学(华东)信息与控制工程学院，山东青岛266580) 摘要：从地震褶积与反褶积过程的等效总系统出发 ，使脉冲响应满足单位脉冲序列形式，导出方差约束下反褶积 输出的峭度最大化准则，作为地震盲反褶积的目标函数 ；针对 目标函数关于反褶积算子的优化问题，采用粒子群 算法进行全局寻优操作，实现地震信号的反褶积。数值模拟和实际资料处理结果表明，该方法不但适用于最小 相位子波的反褶积 ，而且适用于混合相位子波的反褶积。与梯度法优化峭度准则的反褶积结果相比，能够更好 地从地震记录中估计反射系数，拓宽地震资料的频谱，提高地震资料的分辨率。 关键词：地震盲反褶积；峭度 ；反褶积算子}粒子群算法；混合相位子波 DOI．10．3969／j．issn．1000—1441．2012．01．004 中图分类号：P631．4 文献标识码 ：A 反褶积是地震资料常规处理的基本环节之一， 化的检测[9- ；而峭度最大化准则的另一个重要用 其 目的是压缩地震子波，提高地震资料的分辨率。 途，即从地震盲反褶积的角度考虑，可以作为反褶 传统反褶积方法通常假设地震子波为最小相位，反 积输出结果的非高斯性度量准则，且无需对反褶积 射系数为白噪分布，采用地震记录的自相关代替子 输出结果进行概率密度估计，而如何建立地震盲反褶 波的自相关，使用基于二阶统计量的方法(如脉冲反 积与峭度最大化准则的对应关系是需要解决的问题。 褶积)实现反褶积，这在实际应用中取得了一定效 我们从地震褶积和反褶积过程的等效总系统 果。但是，实际地下介质中传播的地震子波通常是 考虑，为达到反褶积的目的，使等效总系统脉冲响 混合相位的，地层反射系数也不呈高斯白噪分布，因 应序列为单位脉冲序列，导出方差约束下反褶积输 而传统反褶积方法的假设条件往往与实际情况不 出的峭度最大化作为独立同分布性准则，建立地震 符，这在一定程度上限制了方法的应用范围。 盲反褶积 目标函数；由于目标函数关于反褶积算子 与传统反褶积方法相比，地震盲反褶积方法采 的高度非线性和梯度法优化峭度准则的局部性，采 用数据的高阶统计信息，不需要 (或者弱化)对地震 用粒子群算 法 (ParticleSwarm Optimization， 子波和反射系数的先验假设，基本可以和实际情况 PSO)对峭度准则进行全局寻优操作，对地震信号 保持一致 是近几年来学者们重点关注和广泛研究 进行反褶积处理。该方法不需要对子波做最小相 的反褶积方法。针对不同的理论基础和适用条件， 位假设，可以避免对反褶积输出结果进行概率密度 地震盲反褶积方法大致可分为两类。一类是基于 估计，能够适用于最小相位子波或者混合相位子波 反射系数的稀疏性假设和 Bernoulli-Gaussian分 的反褶积处理，具有更加广泛的应用范围。 布的数学模型，采用 Bayesian稀疏反演理论的盲 反褶积方法[i-43，该类方法采用的反射系数模型过 1 基本原理 于严格，用稀疏尖脉冲对反射系数进行模拟，而实 1．1 地震褶积与反褶积数学模型 际地层反射系数的变化往往 比稀疏尖脉冲复杂。 假设地层反射系数序列 为零均值、独立同分 另一类是基于反射系数独立同分布和非高斯性的 布的