[发明专利]一种基于时间重采样和同步挤压变换的地震数据处理方法

权利要求书

1.一种基于时间重采样和同步挤压变换的地震数据处理方法，其特征在于：包括以下步骤：

S1.采集记录地震原始数据x(t)；

S2.对原始信号x(t)进行小波变换得到其对应的时频域的复数信号H(t,f)，并由H(t,f)计算得到同步挤压变换的结果P(t,υ)；

S3.利用复数信号H(t,f)对时间t进行傅里叶变换得到其频谱Q(η,f)，计算Q(η,f)的最大动态范围，单位为分贝，找到其对应的有效频带的截止频率，取其最大频率记为G(f)；

S4.确定新的自适应的时间重采样间隔Δτ_g(f)，进而得到新的频谱Q_r(η,f)，对新的频谱Q_r(η,f)进行反傅里叶变换得到新的复数信号H_r(t,f)；

所述步骤S4包括以下子步骤：

S401.确定自适应的时间重采样间隔Δτ_g(f)：

Δt₀(f)＝G(f)^-1·36^-1

S402.计算Q_r(η,f)：

S403.计算H_r(t,f)：

H_r(t,f)＝∫Q_r(η,f)exp(i2πηt)dη；

S5.由新的自适应时间重采样间隔Δτ_g(f)和新的复数信号H_r(t,f)，计算高精度瞬时频率S₁(t,f)，通过H(t,f)和S₁(t,f)计算自适应时间重采样的高分辨率同步挤压变换结果P₁(t,υ)。

2.根据权利要求1所述的一种基于时间重采样和同步挤压变换的地震数据处理方法，其特征在于：所述步骤S2包括以下子步骤：

S201.对原始信号x(t)进行小波变换得到其对应的时频域的复数信号H(t,f)；

S202.设复数信号H(t,f)的实部为h₁(t,f)，虚部为h₂(t,f)，计算瞬时频率S(t,f)：

S203.计算同步挤压变换结果P(t,υ)：

P(t,υ)＝∫H(t,f)df,ifS(t,f)∈[υ,υ+Δυ]；

υ是频率，Δυ是频率的增量，瞬时频率时不进行累加。

3.根据权利要求2所述的一种基于时间重采样和同步挤压变换的地震数据处理方法，其特征在于：所述步骤S3包括以下子步骤：

S301.将H(t,f)对时间t进行傅里叶变换得到其频谱Q(η,f)：

Q(η,f)＝∫H(t,f)exp(-i2πηt)dt；

S302.计算G(f)：

G(f)＝max[C_f(η)]；

|Q(η,f)|是对Q(η,f)求模，表示H(t,f)的振幅谱，B_f(η_p)表示H(t,f)的最大振幅，C_f(η)表示满足振幅谱最大动态范围的频率范围。

4.根据权利要求2所述的一种基于时间重采样和同步挤压变换的地震数据处理方法，其特征在于：所述步骤S5包括以下子步骤：

S501.由Δτ_g(f)和H_r(t,f)，计算关于H_r(t,f)的高精度瞬时频率：

式中，x₁(t,f)和x₂(t,f)分别是H_r(t,f)的实部和虚部；

S502.计算自适应时间重采样的高分辨率同步挤压变换结果P₁(t,υ)：

P₁(t,υ)＝∫H(t,f)df,ifS₁(t,f)∈[υ,υ+Δυ]。