一种适用于AVAF反演的频率域地震道集生成方法与流程

技术特征：

1.一种适用于AVAF反演的频率域地震道集生成方法，其特征在于，包括以下步骤：

获取地震数据；

对所述地震数据进行S变换，得到地震数据的S谱；

针对反演需求设计窗函数，与S谱相乘再积分；

对相乘的结果进行傅里叶反变换，得到频率域道集。

2.根据权利要求1所述的适用于AVAF反演的频率域地震道集生成方法，其特征在于，所述S变换的公式如下：

$S(\mathrm{\τ},v)=\underset{-\∞}{\overset{+\∞}{\∫}}g\left(t\right)\frac{\left|v\right|}{\sqrt{2\mathrm{\π}}}{e}^{-\frac{{(\mathrm{\τ}-t)}^2{v}^2}{2}}{e}^{-i2\mathrm{\πvt}}\mathrm{dt}---\left(1\right)$

其中τ和v分别是S变换之后的时间坐标和频率坐标。

3.根据权利要求1所述的适用于AVAF反演的频率域地震道集生成方法，其特征在于，所述对所述地震数据进行S变换使用有限频率域快速S变换方法，选择地震频带作为需要处理的频率进行分解，则地震数据的S变换可以从α域经过反傅里叶变换得到：

$S(\mathrm{\τ},v)=\underset{-\∞}{\overset{+\∞}{\∫}}\mathrm{\α}(v^{\′},v)\·e^{i2\mathrm{\π}{v}^{\′}\mathrm{\τ}}d{v}^{\′}---\left(2\right).$

4.根据权利要求1所述的适用于AVAF反演的频率域地震道集生成方法，其特征在于，所述针对反演需求设计窗函数，与S谱相乘再积分的步骤中，高斯窗的选择依赖于地震数据，同时与反演子波的高斯窗函数保持一致；

高斯窗函数的个数依赖于反演的频率数，当反演的频率数为N时，窗函数的个数也是N。

5.根据权利要求4所述的适用于AVAF反演的频率域地震道集生成方法，其特征在于，所述针对反演需求设计窗函数，与S谱相乘再积分的步骤，具体公式如下式：

$g_{v_0}\left(t\right)={\∫}_{-\∞}^{\∞}\left\{{\∫}_{-\∞}^{\∞}S(\mathrm{\τ},v)\mathrm{d\τ}\right\}{e}^{i2\mathrm{\πvt}}W(v_0,{\mathrm{\σ}}^2)\mathrm{dv}---\left(3\right)$

其中S(τ,ν)是待分频道集的S谱，W(ν₀,σ²)为窗函数，ν₀为窗函数的均值频率，σ²为窗函数的方差。

6.根据权利要求5所述的适用于AVAF反演的频率域地震道集生成方法，其特征在于，通过公式(3)得到的加窗反变换结果就是最后得到的分频率后的地震道集，该地震道集配合频率域的子波适于进行AVAF反演。

7.根据权利要求2所述的适用于AVAF反演的频率域地震道集生成方法，其特征在于，首先使用公式(1)定义的S变换将地震道集变换到时频域S域。

8.根据权利要求3所述的适用于AVAF反演的频率域地震道集生成方法，其特征在于，在对时频域的数据进行反傅里叶变换时，首先将其在局部时间域上进行全局积分，得到了公式(3)中的然后在频率域进行加窗积分得到