VTI介质qP波广义高阶屏单程传播算子.pdfVIP

维普资讯 2007年12月 石油址球幽理参攘 第42卷 第6期 · 处理方法 · VTI介质qP波广义高阶屏单程传播算子 吴国忱 ① 梁 锴① 王华忠② (① 中国石油大学(华东)地球资源与信息学院；② 同济大学海洋与地球科学学院) 吴国忱 ，梁锴 ，王华忠．VTI介质 qP波广义高阶屏单程传播算子．石油地球物理勘探 ，2007，42(6)：6404650 摘要 本文从VTI介质弹性波波动方程入手，通过本征值方法及声学理论近似建立了VTI介质 qP波频散关 系方程 ，再根据介质扰动理论深入研究了VTI介质 qP波广义高阶屏单程传播算子。将 VTI介质参数场分解 为均匀各 向异性速度场 、速度扰动场 以及各 向异性参数扰动场 ，这样 qP波广义高阶屏单程传播算子就 由背景 介质波场 自由传播项和各 向异性扰动项两部分组成 。背景场 的计算在波数域完成 ，而非均匀各 向异性扰动项 (包括慢度扰动项、各 向异性参数扰动项)在空间域计算。当各 向异性参数横 向变化较为剧烈时，采用高阶屏算 子能提高单程 qP波高角度传播精度 。均匀各 向异性介质和复杂各 向异性介质模型脉冲响应试验表 明：qP波广 义高阶屏算子不仅在大角度传播时能够满足精度要求 ，而且能够适应任意强横 向变速、强横 向各 向异性扰动时 的情况 ，适用于复杂 VT1介质qP波地震深度偏移。 关键词 各 向异性 广义高阶屏 单程传播算子 频散关系 面聚焦效果变差 ]。目前各向异性地震偏移成像 已 1 引言 成为石油勘探领域研究热点问题之一 。 在复杂介质深度偏移成像过程 中，单程波传播 地震反射波是地下地质体信息的综合反映，它 算子 (波场外推算子)是地震偏移成像的核心。近年 包含地质体边界和地质体岩性组分两方面的信息。 来，多域的广义屏波场延拓算子得到了广泛的研究， 确定地质体边界主要采用地震偏移处理技术 ，它是 主要包括拟屏算子、复屏算子、加窗屏算子、高阶屏 研究地下复杂地质体的基础 。波动方程深度偏移方 算子等[7q]。这类算法的基本思路是将速度场分解 法由于能够适应强横向速度变化 ，成像精度高 ，能处 为层 内常速背景和层 内变速扰动 ，因为求解背景场 理振幅和多走时问题，所以获得了广泛而深入的研 相当于求解常速声波方程 ，可通过相移实现 ；求解变 究与应用。波动方程深度偏移主要包括基于射线理 速扰动可认为这种非均匀性相 当于散射源，入射波 论的Kirchhoff积分法、时间一空 间域 的有 限差分 场作用于这些散射源后 ，产生散射波场。算法的具 法 以及频率一波数域算法。总之 ，波动方程深度偏 体实现是先将 问题转换到频率域 ，然后利用快速傅 移技术代表着地震偏移处理的发展方向[1]。 里叶变换 (FFT)在空间域和波数域之间往返转换计 地球介质的各向异性在地震勘探中是广泛存在 算 ，实现多域 (空间域和波数域)操作运算 。这些单 的，大部分沉积岩石都表现为观测上 的弱各 向异 程波传播算子一般包括三项 ：自由传播项 、相位屏项 性[2]。目前研究最为广泛的各向异性介质是具有垂 和高角度补偿项。自由传播是在波数域通过一个具 直对称轴的横向各 向同性介质 (VTI介质)，一般 由 有某种参考速度的均匀介质进行，参考速度可以随 周期性薄互层组成 。在地震成像时忽略各向异性的 深度变化 。其余两项体现波和介质非均匀部分 的相 影响，会导致 目标区的空间归位误差增大 ，在对倾斜 互作用，一般是