保幅共成像道集波动方程偏移.docVIP

保幅共成像道集波动方程偏移 美国斯坦福大学 Paul Sava，Biondo Biondi，劳伦斯伯克利国家实验室 Sergey Fomel 概要 通过向下延拓偏移，我们提出两种方法来计算角度域共成像点道集，而且我们分析他们的反射角振幅响应（AVA）。两种方法的简单实施导致冲突，因此振幅响应显然是不精确的。振幅问题相关的事实是，向下延拓偏移是伴随向上延拓建模的，但这是他的反演差的近似。我们在频率波数域的对角线上得到加权算子，使偏移对建模的反演有一个好的近似。加权后，两种方法计算的角度域共成像道集就一致了。在实际情况下，降低AVA精确性的其他重要因素是有限的采样和偏移范围，和地震数据带限的本质。 引言 传统上，偏移速度分析和AVO用于偏移域共成像点道集，因为零偏移距图像不描述大多数相关信息。但是，用波动方程偏移很难产生这些道集，因为向下延拓数据的偏移尺度随深度而收缩。对这个问题的解决办法是用角度道集来代替偏移道集。通过波动方程偏移获得的角度域共成像点道集在测量速度模型的精确性是非常有力的。ADCLGs也是有吸引力的，因为他们为振幅分析提供更多的直接信息，也就是说，AVA代替普遍的AVO分析。 角度域共成像点道集是由反射点处的入射角排序的地震图像的陈述。角道集可以使用波动方程技术，由 Bruin等提出的激发剖面偏移，或是由Prucha等提出的炮点检波器偏移获得。在这两种情况下，在成像之前，向下延拓波场使用倾斜叠加估计角度道集。成像前分解向下延拓波场产生角度道集作为偏移射线参数代替真正的反射角的函数。角度域道集也能用代替向下延拓波场的倾斜叠加的图像来计算。我们表明这种代替过程直接产生作为反射角函数的角道集。在这两种下，倾斜叠加变换可以很容易的由频率波数域的射线追踪变换变换来执行（Ottolini, 1982）。 AVA分析是重要的，作为一个角度函数，这个过程用于计算ADCIGs反射保幅。它是如此费解以致用这两种方法计算ADCIGs产生相反的振幅，而且偏移也不是相应的向上延拓建模的一个好的近似的方法。为了解决这一难题，我们必须考虑加权函数，它是在偏移过程中由成像过程引进的。这一加权由频率波数域的对角线算子很好的近似。因为计算ADCIGs的两种方法在不同的阶段执行倾斜叠加（一个在成像之前，另一个在成像之后），相应的加权函数是不同的。一旦考虑了加权，这两种方法产生的AVA反应是一致的，偏移能做模型的近似反演。根据反射数据的物理模型，加权能使偏移成为线性反演的一个很好的近 图1 在任意速度介质中的一组反射射线 似。我们采用Stolt and Benson (1986)提出的物理模型，并为计算ADCIGs的两种方法定义合适的加权。如果需要一个迭代估计过程，通过在建模和偏移期间应用加权的平方根也能很容易使其伪统一。 图2例如在图像空间角度和偏移域之间的合成翻转。左图：合成角道集。中图：从角到偏移距的转换。右图：角度域转换。 两种方法计算ADCIGS 在频率域何以很方便的形成角道集。如果我们认为在常速介质中，t是从震源到反射面再返回到地表检波器的旅行时，2h是震源和检波器之间的偏移距，z是反射点的深度，α是地层倾角，γ是反射角(图1)，我们能写方程： 从等式（1）和（2）我们得到： 方程（3）是由常速介质得到的，但在很小意义上它是非常有效的对任意速度介质，如果我们考虑h是反射面深度的有效偏移距，不是表面偏移距（图1）。 在频率波数域，公式（3）取微分形式： 这个等式表明角度道集可以很方便的在频率域偏移算法的帮助下形成 (Stolt, 1978)。此外，波动方程偏移是最适合使用如此一个方法计算角道集，因为反射深度的偏移距精确的描述偏移输出，这是一个模型参数。而不是用表面偏移距描述偏移输出，这是一个数据参数 (Biondi, 1999)。 我们能知道方程（1）只描述在反射率方面波传播的射线参数。使用定义： 因此我们能写一个和（4）相似的关系式来估计在傅里叶域的偏移射线参数： 方程（4）和（6）能用来计算图像道集，通过在傅里叶域德的径向记录道变换。主要的不同时方程（4）是在偏移成像的空间运算的，而方程（6）是在数据空间运算的。这两个方法在其他3个方面也不同： 首先，图像空间的方法与偏移是完全不挂钩的，因此反射角的转换被认为是偏移后的一个处理。这样的后处理是很有趣的，因为它允许从角度域回到偏移域的转换而不用重新偏移（图2），当然，这是不正确的数据空间的方法，其中的转换时数据频率的函数。 其次，从方程（1）可知，如下的偏移射线参数ph也是一个构造倾角α的函数，这是在图像空间不准确估计的反射角。我们用方程（4）得到的角是几何措施，完全独立于构造倾斜。因为AVA的目的，