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

摘要

一种基于频域显著性检测的地震属性分析方法，旨在将显著性检测应用于地震领域，获得新的且有效的地震属性，其方案为读入地震剖面x(t,q),q为地震剖面的道号；对每一道地震数据做傅里叶变换F，获得地震剖面的频域表示X(f,q)；计算地震剖面的频域表示X(f,q)的相位谱P(f,q)；利用相位谱信息P(f,q)进行傅里叶反变换F‑1，初步获得显著性剖面S(t,q)；用低通滤波器g对显著性剖面S(t,q)进行滤波，作为最终的地震显著性属性剖面。

主权项

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)。