本申請(qǐng)屬于地震資料綜合解釋處理領(lǐng)域，尤其涉及一種提取地震數(shù)據(jù)頻散屬性的處理方法、裝置及預(yù)測(cè)系統(tǒng)。

背景技術(shù)：

隨著巖性油氣藏勘探開(kāi)發(fā)的不斷深入，對(duì)于復(fù)雜巖性不僅要求預(yù)測(cè)其分布，而且還通常要求預(yù)測(cè)其含油氣性。

目前主流的流體識(shí)別方法包括利用疊后地震頻率的“吸收衰減”屬性以及多孔彈性介質(zhì)理論指導(dǎo)下的疊前參數(shù)反演預(yù)測(cè)儲(chǔ)層的含油氣性。但是對(duì)于復(fù)雜的儲(chǔ)集體，由于不同流體飽和狀態(tài)下巖石彈性參數(shù)差異很小，流體類(lèi)型識(shí)別困難。同時(shí)，現(xiàn)有常規(guī)方法中依據(jù)含流體填充介質(zhì)常常導(dǎo)致地震波發(fā)生不同程度頻散、衰減的實(shí)驗(yàn)室結(jié)果和理論研究，使得現(xiàn)有的這種方式對(duì)流體的識(shí)別就變得非常的困難。

因此，現(xiàn)有技術(shù)中亟需一種抗噪能力更強(qiáng)、儲(chǔ)層含油氣性預(yù)測(cè)結(jié)果更加可靠和效果更好的方法。

技術(shù)實(shí)現(xiàn)要素：

本發(fā)明目的在于提供一種提取地震數(shù)據(jù)頻散屬性的處理方法、裝置及預(yù)測(cè)系統(tǒng)，可以獲取縱橫波頻散程度參數(shù)，可以用于指示儲(chǔ)層的含油氣性，使得地震波AVO屬性剖面邊界清晰，抗噪能力更強(qiáng)，識(shí)別效果好，提高識(shí)別精度。

本申請(qǐng)?zhí)峁┑囊环N提取地震數(shù)據(jù)頻散屬性的處理方法、裝置及預(yù)測(cè)系統(tǒng)是這樣實(shí)現(xiàn)的：

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

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

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

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

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

所述方法優(yōu)選的實(shí)施例中，所述對(duì)疊前地震數(shù)據(jù)做時(shí)頻變換包括：

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

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

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

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

所述方法優(yōu)選的實(shí)施例中，所述從所述疊前地震提取地震子波包括：

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

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

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

所述方法優(yōu)選的實(shí)施例中，所述對(duì)所述時(shí)頻振幅譜進(jìn)行加權(quán)的時(shí)頻譜處理包括：

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

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

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

所述方法優(yōu)選的實(shí)施例中，所述對(duì)所述處理后頻散時(shí)頻譜進(jìn)行反演包括：

采用下述擴(kuò)展后Aki-Richards的AVO近似公式進(jìn)行反演：

$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è)密度變化率，為橫波與縱波速度比。

所述方法優(yōu)選的實(shí)施例中，從所述疊前地震數(shù)據(jù)中提取的頻散屬性參數(shù)數(shù)據(jù)包括：

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

$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ì)縱波速度差、速度和、橫波速度差、速度和。

所述方法優(yōu)選的實(shí)施例中，所述方法還包括：

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

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

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

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

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

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

所述裝置優(yōu)選的實(shí)施例中，所述時(shí)頻變換模塊對(duì)疊前地震數(shù)據(jù)做時(shí)頻變換包括：

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

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

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

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

所述裝置優(yōu)選的實(shí)施例中，所述時(shí)頻處理模塊對(duì)所述時(shí)頻振幅譜進(jìn)行加權(quán)的時(shí)頻譜處理包括：

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

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

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

所述裝置優(yōu)選的實(shí)施例中，所述頻散屬性提取模塊對(duì)所述處理后頻散時(shí)頻譜進(jìn)行反演包括：

采用下述擴(kuò)展后Aki-Richards的AVO近似公式進(jìn)行反演：

$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è)密度變化率，為橫波與縱波速度比。

所述裝置優(yōu)選的實(shí)施例中，所述頻散屬性提取模塊4從所述疊前地震數(shù)據(jù)中提取的頻散屬性參數(shù)數(shù)據(jù)包括：

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

$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ì)縱波速度差、速度和、橫波速度差、速度和。

所述裝置優(yōu)選的實(shí)施例中，所述裝置還包括：

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

一種基于地震頻散AVO屬性的油氣預(yù)測(cè)系統(tǒng)，所述系統(tǒng)包括：

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

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

所述裝置優(yōu)選的實(shí)施例中，所述系統(tǒng)還包括：

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

