基于表层多次波的地震数据插值方法研究.pdfVIP

第52卷第 6期 石 油 物 探 Vo1．52，No．6 2013年 11月 GE0PHYSICALPR0SPECTING F0RPETR0LEUM Nov．，2013 文章编号-1000—1441(2013)06—0636—06 基于表层多次波的地震数据插值方法研究 郭书娟 ，方伍宝 ，徐兆涛 (1．中国石油化工股份有限公司石油物探技术研究院，江苏南京21]103；2．同济大学海洋与地球 科学学院，上海 200092) 摘要：提出了一种基于原始含表层多次波数据的叠前共炮集地震数据插值新方法。相对于利用相邻道信息变换 或外推插值的缺失地震数据重建，新方法直接利用含表层多次波的原始记录数据，通过互相关构建虚拟一次波， 并在滑动时间一空间窗内对拟地震道进行最小二乘匹配滤波和均方根振幅校正，而后用于缺失近炮检距数据的 插值重建。与已有的利用多次波生成虚拟一次波的方法相 比，新方法无需进行多次波和一次波的提取，所以更 易于实现。模型数据测试和实际资料试处理证明了新方法的有效性和实用性。 关键词：地震数据插值；表层多次波；拟地震道；最小二乘匹配滤波；均方根振幅校正 ID I：10．3969／j．issn．1000—1441．2013．06．012 中图分类号：P631．4 文献标识码 ：A 由于野外采集的条件限制，原始地震数据中往 理中将多次波提取出来并不容易。 往会出现缺道、坏道和数据分布不规则等情况，从 我们提出一种利用本来蕴含在表层多次波数 而严重影响地震资料处理中波动方程偏移等多道 据中的地震信息实现叠前炮集资料近道缺失数 处理算法的精度。因此，如何根据已获取的有用信 据插值的新方法，不用预先提取多次波，而是直 息来恢复缺失的地震信息，成为地震资料处理中急 接利用包含有表层多次波的原始数据，通过互相 需解决的一个问题 。 关构建虚拟一次波，并在滑动时间一空间窗内结 为了解决缺失地震数据的重建问题，有关专家 合最小二乘匹配滤波 (LMF)和均方根 (RMS)振 学者研究了多种地震道插值技术：①波动方程插值 幅校正方法对拟地震道进行校正，而后用于地震 方法11]，包括NMO和反 DMO交替进行的数据 数据插值重建。虽然这种方法构建的准一次波 规则化、AMO方法等；②空间域插值方法[，[1t-x 有一些噪声 ，但是经过校正后用于弥补空缺道是 域预测误差滤波插值等；③频率域插值方法E删 ，如 完全可 以胜任 的。 f-k域道内插值，f-x域预测滤波插值等；④基于 Radon变换的插值方法E书]，如采用抛物线 Radon 1 拟地震道构建 变换、双曲线Radon变换等构建插值算子，进行不 规则采样数据重建等；⑤基于稀疏反演的插值方 采用互相关方法将地震数据多次波路径部分 法[9 ，包括稀疏约束最小二乘反演、非二次模型 包含的地震信息提取出来，构建 “拟地震道”。 权函数的稀疏反演傅里叶重建方法等。这些方法 如图la所示，震源 S点激发的波在地下反射 虽然算法不同，但主要都是利用缺失道的相邻道信 点x反射后传播至接收点A，得到一次波；此一次 息，通过变换或插值外推来填补空 白道。 波在 自由表面反射后继续向下传播，在地下反射点 地震数据中包含有一次反射波和多次反射波， y反射后传播至接收点B，得到表层多次波。那 而通常被看作干扰波的多次波数据 中也蕴含有可 么，A点处地震波的旅行时可表示为 ZSX+TXA，B 利用的地震信息。Shah等[1213]提出通过一次波和 点处地震波的旅行时可表示为TSX+Z-XA+rA+