本發(fā)明涉及天然氣地球物理勘探領(lǐng)域，是一種利用地震數(shù)據(jù)來(lái)提取地震波瞬時(shí)質(zhì)心頻率，進(jìn)而利用地震瞬時(shí)質(zhì)心頻率可靠地指示天然氣的存在。

背景技術(shù)：

在天然氣勘探中，如何從地震信號(hào)的瞬時(shí)屬性中提取更多的有用信息，尤其是利用頻率異常信息，并結(jié)合地質(zhì)、測(cè)井等資料，來(lái)尋找有意義的天然氣儲(chǔ)集帶是石油物探研究人員一直以來(lái)的追求目標(biāo)，同時(shí)也是難點(diǎn)問題。含流體的巖石地層會(huì)造成地震波在傳播過程中發(fā)生能量損失，在含氣層的內(nèi)部及其下部，地震波的能量會(huì)發(fā)生明顯的高頻衰減。譜分解技術(shù)是目前利用地震信號(hào)高頻衰減異常從地震反射數(shù)據(jù)進(jìn)行地質(zhì)解釋及油氣指示的一種常用的有效的含氣性預(yù)測(cè)技術(shù)。地震數(shù)據(jù)體經(jīng)過譜分解可以產(chǎn)生一系列分頻剖面，每個(gè)分頻剖面都是所有地震數(shù)據(jù)體在某個(gè)特定頻率或頻率段處瞬時(shí)振幅的一種反映。不同頻率處的振幅剖面能夠體現(xiàn)不同尺度處的地質(zhì)體的不同的響應(yīng)特征，而含氣性信息可能在某些頻率處的分頻剖面中得到強(qiáng)化反映，更容易被識(shí)別。通常含氣區(qū)域會(huì)體現(xiàn)出“低頻強(qiáng)能量，高頻弱能量”的衰減特征。但是譜分解技術(shù)需要利用一系列地震分頻剖面進(jìn)行分析，然后選取合理的分頻剖面進(jìn)行解釋，工作量較大。

此外，譜分解技術(shù)的核心是利用時(shí)頻分析方法如短時(shí)傅里葉變換、s變換、小波變換等。相比于傳統(tǒng)的地震信號(hào)時(shí)頻分析方法，適用于非線性非平穩(wěn)信號(hào)處理的基于經(jīng)驗(yàn)?zāi)B(tài)分解(empiricalmodedecomposition,emd)的時(shí)頻分析方法具有一些獨(dú)特的優(yōu)勢(shì)，如不受測(cè)不準(zhǔn)原理的限制，時(shí)頻分辨率更高等。目前基于emd的時(shí)頻分析技術(shù)已經(jīng)被廣泛應(yīng)用于天然氣檢測(cè)中，并取得了很好的效果。emd可以把地震信號(hào)從高頻到低頻分解為有限個(gè)單頻或窄帶本征模態(tài)函數(shù)(intrinsicmodefunction,imf)的和。huang等人(1998)認(rèn)為只有本征模態(tài)函數(shù)的瞬時(shí)頻率才有物理意義。雖然基于emd時(shí)頻分析方法的含氣性檢測(cè)技術(shù)較基于傳統(tǒng)時(shí)頻分析方法的含氣性檢測(cè)技術(shù)具有更好的時(shí)空分辨率和精確度，但是由于emd方法本身存在模態(tài)混疊效應(yīng)等問題，會(huì)導(dǎo)致emd分解出的imf缺乏物理意義和地質(zhì)意義，從而降低含氣性檢測(cè)的準(zhǔn)確性。

完備聚合經(jīng)驗(yàn)?zāi)B(tài)分解(completeensembleempiricalmodedecomposition，ceemd)是torres等人(2011)提出的一種克服emd方法中模態(tài)混疊效應(yīng)的改進(jìn)方法，它可以把地震信號(hào)分解的更徹底，誤差在數(shù)值上可以忽略，得到的imf分量更具有物理意義和地質(zhì)意義。

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

