起伏地表叠前波场正演与偏移成像.pdfVIP

40 内蒙古石油化工 2008年第2期 起伏地表叠前波场正演与偏移成像 薛明星，贺振华 (成都理工大学信息工程学院，四川成都 610059) 摘 要：本文借鉴Reshef在偏移中提出的“逐步一累加”法的思路，基于惠更斯原理，实现了起伏地 表条件下的叠前单程波方程的波场正演．对正演数据的偏移结果，表明了该方法的正确性。 关键词：起伏地表；叠前正演；单程波方程；偏移成像 陆上地震采集工作往往在地表起伏的地区进 作为该位置新震源的强度．继续这一过程到最大深 行，特别是西部地区，其地表条件十分复杂，分布着 度．然后，从最大深度开始将新震源向上延拓，最终 广阔的沙漠、戈壁、山地等．如何真实地模拟野外地 到达地面形成记录，此过程相当于自激自收剖面的 震波的传播、有效地反映地下构造的有效波场信息， 正演．该方法对波场延拓算子的选择无限制，可根据 如反射波和绕射波等的特征，已成为当前人们关注 需要任意选用某一算子。 的一个焦点[1]。目前，常规的地震叠前正演方法主要 1．2 起伏地表的处理 有射线追踪法和全程波动方程法．射线追踪法主要 在偏移处理中，以色列学者Reshef首次提出 考虑地震波传播的运动学特征，可以灵活地处理起 “逐步一累加”波场外推的概念，即在采用深度外推 伏的地表条件，所得的地震波的传播时间比较准确， 时进行基准面校正和偏移．具体做法是将波场从一 但计算结果很难保持地震波的动力学特征，而且对 个水平基准面向下延拓，并在每一个深度步长上将 复杂的地质构造会出现盲区。全程波动方程法包含 所截得的地形面上的波场值加到延拓的波场中，这 了地震波场的全部信息，可以处理复杂的地质构造， 样就可以对起伏地形上的记录进行叠前或叠后偏 但是计算效率低，正演记录受多次波干扰，而且假定 移。借鉴这一思路，本文的正演方法如下：首先找到 地表为水平界面，不能处理起伏地表的情况．本文借 起伏地形的最高点和最低点，将这两点所在水平面 鉴Reshef在偏移中提出的“逐步一累加”法的思 之间的部分无速度区域使用常速度填充，该速度应 路D ]，基于惠更斯原理，实现了起伏地表条件下的 接近近地表速度。填充之后的地表不规则区域便可 叠前单程波方程的波场正演。该方法可以根据地质 以使用如前所述的方法进行正演计算。在波场延拓 构造的复杂性，灵活地选取各种单程波场延拓算子 的过程中记录下各个地表检波点处的波场值便得到 进行计算，具有很强的适用性。两组数值模拟实验的 地震记录。 结果表明该方法不仅具有很快的计算速度，而且可 2 数值算例 以获得没有直达波、多次波和其他干扰的高信噪比 设计了两个起伏地表模型进行共炮记录的正演 的共炮点叠前地震记录。 计算以及偏移成像。图1中的地下构造为一平界面。 1 方法原理 模型采用128*128的方式剖分，炮点位于模型起伏 1．1 正演原理 地表的正中。接收方式为全128道接收，道间距2O 本文的正演方法基于惠更斯原理生成共炮道 米。图2为此单炮记录，由该图可以看出，由于地表 集．所谓惠更斯原理，即是认为地下任一反射点在炮 的不规则导致双曲线形的同相轴发生畸变，在隆起 点震源扰动的影响下形成二次震源，所有新震源发 的地表下面地震波的走时明显增加．图3是对图2的 出的球面波向上传播后叠加构成地震记录．首先，将