目前,叠前深度偏移已经广泛应用于生产中,并且对成像的精度要求越来越高。在断坳叠置盆地,工区内断层发育,地层倾角大,从而导致地震波的传播路径变得极为复杂,绕射波和断面波异常发育。此外,构造形态复杂多变,地层和断面之间难以识别和区分,速度横向变化剧烈,局部垂向存在速度反转,目的层速度难以确定,这就给建立合理的速度模型带来了很大的难度,所以如何获得合理的、高精度的速度模型就成为获得高质量复杂构造成像的关键因素。
1 网格层析速度模型建立地震层析反演速度分析就是利用层析反演技术对地震道集中存在剩余时差进行反演处理,进而得到速度模型更新量。基于射线理论的层析反演称为射线层析。在给定初始慢度场的条件下,层析反演速度分析利用射线追踪模拟地震波的传播路径并将共成像道集的剩余时差沿射线路径进行更新,就可以得到慢度更新量和更新后的慢度场,完成速度模型的一次迭代更新。在数学上可以表示为沿着射线路径的Radon变换:
其中s和r分别是激发点和接收点,L(s,r)是从s到r的射线路径,f(x,y,z)是剩余慢度,τ是成像道集的剩余时差。当慢度是离散网格值时,上述积分变为累加求和,公式如下:
ΔIj表示当前射线在j网格内的长度。每一条射线都建立一个方程,所有射线的方程组成一个庞大的稀疏线性方程组:
这样射线层析反问题就变成了迭代的线性问题。通过求解该稀疏方程组,就可以得到慢度或剩余慢度。
通过上述层析速度分析原理的论述,可知进行层析速度反演首先要有一个较好的初始速度模型,进而通过叠前深度偏移得到偏移道集和剖面,再对成像道集拾取剩余深度差,建立层析方程,最后通过射线追踪利用迭代法求解离散的大型稀疏方程组,得到速度修正量,进而更新初始速度模型ÿ