一種利用地震能量差異屬性預(yù)測(cè)儲(chǔ)層的方法
【技術(shù)領(lǐng)域】
[0001]本發(fā)明涉及地球物理勘探方法，是一種利用地震能量差異屬性預(yù)測(cè)儲(chǔ)層的方法。
【背景技術(shù)】
[0002]隨著三維地震勘探技術(shù)的不斷進(jìn)步和三維地震資料在油田勘探開(kāi)發(fā)領(lǐng)域應(yīng)用的不斷深入，地震屬性解釋的重要性日益提高，地震屬性主要有時(shí)間、振幅、頻率、相位、相干和衰減幾大類。在油田勘探開(kāi)發(fā)領(lǐng)域，地震屬性主要應(yīng)用于油藏儲(chǔ)層厚度、孔隙度、滲透率等參數(shù)的預(yù)測(cè)，地震資料的振幅、頻率等單項(xiàng)地震屬性的應(yīng)用比較廣泛，其中，利用振幅屬性預(yù)測(cè)儲(chǔ)層是一種最常用的方法，因?yàn)槠洳僮骱?jiǎn)單而受到地震解釋人員的青睞。但常規(guī)振幅類屬性仍然存在一些問(wèn)題，一方面，振幅類的屬性普遍分辨率較低，且沒(méi)有考慮相位的因素；另一方面，常規(guī)的振幅類屬性沒(méi)有將振幅的差異直觀反映，解釋人員只能通過(guò)反復(fù)調(diào)色標(biāo)來(lái)確定振幅之間的差異，人為干擾因素太強(qiáng)；此外，巖性、含油氣性等因素的影響使不同物理特性的地層地震響應(yīng)具有相似性，常規(guī)振幅類屬性不能將其很好的區(qū)分，因此目前的振幅類屬性往往不能滿足勘探開(kāi)發(fā)的需要。
[0003]本發(fā)明目的在于提供一種可以剔除較弱地震振幅，保留地震強(qiáng)振幅的利用地震能量差異屬性預(yù)測(cè)儲(chǔ)層的方法。

【發(fā)明內(nèi)容】