本發(fā)明提供的一種提取地震數(shù)據(jù)頻散屬性的處理方法、裝置及預(yù)測(cè)系統(tǒng)，以含流體頻散巖石物理理論為基礎(chǔ)，通過(guò)數(shù)值計(jì)算量化含流體介質(zhì)導(dǎo)致地震波發(fā)生頻散的程度，從而指示儲(chǔ)層的含油氣性。利用本申請(qǐng)所述提供的實(shí)施方式計(jì)算得到的縱橫波頻散AVO屬性剖面邊界清晰，含油氣性橫向變化明顯，整體上縱波頻散屬性?xún)?yōu)于橫波頻散屬性，抗噪能力更強(qiáng)。

附圖說(shuō)明

為了更清楚地說(shuō)明本申請(qǐng)實(shí)施例或現(xiàn)有技術(shù)中的技術(shù)方案，下面將對(duì)實(shí)施例或現(xiàn)有技術(shù)描述中所需要使用的附圖作簡(jiǎn)單地介紹，顯而易見(jiàn)地，下面描述中的附圖僅僅是本申請(qǐng)中記載的一些實(shí)施例，對(duì)于本領(lǐng)域普通技術(shù)人員來(lái)講，在不付出創(chuàng)造性勞動(dòng)性的前提下，還可以根據(jù)這些附圖獲得其他的附圖。

圖1是本發(fā)明提供的一種提取地震數(shù)據(jù)頻散屬性的處理方法一種實(shí)施例的方法流程圖；

圖2是本發(fā)明所述一種提取地震數(shù)據(jù)頻散屬性的處理方法在實(shí)際生成應(yīng)用的一種處理流程示意圖；

圖3是利用本發(fā)明一種實(shí)施例預(yù)測(cè)的儲(chǔ)層含油氣性剖面圖；

圖4是分別利用本發(fā)明一種實(shí)施例和常規(guī)高頻衰減方法預(yù)測(cè)的儲(chǔ)層含油氣性平面圖；

圖5是本發(fā)明提供的一種提取地震數(shù)據(jù)頻散屬性的處理裝置一種實(shí)施例的模塊結(jié)構(gòu)示意圖；

圖6是本發(fā)明提供的一種提取地震數(shù)據(jù)頻散屬性的處理裝置另一種實(shí)施例的模塊結(jié)構(gòu)示意圖。

具體實(shí)施方式

為了使本技術(shù)領(lǐng)域的人員更好地理解本申請(qǐng)中的技術(shù)方案，下面將結(jié)合本申請(qǐng)實(shí)施例中的附圖，對(duì)本申請(qǐng)實(shí)施例中的技術(shù)方案進(jìn)行清楚、完整地描述，顯然，所描述的實(shí)施例僅僅是本申請(qǐng)一部分實(shí)施例，而不是全部的實(shí)施例?；诒旧暾?qǐng)中的實(shí)施例，本領(lǐng)域普通技術(shù)人員在沒(méi)有作出創(chuàng)造性勞動(dòng)前提下所獲得的所有其他實(shí)施例，都應(yīng)當(dāng)屬于本申請(qǐng)保護(hù)的范圍。

圖1是本申請(qǐng)所述一種提取地震數(shù)據(jù)頻散屬性的處理方法一種實(shí)施例的方法流程圖。雖然本申請(qǐng)?zhí)峁┝巳缦率鰧?shí)施例或附圖所示的方法操作步驟或裝置結(jié)構(gòu)，但基于常規(guī)或者無(wú)需創(chuàng)造性的勞動(dòng)在所述方法或裝置中可以包括更多或者更少的操作步驟或模塊單元。在邏輯性上不存在必要因果關(guān)系的步驟或結(jié)構(gòu)中，這些步驟的執(zhí)行順序或裝置的模塊結(jié)構(gòu)不限于本申請(qǐng)實(shí)施例或附圖所示的執(zhí)行順序或模塊結(jié)構(gòu)。所述的方法或模塊結(jié)構(gòu)的在實(shí)際中的裝置或終端產(chǎn)品、系統(tǒng)應(yīng)用時(shí)，可以按照實(shí)施例或者附圖所示的方法或模塊結(jié)構(gòu)進(jìn)行順序執(zhí)行或者并行執(zhí)行(例如并行處理器或者多線(xiàn)程處理的環(huán)境、甚至包括分布式處理的實(shí)施環(huán)境，甚至分布式執(zhí)行環(huán)境)。

