[发明专利]一种利用不同模加权稀疏的抛物拉东变换压制多次波的方法

说明书

技术领域

本发明涉及地球物理勘探技术，属于地震资料处理过程，是一种利用不同模稀疏加权的抛物拉东(Radon)变换压制多次波的方法。

背景技术

在地球物理勘探技术资料处理过程中，多次波压制是地震处理中的一个重点和难点。目前，抛物拉东(Radon)变换是用来压制多次波的常用方法之一。拉东变换压制多次波主要是利用一次波和多次波之间的动校正速度差，在Radon变换域中将二者分开，从而识别出多次波并切除。实际应用中，首先对共中心点(CMP)道集做动校正，然后对校正后的CMP道集进行频率域抛物拉东变换得到拉东变换谱(τ-q谱)，接着在拉东变换谱(τ-q谱)中识别出多次波然后进行拉东反变换就可以得到时-空域的多次波，从原始数据中减去得到的多次波就可以获得一次有效波。

拉东变换压制多次波的效果取决于拉东变换谱中多次波和一次波能量的分离程度。但是在常规的最小二乘抛物拉东变换中，由于对不同的抛物线的曲率范围(q值范围)采用相同的加权因子，因此得到的解不够精确，当多次波和一次波动校时差比较小时，在拉东变换谱(τ-q谱)中多次波和一次波的能量不能完全分离，这就会降低压制多次波的效果，从而影响地震资料成像质量。

发明内容

本发明目的是提供一种提高拉东变换的分辨率和地震资料成像质量的利用不同模加权稀疏的抛物拉东变换压制多次波的方法。

本发明通过以下技术方案实现：

1)地震采集，获得动校正后的叠前共中心点道集；

2)对每一个地震道数据做时间参量的富里叶变换，每一频率的所有道的数据，按照炮检距的大小排列，形成数据矩阵；

步骤2)所述的炮检距是炮点到检波点的距离；

步骤2)所述的数据矩阵是列向量，其列数等于数据的道数；

3)按照下述公式计算拉东变换的系数矩阵L和L^II：

L=e-jωq1x12e-jωq2x12Λe-jωqnx12e-jωq1x22e-jωq2x22Λe-jωqnx22MMMMe-jωq1xm2e-jωq2xm2Λe-jωqnxm2---(1)]]>