利用反射波资料实现透射波层析速度反演的方法研究.pdfVIP

利用反射波资料实现透射波 层析速度反演的方法 姚姚刘江平 (中国地质大学，湖北武汉430074) 层析速度反演是地震资料数字处理中的一种重要方法。它可以得到地下速度参数的精细分布。应 用于地震勘探的各个方面。例如：井间层析速度反演的结果用于储层预测，反射波层析速度反演的结果 用于叠前深度偏移，折射波、初至波层析速度反演的结果用于静校正等。上述反演方法的原理基本相同， 不同之处仅是所利用的波型相异。众所周知，理论上已经证明：利用反射波、折射波、初至波的层析反 演都存在多解性问题，而利用透射波的层析反演稳定性强，其解是唯一的。但是，在地震勘探的实践中， 地面反射波资料是大量的，透射波资料很少。本文提出一种利用地面反射波资料提取透射波到达时信息， 进而进行透射波层析速度反演的新方法。该方法既可以利用大量存在的地面反射资料，又可以使用成熟 的利用透射波的层析速度反演算法。理论模型资料和实测资料试验处理的结果得到了令人满意的效果。 利用反射波资料提取透射波到达时信息 从绕射的观点看问题，反射波地面炮集记录是由地下各绕射点所产生的二次绕射波叠加组成的，如 图1-a所示。根据炮点、绕射点的位置，利用波动方程法或其他方法，可以将二次绕射波收敛到其顶点 位置，如图1-b所示。一个炮集可以得到一道由炮点出发到达绕射点的记录。一条测线有多个炮集，由 这些炮集可以得到多个炮点到同一个绕射点的多道记录，如图1-e所示。根据互换原理，这样的记录也 相当于由绕射点出发到达各个检波点(即原炮点)的绕射波(透射波)记录，如图1-d所示。上述过程 表述了由地面炮集记录得到透射波到达时间信息的基本原理。 b c d 图1利用地面反射波资料提取透射波到达时间信息的基本原理 当然，实际计算时不是如此简单，因为要将绕射波收敛到其顶点位置，需要知道绕射点上方介质的 速度，而这个速度正是我们层析反演欲求的解。为了解决这一问题，需要使用迭代的方法计算。首先， 假设绕射点上方介质为均匀介质，给定一个初始速度，按此速度计算时距曲线进行绕射波的收敛运算， 得到如图1-d所示的一个道集。如果给定的速度正确，则该道集上必定出现与给定速度一致的时距曲线； 175 否则，出现的时距曲线与按给定速度计算的时距曲线不一致。出现后一种情况时，求取二时距曲线上的 各道时间差。据此时间差修正按给定速度计算的时距曲线。再按修正后的时距曲线收敛绕射波，再比较 修正后的时距曲线和新得到的道集上的时距曲线。如果一致，则叠代结束；否则继续进行叠代，直至一 致为止。最后得到的修正后的时距曲线上时间值即为要求取的由绕射点到各检波点(即原炮点)的绕射 波(透射波)运行时间值。 透射波层析速度反演 透射波层析速度反演可分为计算上覆介质综合速度和速度结构两类。采用双曲线拟合方法可以反演 上覆介质综合速度。采用透射波射线层析速度反演方法可以获取上覆介质的速度结构。 1．双曲线拟合综合速度反演 当一定范围内横向物性较均匀、速度变化不大时，可以认为由绕射点至各接收点(原炮点)的透射 波到达时所构成的时距曲线满足双曲线方程。采用最小二乘法拟合时距曲线，即可获得上覆介质在绕射 点附近的综合速度。 2．透射波射线层析速度反演 如果速度横向变化较大，绕射点至各接收点(原炮点)的透射波时距曲线不是双曲线，则只能使用 透射波射线层析速度反演算法计算上覆介质的速度结构。目前，在地震层析速度反演算法中，较常用的 等。这些算法各有优缺点。BPT算法重建图像比较粗糙，分辨率比较低，但计算快速简单；SIRT算法 具有较高的成像精度，并能明显克服个别数据误差较大所造成的结果失真和由于射线分布不均所造成的 误差集中等缺点。我们在层析速度反演中综合利用了BPT算法和SIRT算法。首先使用BPT算法进行求 解；将求解的结果作为SIRT算法的迭代初值；并且在迭代重建SIRT过程中利用界面资料进行约束； 从而构成有约束的透射波层析速度反演方法。该方法既解决了数据投影角度不全和数据不均匀覆盖等问 题，又加快了收敛速度，提高了层析速度反演的精度和分辨率