本發(fā)明以含流體頻散巖石物理理論為基礎(chǔ)，通過(guò)數(shù)值計(jì)算量化含流體介質(zhì)導(dǎo)致地震波發(fā)生頻散的程度，從而指示儲(chǔ)層的含油氣性，在此可以稱(chēng)之為基于地震頻散AVO屬性預(yù)測(cè)油氣的方法。例如在一個(gè)實(shí)施例中，可以首先把常規(guī)Aki-Richards理論公式拓展到頻率域，并從理論公式中定義了地震波的頻散程度；其次，在優(yōu)化、分選道集的基礎(chǔ)上，利用廣義S變換計(jì)算出不同角度道集地震數(shù)據(jù)的振幅譜，同時(shí)應(yīng)用子波提取技術(shù)從地震資料中求取時(shí)空變瞬時(shí)子波譜，從而對(duì)地震振幅譜作歸一化處理消除子波譜印記；最后，通過(guò)最優(yōu)化反演方法求取出理論公式中定義的縱橫波頻散程度參數(shù)，用以指示儲(chǔ)層的含油氣性。從理論模型計(jì)算的結(jié)果來(lái)看，縱橫波頻散AVO屬性剖面邊界清晰，含油氣性橫向變化明顯，整體上縱波頻散屬性?xún)?yōu)于橫波頻散屬性，抗噪能力更強(qiáng)。下面就具體的一個(gè)實(shí)施例的實(shí)施過(guò)程進(jìn)行描述，以便更好的理解本發(fā)明方案。如圖1所述，本申請(qǐng)?zhí)峁┑囊环N提取地震數(shù)據(jù)頻散屬性的處理方法的一種實(shí)施例可以包括：

S1：對(duì)疊前地震數(shù)據(jù)做時(shí)頻變換，獲取所述疊前地震數(shù)據(jù)的時(shí)頻振幅譜。

一般的，在地震勘探中，地震屬性通常是指將地震數(shù)據(jù)分解成各種屬性，如振幅、波形、頻率、衰減特性、相位、相關(guān)分析、比率等。地震屬性技術(shù)主要是指提取、存儲(chǔ)、檢驗(yàn)、分析、確認(rèn)、評(píng)估地震屬性以及將地震屬性轉(zhuǎn)換為地質(zhì)特征的一套方法。AVO(Amplitude variation with offset，AVO，振幅隨偏移距的變化)技術(shù)用于研究地震反射振幅隨炮點(diǎn)與接收器之間的距離即炮檢距(或入射角)的變化特征來(lái)探討反射系數(shù)響應(yīng)隨炮檢距(或入射角)的變化，進(jìn)而確定反射界面上覆、下伏介質(zhì)的巖性特征及物性參數(shù)。在利用AVO分析處理時(shí)主要是利用CDP道集(common depth point，CDP，共深度點(diǎn)道集)、CMP道集(common middle point，CMP，共中心點(diǎn)道集)或者角道集資料等，分析反射波振幅隨炮檢距(也即入射角)的變化規(guī)律，估算界面的彈性參數(shù)泊松比，進(jìn)一步推斷地層的巖性和含油氣情況。借助AVO分析，地球物理學(xué)家可以更好地評(píng)估油氣藏巖石屬性，包括孔隙度、密度、巖性與流體含量。

在本申請(qǐng)實(shí)施例中，可以獲取疊前地震數(shù)據(jù)，然后對(duì)疊前地震數(shù)據(jù)做時(shí)頻變換。通過(guò)時(shí)頻分析方法技術(shù)對(duì)地震信號(hào)進(jìn)行分析處理可在獲得地震信號(hào)的瞬時(shí)頻率、瞬時(shí)相位、瞬時(shí)振幅等瞬時(shí)參數(shù)的同時(shí)獲得時(shí)頻譜圖等重要時(shí)頻域信息，實(shí)現(xiàn)對(duì)地震信號(hào)的邊緣檢測(cè)、屬性提取等。在地震數(shù)處理中，常常會(huì)用到S變換，一般的，所述S變換是在小波變換和短時(shí)傅里葉變換發(fā)展起來(lái)的，S變換采用高斯窗函數(shù)且窗寬與頻率的倒數(shù)成正比，免去了窗函數(shù)的選擇和改善了窗寬固定的缺陷，并且時(shí)頻表示中各頻率分量的相位譜與原始信號(hào)保持直接的聯(lián)系。本發(fā)明實(shí)施例中構(gòu)造了一種新廣義S變換方法，它不僅滿(mǎn)足地震信號(hào)低頻段應(yīng)具有很高的頻率分辨率，而且在高頻部分具有很高的時(shí)間分辨率，而且運(yùn)算效率高。在具體實(shí)施時(shí)，例如在優(yōu)化、分選道集的基礎(chǔ)上，可以利用廣義S變換計(jì)算出不同角度道集地震數(shù)據(jù)的振幅譜。

一般的，假設(shè)給定平方可積的信號(hào)x(t)，即x(t)∈L^2(R)，則x(t)的標(biāo)準(zhǔn)S變換可以為：

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

上述中的x(t)為所述疊前AVO地震數(shù)據(jù)采樣點(diǎn)的時(shí)間域振幅表達(dá)式。

而在本發(fā)明的一種實(shí)施方式中定義了一種廣義S變換，通過(guò)引入兩個(gè)參數(shù)對(duì)標(biāo)準(zhǔn)S變換中的窗函數(shù)進(jìn)行控制，其廣義S變換的表達(dá)式可以為：

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

