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

在審

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

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

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

步驟2：對每一道地震數據做a階分數階傅里葉變換Fa，獲得地震剖面的分數域表示Xa(u,q)；

步驟3：計算地震剖面的分數域表示Xa(u,q)的相位譜Pa(u,q)；

步驟4：利用相位譜信息Pa(u,q)進行a階分數階傅里葉反變換F-a，獲得地震顯著性屬性剖面S(t,q)；

所述步驟2涉及分數階傅里葉變換，如下：

a階分數階傅里葉變換算子為Fa，對地震剖面x(t,q)的每一道數據做a階分數階傅里葉變換，獲得分數域表示Xa(u,q)為

式中，q表示地震剖面的道號，a表示分數階傅里葉變換的階次，0＜|a|＜2，其對應的旋轉角度為φ＝aπ/2，Ka(u,t)稱為分數階核函數，u是分數域變量，t是時間變量，j是虛數單位，e是自然常數；

其中所述步驟3涉及分數域相位譜，具體理論如下：

對于分數域表示Xa(u,q)，其相位譜Pa(u,q)為

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

其中所述步驟4涉及利用相位譜信息進行分數階傅里葉反變換，如下：

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

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

式中，-a表示分數階傅里葉變換的階次，0＜|a|＜2，其對應的旋轉角度為φ＝-aπ/2，K-a(u,t)稱為分數階核函數。