一种地震勘探特征子波整形提高分辨率方法与流程

技术特征：

1.一种地震勘探特征子波整形提高分辨率方法，该方法包含以下步骤：先进行子波处理提升高频成分的子波的振幅，再用复数Paul子波替换原来的子波。

2.根据权利要求1所述的一种地震勘探特征子波整形提高分辨率方法，所述提升高频成分的子波的振幅包含以下步骤：

(1)根据实际地震数据计算振幅谱；

(2)计算高频振幅补偿曲线；

(3)对数据进行小波变换得到小子波；

(4)对所述小子波进行振幅补偿，以期望的高频振幅补偿曲线为基准，将每一小子波的最大振幅与到所述期望的高频振幅补偿曲线对应值对应，求出该小子波的振幅补偿因子，然后计算每一小子波补偿后的振幅谱。

3.根据权利要求2所述的一种地震勘探特征子波整形提高分辨率方法，其特征在于所述步骤(1)中的数据是指在选择的分析时窗内，取多道数据：所述振幅谱是多道的傅氏振幅谱；所述步骤(2)的高频振幅补偿曲线计算方式是，首先根据所述傅氏振幅谱求得地层等效吸收衰减曲线，然后根据所述地层等效吸收衰减曲线计算高频振幅补偿曲线。

4.根据权利要求3所述的一种地震勘探特征子波整形提高分辨率方法，其特征在于高频振幅补偿曲线的数值可调。

5.根据权利要求4所述的一种地震勘探特征子波整形提高分辨率方法，其特征在于高频振幅补偿曲线的数值可调降为原来的2/3或1/2。

6.根据权利要求3所述的一种地震勘探特征子波整形提高分辨率方法，其特征在于所述步骤(4)中的计算每一小子波补偿后的振幅谱是通过求得振幅补偿因子来计算的。

7.根据权利要求6所述的一种地震勘探特征子波整形提高分辨率方法，其特征在于所述振幅补偿因子是通过如下步骤获得：

1)进行多道数据的小波变换，得到不同频率的小波振幅谱；

2)在所述小波振幅谱中选取最大主振幅值(Amax(f0))；

3)在fi大于f0的范围内，求取其他小子波的最大主振幅Amax(fi)；

4)每一个小子波的振幅补偿因子k(fi)为，

$k\left(fi\right)=\frac{Amax\left(f0\right)}{Amax\left(fi\right)}\×\mathrm{\α}$

其中，i＝1，2，3，......，0＜α＜1，fi＞f0；

根据上式可计算高频振幅补偿曲线，α值即代表了调降该高频振幅补偿曲线，预设一条期望的高频振幅补偿曲线，将所述每一小波的最大振幅拟合到所述期望高频振幅补偿曲线上，从而求得上述公式中的α值，得到该条高频振幅补偿曲线针对频率fi的该小子波的振幅补偿因子k(fi)；5)在单道数据上，在选择的时窗内，进行小波变换，得到不同频率fi的小子波振幅谱A(fi)；针对频率fi的小子波，用振幅补偿因子k(fi)与该小子波中的所有振幅相乘，得到补偿后小子波的频率为fi的振幅，

A′(fi)＝k×A(fi)，(k＝k(fi))

对每一个频率fi的振幅作如上处理，得到所述小子波补偿后的振幅谱；

6)对每一个小子波做与上述步骤相同的处理，得到每一个小子波的振幅谱。

8.根据权利要求7所述的一种地震勘探特征子波整形提高分辨率方法，其特征在于

将这些新小子波重构，得到高频振幅补偿后的反射波。

9.根据权利要求8所述的一种地震勘探特征子波整形提高分辨率方法，其特征在于对所述重构的反射波进行随机噪声压缩处理。

10.根据权利要求9所述的一种地震勘探特征子波整形提高分辨率方法，其特征在于

所述用主频高-频带宽的复数Paul子波替换原来的子波，采用如下的公式进行复数Paul子波(即_PW)分频重构：

地震反射信号s表示为地震子波r和反射系数b序列的褶积，如下算式：

S_PW(f)＝R(f)·B_PW(f)。