[發明專利]一種基于頻域顯著性檢測的地震屬性分析方法

在審

• 申請號: 201610051252.6
• 申請日: 2016-01-26
• 公開（公告）號: CN105700014A
• 公開（公告）日: 2016-06-22
• 發明（設計）人: 彭真明;王雨青;邢艷;桑湘鑫;韓演;楊維;何艷敏
• 申請（專利權）人: 電子科技大學
• 主分類號: G01V1/30
• 分類號: G01V1/30
• 代理公司: 成都弘毅天承知識產權代理有限公司 51230
• 代理人: 徐金瓊

1.一種基于頻域顯著性檢測的地震屬性分析方法，其特征在于，包括以下步驟：

步驟1：讀入地震剖面x(t,q),q為地震剖面的道號；

步驟2：對每一道地震數據做傅里葉變換F，獲得地震剖面的頻域表示X(f,q)；傅里葉變換如下：

傅里葉變換算子為F，對地震剖面x(t,q)的每一道數據做傅里葉變換，獲得頻域表示X(f,q)

X ( f , q ) = F [ x ( t , q ) ] = ∫ - ∞ ∞ x ( t , q ) e - j 2 π f t d t - - - ( 1 ) ]]>

式中，q表示地震剖面的道號，f是頻率域變量，t是時間變量，j是虛數單位，e是自然常數；

步驟3：計算地震剖面的頻域表示X(f,q)的相位譜P(f,q)；頻域相位譜如下：

對于頻域表示X(f,q)，其相位譜P(f,q)為

P ( f , q ) = a r c t a n ( Im { X ( f , q ) } Re { X ( f , q ) } ) - - - ( 2 ) ]]>

式中，Im{·}表示虛部，Re{·}表示實部，arctan表示反正切函數；

步驟4：利用相位譜信息P(f,q)進行傅里葉反變換F-1，獲得顯著性剖面S(t,q)；利用相位譜信息進行傅里葉反變換如下：

已知相位譜P(f,q)，顯著性剖面S(t,q)表示為

S(t,q)＝|F-1{ejP(f,q)}| (3)

式中，q表示地震剖面的道號，t是時間變量，f是頻率域變量，j是虛數單位，e是自然常數，F-1表示傅里葉反變換，即

F - 1 { e j P ( f , q ) } = ∫ - ∞ ∞ e j P ( f , q ) e j 2 π f t d f - - - ( 4 ) ; ]]>

步驟5：用低通濾波器g對顯著性剖面S(t,q)進行濾波，作為最終的地震顯著性屬性剖面；濾波過程，如下：

對顯著性剖面S(t,q)，q表示地震剖面的道號，采用低通濾波器g，進行濾波后地震顯著性剖面FS(t,q)表示為

FS(t,q)＝g*|S(t,q)|2 (5)。