强横向变速介质中的叠前深度偏移算法的分析.pdf

摘要 九十年代以来，地震勘探逐渐向中深层和复杂构造区域的精细勘探发展。 在这种形势下，作为一种重要的地震数据处理和解释手段，地震偏移迫切需要 新的发展以解决复杂地质体成像问题。实践证明常规的叠后时间偏移在复杂介 质下，很难取得好的成像效果。于是，人们提出使用叠前深度偏移。 常规的叠前深度偏移算法可以分为两类：一类是基于射线追踪的Kirchhoff 积分算法，另一类是基于波场延拓的波动方程解法(如有限差分偏移、相移偏 移和广义屏偏移等)。由于Kirchhoff积分偏移计算效率非常高，它在石油工业 领域得到了广泛的应用。然而，利用程函方程求解初至旅行时的Kirchhoff积分 偏移在成像效果上不如一些波动方程偏移方法。这主要是因为通过焦散区的射 线由于存在多路径而产生的后至波可能携带重要的能量，这部分信息对深度成 像而言是不应当忽略的。尽管基于完全Green函数射线理论(例如在计算中考 虑多至走时及相应的振幅)的Kirchhofr偏移的成像质量有所提高。但不幸的是， 这类算法非常复杂，计算代价特别大，且对焦散区并非完全有效。 基于波动方程的递归偏移方法对Kirchhoff积分偏移有着得天独厚的优势。 首先，它们是从全方程导出的，不是基于高频近似的渐进解。因此，波动方程 偏移方法潜在地比较精确、稳健。其次，当向下延拓方法可用于波场延拓而又 不增加计算维数时(例如零偏移距数据)，这类方法比Kirchhoff积分法的效率 还高。 波场延拓算子是递归偏移方法的关键，它们一般是从单程波方程推导出来 的。通常大家较为关心的单程波方程偏移的典型算法是有限差分方法和付立叶 方法。前者容易处理速度的横向变化，但其缺点在于存在频散和成像倾角限制。 后者不存在频散且对水平层状介质能精确成像，不过，它只适用于水平层状地 口 ，墨。 本文采用共轭梯度法优化傍轴近似方程的系数，提出了一种频率空间域的 有限差分波场延拓方法，并在反射系数估算意义下，推导出了叠前深度偏移的 成像公式。并且发现，当深度偏移算子中的折射项方程采用时移处理，并且和 成像计算中关于坐标变换复原的时移处理合并在一起时，计算可大大简化。 为了利用付立叶方法的优势，许多地球物理学家提出在双域(即频率波数 域和频率空间域)进行地震波成像处理，先后提出了分步付立叶偏移、付立叶 有限差分偏移等算法。论文详细介绍了这两种偏移算法的基本原理，分析了波 场延拓算子的相对误差，而且还对付立叶有限差分偏移算子提出了优化改进的 的思路。 波的传播和成像问题是有着密切联系的。De Hoop(1988)等人为了研究 R 波在随机介质中的传播问题，基于波的散射理论提出了早期的屏方法。Wu L．J．等在近十年来发展了广义屏方法，并且用于叠前深度偏移成像。 S．Huang 论文对比了各种广义屏方法，而且对它们之间的相互关系和稳定性进行了讨论。 递归偏移方法统一起来，并阐述了它们之间的相互联系。借助于脉冲响应和复 杂模型偏移的数值计算结果，文中进一步对比了频率空间域有限差分法、付立 叶有限差分法、分步傅立叶(相屏)法和广义屏方法各自的特点，为我们解决 不同的实际问题选择合适的偏移方法提供了依据。最后对波动方程三维叠前深 度偏移的思路提出了探索性观点。 Abstract Since hasevolvedtowardsrefined of 90’S，seismic exploration processiondeep zoneor bodies．Underthenew complexgeologic data is for methodofseismic and cried processinginterpretation,seismicmigrati