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

权利要求书

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