[發明專利]基于分數階希爾伯特變換的薄砂體平面展布預測方法

有效

• 申請號: 202010419191.0
• 申請日: 2020-05-18
• 公開（公告）號: CN111598340B
• 公開（公告）日: 2021-07-23
• 發明（設計）人: 童明勝;陳杰;陳果;徐明華;方杰;李開建;羅鵬曦
• 申請（專利權）人: 中國石油天然氣集團有限公司;中國石油集團川慶鉆探工程有限公司
• 主分類號: G06Q10/04
• 分類號: G06Q10/04;G06Q50/02;G06F17/14;G01V1/30
• 代理公司: 成都天嘉專利事務所(普通合伙) 51211
• 代理人: 毛光軍

1.一種基于分數階希爾伯特變換的薄砂體平面展布預測方法，其特征在于：先計算10個階數的三維地震資料的分數階希爾伯特變換，再對砂體發育段在縱向上進行多個砂體發育小層段的劃分，然后提取各個砂體發育小層段的屬性平面圖，基于實鉆井在各個砂體發育小層段的發育情況優選最佳的階數對應的振幅屬性平面圖進行薄砂體平面展布預測；

所述方法具體包括如下步驟：

a、輸入疊后三維地震數據體；

b、分別計算10個階數值的分數階希爾伯特變換；

c、根據井震標定分析得出砂體段發育的縱向位置，對10個分數階希爾伯特變換的結果數組在縱向上提取各個砂體發育小層段的10個三維層段數組；

d、統計研究區內的實鉆井在各個砂體發育小層段內的砂體發育情況；

e、對縱向上發育的第一個砂體發育小層段，計算當前層段內的三維數組的10個均方根振幅屬性，并繪制10個均方根振幅平面圖；

f、根據實鉆井在當前層段的砂體的發育情況，計算步驟e中10個均方根振幅平面圖的砂體平面展布預測符合率，選取預測符合率最高的階數對應的均方根振幅平面圖進行薄砂體平面展布預測；

g、對縱向上的其它砂體發育的小層段，重復步驟e和f進行薄砂體平面展布預測。

2.根據權利要求1所述的基于分數階希爾伯特變換的薄砂體平面展布預測方法，其特征在于：所述步驟a中，輸入疊后三維地震數據體w(i,j,n)，i代表Inline測線序號，總測線數記為Nline；j代表Xline測線序號，總測線數記為Nxline；n代表單道地震記錄的采樣點序號，總采樣點數記為N。

3.根據權利要求2所述的基于分數階希爾伯特變換的薄砂體平面展布預測方法，其特征在于：所述步驟b中，令v＝-1.0,-0.8,-0.6,-0.4,-0.2,0.2,0.4,0.6,0.8,1，對三維數組w(i,j,n)的每個地震道記錄進行分數階希爾伯特變換，得到各個階數的三維數組w(i,j,n)的分數階希爾伯特變換的結果數組w_v(i,j,n)，共10個三維數組。

4.根據權利要求3所述的基于分數階希爾伯特變換的薄砂體平面展布預測方法，其特征在于：所述步驟b具體過程如下：

b1、令v＝-1.0；

b2、令i＝1,j＝1提取當前地震道數組s(n)，s(n)為一維數組，數組的長度是N；

b3、對目標一維數組進行s(n)做DFT，得到S(k)(k＝0,1,…,N-1)，其中k＝N/2,…,N-1對應負頻率；

b4、對S(k)進行分頻率計算：

b5、對y(k)進行逆DFT，取其實部即得到s(n)的分數階希爾伯特變換s_v(n)，并寫入三維數組w_v(i,j,n)；

b6、增加i,j的數值，重復步驟b2和步驟b5，直到i＝Nline,j＝Nxline；得到當前階數v的分數階希爾伯特變換的三維數組w_v(i,j,n)；

b7、改變v的數值，重復步驟b2和步驟b6，直到v＝1.0，得到10個階數對應的10個三維數組。

5.根據權利要求4所述的基于分數階希爾伯特變換的薄砂體平面展布預測方法，其特征在于：所述步驟c中，根據井震標定分析得出砂體段發育的縱向位置，對步驟b的分數階希爾伯特變換的結果數組w_v(i,j,n)在縱向上采用等時切片方法提取各個砂體發育小層段的三維層段數組Sand_{l_v}(i,j,nc)，l代表縱向上砂體發育小層段的序號，縱向上砂體發育小層段總數記為Nsand；nc代表當前層段內的地震記錄的采樣點序號，總采樣點數記為Nc。

6.根據權利要求5所述的基于分數階希爾伯特變換的薄砂體平面展布預測方法，其特征在于：所述步驟d中，令研究區內的實鉆井數目是Nwell，統計步驟c中各個層段內的單井砂體發育情況，并寫入二維數組P(a,l)，a代表實鉆井序號，l代表縱向上發育砂體層段的序號，P(a,l)＝0代表該井在該層段砂體不發育，P(a,l)＝1代表該井在該層段砂體發育。