上述中，λ_α、p均為設(shè)定的常數(shù)，可以根據(jù)待分析信號(hào)x(t)的指定數(shù)據(jù)特征進(jìn)行確定，具體可以由作業(yè)人員設(shè)置。

因此，本發(fā)明所述方法的另一種實(shí)施例中，所述對(duì)疊前地震數(shù)據(jù)做時(shí)頻變換可以包括：

S101：通過(guò)引入λ_α、p參數(shù)對(duì)選取的表示S變換中的窗函數(shù)進(jìn)行控制，形成對(duì)疊前地震數(shù)據(jù)做時(shí)頻變換的廣義S變換的表達(dá)式，為：

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

上式中，x(t)為所述疊前地震數(shù)據(jù)采樣點(diǎn)的時(shí)間域振幅表達(dá)式，f為頻率，S(f,τ)為x(t)信號(hào)廣義S變換結(jié)果，λ_α、p為設(shè)定的常數(shù)，根據(jù)x(t)的指定數(shù)據(jù)特征進(jìn)行確定。一般的，λ_α、p取值范圍為正實(shí)數(shù)，例如在分析實(shí)際地震勘探數(shù)據(jù)時(shí)的取值可以分別在2.3、0.9附近。

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

本發(fā)明實(shí)施例中，可以獲取勘探區(qū)域的疊前地震數(shù)據(jù)，然后對(duì)所述疊前地震數(shù)據(jù)做時(shí)頻變換，獲取得到所述疊前地震數(shù)據(jù)的振幅譜。當(dāng)然，地震數(shù)據(jù)時(shí)頻變換的方式有很多，本發(fā)明提供了一種實(shí)施例中可以在優(yōu)化、分選道集的基礎(chǔ)上，利用變換定義的廣義S變換計(jì)算出不同角度道集地震數(shù)據(jù)的振幅譜。在對(duì)疊前地震數(shù)據(jù)進(jìn)行廣義S變換之前可以進(jìn)行道集優(yōu)化、道集分選的處理。

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

本實(shí)施例中可以從所述疊前地震數(shù)據(jù)中提取地震子波，然后利用提取出的地震子波制作子波窗函數(shù)。

地震子波的提取可以采用一些瞬時(shí)子波提取技術(shù)。本發(fā)明提供的一種實(shí)施例中，所述從所述疊前地震提取地震子波可以包括：

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

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

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

通常，在復(fù)賽域的振幅譜上，疊前地震數(shù)據(jù)的地震子波和反射系數(shù)分別位于左右兩端，這樣，本實(shí)施例可以經(jīng)過(guò)上述地震數(shù)據(jù)變換處理提取出地震子波。然后可以利用所述提取的地震子波制作子波窗函數(shù)，這樣我們就可以設(shè)計(jì)一個(gè)濾波器將地震子波的振幅譜分離出來(lái)。

地震振幅譜等于地震子波譜與反射系數(shù)譜乘積，因此分離出地震子波譜后剩下的反射系數(shù)譜即可利用頻散數(shù)學(xué)計(jì)算公式計(jì)算縱橫波頻散屬性。本發(fā)明實(shí)施例中，所述的提取出的地震子波可以用于制作子波窗函數(shù)，可以用于后續(xù)的時(shí)頻譜處理，以消除子波印記。

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

通過(guò)廣義S變換得到的時(shí)頻譜結(jié)果中，反射系數(shù)原有的能量分布特征被地震子波改變，即主頻位置能量最大兩側(cè)隨頻率減小。因此，本發(fā)明的一種實(shí)施例中可以采用加權(quán)的時(shí)頻譜處理方法來(lái)消除“子波疊印”。具體的一種實(shí)施方式中，所述對(duì)所述時(shí)頻振幅譜進(jìn)行加權(quán)的時(shí)頻譜處理可以包括：

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

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

上式中，B_d(t，n，f)為反射系數(shù)譜，即一種處理后頻散時(shí)頻譜；S_d(t，n，f)為疊前地震數(shù)據(jù)的時(shí)頻振幅譜，w_e(f，n)為利用地震子波制作的反子波譜窗函數(shù)。這里的We(f,n)反子波譜窗函數(shù)可以由反子波譜得到，反子波譜是由復(fù)賽譜變換得到。本實(shí)施例中，地震數(shù)據(jù)的時(shí)頻振幅譜與反子波譜窗函數(shù)相乘就可以得到反射系數(shù)譜，窗函數(shù)可以通過(guò)反子波譜進(jìn)一步計(jì)算得到的。

上述加權(quán)的時(shí)頻譜處理后頻譜能量隨著頻率的呈單調(diào)變化，因此本發(fā)明可以夠利用反射系數(shù)譜計(jì)算頻散屬性。

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

