[发明专利]一种提高地震数据初至波质量的初至自动拾取方法

权利要求书

1.一种提高地震数据初至波质量的初至自动拾取方法，其特征在于，包括以下步骤：

(1)、通过时间窗对输入炮集数据进行截取，获得包含初至波的地震数据；

(2)、对于所选取目标炮，其震源位于S位置处，将检波器A和B接收到的参考道和目标道分别通过小波变换映射到小波域；

(3)、在小波域中对参考道和目标道的小波系数进行互相关运算；

(4)、在所有稳相震源位置上对步骤(3)中的互相关运算做积分，得到虚拟道的小波系数，该虚拟道相当于震源位于A且在检波器B接收到的波场记录；

(5)、对创建的虚拟道小波系数谱进行小波域软阈值去噪；

(6)、将去噪后的虚拟道小波系数与参考道小波系数进行褶积运算；

(7)、在所有稳相检波点位置上对步骤(6)中的褶积运算做积分，生成超级虚拟道的小波系数，该超级虚拟道等同于震源位于S且在检波器B接收到的波场记录；

(8)、对生成的超级虚拟道小波系数谱进行小波域软阈值去噪；

(9)、将去噪后的超级虚拟道小波系数逆变换回时空域，得到超级虚拟初至波；

(10)、对目标炮集逐道重复步骤(2)到步骤(9)，得到该炮集的超级虚拟初至波；

(11)、对目标炮集的超级虚拟初至波做反褶积滤波处理；

(12)、通过修改的能量比方法对超级虚拟初至波进行初至拾取；

步骤(12)中修改的能量比计算公式为：

MER_i＝[ER_i·abs(A_i)]³

其中，A_i是道集上第i个样点的振幅，i＝1,2,3,...,N，ER_i是沿着两个滑动的前、后时窗计算出两个时窗内累计能量的比值，定义为：

其中，A_j是道集上第j个样点的振幅，n是一个滑动时窗内样点的数目，前时窗和后时窗长度相同，测试点位于两个时窗的中心。

2.如权利要求1所述的提高地震数据初至波质量的初至自动拾取方法，其特征在于步骤(1)的时间窗长度为预期初至时间之前的一个周期及其之后的一到两个周期。

3.如权利要求1所述的提高地震数据初至波质量的初至自动拾取方法，其特征在于步骤(1)所述时间窗的预估位置可通过道集的偏移距信息除以近地表的近似平均速度实现自动确定。

4.如权利要求1所述的提高地震数据初至波质量的初至自动拾取方法，其特征在于步骤(12)中滑动时窗的长度取决于信号的周期，对于信噪比低的数据，需要通过测试来决定合适的窗口长度。

5.如权利要求1所述的提高地震数据初至波质量的初至自动拾取方法，其特征在于步骤(12)中沿着道集计算出MER_i，确定MER_i属性的峰值，该峰值所对应的时间即为初至波拾取时间。