本發(fā)明的目的在于解決上述現(xiàn)有技術(shù)中存在的難題，提供一種新的地層含氣性信息的頻率異常檢測(cè)方法，采用的是基于ceemd的地震瞬時(shí)質(zhì)心頻率技術(shù)。通過定義一種新的地震屬性：地震瞬時(shí)質(zhì)心頻率，進(jìn)行天然氣儲(chǔ)層含氣性檢測(cè)，提高現(xiàn)有技術(shù)的精確度。

為實(shí)現(xiàn)上述目的，本發(fā)明采用的技術(shù)方案是：一種利用地震瞬時(shí)質(zhì)心頻率進(jìn)行含氣性檢測(cè)的方法，包括如下步驟：

a、綜合利用地質(zhì)、測(cè)井及合成地震記錄等資料準(zhǔn)確標(biāo)定目標(biāo)層；

b、對(duì)各個(gè)地震道進(jìn)行ceemd分解產(chǎn)生不同的imf分量，對(duì)每個(gè)imf分量，結(jié)合希爾伯特變換計(jì)算瞬時(shí)振幅和瞬時(shí)頻率，再利用地震瞬時(shí)質(zhì)心頻率計(jì)算公式計(jì)算各個(gè)imf分量的瞬時(shí)質(zhì)心頻率，通過相關(guān)加權(quán)系數(shù)對(duì)各個(gè)imf分量計(jì)算的瞬時(shí)質(zhì)心頻率進(jìn)行加權(quán)求和，獲得該地震道的瞬時(shí)質(zhì)心頻率值，逐道計(jì)算目標(biāo)層范圍內(nèi)的地震數(shù)據(jù)的瞬時(shí)質(zhì)心頻率，獲得該工區(qū)目標(biāo)層范圍內(nèi)的地震瞬時(shí)質(zhì)心頻率特征體；

c、結(jié)合地質(zhì)、測(cè)井和試油氣信息相關(guān)資料，利用井旁地震瞬時(shí)質(zhì)心頻率特征剖面確定不同巖性、流體等引起的地震瞬時(shí)質(zhì)心頻率屬性體區(qū)別，再將其特征外推到無(wú)井區(qū)域，從而進(jìn)行巖性、含氣性檢測(cè)。

進(jìn)一步的，步驟b具體為：

b1、將各個(gè)地震道進(jìn)行ceemd分解，得到imf分量

1)確定ceemd分解過程中添加的高斯白噪聲、集合成員的個(gè)數(shù)和最佳最大篩選迭代次數(shù)，添加的白噪聲的幅值a和集合成員的個(gè)數(shù)n應(yīng)滿足：

其中，εn是最終的誤差標(biāo)準(zhǔn)偏差，白噪聲的幅值a取值范圍為0.1～0.4，集合成員的個(gè)數(shù)n＞100；

ceemd方法使用特定的噪聲添加在每次分解的過程中，計(jì)算唯一的剩余量以獲得每個(gè)imf分量，最終分解獲得的imf分量是完全的，具有數(shù)值可忽略的誤差；通常，添加的白噪聲幅度不宜過大，否則會(huì)在分解過程中引入虛假的imf分量；添加的白噪聲分量的幅度也不宜太小，否則它會(huì)導(dǎo)致模態(tài)混疊現(xiàn)象不易消除。

2)ceemd分解：在ceemd分解過程中，使用emd經(jīng)過i次分解信號(hào)x(t)+ε0wi(t)(i＝1,2,…,i)以獲得第一個(gè)模態(tài)，并計(jì)算：

其中，wi(t)(i＝1,2,…,i)是不同的白噪聲序列，x(t)是原始地震信號(hào)；

在第一個(gè)階段(k＝1)中，計(jì)算第一個(gè)余量r1(t)：

分解r1(t)+ε1e1(wi(t))(i＝1,2,…,i),直到產(chǎn)生它們的第一個(gè)imf分量，然后定義第二個(gè)imf分量為：

ej(·)表示產(chǎn)生第j個(gè)模態(tài)，

對(duì)于k＝2,…,k,計(jì)算第k次的余量：

