[發(fā)明專利]一種利用疊后地震數(shù)據(jù)進行油氣檢測的方法

有效

• 申請?zhí)枺?/td>: 201410602565.7
• 申請日: 2014-10-31
• 公開（公告）號: CN104360382B
• 公開（公告）日: 2018-12-11
• 發(fā)明（設(shè)計）人: 秦廣勝;蔡其新;姜貽偉;蘇云;張紅霖;晉達;汪露露;常正時
• 申請（專利權(quán)）人: 中國石油化工股份有限公司;中國石油化工股份有限公司中原油田分公司物探研究院
• 主分類號: G01V1/28
• 分類號: G01V1/28
• 代理公司: 鄭州睿信知識產(chǎn)權(quán)代理有限公司 41119
• 代理人: 趙敏

1.一種利用疊后地震數(shù)據(jù)進行油氣檢測的方法，其特征是包括以下步驟:

(1)數(shù)據(jù)準(zhǔn)備：獲取工區(qū)的地震疊加純波數(shù)據(jù)、工區(qū)內(nèi)各井位的坐標(biāo)及各井目的層段含油氣層厚度和解釋的目的層頂、底時間數(shù)據(jù)；

(2)確定地震數(shù)據(jù)待分析時間范圍和滑動時窗長度；

(3)對工區(qū)地震疊后純波數(shù)據(jù)提取能量系數(shù)，構(gòu)建地震波衰減因子,獲得衰減因子數(shù)據(jù)體：

針對地震疊后純波數(shù)據(jù)的每一地震道，在步驟(2)中確定的分析時間范圍內(nèi)對每一時間采樣點以步驟(2)中確定的滑動時窗長度做短時窗傅立葉變換，獲得該時間采樣點的頻譜，確定加權(quán)平均頻率和有效頻帶的低頻端和高頻端；以上述步驟確定的加權(quán)平均頻率為分界點，把該時間采樣點的頻譜劃分為低頻和高頻兩個部分，把頻譜的低頻部分以頻率反向加權(quán)的形式求和，獲得頻譜的低頻部分能量；把頻譜的高頻部分求和，獲得頻譜的高頻部分能量；

利用上述步驟獲取的低頻部分能量與高頻部分能量的比值，獲得該時間采樣點的能量系數(shù)，并通過該時間采樣點的能量系數(shù)和上述步驟確定的加權(quán)平均頻率，構(gòu)建該時間采樣點的地震波衰減因子；

(4)重復(fù)上述步驟(3)，獲得工區(qū)地震疊后純波數(shù)據(jù)每一道在步驟(2)確定分析時間范圍內(nèi)每一時間采樣點的地震波衰減因子，得到工區(qū)地震波衰減因子數(shù)據(jù)體,根據(jù)工區(qū)地震波衰減因子數(shù)據(jù)體數(shù)值的大小，分析、識別工區(qū)內(nèi)有利含油氣目標(biāo)區(qū)；

對工區(qū)地震疊后純波數(shù)據(jù)提取能量系數(shù)，構(gòu)建地震波衰減因子,獲得衰減因子數(shù)據(jù)體，包括以下步驟：

(3.1)針對地震疊后純波數(shù)據(jù)的每一地震道，在確定的分析時間范圍內(nèi)對每一時間采樣點以確定的滑動時窗長度做短時窗傅立葉變換，獲得該時間采樣點的頻譜A(i,f)，其中，i為時間采樣點的序號，f為頻率采樣點的序號，由A(i,f)中最大振幅的1/20為下限確定有效頻帶的低頻端f_L和高頻端f_H，加權(quán)平均頻率f_w通過如下方法確定

(3.2)以上述步驟(3.1)中確定的加權(quán)平均頻率為分界點，把該時間采樣點的頻譜劃分為低頻和高頻兩個部分，把頻譜的低頻部分以頻率反向加權(quán)的形式求和，獲得頻譜的低頻部分能量E_L

把頻譜的高頻部分求和，獲得頻譜的高頻部分能量E_H

(3.3)利用步驟(3.2)獲取的低頻部分能量與高頻部分能量的比值，獲得該時間采樣點的能量系數(shù)C(i)，

通過該時間采樣點的能量系數(shù)和上述步驟(3.1)確定的加權(quán)平均頻率，構(gòu)建該時間采樣點的地震波衰減因子Q(i)

其中，b為常數(shù)，此處取值為1；

(3.4)重復(fù)上述步驟(3.1)-(3.3)，獲得工區(qū)地震疊后純波數(shù)據(jù)每一道在分析時間范圍內(nèi)每一時間采樣點的地震波衰減因子，得到工區(qū)地震波衰減因子數(shù)據(jù)體；

(3.5)對上述步驟(3.4)獲得的工區(qū)地震波衰減因子數(shù)據(jù)體，根據(jù)其數(shù)值的大小，分析、識別工區(qū)內(nèi)有利含油氣目標(biāo)區(qū)。

2.根據(jù)權(quán)利要求1所述的一種利用疊后地震數(shù)據(jù)進行油氣檢測的方法，其特征是：步驟(2)的確定滑動時窗長度包括以下步驟：

(2.1)根據(jù)井位坐標(biāo)在工區(qū)地震疊后純波數(shù)據(jù)中抽取井旁地震道x(i)，i為所抽取井旁地震道時間采樣點的序號；

(2.2)以不同滑動時窗長度，利用傅立葉變換方法獲得地震道x(i)的目的層段每個時間采樣點的頻譜A(i,f)，即

A(i,f)＜--＞FFT(x'(i))

其中，i為時間采樣點的序號，f為頻率采樣點的序號，x'(i)是地震道x(i)加海明窗的結(jié)果，即

海明窗的表達式為

其中，j為海明窗函數(shù)的時間采樣點序號，N為滑動時窗長度的樣點數(shù)，π為圓周率；

(2.3)分析不同時窗長度對頻率域振幅譜形態(tài)及其對油氣層反映的能力，確定滑動時窗長度。