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

摘要

本发明是物探中利用不同模稀疏加权的抛物拉东变换压制多次波的方法。先对动校正后的道集做富里叶变换，在频率域中对每一频率分量做不同模的稀疏加权抛物拉东变换，然后对变换得到的谱做反富里叶变换得到时间域的拉东变换谱，识别并切除一次波，接着做富里叶变换得到频率域中的多次波，然后做富里叶反变换得到时间域的多次波，从原始的动校正后的道集中减去得到的多次波得到一次有效波。本发明通过可以克服常规的最小二乘Radon变换得到的τ-q谱中多次波与一次波的能量混在一起，无法完全压制多次波或损害一次有效波的问题，从而更好的压制多次波。

主权项

1.一种利用不同模加权稀疏的抛物拉东变换压制多次波的方法，特点是通过以下技术方案实现：1)地震采集，获得动校正后的叠前共中心点道集；2)对每一个地震道数据做时间参量的富里叶变换，每一频率的所有道的数据，按照炮检距的大小排列，形成数据矩阵；3)按照下述公式计算拉东变换的系数矩阵L和L^H：L=e-jωq1x12e-jωq2x12Λe-jωqnx12e-jωq1x22e-jωq2x22Λe-jωqnx22MMMMe-jωq1xm2e-jωq2xm2Λe-jωqnxm2---(1)]]>式中，ω是角频率；x_i是每一道的炮检距；q_j是拉东变换参数，即抛物线的曲率；矩阵L^H是矩阵L的共轭转置矩阵；4)按照下述的公式计算不同模的稀疏加权矩阵：WM=b+||Mi-1||maxPb+||Mi-1||P---(2)]]>式中，M_i-1是频率域中拉东变换得到的上一频率的谱；||M_i-1||^P是M_i-1的L_P模，P为常数，或L₁或L₂模；是M_i-1的L_P模的最大值；b是稀疏稳定参数，取的1％；5)按照下述的公式进行L_P模加权稀疏抛物拉东变换：(L^IIL+λW_M)M＝L^HD             (3)上述式中，L和L^H是拉东变换的系数矩阵，由公式(1)计算；W_M是L_P模稀疏加权矩阵，由公式(2)计算；D是数据矩阵，由步骤2)计算得到；λ是稀疏参数，取数据道数的5％-20％；M是拉东变换在频率-空间域中的谱，通过公式(2)(3)得到；6)对所有的频率分量重复步骤2)至步骤5)，抽取每一检波点位置处的所有频率分量，做时间参量的富里叶反变换，得到检波点位置的拉东变换在时间-空间域的拉东变换谱(τ-q谱)；7)在拉东变换谱(τ-q谱)中分别识别出一次波和多次波，并切除一次波；8)对切除一次波后的拉东变换谱(τ-q谱)做时间参量的富里叶变换；9)按照下式对切除一次波后的、频率域中的谱进行拉东反变换，得到频率中的多次波数据；D′＝LM′            (4)式中，L是拉东变换系数，由公式(1)得到；M′是步骤8)得到的多次波数据；10)对所有的频率分量重复步骤9)，抽取每一检波点位置处的所有频率分量，做时间参量的富里叶反变换，得到这一检波点位置的时间-空间域的拉东反变换得到的多次波；11)从原始的动校正后的共中心点道集中减去多次波，得到压制多次波后的一次有效波。