地震波在地層中傳播時(shí)，地震波傳播速度發(fā)生頻散現(xiàn)象，即縱橫波速度變化率隨頻率的變化引起反射系數(shù)隨頻率變化。本發(fā)明中可以選取一種最優(yōu)化反演算法對(duì)所述處理后頻散時(shí)頻譜進(jìn)行反演，求取AVO屬性，可以包括采用非線(xiàn)性數(shù)據(jù)擬合的近似求解方式等。圖2是本發(fā)明所述一種提取地震數(shù)據(jù)頻散屬性的處理方法在實(shí)際生成應(yīng)用的一種處理流程示意圖。

本發(fā)明的一種實(shí)施例提供對(duì)Aki-Richards提出的公式進(jìn)行參數(shù)擴(kuò)展，建立反射系數(shù)(可以表示速率)與頻率的關(guān)系，計(jì)算縱橫波速度隨頻率的變化。Aki-Richards提出的公式為一種最小二乘法，本發(fā)明的一種實(shí)施例對(duì)其進(jìn)行擴(kuò)展和改進(jìn)。因此，本申請(qǐng)的一種實(shí)施例中，所述對(duì)所述處理后頻散時(shí)頻譜進(jìn)行反演可以包括：

S401：采用下述擴(kuò)展后Aki-Richards的AVO近似公式進(jìn)行反演：

$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è)密度變化率，為橫波與縱波速度比。

Aki-Richards近似式：

$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{\ρ}}$

對(duì)Aki-Richards提出的公式進(jìn)行參數(shù)擴(kuò)展，其公式為：

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

擴(kuò)展后Aki-Richards的AVO近似公式反射系數(shù)與入射角和頻率有關(guān)，為計(jì)算地震波速度頻散程度提供了計(jì)算方法。

另外，利用上述擴(kuò)展后Aki-Richards的AVO近似公式得到的為反應(yīng)地震波速率隨頻率變化的量化參數(shù)。具體的一種實(shí)施例中，從所述疊前地震數(shù)據(jù)中提取的頻散屬性參數(shù)數(shù)據(jù)可以包括：

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

$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ì)縱波速度差、速度和、橫波速度差、速度和。

圖2是本發(fā)明提供的一種提取地震數(shù)據(jù)頻散屬性一種實(shí)施過(guò)程的流程示意圖。由圖中可以看出，本發(fā)明的技術(shù)方案可以將儲(chǔ)層的含油氣性通過(guò)流體填充介質(zhì)的頻散屬性進(jìn)行量化體現(xiàn)，可以有效指示儲(chǔ)層油氣，預(yù)測(cè)油氣含量，具有十分明顯的實(shí)際生產(chǎn)應(yīng)用。因此，本發(fā)明所述的一種地震頻散屬性提取方法還可以包括：

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

圖3是利用本發(fā)明一種實(shí)施例預(yù)測(cè)的儲(chǔ)層含油氣性剖面圖，圖4是分別利用本發(fā)明一種實(shí)施例和常規(guī)高頻衰減方法預(yù)測(cè)的儲(chǔ)層含油氣性平面圖，圖4中4-1是采用本發(fā)明方法獲得的結(jié)果，4-2是采用一種高頻衰減方法獲得的結(jié)果。由圖4中的4-1、4-2對(duì)比可以看出，利用本發(fā)明提供的方法使得縱橫波頻散AVO屬性剖面邊界清晰，含油氣性橫向變化明顯，整體上縱波頻散屬性?xún)?yōu)于橫波頻散屬性，抗噪能力更強(qiáng)。從本區(qū)鉆井來(lái)看，河道的南部多為出水井，高頻吸收屬性并不能很好的反應(yīng)這一特征，而頻散AVO屬性對(duì)油氣預(yù)測(cè)更為敏感，預(yù)測(cè)的油氣區(qū)更為合理。

基于本發(fā)明上述提供的提取地震數(shù)據(jù)頻散屬性的處理方法，本發(fā)明還提供一種提取地震數(shù)據(jù)頻散屬性的處理裝置(相應(yīng)的裝置項(xiàng)描述)。圖5是本發(fā)明提供的一種提取地震數(shù)據(jù)頻散屬性的處理裝置一種實(shí)施例的模塊結(jié)構(gòu)示意圖，如圖5所示，所述裝置可以包括：

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

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

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

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

本發(fā)明提供的一種提取地震數(shù)據(jù)頻散屬性的處理裝置的一種實(shí)施例中，所述時(shí)頻變換模塊101對(duì)疊前地震數(shù)據(jù)做時(shí)頻變換可以包括：

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

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

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

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

