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

说明书

本申请公开了地质导向相位解释方法和装置。本申请主要针对传统的瞬时相位属性抗噪性差的问题，提出一种适用于深层复杂构造的新型相位属性。本申请充分利用了构造走向以及多尺度相位信息，能够针对性地突显小尺度异常构造，并避免高频噪声的影响，从而能检测地下更多的细节信息，还能研究河道的相对厚度变化情况，可用于地学领域中沉积构造演化、矿产勘查和油气储层预测等诸多方面的研究。

技术领域

本发明涉及地震勘探中的地震解释领域，尤其涉及一种地质导向相位解释方法和一种地质导向相位解释装置。

背景技术

地震图像中同相轴的不连续性通常与地壳结构和地层异常体有关，例如断层，河道和裂缝群。众所周知，这些地质异常通常有助于圈闭和储层的构造。因此，通过测量不连续性对地下断层，河道和裂缝群的可靠表征在地球科学的许多方面都至关重要。目前已经提出并开发了几种检测方法，例如相干性，瞬时相位，方差，谱分解和熵，测量地下的不连续性，并进一步解释异常体或层位。作为有效的工业化不连续检测器之一，瞬时相位属性具有检测细微变化的优点,但是现有的瞬时相位属性对噪声敏感,抗噪性低,得到的结果往往会给研究人员造成种种的假象。

发明内容

针对传统瞬时相位属性抗噪性差的问题，提出一种适用于深层复杂构造的地震数据解释方法。

根据本发明的一个方面，提供了一种地质导向相位解释方法，所述方法包括：对地震数据进行时频分析，得到多个频率对应的地震数据体瞬时相位θ_ω,k(t)，下标ω表示频率编号，下标k表示复地震道编号，θ_ω,k(t)代表仅包含相位信息的归一化复地震道；针对任意地震道k，对M个时间采样点的瞬时相位进行叠加，得到平均相位道针对任意地震道k，沿空间上的多个方向，分别将从编号为(k-p)至编号为(k+p-1)的2p个平均相位道进行叠加，得到叠加相位道以及分别将从编号为(k-p+1)至编号为(k+p)的2p个平均相位道进行叠加，得到叠加相位道上标du表示空间上方向的编号；计算和的复协方差道并得到复协方差道的相位针对每个频率，选择空间上多个方向对应的多个相位的最大值作为该频率对应的地质导向相位；从多个频率对应的地质导向相位属性中选择最大值作为最终的地质导向相位属性。

在一种可能的实施方式中，所述对地震数据进行时频分析包括：

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

在一种可能的实施方式中，基于下式得到平均相位道

在一种可能的实施方式中，基于下式得到和

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

在一种可能的实施方式中，基于下式计算复协方差道和

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