[发明专利]一种叠后地震波相对波阻抗反演方法及系统

权利要求书

1.一种叠后地震波相对波阻抗反演方法，其特征在于，包括：

获取目标层位的叠后地震波数据；

对叠后地震波数据逐道进行同步挤压小波包变换，以在时频域内提取预设阈值范围内的模极大值曲线；

根据预选的收缩母子波参数，采用同步挤压小波包反变换从模极大值曲线中进行叠后地震波数据的重构；

对重构的叠后地震波数据进行带通积分得到相对波阻抗；

在时频域内提取预设阈值范围内的模极大值曲线的过程为：对经同步挤压小波包变换得到的时频能量分布按列优先原则进行矩阵元素的扫描，得到尺度因子f^j下，区间t∈(t₀-δ,t₀+δ)内，满足|W(f^j,t)|≤|W(f^j,t₀)|的点(f^j,t₀)，则点(f^j,t₀)为同步挤压小波包变换的局部模极大值点，其中，δ为常数；在预设阈值范围内，将局部模极大值点连接后，得到模极大值曲线。

2.如权利要求1所述的一种叠后地震波相对波阻抗反演方法，其特征在于，对叠后地震波数据逐道进行同步挤压小波包变换得到叠后地震波数据的时频能量分布，在时频能量分布中提取预设阈值范围内的模极大值曲线。

3.如权利要求1所述的一种叠后地震波相对波阻抗反演方法，其特征在于，对叠后地震波数据x(t)进行同步挤压小波包变换的过程为：

T_x(v,b)＝∫|W_x(p,b)|²δ(v_x(p,b)-v)dp；

其中，W_x(p,b)是叠后地震波数据x(t)在点(p,b)处的小波包变换，v_x(p,b)为叠后地震波数据x(t)在点(p,b)处的局域小波向量估计。

4.如权利要求3所述的一种叠后地震波相对波阻抗反演方法，其特征在于，叠后地震波数据x(t)的小波包变换W_x(p,b)为：

其中，为母小波的共轭函数，ω_pb(t)＝|p|^sω(|p|^s(t-b))e^2πi(t-b)p，参数s∈(1/2,1)，|p|≥1；p,b∈R²；

叠后地震波数据x(t)的局域小波向量估计v_x(p,b)为：

其中，W_x(p,b)≠0。

5.如权利要求4所述的一种叠后地震波相对波阻抗反演方法，其特征在于，重构的叠后地震波数据为：

其中，为ω_pb(t)的双边框架；

U_k为第k个聚类子集。

6.如权利要求1所述的一种叠后地震波相对波阻抗反演方法，其特征在于，获取目标层位的叠后地震波数据的过程为：综合利用地质信息、测井以及合成地震记录，标定叠后地震波数据的目标层位，确定叠后地震波数据的分析时间范围。