當(dāng)然，地震數(shù)據(jù)時(shí)頻變換的方式有很多，本發(fā)明提供了一種實(shí)施例中可以在優(yōu)化、分選道集的基礎(chǔ)上，利用變換定義的廣義S變換計(jì)算出不同角度道集地震數(shù)據(jù)的振幅譜。在對(duì)疊前地震數(shù)據(jù)進(jìn)行廣義S變換之前可以進(jìn)行道集優(yōu)化、道集分選的處理。

地震子波的提取可以采用一些瞬時(shí)子波提取技術(shù)。通常，在復(fù)賽域的振幅譜上，疊前地震數(shù)據(jù)的地震子波和反射系數(shù)分別位于左右兩端，這樣，本實(shí)施例可以經(jīng)過(guò)上述地震數(shù)據(jù)變換處理提取出地震子波。然后可以利用所述提取的地震子波制作子波窗函數(shù)，這樣我們就可以設(shè)計(jì)一個(gè)濾波器將地震子波的振幅譜分離出來(lái)。地震振幅譜等于地震子波譜與反射系數(shù)譜乘積，因此分離出地震子波譜后剩下的反射系數(shù)譜即可利用頻散數(shù)學(xué)計(jì)算公式計(jì)算縱橫波頻散屬性。本發(fā)明實(shí)施例中，所述的提取出的地震子波可以用于制作子波窗函數(shù)，可以用于后續(xù)的時(shí)頻譜處理，以消除子波印記。

本發(fā)明提供的一種提取地震數(shù)據(jù)頻散屬性的處理裝置的一種實(shí)施例中，所述時(shí)頻處理模塊103對(duì)所述時(shí)頻振幅譜進(jìn)行加權(quán)的時(shí)頻譜處理可以包括：

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

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

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

這里的We(f,n)反子波譜窗函數(shù)可以由反子波譜得到，反子波譜是由復(fù)賽譜變換得到。本實(shí)施例中，地震數(shù)據(jù)的時(shí)頻振幅譜與反子波譜窗函數(shù)相乘就可以得到反射系數(shù)譜，窗函數(shù)可以通過(guò)反子波譜進(jìn)一步計(jì)算得到的。

本發(fā)明的一種實(shí)施例提供對(duì)Aki-Richards提出的公式進(jìn)行參數(shù)擴(kuò)展，建立反射系數(shù)(可以表示速率)與頻率的關(guān)系，計(jì)算縱橫波速度隨頻率的變化。Aki-Richards提出的公式為一種最小二乘法，本發(fā)明的一種實(shí)施例對(duì)其進(jìn)行擴(kuò)展和改進(jìn)。本發(fā)明提供的一種提取地震數(shù)據(jù)頻散屬性的處理裝置的一種實(shí)施例中，所述頻散屬性提取模塊104對(duì)所述處理后頻散時(shí)頻譜進(jìn)行反演可以包括：

采用下述擴(kuò)展后Aki-Richards的AVO近似公式進(jìn)行反演：

$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è)密度變化率，為橫波與縱波速度比。

擴(kuò)展后Aki-Richards的AVO近似公式反射系數(shù)與入射角和頻率有關(guān)，為計(jì)算地震波速度頻散程度提供了計(jì)算方法。

本發(fā)明提供的一種提取地震數(shù)據(jù)頻散屬性的處理裝置的一種實(shí)施例中，所述頻散屬性提取模塊104從所述疊前地震數(shù)據(jù)中提取的頻散屬性參數(shù)數(shù)據(jù)可以包括：

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

$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ì)縱波速度差、速度和、橫波速度差、速度和。。

本發(fā)明的技術(shù)方案可以將儲(chǔ)層的含油氣性通過(guò)流體填充介質(zhì)的頻散屬性進(jìn)行量化體現(xiàn)，可以有效指示儲(chǔ)層油氣，預(yù)測(cè)油氣含量，具有十分明顯的實(shí)際生產(chǎn)應(yīng)用。因此，本發(fā)明提供的一種提取地震數(shù)據(jù)頻散屬性的處理裝置的一種實(shí)施例中，所述裝置還可以包括：

預(yù)測(cè)模塊105，可以用于利用提取出所述疊前地震數(shù)據(jù)的頻散屬性參數(shù)數(shù)據(jù)預(yù)測(cè)儲(chǔ)層的含油氣性。圖6是本發(fā)明提供的一種提取地震數(shù)據(jù)頻散屬性的處理裝置另一種實(shí)施例的模塊結(jié)構(gòu)示意圖。

利用本發(fā)明提供的裝置，可以獲取縱橫波頻散程度參數(shù)，可以用于指示儲(chǔ)層的含油氣性，使得地震波AVO屬性剖面邊界清晰，抗噪能力更強(qiáng)，識(shí)別效果好，提高識(shí)別精度。

