[发明专利]地质导向相位解释方法和装置

权利要求书

1.一种地质导向相位解释方法，其特征在于，所述方法包括：

对地震数据进行时频分析，得到多个频率对应的地震数据体瞬时相位θ_ω,k(t)，下标ω表示频率编号，下标k表示复地震道编号，θ_ω,k(t)代表仅包含相位信息的归一化复地震道；

针对任意地震道k，对M个时间采样点的瞬时相位进行叠加，得到平均相位道

针对任意地震道k，沿空间上的多个方向，分别将从编号为(k-p)至编号为(k+p-1)的2p个平均相位道进行叠加，得到叠加相位道以及分别将从编号为(k-p+1)至编号为(k+p)的2p个平均相位道进行叠加，得到叠加相位道上标du表示空间上方向的编号；

计算和的复协方差道并得到复协方差道的相位

针对每个频率，选择空间上多个方向对应的多个相位中的最大值作为该频率对应的地质导向相位；

从多个频率对应的地质导向相位属性中选择最大值作为最终的地质导向相位属性。

2.根据权利要求1所述的方法，其特征在于，所述对地震数据进行时频分析包括：

采用Gabor变换对地震数据进行时频分析。

3.根据权利要求1所述的方法，其特征在于，基于下式得到平均相位道

4.根据权利要求1所述的方法，其特征在于，基于下式得到和

其中，ω(|l-k|)是加权系数，|*|表示取绝对值；τ_l表示地震道k与地震道l之间的时间延迟。

5.根据权利要求1所述的方法，其中，基于下式计算复协方差道及其相位

其中，[*]^H表示取复共轭。

6.一种地质导向相位解释装置，其特征在于，所述装置包括：

时频分析单元，用于对地震数据进行时频分析，得到多个频率对应的地震数据体瞬时相位θ_ω,k(t)，下标ω表示频率编号，下标k表示复地震道编号，θ_ω,k(t)代表仅包含相位信息的归一化复地震道；

时窗叠加单元，用于针对任意地震道k，对M个时间采样点的瞬时相位进行叠加，得到平均相位道

空间叠加单元，用于针对任意地震道k，沿空间上的多个方向，分别将从编号为(k-p)至编号为(k+p-1)的2p个平均相位道进行叠加，得到叠加相位道以及分别将从编号为(k-p+1)至编号为(k+p)的2p个平均相位道进行叠加，得到叠加相位道上标du表示空间上方向的编号；

协方差计算单元，用于计算和的复协方差道并得到复协方差道的相位

导向相位选择单元，用于针对每个频率，选择空间上多个方向对应的多个相位中的最大值作为该频率对应的地质导向相位；

导向相位确定单元，用于从多个频率对应的地质导向相位属性中选择最大值作为最终的地质导向相位属性。

7.根据权利要求6所述的装置，其特征在于，在所述时频分析单元中，所述对地震数据进行时频分析包括：

采用Gabor变换对地震数据进行时频分析。

8.根据权利要求6所述的装置，其特征在于，在所述时窗叠加单元中，基于下式得到平均相位道

9.根据权利要求6所述的装置，其特征在于，在所述空间叠加单元中，基于下式得到和

其中，ω(|l-k|)是加权系数，|*|表示取绝对值；τ_l表示地震道k与地震道l之间的时间延迟。

10.根据权利要求6所述的装置，其中，在所述协方差计算单元中，基于下式计算复协方差道及其相位

其中，[*]^H表示取复共轭。