一種利用地震能量差異屬性預測儲層的方法
【技術領域】
[0001]本發明涉及地球物理勘探方法，是一種利用地震能量差異屬性預測儲層的方法。
【背景技術】
[0002]隨著三維地震勘探技術的不斷進步和三維地震資料在油田勘探開發領域應用的不斷深入，地震屬性解釋的重要性日益提高，地震屬性主要有時間、振幅、頻率、相位、相干和衰減幾大類。在油田勘探開發領域，地震屬性主要應用于油藏儲層厚度、孔隙度、滲透率等參數的預測，地震資料的振幅、頻率等單項地震屬性的應用比較廣泛，其中，利用振幅屬性預測儲層是一種最常用的方法，因為其操作簡單而受到地震解釋人員的青睞。但常規振幅類屬性仍然存在一些問題，一方面，振幅類的屬性普遍分辨率較低，且沒有考慮相位的因素；另一方面，常規的振幅類屬性沒有將振幅的差異直觀反映，解釋人員只能通過反復調色標來確定振幅之間的差異，人為干擾因素太強；此外，巖性、含油氣性等因素的影響使不同物理特性的地層地震響應具有相似性，常規振幅類屬性不能將其很好的區分，因此目前的振幅類屬性往往不能滿足勘探開發的需要。
[0003]本發明目的在于提供一種可以剔除較弱地震振幅，保留地震強振幅的利用地震能量差異屬性預測儲層的方法。
【發明內容】
[0004]本發明通過以下步驟實現:
[0005]I)激發合成地震并記錄，通過勘探區域已知井的鉆井分層、聲波時差數據、密度數據、地震數據體進行井震時間深度匹配，用地震剖面中目的層段對應的時間范圍作為地震體屬性提取的時間窗口 t ；
[0006]2)在步驟I)確定的時間窗口 t之內，對原始地震數據體提取振幅包絡體屬性，生成目的層段的振幅包絡數據體R ；
[0007]3)在步驟I)確定的時間窗口 t之內，對原始地震數據體提取瞬時相位余弦體屬性,生成目的層段的瞬時相位余弦數據體I ;
[0008]4)應用傅里葉變換提取步驟2)生成的目的層段的振幅包絡數據體R的振幅譜，振幅譜橫坐標為頻率，縱坐標為振幅，振幅譜呈隨頻率增加振幅下降的曲線，曲線由大斜率部分和小斜率部分組成，大斜率和小斜率部分的拐點為突變點，以突變點為界，右側的小斜率部分為高頻段，左側的大斜率部分為低頻段；突變點對應的橫坐標為高頻段和低頻段的門檻頻率f ;
[0009]5)在步驟I)確定的時間窗口 t內，對步驟2)生成的目的層段的振幅包絡數據體R，以步驟4)確定門檻頻率f為界進行高通濾波，生成目的層段振幅包絡數據體R的高通濾波數據體H ；
[0010]所述的高通濾波是將高于f頻率段的振幅保留，低于f頻率段的振幅數據濾除。
[0011]6)對步驟5)生成的目的層段高通濾波數據體H進行負數賦值，將所有負數賦值為O，所有的正值保持不變，得到高通濾波正值數據體H’ ；
[0012]7)將步驟6)的到的高通濾波正值數據體H’與步驟3)生成的目的層段瞬時相位余弦數據體I進行乘法運算，生成能量差異屬性數據體N ；
[0013]所述的乘法運算是H’乘以I。
[0014]8)在原始地震數據體上根據步驟I)的井震時間深度匹配，完成勘探區域內目的層段儲層或其他地質體的頂、底的層位追蹤，生成儲層頂界層位Cl和儲層底界層位C2 ；
[0015]所述的其他地質體是火山巖、礁。
[0016]9)將步驟7)生成的能量差異屬性數據體N以步驟8)中追蹤的層位Cl、C2為間隔時窗提取均方根振幅，生成預測儲層或其他地質體平面分布圖。
[0017]本發明能夠剔除較弱地震振幅，保留地震強振幅，從而能夠預測具有強振幅特征的儲層或其他地質體(火山巖、礁)。
【附圖說明】
[0018]圖1是勘探區域某已知井的合成地震記錄標定剖面。鉆井的目的層的深度范圍是1600-2600米，根據標定剖面，目的層時間范圍為1500-2100毫秒。
[0019]圖2為1500-2100毫秒時窗內的原始地震剖面，剖面中橢圓標識的位置為儲層發育段，砂巖和火山巖不能有效區分。
[0020]圖3為振幅包絡體屬性剖面。是在1500-2100毫秒時窗內對圖2所示的原始地震數據提取的振幅包絡體屬性剖面。
[0021]圖4為瞬時相位余弦體剖面。是在1500-2100毫秒時窗內對圖2所示的原始地震數據提取的瞬時相位余弦體屬性剖面。
[0022]圖5為振幅包絡振幅譜。是對圖3所示的振幅包絡數據體提取的振幅譜，虛線左邊為大斜率部分，虛線右邊為小斜率部分，虛線對應的點為突變點，虛線對應的頻率為門檻頻率，本例中門檻頻率為2HZ。
[0023]圖6是振幅包絡數據體高通濾波剖面。是對圖3所示振幅包絡數據體根據圖5確定的2HZ的門檻頻率進行高通濾波的結果，圖中大于2HZ的成分保留，小于2HZ的成分被濾掉。
[0024]圖7是高通濾波正值體剖面。是對圖6所示的高通濾波數據體進行負數賦值生成的高通濾波正值數據體，原有的負數全部被清零。
[0025]圖8是能量差異屬性剖面。是將圖7所示的高通濾波正值數據體與圖4所示的瞬時相位余弦數據體相乘的結果剖面，剖面中橢圓標識的火山巖和砂巖能夠很好區分。
[0026]圖9為層位解釋剖面。所示的C1、C2是根據圖1所示的標定結果在圖2所示的原始地震數據體上追蹤的儲層頂界和底界的層位。
[0027]圖10為最終儲層分布圖。本例中的儲層為火山巖，圖中強振幅為火山巖發育段。
【具體實施方式】
[0028]以下結合具體實例和附圖詳細說明本發明。
[0029]本發明實例步驟如下:
[0030]I)激發合成地震并記錄，通過勘探區域已知井的鉆井分層、速度數據、地震數據體進行井震時間深度匹配，用地震剖面中目的層段對應的時間范圍作為地震體屬性提取的時間窗口 t，圖1是勘探區域某已知井的合成地震記錄標定剖面，如:鉆井的目的層的深度范圍是1600-2600米，根據標定剖面，目的層時間范圍為1500-2100毫秒，因此將1500-2100毫秒確定為地震體屬性提取的時間窗口 t。
[0031]2)在步驟I)確定的時間窗口 t之內，對原始地震數據體提