1.一种基于多尺度形态学的初至自动拾取方法，其特征在：包括如下步骤：
步骤1、使用能量比法获得初至初始特征值，其计算公式为：
si＝|(b/a)×(b-a)|(1)
其中，
步骤2、对能量比法获得的特征值进行膨胀运算，其计算公式为：
其中，
其中，a、b、c、d可取任意正整数，其大小由初至特征值分布特征而定，其他位置元素均为零；
步骤3、对膨胀后的特征值进行道连通计算，获得各连通区域的道连通数；
步骤4、剔除道连通较少的特征值；
步骤5、保留特征值映射到原始数据中，并局部寻找最优初至。
2.根据权利要求1所述基于多尺度形态学的初至自动拾取方法，其特征在：步骤1中，采用的能量比法中n为计算时窗，该参数越大，抗噪能力越强，但是初至识别精度降低；参数越小，初至识别精度高，但是抗噪能力差，实际应用中，该参数应根据实际初至特征进行选择，也可以经过实验获取。
3.根据权利要求1所述基于多尺度形态学的初至自动拾取方法，其特征在：进行膨胀运算之前，需要首先将特征值剖面进行二值化处理，将每道中最大特征值设置为1，小于0.5的特征值设置为零，大于0.5的特征值，设置为1。
4.根据权利要求1所述基于多尺度形态学的初至自动拾取方法，其特征在：步骤2中，膨胀结构元素a参数a、b、c、d的大小与期望保留的特征值有关，值越大，越多的特征值被连通，保留的特征值越多，噪音较重的时候，算子参数建议使用较小值；信噪比高的时候，算子参数采用较大值，地震剖面采用的a为5*5，即a＝b＝c＝d＝1。
5.根据权利要求1所述基于多尺度形态学的初至自动拾取方法，其特征在：步骤2中,为避免结构元素a的复杂性，在膨胀结构元素应用之前，首先要对特征值进行道方向向下复制的操作，即在当前特征值位置以下若干点都赋值为当前特征值，该操作保证了进行膨胀运算时，道方向的连通性。
6.根据权利要求1所述基于多尺度形态学的初至自动拾取方法，其特征在：步骤3中，进行的是道连通的计算，即经过连通后，需要判断的是，当前连通区域有多少道是连通的，所谓的连通较少的特征值进行剔除，是对连通道较少的区域对应特征值进行剔除。
7.根据权利要求1所述基于多尺度形态学的初至自动拾取方法，其特征在：步骤4中，经过连通处理后，对道连通较少的区域进行剔除后，使用剩余的膨胀结果与原始特征值进行点乘运算，可以剔除掉异常特征值。
8.根据权利要求1所述基于多尺度形态学的初至自动拾取方法，其特征在：步骤5中，剔除后的特征值，并不是最终的初至，需要将这些特征值位置映射到原始的地震剖面上，进行局部寻优，以确定最佳的初至。
9.根据权利要求8所述基于多尺度形态学的初至自动拾取方法，其特征在：所述局部寻优，是在当前特征值位置前后各5个点内寻找正振幅最大值对应位置，该位置对应了初至波峰。
10.根据权利要求9所述基于多尺度形态学的初至自动拾取方法，其特征在：所有初至拾取结束后，可以通过局部拟合插值后再补充部分缺失初至，对大范围连续缺失初至无法进行插值处理。