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

有效

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

技術領域

本發明涉及一種利用地層吸收特性進行分析的油氣預測方法，特別是關于一種基于自適應時窗的吸收分析油氣預測方法。

背景技術

從巖石物理學角度，吸收與速度和密度等參數一樣，表征了巖石的基本物理性質。地球物理油氣預測技術的發展和應用表明，地層吸收特性是識別儲層和油氣藏的有效指標之一。基于地層吸收特征分析的油氣預測技術可以在少井或無井情況下利用地震信息進行油氣預測，較其它油氣預測技術具有其特有的技術優勢，如鶯歌海盆地天然氣預測技術、瓊東南盆地天然氣預測技術和遼東灣巖性氣藏預測技術等都將地層吸收特性作為識別油氣信息有效標志之一。

目前國內外對吸收衰減的研究，主要還是基于固定時窗長度進行吸收分析油氣預測。如Mitchell（Mitchell?J.T.,Derzhi?N.T.,Lichman,E.T.,Energy?absorption?analysis:A?case?study，SEG66-th?Annual?Meeting?Expanded?Abstracts,1996:1785-1788）提出了能量吸收分析的方法，分時窗對地震記錄進行頻譜分析，通過高頻段和低頻段頻率組分衰減速率的異常變化進行油氣預測。Lichman(LichmanE,Goloshubin?G.Unified?approach?to?gas?and?fluid?detection?on?instantaneous?seismic?wavelets[G].SEG73th?Technical?Program?Expanded?Abstracts，2003:1699-1702)提出了瞬時子波吸收分析的油氣預測技術，分時窗對地震記錄進行頻譜分析，通過復賽譜技術提取時窗內的瞬時子波，進行吸收分析油氣預測。Li等(Jiang?Li,Yang?Shaoguo,Zhao?Mingjin,et?al.,Frequency-dependent?seismic?attributes?and?their?application?for?hydrocarbon?detection.80th?Annual?International?Meeting,SEG?Expanded?Abstracts,2010:1566-1569)采用小波變換，利用能量吸收分析方法進行了吸收分析油氣預測。上述吸收分析技術及相關的商業應用軟件都是基于固定時窗長度而實現，在實際應用中這些技術只適用于第三類強振幅“亮點”型油氣藏(張益明,李緒宣,符力耘.瞬時子波吸收分析技術的應用——以中國近海天然氣檢測為例.石油物探,2009,48(4):368-376)，而對于中等強度型油氣藏或弱振幅型油氣藏，這些方法存在著較強的多解性，它的應用效果受到了一定的制約。

發明內容

針對上述問題，本發明的目的是提供一種基于自適應時窗的吸收分析油氣預測方法，該方法能夠降低吸收分析油氣預測技術的多解性，拓展吸收分析油氣預測技術的應用領域，將其技術的地質適用性由強振幅“亮點”型油氣藏推廣到中等強度振幅型油氣藏，乃至弱振幅型油氣藏。

為實現上述目的，本發明采取以下技術方案：一種基于自適應時窗的吸收分析油氣預測方法，其包括以下步驟：1）選擇包含待油氣預測目標儲層地震數據的地震資料，所述的地震數據為經過保幅、保頻率處理后的疊后純波地震數據；2）在步驟1）挑選的地震資料中選取地震數據時間范圍，作為待吸收分析的地震記錄，該地震記錄至少要包括目標儲層的地震數據范圍；3）對地震記錄進行頻譜分析，確定有效頻段范圍的低截頻F₁和高截頻F₂，并設定時頻分析用的能量衰減率α和能量延遲時β；4）計算地震記錄數據體的瞬時頻率，并據此求取相應的自適應時窗；5）對地震記錄逐道進行時間域插值；6）對插值后的地震記錄逐道進行廣義S變換時頻分析，對每個時間采樣點都按照各自的自適應時窗長度計算相應的衰減因子；7）對插值后的地震記錄逐道去除背景吸收值，獲得吸收異常剖面；8）利用已知井油氣層的吸收異常值對吸收異常剖面進行吸收異常值的標定：如果標定顯示實際情況與吸收分析結果相符，就定量地預測潛在油氣；否則，返回步驟3），調整低截頻F₁和高截頻F₂，以及能量衰減率α和能量延遲時β的值，重新進行吸收分析。

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

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