基于高阶双谱的子波整形方法.pdfVIP

第 39卷第5期 物 探 与 化 探 Vo1．39．No．5 2015年 10月 GEOPHYSICAL ＆ GEOCHEMICALEXPLORATION Oct．，2015 d0i：10．11720／wtyht．2015．5．24 杨振，文鹏飞，彭璐，等．基于高阶双谱的子波整形方法 [J]．物探与化探，2015，39(5)：1027—1031．http：／／doi．org／lO．11720／wtyht．2015．5．24 Yangz。WenPF，PengL，eta1．Amethodofwaveletshapingbasedonthehigh—orderbispectmm[J]．GeophysicalandGeochemicalExploration，2015，39 (5)：1027-1031．http：／／doi．org／lO．11720／wtyht．2015．5．24 基于高阶双谱的子波整形方法 杨振 r，文鹏飞 ，彭璐 ，张宝金 (1．国土资源部海底矿产资源重点实验室，广东 广州 510075；2．广州海洋地质调查局，广东 广州 510075；3．中国石化石油物探技术研究院，江苏 南京 211103) 摘 要：地震子波在地震数据处理中扮演着至关重要的角色，传统的子波提取都是假设子波为最小相位的，但是实 际数据 中的子波却是混合相位的。为了能够得到一个较为准确、真实的地震子波 ，笔者采用高阶双谱进行子波重 构．然后利用得到的地震子波对叠后数据进行子波整形处理。仿真实验和实际资料试算结果表明，高阶双谱提取 的地震子波与实际情况吻合度更高 ，利用所提取的子波对叠后数据进行处理 ，可改善剖面的分辨率。 关键词 ：高阶双谱；俞氏子波 ；子波整形；地震资料 ；数据重构 中图分类号 ：P631．4 文献标识码：A 文章编号：1000—8918(2015)05-1027—05 众所周知，在地震勘探中，地震子波反映的是地 1 方法原理 震波的传播特征，在地震数据处理中扮演着至关重 要的角色 。准确的子波估计是地震资料反褶积处 前人研究表明，高阶统计量不但继承了传统二 理、偏移、特征提取、波阻抗反演、地震资料解释和地 阶统计方法的优点，还具有三个性质：非高斯、统计 震正演模型的基础…。如何能够精确地估计子波 独立的随机序列，它的高阶统计量为常数；高斯噪声 是提高地震资料分辨率的重要保证。对于经典的子 的高阶统计量为零；随机信号的高阶统计量能够保 波提取方法(自相关法、希尔伯特变换法)[2-3]，它们 留它的相位信息。由此可见，高阶统计量可以提取 都是基于子波为最小相位假设和噪声为高斯 白噪声 混合相位子波[。 假设。研究表明，实际地震资料的地震子波大部分 1．1 子波的双谱 是混合相位的[4]，这违背了经典算法的基本假设， 实际地震记录Y(n)，由褶积模型可以表示为 所以根据经典算法提取的子波也是不准确的。随着 Y(n)= (n)+ (n)= (n) r(n)+ (n)， (1) 勘探精度的提高，常用的二阶统计方法 (自相关、功 式中：w(n)为地震子波，r(n)为反射系数，(n)为噪 率谱)提取的子波难以满足实际需求。 声。 提取子波的一个重要用途就是进行子波整形处 地震记录Y(n)的三阶相关函数，即三阶统计 理 [5]。以此来提高地震资料的分辨率。如何能够 量_10]可近似表示为 统计一个准确的子波就显得至关重要。基于高阶统 C3(1，r2)=C3(1，丁2)= 计量的基本理论，利用高阶双谱估算子波的振幅谱 1 Ⅳ