[发明专利]一种提高初至波拾取精度和效率的方法

权利要求书

1.一种提高初至波拾取精度和效率的方法，其特征在于采用以下步骤：

1)采集地震数据并得到三维叠前炮集地震数据、以及炮点与接收点关系信息；

2)选择覆盖所有横炮线的纵炮线为种子炮线，从种子炮线中选择信噪比高的一炮为种子炮，得到种子炮的初至时间；

3)以种子炮的初至时间为依据，将种子炮所覆盖的地面观测区域分为方位角较小的子区域，在每个子区域内，按照炮检距大小划分为段，对每一炮检距段内各个道的初至时间进行高程校正，利用最小二乘法对每段内校正后的初至时间炮检距对进行直线拟合，得到直线的斜率，将斜率的倒数作为段初至曲线的视速度，利用距离加权插值计算每道上的瞬时速度，建立视速度场；

4)基于地表一致性关系和初至曲线视速度场，由种子炮初至时间得到非种子炮初至时间；

所述的非种子炮初至时间对于不满足地表一致性关系的未知道，利用已知道上的初至时间，按照折射关系外推未知道的初至时间；

5)利用推算出的非种子炮初至时间作约束，用最大能量初至波拾取方法对非种子炮进行初至拾取；

所述的最大能量初至波拾取方法是将原始地震数据某一地震道记录及进行Hilbert变换后记为：X_i(n)和 

地震数据包络能量为：

时间窗样点长为L+1的包络能量和为：

包络能量和比值函数为：

ME_i，k＝E_i，k/E″_i，k    (4)

搜索ME_i，k的最大值，如果地震数据在k点的值X_i(k)大于零，那么这个最大值的样点为需要拾取的初至时间；

式中：E″_i，k和E_i，k分别代表比值前后L+1个样点的包络能量和；i为地震道序号；n为样点序号；L为以样点数表示计算能量的时窗长度；k为计算能量比值的样点序号；

6)确定非种子炮拾取后的初至时间；

7)把确定非种子炮作为种子炮，按照步骤3)至步骤7)对种子炮附近的非种子炮进行初至拾取，逐渐完成整条炮线的初至拾取；

8)把整条炮线作为种子炮线按照步骤3)至步骤8)完成其它非种子炮线的初至拾取，逐渐递推，完成整个工区的初至时间拾取。

2.根据权利要求1所述的提高初至波拾取精度和效率的方法，其特征在于步骤1)所述的信息包括炮点和接收点的X、Y坐标及炮线号和接收点线号。