基于主分量分析PCA的高分辨率叠加速度分析.pdfVIP

石油天然气学报(江汉石油学院学报) 2006年8月第28卷第4期 JournalofOilandGas Aug．2006V01．28No．4 Technology(J．JPI) 基于主分量分析PCA的高分辨率叠加速度分析 孑!11祥娥 (成都理工大学信息工程学院，。油气藏地质及开发工程’国家重点实验室，I基1)ll成都010050) 钟本善 (成都理工大学信息工程学院。四川成都610059；长江大学电信学院，湖北荆州434023) 周熙襄 (成都理工大学信息工程学院，四JIl成都610059) [摘要]在考虑地震记录的反射同相轴具有一定程度的相关性，而噪声之闯互不相关，同时相邻同相轴与 目的同相轴的相关性存在差异的前提下，以沿着双程旅行时to时间和叠加速度决定的双曲线为中心线截 取数据体，对该敷据体进行主分量分析，将目标同相轴作为主分量，噪声和相邻反射同相轴做为次分量。 并进一步引入加权函数来提高速度扫描时对信号的敏感度，提出了一种基于主分量分量分析的速度检测 准则。通过对人工合成地震记录和实际地震资料的处理，表明本方法能明显地提高叠加速度谱的分辨率。 且具有很强的抗噪能力。 [关键词]叠加速度；主分t分析I高分辨率I速度分析I速度谱 [中围分类号]P631．44 [文献标识码]A 速度是连接地震资料和钻井资料的纽带，也是贯穿地震勘探全过程的一项十分重要的参数，速度选 择的正确与否决定了地震勘探的精度。因此获取准确的速度参数是正确地处理和解释地震资料的中心问 题CZ]。叠加速度分析是各种速度分析中的一个重要环节，其他类型的速度可以通过叠加速度进行估计。 一般的叠加速度分析可以分为叠加类和相关类两大判别准则，叠加类采用叠加能量来进行速度判别，相 关类则采用互相关系数进行速度判别。对低信噪比资料来说，常规速度分析技术得到的速度分辨率不够 高，人们开展了大量的研究工作，提出了许多高分辨率叠加速度分析技术[2叫】。笔者在考虑地震记录的 目的同相轴与相邻同相轴的相关性存在差异，同时噪声互不相关的前提下，提出了一种基于主分量分析 PCA的速度检测准则。 1方法原理 对于给定的数据体z(忌)∈R“，主分量分析(PrimaryComponent 方差矩阵的特征值和特征向量，最大的特征值所对应的特征向量就是随机矢量z的主分量，而其他特征值 所对应的特征向量则是数据体z的次分量，因此主分量对应着协方差矩阵的最大特征值，次分量对应着协 方差矩阵的其他特征值。通过提取数据体z协方差矩阵的特征值，并根据特征值的变化规律，寻找数据体 中的有效信号，找到正确的双曲线轨迹。 假设z是打Xm的样本数据，并且z具有零均值的特性，则z的协方差矩阵C。为： e—z·XT一芝：．；I船 (1) 鬲 式中，．；I；是协方差矩阵C。的特征值，并设|；I。≥|；I。≥…≥．；【。I鳓是对应于特征值．：I；的特征向量，A。对应着 样本z的主分量，．；l。，…，A。对应着样本z的次分量。 对于水平地层，理想CDP道集上各道对同一层的反射波到达时间可视为双曲线分布： t2=t3+事 (2) 式中，z为道偏移距I岛为零偏移距时间It为偏移距为z时反射波的达到时间；口为该层相应的叠加速度。 [收稿日期]2005—11—20 [作者简介】孙祥娥(1970一)，女，1992年大学毕业．副教授。博士生，现主要从事信号处理及其在石油勘探开发中的应用工作。 ·86· 石油天然气学报(江汉石油学院学报) 2006年8月 现在对地震数据沿水平方向做截取，获得数据体xEn，m]，m表示选择的地震道道数，咒表示选择的时 yr-．，m-1进行特征值的提取，提取后的特征值对比曲线如图1所示。 从图1(a)可以看到，对于包含有 双曲线形式的同相轴数据矩阵，其协方 差特征值的分布比较分散，各个特征值 捌 之间的差异也比较小，即主分量和各个 垫