技術(shù)特征：

1.一種提取地震數(shù)據(jù)頻散屬性的處理方法，其特征在于，所述方法包括：

對疊前地震數(shù)據(jù)做時頻變換，獲取所述疊前地震數(shù)據(jù)的時頻振幅譜；

從所述疊前地震數(shù)據(jù)提取地震子波，利用提取出的地震子波制作子波窗函數(shù)；

利用所述子波窗函數(shù)，對所述時頻振幅譜進行加權(quán)的時頻譜處理，得到處理后頻散時頻譜；

對所述處理后頻散時頻譜進行反演，提取出所述疊前地震數(shù)據(jù)的頻散屬性參數(shù)數(shù)據(jù)。

2.如權(quán)利要求1所述的一種提取地震數(shù)據(jù)頻散屬性的處理方法，其特征在于，所述對疊前地震數(shù)據(jù)做時頻變換包括：

通過引入λ_a、p參數(shù)對選取的表示S變換中的窗函數(shù)進行控制，形成對疊前地震數(shù)據(jù)做時頻變換的廣義S變換的表達式，為：

$S(f,\mathrm{\τ})={\∫}_{-\mathrm{\∞}}^{+\mathrm{\∞}}x\left(t\right)\frac{\left|{\mathrm{\λ}}_a\right||f{|}^p}{\sqrt{2\mathrm{\π}}}\exp \[\frac{-{\mathrm{\λ}}_a^2{(\mathrm{\τ}-t)}^2{f}^{2p}}{2}\]\exp (-i2\mathrm{\π}ft)dt$

上式中，x(t)為所述疊前地震數(shù)據(jù)采樣點的時間域振幅表達式，f為頻率，S(f,τ)為x(t)信號廣義S變換結(jié)果，λ_a、p為設(shè)定的常數(shù)，根據(jù)x(t)的指定數(shù)據(jù)特征進行確定；

采用上述廣義S變換對疊前地震數(shù)據(jù)做時頻變換。

3.如權(quán)利要求1所述的一種提取地震數(shù)據(jù)頻散屬性的處理方法，其特征在于，所述從所述疊前地震提取地震子波包括：

獲取所述疊前地震數(shù)據(jù)的瞬時振幅譜，在頻率域求取所述瞬時振幅譜的對數(shù)，得到對數(shù)域的振幅譜：

利用傅立葉反變換將所述對數(shù)域的振幅譜從對數(shù)域變換到復(fù)賽域，得到復(fù)賽域的振幅譜；

在所述復(fù)賽域的振幅譜中確定疊前地震數(shù)據(jù)的地震子波。

4.如權(quán)利要求1所述的一種提取地震數(shù)據(jù)頻散屬性的處理方法，其特征在于，所述對所述時頻振幅譜進行加權(quán)的時頻譜處理包括：

對所述時頻振幅譜中的頻率成分采用如下加權(quán)關(guān)系式進行處理：

B_d(t，n，f)＝S_d(t，n，f)w_e(f，n)

上式中，B_d(t，n，f)為反射系數(shù)譜，即一種處理后頻散時頻譜；S_d(t，n，f)為疊前地震數(shù)據(jù)的時頻振幅譜，w_e(f，n)為利用地震子波制作的反子波譜窗函數(shù)。

5.如權(quán)利要求1所述的一種提取地震數(shù)據(jù)頻散屬性的處理方法，其特征在于，所述對所述處理后頻散時頻譜進行反演包括：

采用下述擴展后Aki-Richards的AVO近似公式進行反演：

$R_{pp}\left(\mathrm{\θ}\right)=\frac{1}{2}(1+{\tan }^2\mathrm{\θ})\frac{{\mathrm{\ΔV}}_p}{V_p}-4\frac{V_s^2}{V_p^2}\frac{{\mathrm{\ΔV}}_s}{V_s}{\sin }^2\mathrm{\θ}+\frac{1}{2}(1-4\frac{V_s^2}{V_p^2}{\sin }^2\mathrm{\θ})\frac{\mathrm{\Δ}\mathrm{\ρ}}{\mathrm{\ρ}}$

