[发明专利]一种基于宽平谱地震子波整形的地震分辨率增强方法

权利要求书

1.一种基于宽平谱地震子波整形的地震分辨率增强方法，其特征在于采用以下具体步骤：

⑴对叠后数据的地震道x(t)进行傅里叶变换求其频谱，并进行归一化处理，获得归一化处理后的频谱X(f)；

⑵采用多项式拟合求取X(f)的拟合频谱F(f)，即最优的初始子波频谱；

⑶选取归一化处理后的频谱的阈值M，确定归一化处理后的频谱的有效频带范围[f_a,f_b]，f_a＜f_b；

[f_a,f_b]＝{X(f)＝M}

⑷根据确定的有效频带范围，利用宽平谱地震子波可自适应调节的特点，设计宽平谱T(f)，其中，建立的宽平谱地震子波的时间域解析式为：

式中，f_a、f_b分别为平谱段的起始频率和截止频率，f_c为终止频率；

对应的频率域解析式为：

式中，a,b分别用于控制地震子波谱低、高截频段的削尖速度，N为频带的宽度；

在该步骤中，为了使宽平谱的高频段窗边变化趋势与拟合频谱的衰减变化趋势更好地吻合，采取以下步骤构建最优的宽平谱T(f):

a)求取拟合频谱F(f)的归一化处理结果G(f)；

b)确定G(f)的主频位置f_m；

c)采用大于主频位置f_m处G(f)的高频段曲线代替宽平谱高频段的衰减下降趋势，构造最优的宽平谱T(f)，其表达式为：

⑸利用宽平谱T(f)与拟合频谱F(f)求取权系数c(f)；

c(f)＝T(f)/F(f)

⑹将归一化处理后的频谱X(f)与权系数c(f)进行相乘，求取整形处理之后的宽频谱Y(f)；

Y(f)＝X(f)·c(f)

⑺将整形处理之后的宽频谱进行反傅里叶变换，即可获得高分辨处理后的地震道y(t)；

y(t)＝IFT[Y(f)]

⑻重复步骤⑴至⑺，直到处理完地震数据中的所有地震道，即可得到整个工区的高分辨地震数据结果。