上述所述的方法或裝置可以用于地震數(shù)據(jù)處理系統(tǒng)中，實(shí)現(xiàn)地震數(shù)據(jù)中頻散屬性的提取，可以用于指示儲(chǔ)層的含油氣性，使得地震波AVO屬性剖面邊界清晰，抗噪能力更強(qiáng)，提高系統(tǒng)識(shí)別精度和識(shí)別效果。因此，本發(fā)明還提供一種基于地震頻散AVO屬性的油氣預(yù)測(cè)系統(tǒng)，一種實(shí)施例中，所述系統(tǒng)可以包括：

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

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

當(dāng)然，本發(fā)明提供的所述基于地震頻散AVO屬性的油氣預(yù)測(cè)系統(tǒng)的其他實(shí)施例中，如具體的對(duì)疊前地震數(shù)據(jù)做時(shí)頻變換、從所述疊前地震提取地震子波、對(duì)所述時(shí)頻振幅譜進(jìn)行加權(quán)的時(shí)頻譜處理、對(duì)所述處理后頻散時(shí)頻譜進(jìn)行反演以及所述疊前地震數(shù)據(jù)中提取的頻散屬性參數(shù)數(shù)據(jù)等的實(shí)施方式可以參照上述方法或裝置的相關(guān)描述，在此不做贅述。

進(jìn)一步的，所述系統(tǒng)得到疊前地震數(shù)據(jù)的頻散屬性參數(shù)數(shù)據(jù)后，如Ia、Ib等，可以將這些頻散屬性參數(shù)數(shù)據(jù)輸出到其他地震數(shù)據(jù)處理系統(tǒng)，或者所述系統(tǒng)的其他功能模塊中，進(jìn)一步利用這些頻散屬性參數(shù)數(shù)據(jù)或者結(jié)合其他參數(shù)預(yù)測(cè)、分析儲(chǔ)層的含油氣性。因此，本發(fā)明的所述系統(tǒng)的另一種實(shí)施例中，所述系統(tǒng)還可以包括：

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

上述實(shí)施例提供的基于地震頻散AVO屬性的油氣預(yù)測(cè)系統(tǒng)，可以獲取縱橫波頻散程度參數(shù)，可以用于指示儲(chǔ)層的含油氣性，使得地震波AVO屬性剖面邊界清晰，抗噪能力更強(qiáng)，提高系統(tǒng)儲(chǔ)層含油氣性的識(shí)別精度和識(shí)別效果。

盡管本申請(qǐng)內(nèi)容中提到添加了λ_a、p參數(shù)的廣義S變換、擴(kuò)展后ki-Richards的AVO近似公式、地震子波提取、窗函數(shù)構(gòu)建等的地震數(shù)據(jù)處理、不同域變換、常規(guī)方法變形的處理/交互/判斷方式等的描述，但是，本申請(qǐng)并不局限于必須是符合標(biāo)準(zhǔn)定義或?qū)嵤├枋龅那闆rS變換、公式擴(kuò)展、數(shù)據(jù)處理等，某些行業(yè)標(biāo)準(zhǔn)、通用數(shù)據(jù)處理方法或者使用自定義方式或?qū)嵤├枋龅膶?shí)施基礎(chǔ)上略加修改后的實(shí)施方案也可以實(shí)現(xiàn)上述實(shí)施例相同、等同或相近、或變形后可預(yù)料的實(shí)施效果。應(yīng)用這些修改或變形后的數(shù)據(jù)獲取、存儲(chǔ)、判斷、處理方式等獲取的實(shí)施例，仍然可以屬于本申請(qǐng)的可選實(shí)施方案范圍之內(nèi)。

雖然本申請(qǐng)?zhí)峁┝巳鐚?shí)施例或流程圖所述的方法操作步驟，但基于常規(guī)或者無(wú)創(chuàng)造性的手段可以包括更多或者更少的操作步驟。實(shí)施例中列舉的步驟順序僅僅為眾多步驟執(zhí)行順序中的一種方式，不代表唯一的執(zhí)行順序。在實(shí)際中的裝置或客戶(hù)端產(chǎn)品執(zhí)行時(shí)，可以按照實(shí)施例或者附圖所示的方法順序執(zhí)行或者并行執(zhí)行(例如并行處理器或者多線(xiàn)程處理的環(huán)境，甚至為分布式數(shù)據(jù)處理環(huán)境)。術(shù)語(yǔ)“包括”、“包含”或者其任何其他變體意在涵蓋非排他性的包含，從而使得包括一系列要素的過(guò)程、方法、產(chǎn)品或者設(shè)備不僅包括那些要素，而且還包括沒(méi)有明確列出的其他要素，或者是還包括為這種過(guò)程、方法、產(chǎn)品或者設(shè)備所固有的要素。在沒(méi)有更多限制的情況下，并不排除在包括所述要素的過(guò)程、方法、產(chǎn)品或者設(shè)備中還存在另外的相同或等同要素。