分解rk(t)+εkek(wi(t))(i＝1,2,…,i)直到獲得它們的第一個(gè)emd模態(tài)，然后定義第(k+1)個(gè)模態(tài)為：

重復(fù)前述步驟計(jì)算第k個(gè)模態(tài)；

b2、計(jì)算每個(gè)地震道分解產(chǎn)生的imf分量的瞬時(shí)屬性

對(duì)每一條地震道提取出來(lái)的各個(gè)模態(tài)分量c(t)，分別利用希爾伯特變換用下式提取瞬時(shí)幅度a(t)和瞬時(shí)頻率ω(t)：

其中，h[·]表示hilbert變換，p為柯西主值；為了避免(7)式中計(jì)算瞬時(shí)頻率時(shí)相位解卷繞導(dǎo)致的模糊度，采用下式計(jì)算瞬時(shí)頻率ω(t)：

b3、計(jì)算各條地震道的瞬時(shí)質(zhì)心頻率

瞬時(shí)質(zhì)心頻率f定義為：

其中，ωk為每個(gè)imf分量在每個(gè)時(shí)間采樣點(diǎn)處的瞬時(shí)頻率；ak為每個(gè)imf分量在每個(gè)時(shí)間采樣點(diǎn)處的瞬時(shí)振幅；

本發(fā)明中，瞬時(shí)質(zhì)心頻率f可以更有效的反映頻率異常信息，這里，首先計(jì)算每條地震道各個(gè)imf分量計(jì)算獲得的瞬時(shí)質(zhì)心頻率，然后，采用相關(guān)加權(quán)系數(shù)對(duì)該條地震道計(jì)算的各個(gè)imf分量對(duì)應(yīng)的瞬時(shí)質(zhì)心頻率進(jìn)行加權(quán)求和，結(jié)果作為該地震道的瞬時(shí)質(zhì)心頻率，相關(guān)加權(quán)方案如下：

利用各條地震道分解后生成的imf分量與原始地震道的相關(guān)系數(shù)r的大小對(duì)各個(gè)imf分量獲得的瞬時(shí)質(zhì)心頻率進(jìn)行加權(quán)，相關(guān)加權(quán)系數(shù)rc定義為：

式(10)對(duì)具有強(qiáng)相關(guān)的imf分量計(jì)算得到的瞬時(shí)質(zhì)心頻率保持不變，對(duì)具有較強(qiáng)相關(guān)的imf分量計(jì)算得到的瞬時(shí)質(zhì)心頻率進(jìn)行10^-1的衰減，對(duì)具有弱相關(guān)的imf分量計(jì)算得到的瞬時(shí)質(zhì)心頻率進(jìn)行10^-2的衰減；該操作可以加強(qiáng)該地震道的主要貢獻(xiàn)成分，同時(shí)減弱次要貢獻(xiàn)分量的信息成分。

對(duì)計(jì)算出的最終該條地震道的瞬時(shí)質(zhì)心頻率進(jìn)行歸一化，采用如下歸一化計(jì)算公式：

x2＝(x-min(x))/(max(x)-min(x))……(11)

其中，x為該條地震道的瞬時(shí)質(zhì)心頻率，x2為該條地震道的歸一化瞬時(shí)質(zhì)心頻率，將結(jié)果歸一化到區(qū)間[0,1]范圍內(nèi)，min(·)表示取數(shù)據(jù)的最小值，max(·)表示取數(shù)據(jù)的最大值。

本發(fā)明的有益技術(shù)效果是：本發(fā)明的基于ceemd的地震瞬時(shí)質(zhì)心頻率提取方法，具有如下特點(diǎn)，主要表現(xiàn)為：

(1)使用了克服模態(tài)混疊效應(yīng)的ceemd方法進(jìn)行地震道的分解。通過ceemd方法獲得的各個(gè)imf分量是具有不同頻帶寬度的窄帶信號(hào)，且較常規(guī)emd方法更具物理意義和地質(zhì)意義。準(zhǔn)確的imf分量保證了通過它計(jì)算獲得的瞬時(shí)頻率具有物理意義和更為準(zhǔn)確、明確的地質(zhì)意義。

