基于高阶统计约束实现同态盲地震反褶积.pdfVIP

2013年2月 石Ii由址球幽理参攘 第48卷 第1期 · 处理技术 · 文章编号 ：1000—7210(2013)01—0031—06 基于高阶统计约束实现同态盲地震反褶积 孙建业④ 高 伟 ② 刘怀山① 刘喜武③ (①中国海洋大学海洋地球科学学院，山东青岛 266100；② 国家深海基地管理 中心 ， 山东青 岛 266061；③中N石油化工股份有限公司勘探开发研究院，北京 100083) 孙建业 ，高伟 。刘怀 山，刘喜武．基于高阶统计约束实现 同态盲地震反褶积．石油地球物理勘探 ，2013，48(1)： 31～ 36 摘要 本文通过引入高阶统计学的独立变量分析算法(IcA)，利用双谱估计出的地震子波作为约束条件来实现 同态盲地震反褶积 。首先在无噪声的条件下，将地震记录 由时间域变换到 同态域 ，利用相邻两道地震道作为组 合，将地震褶积模 型转换为线性混合 ICA模 型；然后 ，利用高阶统计学 中的双谱估计地震子波 ，利用估计 出的地 震子波信息作为约束条件应用于快速独立变量分析算法(FastICA)中，实现地震子波和反射系数在 同态域的分 离；最后将分离的地震子波和反射系数再反变换到时间域，即可得到反射系数和地震子波。模拟和实际地震数 据 的数值算例均表 明，在不对反射系数作高斯 白噪假设及子波最小相位假设的条件下，基于高阶统计约束的同 态盲反褶积方法能够取得较好 的反褶积效果。 关键词 高阶统计 独立变量分析 双谱 盲反褶积 同态域 中图分类号：P631 文献标识码 ：A 足而引起的误差_4]。同态反褶积虽然对地震子波不 1 引言 做最小相位假设 ，利用对数谱平均提取子波直接分 离出地震子波和反射系数_5]，但是在用同态反褶积 地震反褶积基本上是一个盲过程口]，因为反射 方法估算 出的若干个地震子波 中，如何选取一个最 系数未知，通常地表爆炸激发的地震子波也未知 ，但 佳的地震子波却是 比较困难 的问题 。因为人工挑选 在地震记录上呈现的反射地震波又是重叠 的，而地 不仅需要花费大量的时间，而且还需要一定的经验， 下地层的厚度通常要小于地震子波 的波长 ，也无法 挑选结果往往 因人而异 。Oppenheim[6首先提出处 从地震记录中分离出地震子波。在地震数据处理中 理非线性系统的同态反褶积方法 ，在 图像信号和语 应用 的传统反褶积方法通常依靠线性滤波 (如 ，维纳 音信号的处理方面已获得了一定的成功 。在解决地 滤波、预测反褶积)或最优滤波来估计子波特性 ，往 震信号处理方面 ，国内外的一些学者也进行了研究 。 往要做统计性假设 。例如，假设反射系数是 白噪序 Schaferl7提 出了相位 曲线展开法 ，但该方法 比较复 列 ，地震子波为最小相位等 ]。实 际应用 中，无法 杂 ，而且 由于用此法展开的相位 曲线不规则 ，在最后 保证这些假设条件总是正确的。因此，在无上述假设 分离 出的反射系数 中会产生 明显 的附加噪声 。 条件下进行地震反褶积或盲反褶积就显得十分重要。 Stoffa等 ．9j还提出了先求相位 曲线 的微分值而后 长期 以来 ，人们一直在研究对地震子波不做最 积分 的方法 ，将 z变换对数谱展开成泰勒级数 ，但 小相位假设的反褶积方法 ，推出了一些混合相位地 该法假设反射系数序列是最小相位 的。俞寿朋_l】 震子波反褶积方法 ，同态反褶积就是其 中一种。与 对子波相位的展开进行了探讨 ，周兴元_5为克服实 传统 的反褶