上述實(shí)施例闡明的單元、裝置或模塊等，具體可以由計(jì)算機(jī)芯片或?qū)嶓w實(shí)現(xiàn)，或者由具有某種功能的產(chǎn)品來(lái)實(shí)現(xiàn)。為了描述的方便，描述以上裝置時(shí)以功能分為各種模塊分別描述。當(dāng)然，在實(shí)施本申請(qǐng)時(shí)可以把各模塊的功能在同一個(gè)或多個(gè)軟件和/或硬件中實(shí)現(xiàn)，也可以將實(shí)現(xiàn)同一功能的模塊由多個(gè)子模塊或子單元的組合實(shí)現(xiàn)等。

本領(lǐng)域技術(shù)人員也知道，除了以純計(jì)算機(jī)可讀程序代碼方式實(shí)現(xiàn)控制器以外，完全可以通過(guò)將方法步驟進(jìn)行邏輯編程來(lái)使得控制器以邏輯門(mén)、開(kāi)關(guān)、專(zhuān)用集成電路、可編程邏輯控制器和嵌入微控制器等的形式來(lái)實(shí)現(xiàn)相同功能。因此這種控制器可以被認(rèn)為是一種硬件部件，而對(duì)其內(nèi)部包括的用于實(shí)現(xiàn)各種功能的裝置也可以視為硬件部件內(nèi)的結(jié)構(gòu)?；蛘呱踔?，可以將用于實(shí)現(xiàn)各種功能的裝置視為既可以是實(shí)現(xiàn)方法的軟件模塊又可以是硬件部件內(nèi)的結(jié)構(gòu)。

本申請(qǐng)可以在由計(jì)算機(jī)執(zhí)行的計(jì)算機(jī)可執(zhí)行指令的一般上下文中描述，例如程序模塊。一般地，程序模塊包括執(zhí)行特定任務(wù)或?qū)崿F(xiàn)特定抽象數(shù)據(jù)類(lèi)型的例程、程序、對(duì)象、組件、數(shù)據(jù)結(jié)構(gòu)、類(lèi)等等。也可以在分布式計(jì)算環(huán)境中實(shí)踐本申請(qǐng)，在這些分布式計(jì)算環(huán)境中，由通過(guò)通信網(wǎng)絡(luò)而被連接的遠(yuǎn)程處理設(shè)備來(lái)執(zhí)行任務(wù)。在分布式計(jì)算環(huán)境中，程序模塊可以位于包括存儲(chǔ)設(shè)備在內(nèi)的本地和遠(yuǎn)程計(jì)算機(jī)存儲(chǔ)介質(zhì)中。

通過(guò)以上的實(shí)施方式的描述可知，本領(lǐng)域的技術(shù)人員可以清楚地了解到本申請(qǐng)可借助軟件加必需的通用硬件平臺(tái)的方式來(lái)實(shí)現(xiàn)?；谶@樣的理解，本申請(qǐng)的技術(shù)方案本質(zhì)上或者說(shuō)對(duì)現(xiàn)有技術(shù)做出貢獻(xiàn)的部分可以以軟件產(chǎn)品的形式體現(xiàn)出來(lái)，該計(jì)算機(jī)軟件產(chǎn)品可以存儲(chǔ)在存儲(chǔ)介質(zhì)中，如ROM/RAM、磁碟、光盤(pán)等，包括若干指令用以使得一臺(tái)計(jì)算機(jī)設(shè)備(可以是個(gè)人計(jì)算機(jī)，移動(dòng)終端，服務(wù)器，或者網(wǎng)絡(luò)設(shè)備等)執(zhí)行本申請(qǐng)各個(gè)實(shí)施例或者實(shí)施例的某些部分所述的方法。

本說(shuō)明書(shū)中的各個(gè)實(shí)施例采用遞進(jìn)的方式描述，各個(gè)實(shí)施例之間相同或相似的部分互相參見(jiàn)即可，每個(gè)實(shí)施例重點(diǎn)說(shuō)明的都是與其他實(shí)施例的不同之處。本申請(qǐng)可用于眾多通用或?qū)Ｓ玫挠?jì)算機(jī)系統(tǒng)環(huán)境或配置中。例如：個(gè)人計(jì)算機(jī)、服務(wù)器計(jì)算機(jī)、手持設(shè)備或便攜式設(shè)備、平板型設(shè)備、多處理器系統(tǒng)、基于微處理器的系統(tǒng)、置頂盒、可編程的電子設(shè)備、網(wǎng)絡(luò)PC、小型計(jì)算機(jī)、大型計(jì)算機(jī)、包括以上任何系統(tǒng)或設(shè)備的分布式計(jì)算環(huán)境等等。

雖然通過(guò)實(shí)施例描繪了本申請(qǐng)，本領(lǐng)域普通技術(shù)人員知道，本申請(qǐng)有許多變形和變化而不脫離本申請(qǐng)的精神，希望所附的權(quán)利要求包括這些變形和變化而不脫離本申請(qǐng)的精神。