相似系数加权拉冬变换及其在多次波压制中的应用.pdfVIP

相似系数加权拉冬变换及其在多次波压制中的应用 Bancroft 曹志红。John CREWES T2NIN4 Project，UniversityofCalgary,Calgary，AB，Canada 1．引言 海洋勘探中多次波的存在严重地困扰着地震勘探工作。近些年拉冬变换在多次波压制中得到广泛的应用。 但拉冬变换的分辨率问题一直难以解决，地球物理学家们为此做出了不屑的努力。本文介绍了以相似系数加权 的拉冬变换，大大提高了算法的分辨率。 2．方法 传统的拉冬变换定义为沿着一特定的曲线将某种物理属性的叠加。地震勘探中，在反射界面水平情况下， 反射波在炮集或者CMP道集上表现为双曲线，因此拉冬变换公式可以如下： u(r，p)=∑d(f2=r2+px2，X) (1) 其中烈f’曲是炮集或者CMP道集，H(瓦P)是拉冬变换，t是双程旅行时，x是偏移距，f是零偏移距双程旅 行时，P=l／P：，是曲线斜率，‰是RMS速度。在理想情况下，反射波在拉冬变换里为一个点。同一时间 到达的多次波速度通常比一次反射速度小，其拉冬变换里的P值不同，所以可以互相分开。但由于传统拉冬变 换分辨率低，当多次波与一次反射的速度相差不大时，很难将二者分开。图l所示为包括两个一次反射和两个 多次波的CMP道集，图2为其由传统方法计算的拉冬变换。 Ng(2004)提出了利用高斯迭代法计算的相似系数加权的拉冬变换： 2=f2+px2 (2) “(『，p)=sp，p废d0 X) f∑d02=f2+px2，训 相似系数刮‘p)定义为：s(r，p)2手专支i≯虿‘葫 ‘3) Ⅳ^是计算相似系数时所应用的地震道数，，是平滑窗口长度。相似系数本身是对反射波同相轴曲线的一种估算， 能加强那些与同相轴相似度较大的双曲线对应的拉冬变换，而减弱与同相轴相似度小的双曲线对应的拉冬变换， 从而提高拉冬变换的分辨率。图3所示为相似系数加权法计算的拉冬变换。与图2相比，分辨率大大提高。 借助Sacchi(1995)提出的以先前的拉冬变换计算结果为先验条件再一次计算拉冬变换的思路，本文提出 了以相似系数为先验条件来计算拉冬变换。根据相似系数搜索与同相轴最吻合的双曲线所对应的P值，计算其 拉冬变换，然后将该P值对应的能量从CMP道集中减去；然后再计算下一个最吻合双曲线对应P值的拉冬变换。 利用先到先得的原理提高拉冬变换的分辨率。图4是该算法计算出的拉冬变换，其分辨率明显好于前：者。本 文提出的方法与Sacchi(1995)的方法相比，优点在于相似系数不受振幅影响，而只与反射波同相轴形状有关。 3．拉冬变换在多次波压制中的应用 为了便于在多次波压制中的应用，该算法更适于应用在动校正后的CMP道集上(图5)。此时P与RMS ·329· 速度的关系发生了变化，为P=l／‰一l／矿P，‰为多次波的RMS速度，％为一次反射的RMS速度。一次反 射的P值都为零，而多次波P值都不为零。其拉冬变换为图6所示。将一次反射能量切除，保留多次波专在时 间域内恢复多次波，图7专最后从输入剖面中压制多次波，得到一次反射如图8所示。 4．结论 相似系数加权及先验条件的应用大大提高了拉冬变换的分辨率。同时因为是在时间域内实现拉冬变换，避 免了在频率域内因采用双曲线作为叠加曲线而大大增加的计算费用。 5．参考文献 resolutionRadontransforminthet-xdomain：CSEGabstract．2004． Ng，M．，2004，High and Sacchi，M．D．，and offset reconstructio