起伏地表的自由基准面波场成像方法研究.pdfVIP

维普资讯 石油天然气学报 (江汉石油学院学报) 2008年6月 第30卷 第3期 JournalofOilandGasTechnology(J．JPI) Jun．2008 Vo1．30 No．3 起伏地表的 自由基准面波场成像方法研究 熊晓军 ，贺振华 ，黄德济 (成都理工大学信息工程学院，四JII成都610059) [摘要]基于波动方程基准面延拓的偏移方法理论 ，从地震波场正演模拟 的角度 出发，提 出了一种新的适 用于起伏地形 的地震叠后正演方法。该方法综合利用 了 “地形淹没”和 “零速层偏移” 的思想 ，通过基 准面变换和设置零波场传播层 ，将原始模型的地表作为新模型的一部分，有效地消除了起伏地形的影响， 并采用常规 的波场延拓算子进行波场延拓与波场成像计算，实现 了起伏地表模型的地震叠后数值模拟 。 基于该方法理论 ，还发展 了一套完整的 自由基准面波场延拓流程，可 以将地表观测 的地震波场延拓至任 意的基准面之上 ，对研究野外地震波的传播规律、后续的偏移处理 以及地震资料 的解释具有非常重要 的 意义。仿真实验 的结果均验证 了地震正演模拟方法和 自由基准面波场延拓流程的正确性和有效性。 [关键词]波场成像 ； 自由基准面；起伏地表 ；叠后正演 ；波场延拓 [中图分类号]P631．44 [文献标识码]A [文章编号]i000—9752 (2008)O3—0075—03 传统的静校正方法通过垂向拉伸把地震资料校正到一个假想的水平观测面上，由于以地表一致性假设 为基础，仅考虑地震波的垂向运动特征 ，对地表复杂地区会产生新 的畸变 。近年来 ，随着人们对地震波场 传播机理研究的不断深入 ，提 出了一些基于波动方程 的基准面延拓方法 ，例如 Berryhill(1979)口]基于 Kirchhoff积分解的基准面延拓算法 ，Reshef(1991)2【基于相移法 的 “逐步一累加”波场延拓方法，Bolte (2001)[3不依赖于速度模型的基准面延拓方法 ，以及江凡 (2006)E43等从波场传播算子的角度对 Reshef方 法的改进等等。这些方法从地震偏移的角度 出发，采用波动方程进行计算，兼顾 了地震波传播的运动学和 动力学特征 ，较好地消除了地表起伏的影响。笔者反其道而行之，从地震波场正演模拟的角度 出发，重点 研究了基准面变换与波场成像的机理，提出了一种新的适用于起伏地形的地震叠后正演模拟方法，并以此 为基础，发展了一套波场延拓的计算流程，其可以将地表观测的地震波场向上或下延拓至任意基准面之 卜，为深入研究 目标地质体的反射特征以及后续偏移计算提供了有效的数据支持。 1 理论原理 1．1 起伏地表的地震叠后正演 常规的基于波动方程的地震叠后正演方法采用数值方法直接求解波动方程，包含了地震波传播的动 力学特征 ，波场信息丰富，模拟结果较为准确 ，但是只能处理水平地表的情况 ，波场延拓的基准面必须 为一个水平界面。为了满足实际起伏地形的要求 ，采用矩形框对 目标地质模型进行整体封装 ，矩形框的 顶面即为输出的基准面，其为一个水平界面，且位于起伏地表的最高点之上 。然后，对封装后的地质模 型进行规则的网格剖分，起伏的地表形态被离散成为规则网格的一部分，起伏地表之上的网格采用任意 的速度值和密度值进行充填，于是得到了新的地质模型。 新模型包含了 “地形淹没”的思想 (Bevc，1997)]，其将地表作为模型的一部分参与计算。因此， 在后续的波场延拓与成像计算过程中，需要消除地表的影响。研究借用 “零速度层”的思想 (Beasley， l992)=[6]，通过在输人基准面与输出基准面之间设置一个零波场层参与波场延拓计算 ，不仅可以很好地 处理地形的起伏变化，而且十分灵活方便，其输人基准面与输 出基准面之间的速度和密度可以取任意 [收稿日期]2008—01—19 [基金项目]国家自然科学基金项 目 ；国家 863 计划项目 (2006AA09A102—12)。 [作者简介]熊晓军 (1980一)，勇，2001年大学毕业，博士 ，讲师，现主要从事地震正、反演方法方面的研究工作。 维普资讯 石油天然气学报 (江汉石油学院学报) 2008年 6月 值