(2)對(duì)地震道進(jìn)行ceemd分解產(chǎn)生的各個(gè)不同imf分量分別結(jié)合希爾伯特變換提取瞬時(shí)振幅和瞬時(shí)頻率，再計(jì)算瞬時(shí)質(zhì)心頻率，通過相關(guān)加權(quán)系數(shù)對(duì)各個(gè)imf分量計(jì)算的瞬時(shí)質(zhì)心頻率進(jìn)行加權(quán)求和，弱化了地震信號(hào)中的次要貢獻(xiàn)成分，突出了主要貢獻(xiàn)成分，能夠更有效地反映地震信號(hào)的微弱變化，加強(qiáng)了天然氣地震響應(yīng)特征，有利于含氣性檢測(cè)。

(3)地震瞬時(shí)質(zhì)心頻率屬性體的計(jì)算可適用于二維或三維數(shù)據(jù)的剖面分析、沿層切片或者時(shí)間切片等分析，計(jì)算方式靈活多樣。

(4)基于ceemd的地震瞬時(shí)質(zhì)心頻率屬性體估計(jì)方法較常規(guī)譜分解技術(shù)和衰減分析技術(shù)更為簡(jiǎn)單方便，同時(shí)具有較高的分辨率。

附圖說明

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

圖1為含氣地質(zhì)衰減模型；

圖2為模型的地震響應(yīng)；

圖3為模型的基于ceemd的地震瞬時(shí)質(zhì)心頻率剖面；

圖4為四川盆地某地的須家河組一條二維疊后偏移過含氣井剖面；

圖5為利用本技術(shù)計(jì)算的該地震剖面基于ceemd的地震瞬時(shí)質(zhì)心頻率屬性剖面；

圖6為利用基于短時(shí)傅里葉變換的譜分解方法計(jì)算的該地震剖面的分頻剖面的對(duì)比結(jié)果；

圖7為利用基于小波變換的譜分解方法計(jì)算的該地震剖面的分頻剖面的對(duì)比結(jié)果；

圖8為四川盆地某地的須家河組另一條二維疊后偏移過氣水井剖面；

圖9為利用本技術(shù)計(jì)算的該疊后偏移過氣水井剖面的地震瞬時(shí)質(zhì)心頻率屬性剖面。

具體實(shí)施方式

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

實(shí)施例1

一種利用地震瞬時(shí)質(zhì)心頻率進(jìn)行含氣性檢測(cè)的方法，包括如下步驟：

a、綜合利用地質(zhì)、測(cè)井及合成地震記錄等資料準(zhǔn)確標(biāo)定目標(biāo)層；

b、對(duì)各個(gè)地震道進(jìn)行ceemd分解產(chǎn)生不同的imf分量，對(duì)每個(gè)imf分量，結(jié)合希爾伯特變換計(jì)算瞬時(shí)振幅和瞬時(shí)頻率，再利用地震瞬時(shí)質(zhì)心頻率計(jì)算公式計(jì)算各個(gè)imf分量的瞬時(shí)質(zhì)心頻率，通過相關(guān)加權(quán)系數(shù)對(duì)各個(gè)imf分量計(jì)算的瞬時(shí)質(zhì)心頻率進(jìn)行加權(quán)求和，獲得該地震道的瞬時(shí)質(zhì)心頻率值，逐道計(jì)算目標(biāo)層范圍內(nèi)的地震數(shù)據(jù)的瞬時(shí)質(zhì)心頻率，獲得該工區(qū)目標(biāo)層范圍內(nèi)的地震瞬時(shí)質(zhì)心頻率特征體；具體為：

b1、將各個(gè)地震道進(jìn)行ceemd分解，得到imf分量

1)確定ceemd分解過程中添加的高斯白噪聲、集合成員的個(gè)數(shù)和最佳最大篩選迭代次數(shù)，添加的白噪聲的幅值a和集合成員的個(gè)數(shù)n應(yīng)滿足：

