基于弹性矢量波叠前逆时偏移方法.pdfVIP

l ·577· 中国地球物理201 (17)油气田与煤田地球物理勘探 基于弹性矢量波的叠前逆时偏移方法 on Prestackreversetime methodbasedelasticvectorwave migration 张慧‘李振春 HuiLiZhenchun Zhang 中国石油大学(华东)地球资源与信息学院青岛266555 1．前言 近年来，随着多波多分量地震勘探的不断深入和发展，人们对弹性波叠前偏移理论与技术给予了极大 的重视与关注，尤其是弹性波叠前逆时偏移。弹性波逆时偏移能够对多波数据进行真振幅成像，获取多个 独立的弹性波运动学和动力学属性，并有可能使成像结果为弹性波振幅随炮检距变化规律(AVO)分析和 弹性波叠前反演所应用，因此弹性波叠前逆时偏移是目前主要发展的地震叠前偏移成像技术。 2．方法原理 对弹性波波场偏移的最佳方法是在波场重构和成像时使用全矢量波场。弹性波逆时偏移算法就是以矢 量波理论为基础的深度域偏移成像方法，它使用弹性波方程的数值解，并且将垂直分量和水平分量接收到 的波场值同时作为边界条件的输入。基于Whitmore提出的逆时偏移的基本物理思想，弹性波逆时偏移可 以视为弹性波正演的逆过程，主要包括两部分内容：多分量波场记录的一致性外推和成像条件的运用。用 于正演模拟的有限差分算法同样可使用于逆时外推，因此弹性逆时偏移算法也面临着计算量大和计算效率 低的问题。针对这一问题，本文采用可变网格进行逆时延拓。利用可变网格进行逆时延拓，所采用的是与 传统交错网格算法一样的过程，只是在波场外推的过程中添加波场在全局粗网格和局部精细网格之间的转 移。具体步骤如下：第一步，进行震源波场的正向外推。对每一全局时刻，f都判断过渡区域的波场值，若 此时刻过渡区域波场值控制函数为零，则按照常规网格的外推算法进行计算。若此时刻过渡区域波场值控 制函数不为零，则将之后的每一全局时间波场进行精细化处理，在空间网格步长变化处采用变网格算法求 解两个全局时间层平面的边界部分，并利用插值公式计算两个全局时间层之间的精细时间采样平面的边界 值。然后采用高阶有限差分算法求解每一局部时间层内的精细网格剖分区域的内部网格点波场值。最后在 计算结束后，将局部时间层的精细网格剖分区域波场值返回到全局时间层内。第二步，进行检波点波场的 逆向外推，作为正演的逆过程，过程同第一步算法一致。 3．成像条件 目前，叠前逆时偏移所用的成像条件主要激发时间成像条件和零延迟互相关成像条件。Kselin和 以消除由震源引起的一部分成像假象，而且在成像值中消除了震源的影响，更接近与真实的反射系数。震 源归一化互相关成像条件的数学表达式为： mage(x，z)一[一s(z，z，，)月(工，z，t)dt／[…S2(x，z，t)dt 4。认识和结论 弹性逆时偏移采用双程波动方程．适用于任意复杂的速度模型，无倾角限制，是目前解决复杂构造成 像问题最有效的手段。在弹性波场外推过程中引入时空双变网格算法，提高了逆时偏移的计算效率。利用 震源归一化互相关成像条件，补偿了地震波逆时延拓过程中的能量损失，提高了对深层构造的成像能力。 从模型测试结果和模拟结果看出，本文的弹性逆时偏移算法能够准确的成像出地下反射界面的真实位置和 形态，是一种正确的高效算法。 参考文献 A．．Elastic Q reverse-time 【l】Changw．F．，McMechan migration[J]．Geophysics，1987b，52：1365～1375． 【2】张慧，李振春．基于双变网格算法的地震波正演模拟【J】．地球物理学报，2011。54(0：77—86． 【3】何兵寿