一种适用于复杂构造成象的平面波叠前深度偏移方法.pdfVIP

一种适用于复杂构造成象的平面波 叠前深度偏移方法 孙沛勇 张叔伦 (石油地球物理勘探局研究院)(大连理工大学数学系) 摘 要 本文提出了一种斯的利用最大能量震源渡场的平面菠叠前深度儡穆方涪。试方法保留了 平面泣儡移方法的优点。在成泉过程中，通过同时廷拓额域稀疏呆样的平面渡蔓豫渡场和最 大巍量旅行时，求取震潭波场的蕞大能量部分．从而显著提高了平面波叠前深度偏穆方法的效 率。利用震源渡埽的最大能量部分成策，保证了一定的成泉精度。Ma脚ou《模型的试算结果 表明，在复杂构造情况下．车文提出的方法可获褥同计算全波场的平面渡叠前深度偏移方法类 似的良好戚泉结果。 引 言 基于波场直接延拓的非Kirehhoff积分类型的叠前深度偏移方法正越来越受到重视。这 类方法具有理论简单、最有高频近似、成象精度高的优点。但同一般积分法相比存在着计算 成本偏高的明显缺点。 Berkhom首先提出的“面炮”偏移方法，利用面积震源和合成记录的向下延拓完成成象 过程。同单炮偏移方法相比，由于减少了叠前数据量，因而显著降低了叠前深度偏移的计算 量，实质上是一种基于平面波台成的偏移方法。利用优化系数的傅里叶有限差分延拓算子 实现震源波场可控的平面波叠前深度偏移。本文的目的是在保证一定的平面渡叠前深度偏 移精度的前提下，进一步提高算法效率。 我们采取的思路是：在直角坐标系下，通过同时延拓频域稀疏采样的平面波震源波场和 最大能量旅行时，计算出平面波震源渡场中的最大能量部分，并在成象过程中用它代替震源 全波场。运用此种方法对Marmousi模型试算，可获得同计算全波场的平面渡偏移方法可比 拟的成象结果。 一、基本原理 1．平面波成象 平面波偏移中利用合成算子r(南)在地表却利用地袭的震源矩阵s(zo)和原始记录矩 95 算子，成象时对s(“)，P(矗)应用成象准则。 2．傅里叶有限差分延拓算子 本文波场延拓中采用了傅里叶有限差分算子。傅里叶有限差分延拓算子是一种混合算 子，兼具有相移和有限差分算子的优点，由3部分组成：包括传统的相移算子，筒单的相移 修正算子、类似45度方程的差分算子。利用连分式展开得到傅里叶有限差分延拓算子．并 对算子作了必要的优化。 3．最大能量平面波震诹波场的求取 频域重采样后，在时间域震源波场将发生混迭现象。如果波场混迭的周期大于原震源 波场的时间长度，在混迭后的波场中可以保留原波场。因此为了估计最大能量震源渡场，应 当在延拓频域稀疏采样的震源波场的同时，延拓最大能量旅行时，以便估计出最大能量波场 的大致分布范围。 为了从混迭的平面渡震源波场中估计出最大能量部分，本文采取如下简单的步骤： 采样的震源波场s(z。)。 (2)对深度五处横向位置xI依据深度Z。。；处的旅行时扫，估计出最大能量旅行时々。 (3)对深度z。横向位置≈处震源渡场s，(z．)作反付氏变换。求取以‘f为中心的时窗内 的波场。 (4)在求取的时窗内计算出最大能量旅行时“。 ^h (5)以‘，为中心选取成象用的最大能量震源波场sJ(五)，再将s，(瓦)变换至频域。在 A 获得了震源波场s(厶)后，利用延拓后的记录波场P(z。)即可以成象。 二、Marmousi模型试算 度偏移算法的精度。 图2是Marmousi模型利用全渡场震源的平面波偏移结果。 图3是仅利用相同频带的20个等间隔稀疏频率采样的平面波最大能量偏移结果。 个平面渡震源，在深度1450m和2200m合成平面波震源。成象中均使用了最小平方准则。 两图的偏移结果很相似，说明在较好地估计出平面渡最大能量震源渡场的前提下，可以仅延 拓少量频率的震源波场，以较高的效率，获得良好的成象结果。 ≯j兰；?