其中，εn是最終的誤差標(biāo)準(zhǔn)偏差，白噪聲的幅值a取值范圍為0.1～0.4，集合成員的個(gè)數(shù)n＞100；

ceemd方法使用特定的噪聲添加在每次分解的過程中，計(jì)算唯一的剩余量以獲得每個(gè)imf分量，最終分解獲得的imf分量是完全的，具有數(shù)值可忽略的誤差；通常，添加的白噪聲幅度不宜過大，否則會(huì)在分解過程中引入虛假的imf分量；添加的白噪聲分量的幅度也不宜太小，否則它會(huì)導(dǎo)致模態(tài)混疊現(xiàn)象不易消除。

2)ceemd分解：在ceemd分解過程中，使用emd經(jīng)過i次分解信號(hào)x(t)+ε0wi(t)(i＝1,2,…,i)以獲得第一個(gè)模態(tài)，并計(jì)算：

其中，wi(t)(i＝1,2,…,i)是不同的白噪聲序列，x(t)是原始地震信號(hào)；

在第一個(gè)階段(k＝1)中，計(jì)算第一個(gè)余量r1(t)：

分解r1(t)+ε1e1(wi(t))(i＝1,2,…,i),直到產(chǎn)生它們的第一個(gè)imf分量，然后定義第二個(gè)imf分量為：

ej(·)表示產(chǎn)生第j個(gè)模態(tài)，

對(duì)于k＝2,…,k,計(jì)算第k次的余量：

分解rk(t)+εkek(wi(t))(i＝1,2,…,i)直到獲得它們的第一個(gè)emd模態(tài)，然后定義第(k+1)個(gè)模態(tài)為：

重復(fù)前述步驟計(jì)算第k個(gè)模態(tài)；

b2、計(jì)算每個(gè)地震道分解產(chǎn)生的imf分量的瞬時(shí)屬性

對(duì)每一條地震道提取出來(lái)的各個(gè)模態(tài)分量c(t)，分別利用希爾伯特變換用下式提取瞬時(shí)幅度a(t)和瞬時(shí)頻率ω(t)：

其中，h[·]表示hilbert變換，p為柯西主值；為了避免(7)式中計(jì)算瞬時(shí)頻率時(shí)相位解卷繞導(dǎo)致的模糊度，采用下式計(jì)算瞬時(shí)頻率ω(t)：

b3、計(jì)算各條地震道的瞬時(shí)質(zhì)心頻率

瞬時(shí)質(zhì)心頻率f定義為：

其中，ωk為每個(gè)imf分量在每個(gè)時(shí)間采樣點(diǎn)處的瞬時(shí)頻率；ak為每個(gè)imf分量在每個(gè)時(shí)間采樣點(diǎn)處的瞬時(shí)振幅；

本發(fā)明中，瞬時(shí)質(zhì)心頻率f可以更有效的反映頻率異常信息，這里，首先計(jì)算每條地震道各個(gè)imf分量計(jì)算獲得的瞬時(shí)質(zhì)心頻率，然后，采用相關(guān)加權(quán)系數(shù)對(duì)該條地震道計(jì)算的各個(gè)imf分量對(duì)應(yīng)的瞬時(shí)質(zhì)心頻率進(jìn)行加權(quán)求和，結(jié)果作為該地震道的瞬時(shí)質(zhì)心頻率，相關(guān)加權(quán)方案如下：

利用各條地震道分解后生成的imf分量與原始地震道的相關(guān)系數(shù)r的大小對(duì)各個(gè)imf分量獲得的瞬時(shí)質(zhì)心頻率進(jìn)行加權(quán)，相關(guān)加權(quán)系數(shù)rc定義為：

式(10)對(duì)具有強(qiáng)相關(guān)的imf分量計(jì)算得到的瞬時(shí)質(zhì)心頻率保持不變，對(duì)具有較強(qiáng)相關(guān)的imf分量計(jì)算得到的瞬時(shí)質(zhì)心頻率進(jìn)行10^-1的衰減，對(duì)具有弱相關(guān)的imf分量計(jì)算得到的瞬時(shí)質(zhì)心頻率進(jìn)行10^-2的衰減；該操作可以加強(qiáng)該地震道的主要貢獻(xiàn)成分，同時(shí)減弱次要貢獻(xiàn)分量的信息成分。

