[发明专利]地震数据反射系数标签的自动标定方法

权利要求书

1.一种地震数据反射系数标签的自动标定方法，其特征在于，包括以下步骤：

根据褶积模型重构变子波的地震记录；

将地震信号自适应地划分为多个平稳的信号片段；

分别在每个平稳的片段内提取瞬时子波；

以反射界面位置为约束，通过提取到的瞬时子波对所述地震记录进行反演计算获得反射系数，将所述反射系数作为所述地震记录的标签。

2.根据权利要求1所述的地震数据反射系数标签的自动标定方法，其特征在于，在某个地层片段内子波是常子波，则地震记录S(t)满足下列表达式：

其中，W^j(t)为第j个地层片段内的子波；

L为地层片段的个数；

为反射系数；

t为时间；

Δt为采样时间间隔；

kΔt为第k个采样时刻；

c_k为第k时刻反射系数的值；

δ为脉冲函数符号；

k，M₁和M₂均为自然数。

3.根据权利要求2所述的地震数据反射系数标签的自动标定方法，其特征在于，对地震记录S(t)进行傅里叶域变换，获得变子波地震记录的傅里叶域表达式：

其中，W^j(w)为第j个地层片段内的频率域瞬时子波；

t_k为第k个采样上的时刻，t_k＝kΔt；

i为虚数单位。

4.根据权利要求3所述的地震数据反射系数标签的自动标定方法，其特征在于，第j个地层片段内，瞬时子波的振幅谱s^j(w)满足下列表达式：

5.根据权利要求1-4中任一项所述的地震数据反射系数标签的自动标定方法，其特征在于，通过构造自适应Gabor标架，以使得覆盖在此区域的地震子波是平稳的来将地震信号自适应地划分为多个平稳的信号片段。

6.根据权利要求1-4中任一项所述的地震数据反射系数标签的自动标定方法，其特征在于，将地震信号自适应地划分为多个平稳的信号片段包括以下子步骤：

在各个采样点均给定一个相同的Gabor小原子，其原子窗为高斯窗；

采用希尔伯特变换的方式求取地震子波的瞬时频率；

根据所述地震子波的瞬时频率获得地震子波的振幅包络加权频率，对所述振幅包络加权频率曲线采取保边平滑滤波处理；

将所述高斯窗作用于保边平滑滤波处理后的结果上，以将地震信号自适应地划分为多个平稳的信号片段。

7.根据权利要求6所述的地震数据反射系数标签的自动标定方法，其特征在于，地震子波的瞬时频率f(t)满足下列表达式：

其中，s(t)是地震信号的实部；

s^*(t)是地震信号的虚部；

a²(t)是瞬时包络振幅；

a²(t)满足下列表达式：

a²(t)＝s²(t)+s²*(t) (13)。

8.根据权利要求6所述的地震数据反射系数标签的自动标定方法，其特征在于，地震子波的振幅包络加权频率f'(t)满足下列表达式：

其中，W(t)是瞬时振幅包络。

9.根据权利要求6所述的地震数据反射系数标签的自动标定方法，其特征在于，原子窗满足以下高斯函数：

其中，T为小原子窗的半周期。

10.根据权利要求6所述的地震数据反射系数标签的自动标定方法，其特征在于，在每个平稳的信号片段内子波均为0相位。