[发明专利]多分量地震数据的矢量去噪方法

摘要

本发明公开了一种多分量地震数据的矢量去噪方法，包括：取得多分量地震数据的波矢量，其中多分量地震数据的波矢量包括有效信号、地滚波和随机噪音的波矢量。利用第一移动时窗，对多分量地震数据的波矢量进行计算出第一平均波矢量、并将第一平均波矢量进行中值滤波，以取得地滚波的真实模量，且将多分量地震数据的波矢量与地滚波的波矢量相减，以取得矢量时间序列。利用第二移动时窗，对所述矢量时间序列进行均值滤波，以取得第二平均波矢量，并对所述第二平均波矢量进行中值滤波以及在同一接收线内进行中值滤波，以压制矢量时间序列中的随机噪音的波矢量，进而取得净化后的有效信号的波矢量。

主权项

1.一种多分量地震数据的矢量去噪方法，其特征在于，包括：取得多分量地震数据的波矢量，其中所述多分量地震数据的波矢量包括有效信号、地滚波和随机噪音的波矢量；利用第一移动时窗，对多分量地震数据的波矢量进行计算出第一平均波矢量，并将所述第一平均波矢量进行中值滤波，以取得所述地滚波的真实模量，且将多分量地震数据的波矢量与所述地滚波的波矢量相减，以取得矢量时间序列，所述第一移动时窗的长度的下限为所述有效信号的最大视周期的一半，且所述第一移动时窗的长度小于所述地滚波的最小视周期的一半；利用第二移动时窗，对所述矢量时间序列进行均值滤波，以取得第二平均波矢量，并对所述第二平均波矢量进行中值滤波以及在同一接收线内进行中值滤波，以压制矢量时间序列中的随机噪音的波矢量，进而取得净化后的有效信号的波矢量，所述第二移动时窗的长度设置为地震子波旁瓣的1/4；所述多分量地震数据的波矢量满足如下公式：A(t)＝B(t)+G(t)+N(t)，其中，t表示时间，B(t)、G(t)和N(t)分别表示在t时刻的有效信号、地面滚动和随机噪音的波矢量；所述第一平均波矢量满足如下公式：其中，Kj表示第一平均波矢量，T1表示第一移动时窗，△t表示地震数据的时间采样率，将所述第一平均波矢量进行中值滤波满足如下公式：M(t)＝argmin(D(Kj))，其中，M(t)表示中值矢量，D(Kj)表示所述第一平均波矢量中一个矢量与所有其它矢量的距离的总和且满足如下公式：其中，所述地滚波的真实模量满足如下公式：G(t)＝γ1(t)M(t)，其中，γ1(t)表示第一比例因子；所述第二平均波矢量满足如下公式：其中，M1(t)表示第二平均波矢量，T2表示第二移动时窗，C(t)表示所述矢量时间序列，对所述第二平均波矢量进行中值滤波以及在同一接收线内进行中值滤波满足如下公式：其中，M2(t)表示压制随机噪声矢量后的波矢量，M3(t)表示压制其它非地滚波的相干噪声的波矢量，k为中值滤波的输出道，K为中值滤波的移动道窗长度，所述净化后的有效信号的波矢量满足如下公式：B(t)＝γ2(t)M3(t),其中，γ2(t)表示第二比例因子。