[发明专利]一种时频域可控震源近地表吸收补偿方法

权利要求书

1.一种时频域可控震源近地表吸收补偿的方法，特点是通过以下步骤实现：

1）使用可控震源激发并记录地震波数据；

2）将地震波数据在时频域进行分解；所述的分解是在频率域采用余弦函数构造的分频滤波函数对地震波数据进行分频处理；

3）对步骤2）分解后的各个频段的分频数据进行共炮集频率域中值滤波，产生每一个频段的炮集模型道；

4）采用高阶e指数函数作为近地表吸收衰减的拟合函数，利用炮集模型道数据拟合求取每一炮各个频段的统计吸收衰减函数；

5）将扫描信号在时频域进行分解为不同的频段，利用吸收衰减函数在相应频段形成全频带的扫描信号，处理得到含近地表吸收的扫描信号，再与原始扫描信号相关，得到含近地表吸收的可控震源模拟子波；

6）将步骤5）得到的每一炮的含近地表吸收的可控震源模拟子波，应用最小平方算法求出反褶积因子，应用至对应的炮集数据；

7）对所有炮数据进行步骤1）-6）的处理，实现时频域可控震源近地表吸收补偿。

2.根据权利要求1的方法，特点是步骤2）所述的分频滤波函数是在频率域将0-f_N折叠频率平均分为17～25个频段，在每个频段用余弦函数产生一个滤波器，相邻每个滤波器重叠一半，构造覆盖全部频率的多个滤波器。

3.根据权利要求1的方法，特点是步骤3）所述的中值滤波是将某一炮点或检波点各个地震道按采样点排列后取中值，获得炮集模型道。

4.根据权利要求1的方法，特点是步骤4）所述的拟合计算公式如下：

ϵ=Σt=1tN{lnMIDj=1,2,···MA[Xij(t,fk)]-αi(t,fk)}2---(2)]]>

式中:A[·]为振幅；—中值统计滤波；X_ij(t,f_k)为分频后的炮集数据；f_k为第k个频带；t_n为采样时间；i—炮号；j—道号；

其中:α_i(t,f_k)＝α₀(f_k)+α₁(f_k)t+…+α_p(f_k)t^p    （3）

式中：α_i(t,f_k)为第i炮的拟合函数；p为拟合阶数，为3～5；α_p为拟合系数；f_k为第k个频段；k=1，2，…，N，k是指步骤2）所述分频段的序号，N是频段的个数。

5.根据权利要求1的方法，特点是步骤5）所述的分解是将扫描信号按步骤2）的方法分解为不同的频段。