[发明专利]基于三角滤波的局部倾斜叠加波束形成偏移假频压制方法

摘要

基于三角滤波的局部倾斜叠加波束形成偏移假频压制方法。该三角反假频滤波算子直接在时域实时完成，在对地震数据进行因果和反因果积分后，可以直接与时域倾斜叠加波束形成结合完成，避免了将数据变换至频域进行滤波操作。同时，本发明中提出在反假频公式中加入权重系数，以对反假频程度进行控制。本发明提出的反假频方法是在叠前数据域完成，根据地表速度决定的慢度分量在数据域滤波，反假频后数据可进行压缩存储后重复使用，不受速度模型更新影响。

主权项

1.基于三角滤波的局部倾斜叠加波束形成偏移假频压制方法，其特征在于，包括以下步骤：(1)对给定的三维叠前反射数据u(x_s,y_s,x_r,y_r,t)，对数据划分形成超道集其中，x_s,y_s分别是炮点在联络线和主测线方向上的位置坐标，x_r,y_r分别是检波点在相应方向上的位置坐标，t是时间坐标，分别是每一道数据相对于超道集中心位置的局部相对坐标，分别是超道集炮点中心位置坐标，分别是超道集检波点中心位置坐标；(2)对形成的超道集进行坐标变换，变换后的超道集利用局部中心点局部偏移距坐标表示为其中，Δh_x和Δh_y分别表示局部偏移距，Δm_x和Δm_y分别表示局部中心点坐标；(3)对坐标变换后的超道集进行部分正常时差校正，消除局部偏移距Δh_x和Δh_y的影响，得到部分NMO后的超道集数据(4)对NMO后超道集进行因果和反因果积分后，进行时移叠加即可实现带反假频效果的倾斜叠加波束形成结果u(p_mx,p_my,τ)；p_mx,p_my为中心点射线参数(p_mx,p_my)的x分量和y分量，τ是倾斜叠加的截距时间；具体过程为：(1)对部分正常时差校正后的超道集数据进行因果积分和反因果积分，得到(2)对给定的中心点射线参数(p_mx,p_my)，计算倾斜叠加时移量t_shift＝p_mxΔm_x+p_myΔm_y；(3)对给定的中心点射线参数(p_mx,p_my)，计算反假频滤波算子时移量以及归一化因子α＝(k+1)²，其中，最大频率dt是时间采样间隔，dx和dy分别是联络线和主测线方向上的空间采样间隔，w_x和w_y是引入的反假频权重系数，分别控制联络线和主测线方向上反假频量；(4)对利用偏移倾斜叠加时移量和反假频滤波算子时移量，累加至倾斜叠加结果其中，τ是倾斜叠加的截距时间。