[发明专利]三维高密度地震勘探记录叠前面波压制方法

权利要求书

1.一种三维高密度地震勘探记录叠前面波压制方法，其特征是采用以下技术步骤实现：

1)激发地震波，野外采集三维高密度地震数据，解编后按照炮线和检波线进行分选，得到多个正交子集记录；

2)根据子集数据的炮点间的距离Δx和道间距Δy，根据Nyquist定理得到空间频率的范围kx=0~12*Δx,]]>ky=0~12*Δy,;]]>

3)根据子集数据的时间采样率Δt，根据Nyquist定理得到时间方向的频率范围w=12*Δt;]]>

4)然后根据炮点数量m、检波点数量n和每一道的样点数p，得到子集数据中每一个点的空间波数值和频率值i为该点所在的炮数，j为该点所在的道数，k为该点所在的样点数；

5)通过叠前记录中面波传播的距离X和传播该距离所需时间t的比值得到面波的视速度V_L；

6)从第1个正交子集开始，利用以下傅氏变换的公式进行三维傅氏变换，得到该正交子集记录的三维频率波数谱；

F(w,kx,ky)=∫-∞∞∫-∞∞∫-∞∞X(t,x,y)e-i(wt+kxx+kyy)dtdxdy---(1)]]>

式中F(w，k_x，k_y)是正交子集记录的三维频率波数谱，X(t，x，y)是正交子集地震记录，w是时间方向的频率，k_x是沿炮线方向的空间频率，k_y是沿检波线方向的空间频率；

7)根据上述4)和5)得到的数值，带入以下公式，得到FKK域压制面波的滤波器响应F′(w，k_x，k_y)；

F′(w,kx,ky)=1w/(kx2+ky2)1/2>VL0w/(kx2+ky2)1/2≤VL---(2)]]>

8)利用以下公式，将第6)和7)得到的数据相乘，得到压制面波后的三维频率波数谱；

F(w,kx,ky)=F′(w,kx,ky)∫-∞∞∫-∞∞∫-∞∞X(t,x,y)e-i(wt+kxx+kyy)dtdxdy---(3)]]>

9)将得到的三维频率波数谱利用傅氏反变换的公式进行三维傅氏反变换，得到压制面波后的时间域子集数据；

10)换一个子集数据，重复上述步骤6)～9)，完成整个三维高密度空间采样数据的叠前面波压制。