對(duì)計(jì)算出的最終該條地震道的瞬時(shí)質(zhì)心頻率進(jìn)行歸一化，采用如下歸一化計(jì)算公式：

x2＝(x-min(x))/(max(x)-min(x))……(11)

其中，x為該條地震道的瞬時(shí)質(zhì)心頻率，x2為該條地震道的歸一化瞬時(shí)質(zhì)心頻率，將結(jié)果歸一化到區(qū)間[0,1]范圍內(nèi)，min(·)表示取數(shù)據(jù)的最小值，max(·)表示取數(shù)據(jù)的最大值。

c、結(jié)合地質(zhì)、測(cè)井和試油氣信息相關(guān)資料，利用井旁地震瞬時(shí)質(zhì)心頻率特征剖面確定不同巖性、流體等引起的地震瞬時(shí)質(zhì)心頻率屬性體區(qū)別，再將其特征外推到無(wú)井區(qū)域，從而進(jìn)行巖性、含氣性檢測(cè)。

應(yīng)用實(shí)施例1

(1)利用地質(zhì)、測(cè)井、速度等資料結(jié)合彌散黏滯方程一個(gè)建立地質(zhì)衰減模型。表1是依據(jù)鄂爾多斯盆地某氣田的資料建立的含氣衰減地質(zhì)模型的參數(shù)。其中，層④為含氣層，層③為干層(地質(zhì)模型如圖1所示)。

表1模型的地質(zhì)參數(shù)

其中，vp是縱波速度，ρ是密度，ζ是彌散系數(shù)，η是黏滯系數(shù)。q是衰減系數(shù)。

模型含氣層厚度為75m，子波頻率為40hz，采樣頻率為512hz。

(2)圖2為模型的地質(zhì)響應(yīng)。

(3)圖3為模型對(duì)應(yīng)的基于ceemd的瞬時(shí)質(zhì)心頻率剖面。如圖所示，本技術(shù)很好的檢測(cè)到了含氣區(qū)域。

(4)圖4為四川盆地某地的須家河組一條二維疊后偏移過井剖面。該區(qū)域?yàn)橹旅苌皫r儲(chǔ)層。圖中上部的橢圓所示區(qū)域?yàn)楹瑲鈪^(qū)域，下部的橢圓所示區(qū)域?yàn)楹畢^(qū)域。含氣井命名為well。

(5)利用本技術(shù)計(jì)算的基于ceemd的地震瞬時(shí)質(zhì)心頻率剖面如圖5所示。從圖中可以看到，在上部分的橢圓所示的含氣區(qū)域中存在強(qiáng)烈的瞬時(shí)質(zhì)心頻率異常值。在下部分的橢圓所示的含水區(qū)域有較強(qiáng)的瞬時(shí)質(zhì)心頻率異常值。排除巖性等其他地層影響因素，所提方法給出了一個(gè)含烴類的統(tǒng)計(jì)性解釋結(jié)果。同時(shí)，所提方法對(duì)于氣、水的衰減程度在檢測(cè)上有所差異，含氣區(qū)域的異常值更大，水層異常值較小。氣層、水層區(qū)分較大。圖像時(shí)間分辨率和空間分辨率較強(qiáng)。

(6)圖6為地震剖面的利用基于短時(shí)傅里葉變換的譜分解方法計(jì)算的分頻剖面的對(duì)比結(jié)果。圖6a和6b所示為利用短時(shí)傅里葉變換分別提取的低頻分頻剖面和高頻分頻剖面。該數(shù)據(jù)主頻為50hz。這里，低頻分頻剖面為45hz，高頻分頻剖面為65hz。從圖中可以看出，基于短時(shí)傅里葉變換的低頻分頻剖面在含氣區(qū)域和含水區(qū)域都檢測(cè)到了較強(qiáng)的強(qiáng)振幅異常，在高頻分頻剖面中在含氣區(qū)域和含水區(qū)域強(qiáng)能量減弱，都表現(xiàn)出了“低頻強(qiáng)能量、高頻弱能量”的特征，排除巖性等其他地層影響因素，都給出了一個(gè)含烴類的統(tǒng)計(jì)解釋結(jié)果。但是分頻剖面時(shí)空分辨率較低，對(duì)氣、水區(qū)域差異不大。

