[發(fā)明專利]一種基于譜分解和相干屬性融合流程的水道復合體沉積要素精細識別方法

在審

• 申請?zhí)枺?/td>: 201810995608.0
• 申請日: 2018-08-29
• 公開（公告）號: CN109116414A
• 公開（公告）日: 2019-01-01
• 發(fā)明（設計）人: 李全;康洪全;李欣;楊婷;逄林安;曹向陽;常邁;李丹;解東寧;段俊;劉勵云;姜慧超
• 申請（專利權）人: 中國海洋石油集團有限公司;中海油研究總院有限責任公司
• 主分類號: G01V1/28
• 分類號: G01V1/28;G01V1/30
• 代理公司: 北京紀凱知識產(chǎn)權代理有限公司 11245
• 代理人: 關暢;王春霞

1.一種基于譜分解和相干屬性融合流程的水道復合體沉積要素精細識別方法，它包括以下步驟：

(1)數(shù)據(jù)體頻率分析：根據(jù)目標水道復合體發(fā)育的位置分析地震資料頻帶，得到所述地震資料的最低頻率、最高頻率和主頻，并計算得到能識別的沉積要素最大和最小厚度；

(2)形成調(diào)諧數(shù)據(jù)體：結(jié)合步驟(1)得到的數(shù)據(jù)對目標水道復合體進行S變換，S變換的時窗通過沿層開取時窗或在水道頂?shù)讓又g進行；通過計算，在垂向上生成頻率連續(xù)變化的振幅和相位數(shù)據(jù)體，得到目標水道復合體內(nèi)沉積要素的調(diào)諧頻率；然后通過觀察不同頻率下對應的沉積要素的成像特征，確定不同沉積要素對應的有效頻率；

(3)建立離散時頻率體：將步驟(1)中所述地震資料分解成一系列時間域離散頻率振幅數(shù)據(jù)體，該數(shù)據(jù)體在垂向上為時間，每個數(shù)據(jù)體只包含需要表征的沉積要素對應的頻率成分，反映了特定頻率振幅的空間變化特征；

(4)進行相干體運算：通過地震上波形和振幅側(cè)向上的變化與斷層或者沉積要素的邊界進行相干體運算，得到沉積要素的邊界信息；

(5)體融合運算：將步驟(3)中得到的不同沉積要素的對應的離散時頻體分別與步驟(4)運算的相干體疊加、融合顯示，得到不同沉積要素對應下的相干時頻體；

(6)應用等時切片或地層切片技術對所述相干視頻體進行可視化顯示，重構(gòu)所述目標水道復合體單一沉積要素垂向演化過程和平面分布特征，以精細識別沉積要素。

2.根據(jù)權利要求1所述的方法，其特征在于：步驟(1)中，所述沉積要素包括天然堤、點壩、單期水道、決口扇和沉積波；

計算得到能識別的沉積要素最大和最小厚度的公式如式Ⅰ和式Ⅰ-1所示：

H_max＝V/4f_max 式Ⅰ

H_min＝V/4f_min 式Ⅰ-1

式Ⅰ和式Ⅰ-1中，f_max為步驟1中目標水道復合體的最大頻率；f_min為步驟1中目標水道復合體的最小頻率；V為地震資料速度；H_max和H_min分別為可分辨沉積要素的最大和最小厚度。

3.根據(jù)權利要求1或2所述的方法，其特征在于：所述步驟(2)中的S變換是時間和頻率的聯(lián)合函數(shù)，以時間和頻率為變量來描述信號的能量密度或信號強度，其根據(jù)如下式Ⅱ所示公式進行計算：

式Ⅱ中，t表示時間，f表示頻率，τ為高斯窗在時間軸上的位置，時窗寬度隨頻率的變化而改變，頻率的倒數(shù)決定了時窗的大小，s(t)表示時間域輸入信號，exp為自然對數(shù)函數(shù)，π為圓周率，i表示虛數(shù)單位，dt表示積分函數(shù)。

4.根據(jù)權利要求1-3中任一項所述的方法，其特征在于：步驟(3)中根據(jù)步驟(2)中確定的目標水道復合體的沉積要素對應頻率進行地震資料分解；

所述步驟3中的時頻率體，包括不同厚度對應下的沉積要素頻率15HZ、20HZ、25Hz、30Hz、40Hz和50Hz頻率體。

5.根據(jù)權利要求1-4中任一項所述的方法，其特征在于：所述步驟(4)中，所示相干體的運算采取BG Fast算法，該算法在所示地震數(shù)據(jù)中傾角定義為水平時間梯度與垂直時間梯度平方和的均方根值，而方位角定義為垂直時間梯度與水平時間梯度的比值的反正切值。

6.根據(jù)權利要求1-5中任一項所述的方法，其特征在于：所述步驟(5)中，所述體融合運算通過蘭德馬克公司的Geoprobe軟件陰影疊置功能實現(xiàn)。