上式中，R_pp(θ)為縱波反射系數(shù)；θ為入射角和透射角的平均值，為縱波速度變化率，為橫波速度變化率，為反射界面兩側(cè)密度變化率，為橫波與縱波速度比。

6.如權(quán)利要求5所述的一種提取地震數(shù)據(jù)頻散屬性的處理方法，其特征在于，從所述疊前地震數(shù)據(jù)中提取的頻散屬性參數(shù)數(shù)據(jù)包括：

采用最小二乘法計算獲得的所述疊前地震數(shù)據(jù)中縱橫波速度變化率的導(dǎo)數(shù)，其表達式為：

$I_a=\frac{d}{df}\left(\frac{{\mathrm{\Δv}}_p}{v_p}\right),I_b=\frac{d}{df}\left(\frac{{\mathrm{\Δv}}_s}{v_s}\right)$

上式中，I_a表示為縱波速度隨頻率變化率定義的縱波頻散屬性，I_b表示為橫波速度隨頻率變化率定義的橫波頻散屬性，Δv_p、v_p、Δv_s、v_p分別為上下層介質(zhì)縱波速度差、速度和、橫波速度差、速度和。

7.如權(quán)利要求1至6中任意一項所述的一種提取地震數(shù)據(jù)頻散屬性的處理方法，其特征在于，所述方法還包括：

利用提取出所述疊前地震數(shù)據(jù)的頻散屬性參數(shù)數(shù)據(jù)預(yù)測儲層的含油氣性。

8.一種提取地震數(shù)據(jù)頻散屬性的處理裝置，其特征在于，所述裝置包括：

時頻變換模塊，用于對疊前地震數(shù)據(jù)做時頻變換，獲取所述疊前地震數(shù)據(jù)的時頻振幅譜；

子波窗函數(shù)模塊，用于從所述疊前地震數(shù)據(jù)提取地震子波，利用提取出的地震子波制作子波窗函數(shù)；

時頻處理模塊，用于利用所述子波窗函數(shù)，對所述時頻振幅譜進行加權(quán)的時頻譜處理，得到處理后頻散時頻譜；

頻散屬性提取模塊，用于對所述處理后頻散時頻譜進行反演，提取出所述疊前地震數(shù)據(jù)的頻散屬性參數(shù)數(shù)據(jù)。

9.如權(quán)利要求8所述的一種提取地震數(shù)據(jù)頻散屬性的處理裝置，其特征在于，所述時頻變換模塊對疊前地震數(shù)據(jù)做時頻變換包括：

通過引入λ_a、p參數(shù)對選取的表示S變換中的窗函數(shù)進行控制，形成對疊前地震數(shù)據(jù)做時頻變換的廣義S變換的表達式，為：

$S(f,\mathrm{\τ})={\∫}_{-\mathrm{\∞}}^{+\mathrm{\∞}}x\left(t\right)\frac{\left|{\mathrm{\λ}}_a\right||f{|}^p}{\sqrt{2\mathrm{\π}}}\exp \[\frac{-{\mathrm{\λ}}_a^2{(\mathrm{\τ}-t)}^2{f}^{2p}}{2}\]\exp (-i2\mathrm{\π}ft)dt$

上式中，x(t)為所述疊前地震數(shù)據(jù)采樣點的時間域振幅表達式，f為頻率，S(f,τ)為x(t)信號廣義S變換結(jié)果，λ_a、p為設(shè)定的常數(shù)，根據(jù)x(t)的指定數(shù)據(jù)特征進行確定；

采用上述廣義S變換對疊前地震數(shù)據(jù)做時頻變換。

10.如權(quán)利要求8所述的一種提取地震數(shù)據(jù)頻散屬性的處理裝置，其特征在于，所述時頻處理模塊對所述時頻振幅譜進行加權(quán)的時頻譜處理包括：

對所述時頻振幅譜中的頻率成分采用如下加權(quán)關(guān)系式進行處理：

