[发明专利]宽频时空变地震子波提取方法

权利要求书

1.宽频时空变地震子波提取方法，其特征在于，该宽频时空变地震子波提取方法包括：

步骤1，获取地震数据，并对地震数据进行预处理；

步骤2，从被预处理后的地震数据中获取地震解释层位，并在地震解释层位中截取层位时窗，层位时窗表示目标区域；

步骤3，在目标区域中，划定一个或多个子波空变区域；

步骤4，在一个或多个子波空变区域中，划定一个或多个时变滑动延拓时窗；

步骤5，将所有子波空变区域中，将所有时变滑动延拓时窗得到的估算地震子波进行组合，得到宽频时空变地震子波。

2.根据权利要求1所述的宽频时空变地震子波提取方法，其特征在于，在步骤1中，地震数据为叠后地震数据，可以为全叠加叠后地震数据，也可为分偏移距叠加地震数据、分角度叠加地震数据或分方位角叠加地震数据；预处理包括去噪、提取道头信息。

3.根据权利要求1所述的宽频时空变地震子波提取方法，其特征在于，在步骤3中，子波空变区域，是指根据地质特征、地震资料采集处理情况和研究需要划分子波空变区域，若仅提取时变子波，则将目标区划分为单一的空变区域；在不同的空变区域中，分别提取时变子波，得到整个目标区的时空变子波。

4.根据权利要求1所述的宽频时空变地震子波提取方法，其特征在于，在步骤4中，时变滑动延拓时窗是指对于每一个空变区域，假设起始时变采样点为τ₀，时变滑动间隔为Δt，采样次数为N，终止时变采样点为τ₀+N·Δt，延拓时窗长度为2ΔT，则时变采样点τ的集合表示为：τ∈{τ₀+i·Δt|i＝0，1，2…N}，时变滑动延拓时窗表示为[τ-ΔT，τ+ΔT]，通过滑动时变采样点，在时变滑动延拓时窗范围内提取该时变采样点的地震子波，最终组成时变子波。

5.根据权利要求4所述的宽频时空变地震子波提取方法，其特征在于，在步骤4中，对于一个或多个时变滑动延拓时窗中的至少一个，执行以下步骤：

步骤4.1，使用窗函数，从一个或多个时变滑动延拓时窗中的至少一个中提取多道随机短时窗地震数据；

步骤4.2，对于多道随机短时窗地震数据，进行频域滤波处理和同态域平均处理，以获得同态域子波数据；

步骤4.3，变换同态域子波数据，以获得时域子波数据；

步骤4.4，对时域子波数据，进行子波整形处理，以获得估算地震子波。

6.根据权利要求5所述的宽频时空变地震子波提取方法，其特征在于，在步骤4.1中，窗函数是指汉宁窗函数Win(t)，对于时变采样点τ_i，其表达式如下式给出：

式中，t₁、t₂分别为随机短时窗起、止时间，t为采样点时间，

7.根据权利要求6所述的宽频时空变地震子波提取方法，其特征在于，在步骤4.1中，多道随机短时窗，是指在该空变区域内随机选取若干道数据，在时变滑动延拓时窗范围内随机选取若干个短时窗，使用汉宁窗函数得到短时窗地震数据，其中随机的道数和短时窗数由用户根据实际需要给定，随机道数不少于总道数的10％，单道随机短时窗数不少于25个。

8.根据权利要求5所述的宽频时空变地震子波提取方法，其特征在于，在步骤4.2中，同态域子波数据，是指对于某一时变采样点τ_i，s_j(t，τ_i)为步骤4.1中得到的第j个随机短时窗地震数据，s_j(t，τ_i)由傅里叶变换并进行降噪滤波得到频率域的随机短时窗地震数据表示为S_j(ω，τ_i)，则同态域子波由下式获得：

式中，为同态域子波，N为随机短时窗总数。