(7)圖7為利用基于小波變換的譜分解方法計(jì)算的地震剖面的分頻剖面的對(duì)比結(jié)果。圖7a和7b所示為利用小波變換分別提取的低頻分頻剖面和高頻分頻剖面。這里，低頻分頻剖面為45hz，高頻分頻剖面為65hz。從圖中可以看出，基于小波變換的低頻分頻剖面和高頻分頻剖面在含氣區(qū)域和含水區(qū)域都呈現(xiàn)出了“低頻強(qiáng)能量、高頻弱能量”的特征，都給出了一個(gè)含烴類的統(tǒng)計(jì)解釋結(jié)果。與圖6基于短時(shí)傅里葉變換提取的分頻剖面相比，基于小波變換的分頻剖面的時(shí)空分辨率較高，但是從圖6a可以看出，氣、水區(qū)域差異不大。對(duì)比圖6、圖7利用常規(guī)譜分解技術(shù)與圖5所提技術(shù)得到的結(jié)果，本技術(shù)給出的結(jié)果時(shí)空分辨率較好，對(duì)儲(chǔ)層信息體現(xiàn)較好，對(duì)氣、水的區(qū)分能力更好。同時(shí)，本技術(shù)不需要類似常規(guī)譜分解技術(shù)需要分析大量的分頻剖面。

(8)圖8為四川盆地某地的須家河組另一條二維疊后偏移過氣水井剖面。圖中的橢圓所示區(qū)域?yàn)楹瑲馑畢^(qū)域。

(9)圖9為利用本技術(shù)計(jì)算的該疊后偏移過氣水井剖面的地震瞬時(shí)質(zhì)心頻率屬性剖面。從圖中可以看到，在含氣水區(qū)域存在較強(qiáng)的瞬時(shí)質(zhì)心頻率異常值。排除巖性等其他地層影響因素，所提方法給出了一個(gè)含烴類的統(tǒng)計(jì)性解釋結(jié)果。

本技術(shù)領(lǐng)域技術(shù)人員可以理解，除非另外定義，這里使用的所有術(shù)語(yǔ)(包括技術(shù)術(shù)語(yǔ)和科學(xué)術(shù)語(yǔ))具有與本發(fā)明所屬領(lǐng)域中的普通技術(shù)人員的一般理解相同的意義。還應(yīng)該理解的是，諸如通用字典中定義的那些術(shù)語(yǔ)應(yīng)該被理解為具有與現(xiàn)有技術(shù)的上下文中的意義一致的意義，并且除非像這里一樣定義，不會(huì)用理想化或過于正式的含義來(lái)解釋。

本領(lǐng)域普通技術(shù)人員可以理解：實(shí)施例中的裝置中的部件可以按照實(shí)施例描述分布于實(shí)施例的裝置中，也可以進(jìn)行相應(yīng)變化位于不同于本實(shí)施例的一個(gè)或多個(gè)裝置中。上述實(shí)施例的部件可以合并為一個(gè)部件，也可以進(jìn)一步拆分成多個(gè)子部件。

最后所應(yīng)說明的是：以上實(shí)施例僅用以說明而非限制本發(fā)明的技術(shù)方案，盡管參照上述實(shí)施例對(duì)本發(fā)明進(jìn)行了詳細(xì)說明，本領(lǐng)域的普通技術(shù)人員應(yīng)該理解：依然可以對(duì)本發(fā)明進(jìn)行修改或者等同替換，而不脫離本發(fā)明的精神和范圍的任何修改或局部替換，其均應(yīng)涵蓋在本發(fā)明的權(quán)利要求范圍當(dāng)中。