[发明专利]初至时间拾取方法

权利要求书

1.一种初至时间拾取方法，包括：

基于原始地震道，计算初至搜索中心时间，获得搜索窗口半径，进而确定初至搜索范围；

基于所述初至搜索中心时间与所述原始地震道的数据，计算所述原始地震道的稳定因子；

在所述原始地震道中的所述初至搜索范围内，针对每个样点，分别在所述样点的前后划分一个目标时窗，进而计算所述样点的地震属性值；

基于每个样点的地震属性值与所述稳定因子，计算每个样点对应的能量比；

提取具有最大的能量比的样点所对应的样点时间作为初至时间；

其中，所述初至搜索中心时间为：

其中，m表示初至搜索中心时间，x表示原始地震道的偏移距，v表示地下近地表速度的估计值；

所述稳定因子为：

其中，β表示稳定因子，s_i表示第i个样点的样点值，i表示样点序号，N表示样点个数；

地震属性值为前后目标时窗的互相关的最大值，或前后两个目标时窗的余弦距离或相关距离。

2.根据权利要求1所述的初至时间拾取方法，其中，所述能量比为：

其中，ER(i)表示第i个样点对应的能量比，pdist(j)表示第j个样点的地震属性值，i1表示第i个样点的前一个目标时窗的起始样点的样点序号，i2表示第i个样点的后一个目标时窗的终止样点的样点序号。

3.根据权利要求2所述的初至时间拾取方法，其中，所述地震属性值为：

其中，pdist1(i)表示第i个样点的地震属性值，N_c表示每个目标时窗包含的样点数，s_l表示第l个样点的样点值，s_l+k表示第l+k个样点的样点值，k为整数。

4.根据权利要求2所述的初至时间拾取方法，其中，所述地震属性值为：

其中，pdist2(i)表示第i个样点的地震属性值，N_c表示每个目标时窗包含的样点数，s_l表示第l个样点的样点值，s_l+k表示第l+k个样点的样点值，k为整数，表示第i个样点的前一个目标时窗的样点值的平均值，表示第i个样点的后一个目标时窗的样点值的平均值。

5.根据权利要求3或4所述的初至时间拾取方法，其中，所述k为1或2。