引 言
在地面记录的未经偏移的反射地震能量(例如共偏移距记录)通常可提供地下构造的强相干图像。尽管是强相干的,但该图象是不正确的,原因在于它受到几种畸变作用的影响,最显著的影响是来自地质界面截断的绕射波和倾斜界面上反射点与地面位置间能量的横向移动。地震成像就是校正这些畸变影响的地震处理手段。
在常规处理中,地震成像是众多处理手段的最终阶段。全部地震处理过程包括能量增益化、反褶积、静校正和速度分析、NMO和DMO、叠加,最后到主要用于构造成像的偏移。根据偏移数据体绘制地震构造图,然后在构造图上布置钻探目的层的勘探井位。这种做法一直延续至今,但是偏移后的数据已逐渐被用于后续进一步处理的输入数据。例如,当地质学家和工程师们联合攻关已做过3D地震工作的油田的油气开发问题时,他们把偏移后的地震数据体用作地震属性分析处理的输入数据,诸如瞬时相位或振幅随偏移距变化的处理中。通过已知井与建议井位处地震属性特征的对比,我们就能预测出建议井位处的岩性特征。例如,倘若在几口井中某一属性可能都与孔隙度观测值的相关性很好,那么这个属性就可用于估算其它井位处的孔隙度。多年来,地球物理学家都在用地震数据在比地震波长小得多的尺度上估算岩石物性特征,但直到最近他们才对偏移数据实施这种处理,它代表了偏移应用的一个全新领域。这也促使我们在偏移中要非常重视象振幅这样的参数的应用。
时间和深度、叠后和叠前、2D和3D
地球物理学家习惯使用时间域记录的地震数据。以往,他们乐于解释时间记录,甚至是偏移后的时间记录。这就与在深度域工作的地质学家和油藏工程师产生了根本性的矛盾,同时也给许多通过设计偏移算法来摸索地球物理学规律的数学家和物理学家带来困惑。尽管深度域的解释结果明显会更正确,但迫使地球物理学家这样做会妨碍他们解释地震数据的能力,例如在资料质量差的地区,在偏移前后剖面上难以确定一一对应的反射波同相轴。所以,就有了时间域偏移(时间偏移)和深度域偏移(深度偏移),它们之间的差异远不是表面上的。通常时间偏移比深度偏移更容易实现,且已证实时间偏移对偏移速度的精度要求比深度偏移要低得多。这样,在精度要求和投入成本低的时候,我们就可以采用时间域偏移而非深度域偏移方法。另一方面,当我们试图利用偏移不仅生成地下图像,而且还要估算地层速度时,时间偏移就不如深度偏移有效了。
理想情况下,应该对采集到的全部地震道实施偏移,而且我们也正朝着这个目标快步迈进。但是,即使在几年之前,全部采集道的偏移——叠前偏移——几乎是一个梦想,很简单的原因在于成本太高难以实现。直到近来提出了一个标准的处理步骤,即在偏移前通过叠加减少数据体的大小,那就是对大量的道实施动校正和倾角时差校正,然后将这些道相加成一个以地表某一位置为参考的单道。许多道叠加成一道的处理过程模似了这样一种试验,即用零偏移距替代多偏移距,由此获得偏移成本大为降低的数据体。多少年来的实践已经证明叠后偏移处理是很成功的,尤其是在构造适度复杂地区的成像。另一方面,我们早已认识到在一些射线弯曲严重的地区并不服从叠加处理中的假设条件,以致地震成像质量很差。换句话说,横向速度变化很大的地区,在叠前就应该进行偏移处理。需作叠前偏移的典型例子是墨西哥湾的盐体之下的构造,由于墨西哥湾地区的速度变化极大,在叠后偏移处理出的结果中几乎见不到盐下反射同相轴。
需作叠前偏移的这些盐体还须做三维(3D)的成像。如果能采集到穿过径向形状盐丘的一条单倾向地震测线,通过二维(2D)偏移就能正确成像;但是如果未能较好地确定倾向或走向方向,则穿过沉积围岩的扭曲流状深根型盐体往往被偏移成一扁平走向的盐体。利用2D偏移方法时,盐下构造成像不能完全消除来自采集面以外的能量,且由于波传播到采集面之外而使得采集面内缺少成像所需的反射能量。
叠前时间或深度偏移代表了当今最新的应用水平。海上地震采集主要是3D采集,可以说目前开展三维叠前偏移的大多数工作量是作海洋资料偏移。与陆上采集条件正好相反,海上采集主要受地下成像目标的控制,而受海面条件限制较少。所以,我们一般对陆地数据是作2D成像,特别是山区数据;当地下变化剧烈时,则应作3D成像,因为2D叠前偏移不能解决问题。但采集3D陆上数据时,往往因为经费紧张而使得采集的数据集没有足够的覆盖次数来完成合适的静校正或速度分析。其结果是我们从陆上地震数据得到高质量成像的能力通常落后于海洋数据。
偏移方法
大多数勘探地球物理学家对偏移做什么都有所了解,且许多地球物理学家知道哪几种方法可用于偏移,但只有很少的地球物理学家熟悉偏移方法的本质特征。描述所有实用的偏移方法已超出了本文的讨论范围,所以,我们将仅介绍当今对地震数据成像有较大影响的几种代表性的方法。Gardner(1985)曾详细论述了所有这些方法;在其论著的再版卷中讨论的大多数方法为叠后偏移。叠前偏移方法远非仅仅是叠后方法的“拓展”,而且涉及到其自身概念上的重大难度。在叠前偏移中这种难度与需要处理两个波场有关(对于某个炮记录,一个波场来自炮点位置,另一个波场来自记录排列的所有检波点位置),而叠后偏移仅处理单个波场(零偏移距未偏移剖面)。
为了描述这些方法,我们需要指出支撑偏移的模型的线性性,即仅考虑一次波的声波或弹性波方程。假定没有发生多次散射(多次散射为地震能量从震源传播到一个地下绕射点,再到第二个绕射点,最后到检波器),在下面的两种方式中就可以观察到这种线性性。第一种方式是将地下看作为一组分离的散射点,而连续的反射层是由这些反射点的集合构建起来的。在这种方式下,很容易看到如果另外增加一个散射点到地下层位上,那么在未偏移记录中就会另外增加一个确定的散射同相轴。这样,在考虑散射成像时,一个时间只用一个散射点。第二种可观察到线性性的方式是将未偏移的地震记录看作为一组分离的尖脉冲信号,而连续反射同相轴则由这些脉冲的集合所构成。每一个脉冲都来自地下一个连续反射层,根据从震源到反射层的任意点再到检波点的双程旅行时等于脉冲的记录时间来确定成像位置(例如,在常速介质中,产生未偏移响应脉冲的点的集合是一个旋转的椭球,其焦点在震源和检波点位置)。增添一个尖脉冲到未偏移记录中,将导致另一个偏移“脉冲响应”被加到偏移剖面中。令人吃惊的是,波动方程提供的数学手段能消除所有偏移画弧的噪声,仅留下未扭曲的地质图像。
最熟悉或许是最容易理解的方法是Kirchhoff偏移。这种方法是在原有的空间—时间域内处理地震数据,且在作一些初始的滤波后,直接地实施前述两种方式的任意一种。Kirchhoff偏移能用于时间偏移或深度偏移,且它是非常灵活的,即能够将未偏记录的任意一组采样偏移到偏移记录上的任意一组采样。由于3D叠前偏移非常之依赖于这种能力,因此在现行的多数大规模偏移中,它仍是优选的偏移方法,但这并不意味着它是叠后偏移中最精确、最快捷的方法。Kirchhoff偏移能用于高陡倾角构造的成像,但是在三维偏移中选用足以对陡倾角构造成像的大偏移孔径时,其高昂的处理成本也是无法接受的。
在叠后偏移中仍旧流行的早期数字偏移,产生于对波动方程的有限差分近似,即只允许波场向下传播(相对于可向上、向下传播的全波动方程而言)。现有几种采用显式和隐式数值计算的有限差分偏移方法,它们的共同特点是一个完整的波场向下延拓(并不象Kirchhoff偏移那样选道偏移),即对单程(向下传播)波动方程采用有限差分近似,波从记录面穿过一个中间厚度层再到成像深度点。这些方法能在时间域或深度域中实施。当倾角适度时,这些方法相当精确,但当倾角很陡时则就不那么精确了。
最快捷的偏移方法是在f—k(频率—波数)域运行的。这些方法将记录波场分解成平面波分量,并分别偏移每个分量。象有限差分法一样,f—k方法也是从记录面到成像深度向下延拓波场。这些方法对波动方程不采用有限差分近似,而是采用相移法实现向下延拓。即使在倾角很陡的情况下这些方法也很精确,但是它们不能很好地适合横向速度变化很大的情况。而将其扩展到适应横向速度变化情况是可能的,但此时则不再是最快而变成了最慢的方法之一。
最后,一种采用全波场双程波动方程的偏移方法称作逆时偏移。从记录到的波场开始,将某个时间的时间切片逆时传播进入地下,即依据波动方程,波场从切片时间点反向传播。通过这种方式,逆时偏移构筑了地球内部每个时间步的完整波场,以至于我们能看到能量快照从记录面向反射点移动。相对其它方法而言,逆时偏移对地震波传播方程的近似较少,因而它是现有最精确的偏移方法。很遗憾的是,通常它也是最慢的方法。然而,它已经成功地用作为标定其它偏移技术的概念模型,以及作为具有很好速度控制情况下极度复杂地区的最终偏移手段。
我们看到似乎有太多的偏移方法。为什么有如此之多的方法呢?在这里我们仅介绍了其中的几种,此外还有其它许多常用的方法。如此多的偏移方法在迷惑了许多地球物理学家的同时,也已证实在解决我们遇到的各种成像问题时是很有用的。一些方法适合于高陡倾角复杂地质构造的成像,而其它方法更适合于低成本的高分辨率地层成像。由于Kirchhoff偏移的灵活性,它适合于处理很多问题,但是对于许多偏移问题来说,它既不是最快也不是最精确的方法。我们可能会提出这样的问题:“为什么有如此多的偏移方法?”或换句方式问“较少的方法是否会更好些呢?”考虑到我们面对着大型成像项目的高额经费投入,和由于成像失败而导致干井的巨大资金浪费,对第二个问题的回答必须是强有力的“不!”
速度估计
为了实现地震数据的成像,偏移处理需要已知从震源和检波点到地下位置的旅行时,而要确定旅行时就需要速度。通常时间偏移是利用地面位置和成像位置间的一种平均负载速度,也可称之为成像速度;而深度偏移则是利用一种更详尽的层速度函数。对于叠后偏移,这些函数通常完全能够从叠加速度分析得到的速度可靠地估计出来。对于叠前偏移就不能采用叠加速度了,而必须利用其偏移处理本身来估算速度。对于这种处理,深度偏移比时间偏移难度要大得多。已经提出了几种用于叠前深度偏移的层速度估算方法。这些方法的最简单方式是在每个水平位置按垂直方向修改速度,即假定某个地下位置的偏移后数据仅与其位置正上方的速度有关。常用的最完善的速度估计技术为地震层析成像,即在每个成像位置求偏移后记录间的代数最小误差。象偏移方法一样,每种求取速度的方法都有其用途,即简单的方法适合于弱至适度横向速度变化的地区。然而,没有一种方法能非常适合于数据质量差和速度剧烈变化的地区,而恰恰我们最需要很适应这些复杂地区的速度估算方法。
面向未来
地震成像比其它大多数地震处理的计算量要大得多,故它一直紧随着计算机硬件技术的发展而进步。至今我们仍在对大量的叠后数据作偏移,这就意味着叠前偏移的计算成本仍然太高,因而我们还需要具有更大内存和更快磁盘存取的高速计算机来降低处理成本。然而,我们已经从在特制的阵列计算机上运行2D叠后偏移就需要花费数天的时代前进了一大步;我们永远告别了第一代偏移,那时的偏移是用量尺和罗盘来机械地实现的。
过去几年里,地震成像已在主体标量计算机、向量阵列处理机、向量超级计算机、并行超级计算机和工作站上运行和应用。再过几年,我们将会看到在个人计算机上完成的偏移处理工作量将与日俱增,甚至可作大规模的3D叠前偏移,因为现在的个人计算机具有10年前超级计算机的运行速度,且内存更大得多。偏移处理能力的提高应促使我们更有效地解决偏移速度的估算问题(常常为各向异性),以避免因偏移迭代次数多而增加成本。随着偏移处理能力的提高和速度估算技术的改进,3D叠前偏移将由几乎不可能实现的王国(10年前)过渡到常规处理王国(几年后)。
正如前文所述,地震偏移成像正逐渐被用作后续处理的输入数据,即用其估算3D数据体中基于地震属性的岩石性质。随着对这项技术的进一步了解,我们就会意识到成像方法中保持地震振幅的重要性。这将使我们既能分析诸如反射系数随偏移距变化这类的简单属性,又能分析比现在更可靠、更复杂的新属性参数。
可以预计,弹性波成像的应用将会与日俱增。朝这一方向努力的第一步是实现转换波成像。有关的工作已获得了一些成功,例如在那些P波记录上很模糊的气云区中用转换波数据成像获得了成功,但我们还需要开发满足全波动方程的矢量波场方法。未来必将证实,当今最慢的声波(单分量)偏移方法——逆时偏移,将成为矢量波场成像的最简单、最快速的方法。