一种基于时频分析技术的井控有色拓频方法与流程

技术特征：

1.一种基于时频分析技术的井控有色拓频方法，其特征在于包括如下步骤：

步骤S1：首先对深度域测井数据的声波数据求取反射系数，通过时深转换转换到时间域，对所获得的反射系数做傅里叶变换，然后获得反射系数的振幅谱，统计分析研究区测井资料的反射系数振幅谱，利用最优化拟合方法获得修正因子C(f)＝a+b*arctan(f)中的未知参数a和b；

步骤S2：对地震数据进行时频分析，对时域地震数据x(t)进行S变换，选用的窗函数为可得时频域处理结果：

$X(\mathrm{\τ},f)={\∫}_{-\mathrm{\∞}}^{+\mathrm{\∞}}x\left(t\right)g_f(\mathrm{\τ}-t)\exp (-i2\mathrm{\π}ft)dt$

然后获得时间t，频率f下的地震数据对应的振幅谱Amp(t,f)和相位谱P(t,f)：

$Amp(t,f)=\sqrt{{\left(real\right(X\left(t,f\right)\left)\right)}^2+{\left(imag\right(X\left(t,f\right)\left)\right)}^2}$

$P(t,f)={\tan }^{-1}\left(\frac{imag\left(X\right(t,f\left)\right)}{real\left(X\right(t,f\left)\right)}\right)$

然后对地震数据振幅谱有色修正，对S变换获得的振幅谱采用有色因子进行校正，获得校正后的有色振幅谱CF(t,f)；

步骤S3：对获得的有色振幅谱R_HC和原始的相位谱I_HC相结合：

R_HC＝CF(t,f)*cos(P(t,f))

I_HC＝CF(t,f)*sin(P(t,f))

HC(t,f)＝complex(R_HC,I_HC)；

步骤S4：最后完成地震道重建，利用逆S变换实现从时频域到时域的转换，得到其拓频处理后的地震道。