各向异性介质叠前逆时偏移.pdfVIP

·744· 中国地球物理2011 (20)地震波传播与成像探查 各向异性介质叠前逆时偏移 ReverseTime Pre-stack in Medium Anisotropy Migration 马婷1李庆春’1周学明2 Ma Li Zhou Ting QingchunXueming 1长安大学地质工程与测绘学院 西安710054；2铁道第三勘察设计院集团有限公司天津300251 偏移成像问题历来是地震勘探的核心问题，随着成像目标越来越复杂，特别是在勘探中遇到金属矿体、 盐体及盐下成像等，基于射线及基于单程波动方程的偏移成像方法已不能很好地对这种构造成像，而逆时 偏移则可以比较好的解决这类问题。基于双程波动方程的逆时偏移完全遵守波动方程，可以对所有传播路 径的波进行成像，被认为是最为精确的偏移成像方法。虽然P．S波场分离是合理的，但是在P．S波干涉严 重的情况下，多分量数据的处理以矢量形式将更为合理，而不是将每种波场按类似P波的处理方法进行处 理。本文采用交错网格高阶差分格式完成矢量形式的VTI及TTI各向异性介质弹性波叠前逆时偏移，对逆 时偏移中假象的消除进行了分析研究。 1．叠前逆时偏移方法原理及实现过程 逆时偏移的成像理论基础是时间一致性成像原理，即反射面存在于地层内下行波波至时间与某一上行 波波至时间相符一致之处。叠前逆时偏移的成像就是在震源波场和接收波场外推到相同时刻时应用成像条 件即可，分为三个部分：一是震源波场的正向模拟：二是接收波场的逆时外推；三是成像条件的应用。分 别在正演模拟及逆时传播计算中使用各向异性介质中的弹性波方程即完成了各向异性介质叠前逆时偏移。 震源的正演模拟根据使用的成像条件而不同，这里采用互相关成像条件，这时需要计算并保存每个时间步 长的震源波场。接收波场的逆时外推是波的逆时传播过程，实质上是求解波动方程的边值问题。接收波场逆时 外推的具体实现从接收的最大时间开始，每逆时外推一步就将该时刻接收点记录到的波场值插入到接收点位置 (当记录的采样率与外推时间步长不一致时可以进行插值或者在外推到一个采样率时间时才插入波场值)，作为 这个时刻接收点位置的波场值，再进行下一时刻的外推，一直进行到t=0为止，完成整个接收波场的外推过程。 对于互相关成像条件，每逆时外推一步，则将该时刻的波场与正演模拟中对应时刻的波场进行互相关 成像，然后将成像值累加至成像剖面上，直到接收波场外推完成得到最终的偏移成像剖面。 2．逆时偏移假象的消除 逆时偏移的一个重要缺陷是在成像剖面中会出现假象，这主要是由于逆时偏移使用双程波动方程。波 在正演模拟和逆时传播过程中遇到分界面都会产生反射。在应用成像条件时就会将这些反射也作为成像值 加入到成像剖面中，从而产生假象，这些假象通常表现为低频干扰。特别是在浅层出现强反射界面的时候， 产生的干扰尤为明显，常常会掩盖干扰区域的成像。通常为消除这些假象干扰，可以通过在偏移前的预处 理，成像过程中的优化成像条件和在成像后对剖面进行滤波或波场重建来完成。 本文根据能量密度矢量的方向性判断波场的传播方向进行优化成像。通过结合能流矢量密度的互相关成像 条件可以根据地下构造的情况选择不同传播方向的波进行成像，可以对复杂构造体进行理想的成像，如盐丘侧 面，盐下构造成像等，同时能极大地消除逆时偏移中由于波的正演模拟和逆向传播在反射面的反射波和透射波 形成的干扰。本文采用的是双程波动方程进行正向和逆向传播，然后再根据能流密度矢量判断波的传播方向，保 证了所有路径波的传播，并且能流密度矢量在波场正演模拟传播过程中可以直接计算，不会增加过多的计算时间。 3．结论 论文研究了VTI介质与TTI介质中弹性波的叠前逆时偏移成像及假象干扰的消除等相关问题，采用结 合能流矢量密度的互相关成像条件可以极大的消除假象干扰。用简绕射点模型、水平层状模型及比较复杂 的断层模型及金属矿模型进行弹性波叠前逆时偏移成像试算，均取得了很好的成像效果，并用断