[发明专利]一种提高边缘检测能力的地震数据处理方法

说明书

本发明公开了一种提高边缘检测能力的地震数据处理方法，包括以下步骤：S1.确定待处理的地震数据x(t)；S2.对地震数据x(t)进行解析，得到其解析信号z(t)；S3.对解析信号z(t)求时频谱得到Z₁(t,f)，将时频谱的实部和虚部表示为G₁(t,f)和G₂(t,f)；S4.引入阶数p和符号函数sign( )保护相位，对时频谱Z₁(t,f)的实部G₁(t,f)和虚部G₂(t,f)进行保相位的广义希尔伯特变换求得g₁(t)和g₂(t)；S5.计算得到保相位的广义希尔伯特变换信号的实部x_p(t)和虚部y_p(t)，进而得到保相位广义希尔伯特变换后的信号z_p(t)。本发明由于引入了符号函数保护相位，故避免了广义希尔伯特变换会更改实部或虚部的符号，导致错误输出的问题。

技术领域

本发明涉及地震数据处理，特别是涉及一种提高边缘检测能力的地震数据处理方法。

背景技术

提高地震资料中地质体边缘检测能力一直以来都是地震资料处理人员追求的目标，目前用来提高地震资料边缘检测能力的方法有希尔伯特变换、广义希尔伯特变换等已广泛应用于信号处理中。

广义希尔伯特变换通过在时频域引入阶数n来提高希尔伯特变换的分辨率和抗噪性。但是，当时频谱的虚部或实部小于零并且选择阶数n作为偶数与奇数之比时，广义希尔伯特变换会更改实部或虚部的符号，这将导致错误的输出，无法产生理想的结果。

发明内容

本发明的目的在于克服现有技术的不足，提供一种提高边缘检测能力的地震数据处理方法，先对地震数据进行解析，然后求解析信号的时频谱，再结合符号函数进行保相位的广义希尔伯特变换，由于引入了符号函数保护相位，故避免了广义希尔伯特变换会更改实部或虚部的符号，导致错误输出的问题。

本发明的目的是通过以下技术方案来实现的：一种提高边缘检测能力的地震数据处理方法，包括以下步骤：

S1.确定待处理的地震数据x(t)；

S2.对地震数据x(t)进行希尔伯特变换，得到其解析信号z(t)；

S3.对解析信号z(t)进行改进的短时傅里叶变换，得到时频谱得到Z₁(t,f)，将时频谱的实部和虚部表示为G₁(t,f)和G₂(t,f)；

所述的G₁(t,f)和G₂(t,f)表示如下：

其中，Re[]表示取实部，Im[]表示取虚部。

S4.引入阶数p和符号函数sign()保护相位，对时频谱Z₁(t,f)的实部G₁(t,f)和虚部G₂(t,f)进行保相位的广义希尔伯特变换求得g₁(t)和g₂(t)，具体方式如下：

S5.计算得到保相位的广义希尔伯特变换信号的实部x_p(t)和虚部y_p(t)，进而得到保相位的广义希尔伯特变换信号z_p(t)，具体包括以下步骤：

计算保相位的广义希尔伯特变换信号的实部x_p(t)和虚部y_p(t)：

根据变换信号的实部和虚部，确定保相位的广义希尔伯特变换信号z_p(t)：