[0004]本發(fā)明通過(guò)以下步驟實(shí)現(xiàn):
[0005]I)激發(fā)合成地震并記錄，通過(guò)勘探區(qū)域已知井的鉆井分層、聲波時(shí)差數(shù)據(jù)、密度數(shù)據(jù)、地震數(shù)據(jù)體進(jìn)行井震時(shí)間深度匹配，用地震剖面中目的層段對(duì)應(yīng)的時(shí)間范圍作為地震體屬性提取的時(shí)間窗口 t ；
[0006]2)在步驟I)確定的時(shí)間窗口 t之內(nèi)，對(duì)原始地震數(shù)據(jù)體提取振幅包絡(luò)體屬性，生成目的層段的振幅包絡(luò)數(shù)據(jù)體R ；
[0007]3)在步驟I)確定的時(shí)間窗口 t之內(nèi)，對(duì)原始地震數(shù)據(jù)體提取瞬時(shí)相位余弦體屬性,生成目的層段的瞬時(shí)相位余弦數(shù)據(jù)體I ;
[0008]4)應(yīng)用傅里葉變換提取步驟2)生成的目的層段的振幅包絡(luò)數(shù)據(jù)體R的振幅譜，振幅譜橫坐標(biāo)為頻率，縱坐標(biāo)為振幅，振幅譜呈隨頻率增加振幅下降的曲線，曲線由大斜率部分和小斜率部分組成，大斜率和小斜率部分的拐點(diǎn)為突變點(diǎn)，以突變點(diǎn)為界，右側(cè)的小斜率部分為高頻段，左側(cè)的大斜率部分為低頻段；突變點(diǎn)對(duì)應(yīng)的橫坐標(biāo)為高頻段和低頻段的門(mén)檻頻率f ;
[0009]5)在步驟I)確定的時(shí)間窗口 t內(nèi)，對(duì)步驟2)生成的目的層段的振幅包絡(luò)數(shù)據(jù)體R，以步驟4)確定門(mén)檻頻率f為界進(jìn)行高通濾波，生成目的層段振幅包絡(luò)數(shù)據(jù)體R的高通濾波數(shù)據(jù)體H ；
[0010]所述的高通濾波是將高于f頻率段的振幅保留，低于f頻率段的振幅數(shù)據(jù)濾除。
[0011]6)對(duì)步驟5)生成的目的層段高通濾波數(shù)據(jù)體H進(jìn)行負(fù)數(shù)賦值，將所有負(fù)數(shù)賦值為O，所有的正值保持不變，得到高通濾波正值數(shù)據(jù)體H’ ；
[0012]7)將步驟6)的到的高通濾波正值數(shù)據(jù)體H’與步驟3)生成的目的層段瞬時(shí)相位余弦數(shù)據(jù)體I進(jìn)行乘法運(yùn)算，生成能量差異屬性數(shù)據(jù)體N ；
[0013]所述的乘法運(yùn)算是H’乘以I。
[0014]8)在原始地震數(shù)據(jù)體上根據(jù)步驟I)的井震時(shí)間深度匹配，完成勘探區(qū)域內(nèi)目的層段儲(chǔ)層或其他地質(zhì)體的頂、底的層位追蹤，生成儲(chǔ)層頂界層位Cl和儲(chǔ)層底界層位C2 ；
[0015]所述的其他地質(zhì)體是火山巖、礁。
[0016]9)將步驟7)生成的能量差異屬性數(shù)據(jù)體N以步驟8)中追蹤的層位Cl、C2為間隔時(shí)窗提取均方根振幅，生成預(yù)測(cè)儲(chǔ)層或其他地質(zhì)體平面分布圖。
[0017]本發(fā)明能夠剔除較弱地震振幅，保留地震強(qiáng)振幅，從而能夠預(yù)測(cè)具有強(qiáng)振幅特征的儲(chǔ)層或其他地質(zhì)體(火山巖、礁)。
【附圖說(shuō)明】
[0018]圖1是勘探區(qū)域某已知井的合成地震記錄標(biāo)定剖面。鉆井的目的層的深度范圍是1600-2600米，根據(jù)標(biāo)定剖面，目的層時(shí)間范圍為1500-2100毫秒。
[0019]圖2為1500-2100毫秒時(shí)窗內(nèi)的原始地震剖面，剖面中橢圓標(biāo)識(shí)的位置為儲(chǔ)層發(fā)育段，砂巖和火山巖不能有效區(qū)分。
[0020]圖3為振幅包絡(luò)體屬性剖面。是在1500-2100毫秒時(shí)窗內(nèi)對(duì)圖2所示的原始地震數(shù)據(jù)提取的振幅包絡(luò)體屬性剖面。
[0021]圖4為瞬時(shí)相位余弦體剖面。是在1500-2100毫秒時(shí)窗內(nèi)對(duì)圖2所示的原始地震數(shù)據(jù)提取的瞬時(shí)相位余弦體屬性剖面。
[0022]圖5為振幅包絡(luò)振幅譜。是對(duì)圖3所示的振幅包絡(luò)數(shù)據(jù)體提取的振幅譜，虛線左邊為大斜率部分，虛線右邊為小斜率部分，虛線對(duì)應(yīng)的點(diǎn)為突變點(diǎn)，虛線對(duì)應(yīng)的頻率為門(mén)檻頻率，本例中門(mén)檻頻率為2HZ。
[0023]圖6是振幅包絡(luò)數(shù)據(jù)體高通濾波剖面。是對(duì)圖3所示振幅包絡(luò)數(shù)據(jù)體根據(jù)圖5確定的2HZ的門(mén)檻頻率進(jìn)行高通濾波的結(jié)果，圖中大于2HZ的成分保留，小于2HZ的成分被濾掉。
[0024]圖7是高通濾波正值體剖面。是對(duì)圖6所示的高通濾波數(shù)據(jù)體進(jìn)行負(fù)數(shù)賦值生成的高通濾波正值數(shù)據(jù)體，原有的負(fù)數(shù)全部被清零。
[0025]圖8是能量差異屬性剖面。是將圖7所示的高通濾波正值數(shù)據(jù)體與圖4所示的瞬時(shí)相位余弦數(shù)據(jù)體相乘的結(jié)果剖面，剖面中橢圓標(biāo)識(shí)的火山巖和砂巖能夠很好區(qū)分。
[0026]圖9為層位解釋剖面。所示的C1、C2是根據(jù)圖1所示的標(biāo)定結(jié)果在圖2所示的原始地震數(shù)據(jù)體上追蹤的儲(chǔ)層頂界和底界的層位。
[0027]圖10為最終儲(chǔ)層分布圖。本例中的儲(chǔ)層為火山巖，圖中強(qiáng)振幅為火山巖發(fā)育段。
【具體實(shí)施方式】
[0028]以下結(jié)合具體實(shí)例和附圖詳細(xì)說(shuō)明本發(fā)明。
[0029]本發(fā)明實(shí)例步驟如下:
[0030]I)激發(fā)合成地震并記錄，通過(guò)勘探區(qū)域已知井的鉆井分層、速度數(shù)據(jù)、地震數(shù)據(jù)體進(jìn)行井震時(shí)間深度匹配，用地震剖面中目的層段對(duì)應(yīng)的時(shí)間范圍作為地震體屬性提取的時(shí)間窗口 t，圖1是勘探區(qū)域某已知井的合成地震記錄標(biāo)定剖面，如:鉆井的目的層的深度范圍是1600-2600米，根據(jù)標(biāo)定剖面，目的層時(shí)間范圍為1500-2100毫秒，因此將1500-2100毫秒確定為地震體屬性提取的時(shí)間窗口 t。
[0031]2)在步驟I)確定的時(shí)間窗口 t之內(nèi)，對(duì)原始地震數(shù)據(jù)體提