B_d(t，n，f)＝S_d(t，n，f)w_e(f，n)

上式中，B_d(t，n，f)為反射系數(shù)譜，即一種處理后頻散時頻譜；S_d(t，n，f)為疊前地震數(shù)據(jù)的時頻振幅譜，w_e(f，n)為利用地震子波制作的反子波譜窗函數(shù)。

11.如權(quán)利要求8所述的一種提取地震數(shù)據(jù)頻散屬性的處理裝置，其特征在于，所述頻散屬性提取模塊對所述處理后頻散時頻譜進行反演包括：

采用下述擴展后Aki-Richards的AVO近似公式進行反演：

$R_{pp}\left(\mathrm{\θ}\right)=\frac{1}{2}(1+{\tan }^2\mathrm{\θ})\frac{{\mathrm{\ΔV}}_p}{V_p}-4\frac{V_s^2}{V_p^2}\frac{{\mathrm{\ΔV}}_s}{V_s}{\sin }^2\mathrm{\θ}+\frac{1}{2}(1-4\frac{V_s^2}{V_p^2}{\sin }^2\mathrm{\θ})\frac{\mathrm{\Δ}\mathrm{\ρ}}{\mathrm{\ρ}}$

上式中，R_pp(θ)為縱波反射系數(shù)；θ為入射角和透射角的平均值，為縱波速度變化率，為橫波速度變化率，為反射界面兩側(cè)密度變化率，為橫波與縱波速度比。

12.如權(quán)利要求11所述的一種提取地震數(shù)據(jù)頻散屬性的處理裝置，其特征在于，所述頻散屬性提取模塊從所述疊前地震數(shù)據(jù)中提取的頻散屬性參數(shù)數(shù)據(jù)包括：

采用最小二乘法計算獲得的所述疊前地震數(shù)據(jù)中縱橫波速度變化率的導(dǎo)數(shù)，其表達式為：

$I_a=\frac{d}{df}\left(\frac{{\mathrm{\Δv}}_p}{v_p}\right),I_b=\frac{d}{df}\left(\frac{{\mathrm{\Δv}}_s}{v_s}\right)$

上式中，I_a表示為縱波速度隨頻率變化率定義的縱波頻散屬性，I_b表示為橫波速度隨頻率變化率定義的橫波頻散屬性，Δv_p、v_p、Δv_s、v_p分別為上下層介質(zhì)縱波速度差、速度和、橫波速度差、速度和。

13.如權(quán)利要求8至12中任意一項所述的一種提取地震數(shù)據(jù)頻散屬性的處理裝置，其特征在于，所述裝置還包括：

預(yù)測模塊，用于利用提取出所述疊前地震數(shù)據(jù)的頻散屬性參數(shù)數(shù)據(jù)預(yù)測儲層的含油氣性。

14.一種基于地震頻散AVO屬性的油氣預(yù)測系統(tǒng)，其特征在于，所述系統(tǒng)包括：

I/O接口，用于獲取疊前地震數(shù)據(jù)和輸出計算得到的頻散屬性參數(shù)數(shù)據(jù)；

處理器，用于對疊前地震數(shù)據(jù)做時頻變換，獲取所述疊前地震數(shù)據(jù)的時頻振幅譜；還用于從所述疊前地震數(shù)據(jù)提取地震子波，利用提取出的地震子波制作子波窗函數(shù)；還用于利用所述子波窗函數(shù)，對所述時頻振幅譜進行加權(quán)的時頻譜處理，得到處理后頻散時頻譜；還用于對所述處理后頻散時頻譜進行反演，提取出所述疊前地震數(shù)據(jù)的頻散屬性參數(shù)數(shù)據(jù)。

15.如權(quán)利要求14所述的一種基于地震頻散AVO屬性的油氣預(yù)測系統(tǒng)，其特征在于，所述系統(tǒng)還包括：

預(yù)測處理單元，用于利用提取出所述疊前地震數(shù)據(jù)的頻散屬性參數(shù)數(shù)據(jù)預(yù)測儲層的含油氣性。