[发明专利]多分量地震数据的t-f-k域极化滤波方法

摘要

本发明公开一种多分量地震数据的t‑f‑k域极化滤波方法，包括：输入多分量地震炮集数据，取得多分量地震炮集数据的垂直与径向分量；以垂直分量为虚部、径向分量为实部，产生复矢量炮集地震数据；给定时窗长度、椭圆极化率滤波门限值与频率门限值；对复矢量炮集地震数据进行t‑f‑k变换，得到三维时变频波谱；根据三维时变频波谱的正频率与负频率，取得椭圆极化率；根据椭圆极化率、椭圆极化率滤波门限值与频率门限值，取得滤波函数；对滤波函数进行三维中值滤波；以滤波函数，对三维时变频波谱进行滤波；将滤波后的三维时变频波谱进行逆t‑f‑k变换，取逆变换的实部为滤波后的径向分量，虚部为滤波后的垂直分量。通过本发明，以解决现有技术存在的由于横波与瑞雷面波的频率和视速度差异较小而难以有效分离瑞雷面波与横波的问题。

主权项

1.一种多分量地震数据的t‑f‑k域极化滤波方法，其特征在于，包括：输入多分量地震炮集数据，以取得所述多分量地震炮集数据的垂直分量与径向分量；以所述垂直分量为虚部、所述径向分量为实部，产生复矢量炮集地震数据；给定时窗长度、椭圆极化率滤波门限值与频率门限值；对所述复矢量炮集地震数据进行t‑f‑k变换，以得到三维时变频波谱；根据所述三维时变频波谱的正频率与负频率，取得椭圆极化率；根据所述椭圆极化率、所述椭圆极化率滤波门限值与所述频率门限值，取得滤波函数；对滤波函数进行三维中值滤波；以所述滤波函数，对三维时变频波谱进行滤波；将滤波后的三维时变频波谱进行逆t‑f‑k变换，取逆变换的实部为滤波后的径向分量，虚部为滤波后的垂直分量；其中，所述复矢量炮集地震数据满足如下公式：s(τ,x)＝sr(τ,x)+j*sz(τ,x)，其中，x为偏移距参数，τ为时间参数，j为虚数单位，sr(τ,x)为径向分量，sz(τ,x)为垂直分量；其中，所述三维时变频波谱满足如下公式：其中，g(τ)为时间窗函数，f表示频率，k表示波数，t表示时间；其中，所述椭圆极化率满足如下公式：其中，椭圆极化率ρ(t,f,k)是关于原点对称的。