[發(fā)明專利]一種基于自適應時窗的吸收分析油氣預測方法

有效

• 申請?zhí)枺?/td>: 201310054085.7
• 申請日: 2013-02-20
• 公開（公告）號: CN103197347A
• 公開（公告）日: 2013-07-10
• 發(fā)明（設計）人: 張顯文;李維新;姚振興;趙偉;李緒宣;胡光義;范廷恩;桑淑云;朱振宇;王建花;印海燕;王清振
• 申請（專利權）人: 中國海洋石油總公司;中海油研究總院
• 主分類號: G01V1/30
• 分類號: G01V1/30
• 代理公司: 北京紀凱知識產(chǎn)權代理有限公司 11245
• 代理人: 徐寧;關暢

1.一種基于自適應時窗的吸收分析油氣預測方法，其包括以下步驟：

1）選擇包含待油氣預測目標儲層地震數(shù)據(jù)的地震資料，所述的地震數(shù)據(jù)為經(jīng)過保幅、保頻率處理后的疊后純波地震數(shù)據(jù)；

2）在步驟1）挑選的地震資料中選取地震數(shù)據(jù)時間范圍，作為待吸收分析的地震記錄，該地震記錄至少要包括目標儲層的地震數(shù)據(jù)范圍；

3）對地震記錄進行頻譜分析，確定有效頻段范圍的低截頻F₁和高截頻F₂，并設定時頻分析用的能量衰減率α和能量延遲時β；

4）計算地震記錄數(shù)據(jù)體的瞬時頻率，并據(jù)此求取相應的自適應時窗；

5）對地震記錄逐道進行時間域插值；

6）對插值后的地震記錄逐道進行廣義S變換時頻分析，對每個時間采樣點都按照各自的自適應時窗長度計算相應的衰減因子；

7）對插值后的地震記錄逐道去除背景吸收值，獲得吸收異常剖面；

8）利用已知井油氣層的吸收異常值對吸收異常剖面進行吸收異常值的標定：

如果標定顯示實際情況與吸收分析結果相符，就定量地預測潛在油氣；

否則，返回步驟3），調(diào)整低截頻F₁和高截頻F₂，以及能量衰減率α和能量延遲時β的值，重新進行吸收分析。

2.如權利要求1所述的一種基于自適應時窗的吸收分析油氣預測方法，其特征在于：

所述步驟6）包括以下小步驟：

①選取插值后的地震記錄中的一道記錄Trace(i)，i＝1～N_trace，進行廣義S變換計算此道記錄的時頻譜TF_Spectrum(t，f)；

②選取此道記錄中的一個時間采樣點j，j＝1～NT，其中，NT是此道記錄中采樣點個數(shù)，以廣義S變換后得到的時頻譜的瞬時點譜TF_Spectrum(t_j，f)作為該時間采樣點的頻譜F_Spectrum(f)；

③利用復賽譜技術計算該時間采樣點的瞬時子波振幅譜IW_FSpectrum(f)；

④計算該時間采樣點的瞬時子波振幅譜IW_FSpectrum(f)的主頻f₀；

⑤對瞬時子波振幅譜IW_FSpectrum(f)在低頻段范圍F₁～f₀內(nèi)進行擬合，計算該時間采樣點的流體衰減因子Q_f(j)；

⑥對瞬時子波振幅譜IW_FSpectrum(f)在高頻段范圍f₀～F₂內(nèi)進行擬合，計算該時間采樣點的氣體衰減因子Q_g(j)；

⑦檢查此道記錄中是否還有未計算衰減因子的時間采樣點：

如果是，返回步驟②；

否則，進入步驟⑧；

⑧檢查插值后的地震記錄中是否還有未進行吸收分析的記錄：

如果是，返回步驟①；

否則，進入步驟7）。

3.如權利要求1所述的一種基于自適應時窗的吸收分析油氣預測方法，其特征在于：所述步驟7）中，對步驟6）獲得的衰減因子逐道進行二階多項式曲線擬合，得到每道記錄的背景吸收值，將各衰減因子減去自己的背景吸收值，從而得到吸收異常剖面。

4.如權利要求2所述的一種基于自適應時窗的吸收分析油氣預測方法，其特征在于：所述步驟7）中，對步驟6）獲得的衰減因子逐道進行二階多項式曲線擬合，得到每道記錄的背景吸收值，將各衰減因子減去自己的背景吸收值，從而得到吸收異常剖面。