基于起伏地表的频率空间域有限差分法叠前深度偏移研究.pdfVIP

基于起伏地表的频率空问域 有限差分法叠前深度偏移 张景涛 李振春 张辉 (中国石油大学地球资源与信息学院，山东东营257061) 摘 要 起伏地形勘探越来越被人们所关注。直接从起伏地表开始的叠前深度偏移方法是对付复杂地表和复 杂地质构造成像的有效手段。本文实现了直接从起伏地表向下延拓的波动方程频率空间域有限差分法叠前深 度偏移，从而达到既能使复杂构造精确成像，又能处理任意起伏地表的目的。数值试算表明，与“波场上延” 法相比，该方法对虚拟层速度要求很低，同时能达到很好的成像效果。 关键词起伏地表频率空间域有限差分法叠前深度偏移 “波场上延”法 随着地震勘探的不断发展，油气勘探的重点正转向起伏地表条件和复杂地质条件的区域。熊翥(2000) 在“中国西部地区物探工作的思考”一文中指出：“现有的偏移程序，大都建立在激发点和接收点位于 同一个水平面上，这与我们需要进行叠前深度偏移处理地区的实际观测条件不相符。过去我们用静校正 技术来解决这个问题，从波场延拓的角度讲，静校正使波场产生了畸变，再作深度偏移时就会产生一系 列的误差，严重影响深度偏移的效果。” 基于起伏地表的叠前深度偏移可以同时解决这种地质上的双重复杂性(复杂地表和复杂地下构造)。 何英等(2002)提出了“波场上延”法，但对虚拟层速度有较高的要求，如果不能准确估算近地表速度， 很难保证成像效果。本文首先介绍“逐步一累加”波场外推的概念和基本思想，然后用波动方程频率空 间域有限差分法实现从复杂地表开始向下延拓的叠前深度偏移，并给出相应的偏移结果。频率空间域有 限差分法叠前深度偏移可以很好地适应横向变速剧烈的介质。 基本理论 1． “逐步一累加”法 面校正和偏移的方法，具体做法是将波场从一个水平基准面向下延拓，并在每一个深度步长上将所截得 的地形面上的波场值加到延拓的波场中，这样就可以对起伏地形上的记录进行叠前或叠后偏移。 2．频率空间域有限差分法叠前深度偏移 从三维声波方程 1 a2Ua2Ua2U a2U ，．、 … 1，2 Ot2 孤2。a、，2。az2 出发，经浮动坐标变换得到： 要=z争0)2一《∥nivco万 协a， 堕iO(1一!诃 (2-b) 一Oz=-iO(--1，一～c皿 211 称(2一a)式为绕射项方程，它可使绕射波收敛；(2-b)为折射项方程，它的作用是校正由于横向变速引 起的时差。为用有限差分法求解(2-a)式，采用有理分式逼近，联立(2-b)式，得到频率一空间域任 意变速情况下的下行波的波场深度外推方程。同理可得到上行波的波场深度外推方程。 模型试算 该模型模拟的是山包沟谷地形(见图1)。模拟的野外放炮方式为：中间放炮两边接收，接收道数 时间采样率At=lms；总炮数为 80炮。 在图2中给出了一个模拟单 炮记录。使用本文提出的方法， 对该模型记录直接进行叠前深 度偏移，结果显示在图3中。从 图中我们可以看出，无论是浅层 还是中深层，地下构造都能被很 好的成像，而且地表形态较为准 确。这个结果与“波场上延”法 的结果相当，而该方法不需要向 上延拓，节省了机时，同时提高 了稳定性。 改变虚拟层速度，比较每次 成像的结果，我们发现该方法对 虚拟层速度要求很低，在虚拟层 速度为500m／s的时候仍能较好地 成像(如图4所示)。 图1速度一深度模型 图2模拟单炮纪录