[发明专利]衰减近地表散射波的方法

权利要求书

1.一种衰减近地表散射波的方法，其特征在于，包括以下步骤：

输入能量均衡的炮集波场和直达波场；

根据直达波地震干涉测量技术对所述炮集波场和所述直达波场进行处理，获得时域近地表散射波场；

用多道自适应匹配滤波技术将所述时域近地表散射波场从所述炮集波场中滤除。

2.根据权利要求1所述的衰减近地表散射波的方法，其特征在于，所述根据直达波地震干涉测量技术对所述炮集波场和所述直达波场进行处理，获得时域近地表散射波场，具体包括以下步骤：

输入所述炮集波场和所述直达波场，分别进行傅里叶变换，从而得到频域炮集波场和频域直达波场；

根据公式ps1(rA,rB,ω)=1NSΣrSp0(rA,rS,ω)p*(rB,rS,ω)p(rB,rS,ω)p*(rB,rS,ω)]]>对所述频域炮集波场及所述频域直达波场的近源道进行反褶积型地震干涉测量计算，其中，p_s1(r_A,r_B,ω)为近源道频域近地表散射波场，p₀(r_A,r_S,ω)为频域直达波场，p(r_B,r_S,ω)为频域炮集波场，p^*(r_B,r_S,ω)为p(r_B,r_S,ω)的复共轭，ω为圆频率，r_S为震源位置，r_A和r_B为接收点位置，NS为叠加所用的炮数，r_B为离r_S最近的接收点位置，r_S在r_B和r_A的单边位置上；

根据公式ps2(rA,rB,ω)=1NSΣrSp0(rA,rS,ω)p(rB,rS,ω)p(rB,rS,ω)p*(rB,rS,ω)]]>对所述频域炮集波场及所述频域直达波场的远源道进行褶积型地震干涉测量计算，其中，p_s2(r_A,r_B,ω)为远源道频域近地表散射波场，其与所述p_s1(r_A,r_B,ω)共同构成一个完整的频域近地表散射波场，p₀(r_A,r_S,ω)为频域直达波场，p(r_B,r_S,ω)为频域炮集波场，p^*(r_B,r_S,ω)为p(r_B,r_S,ω)的复共轭，ω为圆频率，r_S为震源位置，r_A和r_B为接收点位置，NS为叠加所用的炮数，r_B为离r_S最近的接收点位置，r_S在r_B和r_A之间；

将所述频域近地表散射波场进行傅里叶反变换，获得对应的时域近地表散射波场。