[发明专利]一种基于分数域显著性检测的地震属性分析方法

摘要

一种基于分数域显著性检测的地震属性分析方法，旨在提出一种基于分数域显著性检测模型，并将其应用于地震领域，获得新的地震显著性属性。其方案为读入地震剖面x(t,q),q为地震剖面的道号；对每一道地震数据做a阶分数阶傅里叶变换Fa，获得地震剖面的分数域表示Xa(u,q)；计算地震剖面的分数域表示Xa(u,q)的相位谱Pa(u,q)；利用相位谱信息Pa(u,q)进行a阶分数阶傅里叶反变换F‑a，获得地震显著性属性剖面S(t,q)。

主权项

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)称为分数阶核函数。