[发明专利]基于反射槽波频散相似度的断层成像方法

摘要

本发明公开了一种基于反射槽波频散相似度的断层成像方法，通过以下步骤实现：步骤一：时频分解；步骤二：求直达槽波频散曲线；步骤三：消除直达槽波；步骤四：成像。本发明充分利用了槽波的频散性质，依据频散性减少了其他噪声对成像的影响，减少了画弧现象，大大提高了断层成像的精度。

主权项

1.一种基于反射槽波频散相似度的断层成像方法，其特征在于：通过以下步骤实现：步骤一：时频分解，用广义S变换求出各道地震数据的时频谱其中sn(t)为第n道地震数据的时域形式，Sn(τ,f)为其时频谱，n∈[1,N]，N为总道数；步骤二：求直达槽波频散曲线，依据各道的炮点位置(xsn,ysn)与检波点位置(xrn,yrn)，求出直达槽波的传播距离Ln再由Ln对Sn(τ,f)进行坐标转换，可求出各道直达槽波的频散曲线，叠加后获得总体数据中直达槽波的频散曲线C(v,f)步骤三：消除直达槽波，根据Ln将C(v,f)从“速度‑频率”域转换回“时间‑频率”域后，与Sn(τ,f)相减，就可以得到消除了直达槽波后的各道时频谱S’n(τ,f)步骤四：成像，首先选定一个成像区域，区域的一个点P，坐标为(xp,yp)，假设P点是一个反射点，则槽波从第n道数据的炮点(xsn,ysn)出发，在P点反射后，被(xrn,yrn)处的检波器接收到，则传播路径的总长度为如果P点是反射点，则根据Lpn求出的频散曲线Cpn(v,f)应该与C(v,f)一致，因此，将所有道Cpn(v,f)与C(v,f)的相关系数的和Rp作为该成像点的值，当该点确实为反射点时，Rp较大，而当它不是反射点时，Rp较小，对成像区域内所有成像点进行上述处理，得到整个区域的成像结果。