[发明专利]一种时间域扇形滤波器及其应用

权利要求书

1.一种时间域扇形滤波器，其特征在于：所述时间域扇形滤波器的时间 频率范围为f₁～f₂f₁为起始时间频率，f₂为截止时间频率，且f₁＜f₂；期望输出为 y(T_n，X_m)，其滤波算子如下：

假设m为空间采样点数，n为时间采样点数，采样间隔为Δt，c为整数，

若m≠n≠0，则

y(Tn,Xm)=12(c2m3-mn2)(cm+n)(cos(2Δtf1(cm-n)π)-cos(2Δtf2(cm-n)π))+(cm-n)(cos(2Δtf1(cm+n)π)-cos(2Δtf2(cm+n)π))---(4)]]>

若c·m＝n≠0，则

y(Tn,Xm)=c(cos(4Δtf1nπ)-cos(4Δtf2nπ))4n2π2---(5)]]>

若m＝n＝0，则

y(Tn,Xm)=2cΔt2f22-2cΔt2f12---(6)]]>

若m＝0，n≠0，则

y(Tn,Xm)=c(cos(2Δtf2nπ)-cos(2Δtf1nπ)+2Δtnπ(f2sin(2Δtf2nπ)-f1sin(2Δtf1nπ)))n2π2---(7).]]>

2.一种利用权利要求1所述时间域扇形滤波器实现的滤波方法，其特征 在于：所述方法包括以下步骤：

(A)输入二维地震数据，利用相对时差校正量将有效信号拉平，形成输入 数据；

(B)选择滤波器算子维度，利用所述公式(4)-(7)得到时间域扇形滤 波器的滤波算子；

(C)将滤波算子作为输入得到褶积算子，然后利用褶积算子将所述时间 域扇形滤波器应用到步骤(A)产生的输入数据上，得到去噪后的地震数据。